source: LMDZ6/branches/LMDZ-tracers/libf/phylmd/inlandsis/VARySV.F90 @ 3871

Last change on this file since 3871 was 3792, checked in by evignon, 4 years ago

Ajout de INLANDSIS, nouvelle interface entre LMDZ et la neige de SISVAT
Etienne, 04/01/2021
File size: 5.9 KB
Line 
1MODULE VARySV
2
3USE VAR_SV, only : klonv, nsol, nsno
4
5
6IMPLICIT NONE
7
8! +
9      INTEGER, DIMENSION(:),SAVE,ALLOCATABLE ::   NLaysv  ! New   Snow     Layer   Switch   
10!$OMP THREADPRIVATE(NLaysv)
11      INTEGER, DIMENSION(:),SAVE,ALLOCATABLE ::   i_thin  ! Index of the thinest Layer     
12!$OMP THREADPRIVATE(i_thin)
13      INTEGER, DIMENSION(:),SAVE,ALLOCATABLE ::   LIndsv  ! Contiguous Layer relative Index
14!$OMP THREADPRIVATE(LIndsv)
15      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   albisv  ! Integrated Surface Albedo
16!$OMP THREADPRIVATE(albisv)
17      !hj prelim                                           
18      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   alb1sv  ! Surface Albedo VIS     
19!$OMP THREADPRIVATE(alb1sv)
20      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   alb2sv  ! Surface Albedo NIR     
21!$OMP THREADPRIVATE(alb2sv)
22      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   alb3sv  ! Surface Albedo FIR     
23!$OMP THREADPRIVATE(alb3sv)
24
25      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   albssv  ! Soil               Albedo [-]   
26!$OMP THREADPRIVATE(albssv)
27      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   SoSosv  !? Surface Absorbed Solar Radiat. 
28!$OMP THREADPRIVATE(SoSosv)
29      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   Eso_sv  !? Emissivity of       Soil+Snow   
30!$OMP THREADPRIVATE(Eso_sv)
31      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   HSv_sv  ! Sensible Heat Flux  [W/m2]     
32!$OMP THREADPRIVATE(HSv_sv)
33      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   HLv_sv  ! Latent   Heat Flux  [W/m2]     
34!$OMP THREADPRIVATE(HLv_sv)
35      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   HSs_sv  !? Sensible Heat Flux (t)         
36!$OMP THREADPRIVATE(HSs_sv)
37      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   HLs_sv  !? Latent   Heat Flux (t)         
38!$OMP THREADPRIVATE(HLs_sv)
39      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   sqrCm0  ! in Neutral Drag Coef.Moment.   
40!$OMP THREADPRIVATE(sqrCm0)
41      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   sqrCh0  ! in Neutral Drag Coef.Heat       
42!$OMP THREADPRIVATE(sqrCh0)
43      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   Lx_H2O  ! Latent Heat of Vaporiz./Sublim.
44!$OMP THREADPRIVATE(Lx_H2O)
45      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   ram_sv  ! Aerodyn.Resistance (Moment.)   
46!$OMP THREADPRIVATE(ram_sv)
47      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   rah_sv  ! Aerodyn.Resistance (Heat)       
48!$OMP THREADPRIVATE(rah_sv)
49      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   Fh__sv  ! Stability Function             
50!$OMP THREADPRIVATE(Fh__sv)
51      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   dFh_sv  ! Stability Function (Deriv.)     
52!$OMP THREADPRIVATE(dFh_sv)
53      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   Evp_sv  !x Evaporation        [kg/m2]     
54!$OMP THREADPRIVATE(Evp_sv)
55      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   EvT_sv  !x Evapotranspiration [kg/m2]     
56!$OMP THREADPRIVATE(EvT_sv)
57      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   LSdzsv  ! Land/Sea Vert. Discretiz. Fact.
58!$OMP THREADPRIVATE(LSdzsv)
59      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   Tsrfsv  ! Surface    Temperature         
60!$OMP THREADPRIVATE(Tsrfsv)
61      REAL,SAVE,ALLOCATABLE  ::   sEX_sv(:,:)  ! Verticaly Integr.Extinct.Coef. 
62!$OMP THREADPRIVATE(sEX_sv)
63      REAL,SAVE,ALLOCATABLE  ::   zzsnsv(:,:)    ! Snow  Pack Thickness      [m]   
64!$OMP THREADPRIVATE(zzsnsv)
65      REAL,SAVE,ALLOCATABLE  ::   psi_sv(:,:)    ! Soil   Water        Potential   
66!$OMP THREADPRIVATE(psi_sv)
67      REAL,SAVE,ALLOCATABLE  ::   Khydsv(:,:)    ! Soil   Hydraulic    Conductiv. 
68!$OMP THREADPRIVATE(Khydsv)
69      REAL, DIMENSION(:),SAVE,ALLOCATABLE    ::   EExcsv  ! Energy in Excess, current       
70!$OMP THREADPRIVATE(EExcsv)
71
72CONTAINS
73  SUBROUTINE INIT_VARySV
74  IMPLICIT NONE
75      INTEGER                ::   ikl
76
77      ALLOCATE(NLaysv(klonv))  ! New   Snow     Layer   Switch   
78      ALLOCATE(i_thin(klonv))  ! Index of the thinest Layer     
79      ALLOCATE(LIndsv(klonv))  ! Contiguous Layer relative Index
80      ALLOCATE(albisv(klonv))  ! Integrated Surface Albedo
81      !hj prelim                                           
82      ALLOCATE(alb1sv(klonv))  ! Surface Albedo VIS     
83      ALLOCATE(alb2sv(klonv))  ! Surface Albedo NIR     
84      ALLOCATE(alb3sv(klonv))  ! Surface Albedo FIR       
85
86      !
87      ALLOCATE(albssv(klonv))  ! Soil               Albedo [-]   
88      ALLOCATE(SoSosv(klonv))  !? Surface Absorbed Solar Radiat. 
89      ALLOCATE(Eso_sv(klonv))  !? Emissivity of       Soil+Snow   
90      ALLOCATE(HSv_sv(klonv))  ! Sensible Heat Flux  [W/m2]     
91      ALLOCATE(HLv_sv(klonv))  ! Latent   Heat Flux  [W/m2]     
92      ALLOCATE(HSs_sv(klonv))  !? Sensible Heat Flux (t)         
93      ALLOCATE(HLs_sv(klonv))  !? Latent   Heat Flux (t)     
94      ALLOCATE(sqrCm0(klonv))  ! in Neutral Drag Coef.Moment.   
95      ALLOCATE(sqrCh0(klonv))  ! in Neutral Drag Coef.Heat       
96      ALLOCATE(Lx_H2O(klonv))  ! Latent Heat of Vaporiz./Sublim.
97      ALLOCATE(ram_sv(klonv))  ! Aerodyn.Resistance (Moment.)   
98      ALLOCATE(rah_sv(klonv))  ! Aerodyn.Resistance (Heat)       
99      ALLOCATE(Fh__sv(klonv))  ! Stability Function             
100      ALLOCATE(dFh_sv(klonv))  ! Stability Function (Deriv.)     
101      ALLOCATE(Evp_sv(klonv))  !x Evaporation        [kg/m2]     
102      ALLOCATE(EvT_sv(klonv))  !x Evapotranspiration [kg/m2]     
103      ALLOCATE(LSdzsv(klonv))  ! Land/Sea Vert. Discretiz. Fact.
104      ALLOCATE(Tsrfsv(klonv))  ! Surface    Temperature         
105      ALLOCATE(sEX_sv(klonv,-nsol:nsno+1))  ! Verticaly Integr.Extinct.Coef. 
106      ALLOCATE(zzsnsv(klonv,    0:nsno))    ! Snow  Pack Thickness      [m]   
107      ALLOCATE(psi_sv(klonv,-nsol:0   ))    ! Soil   Water        Potential   
108      ALLOCATE(Khydsv(klonv,-nsol:0   ))    ! Soil   Hydraulic    Conductiv. 
109      ALLOCATE(EExcsv(klonv))  ! Energy in Excess, current         
110
111      DO ikl=1,klonv
112        NLaysv(ikl)   = 0
113        i_thin(ikl)   = 0
114        LIndsv(ikl)   = 0
115      END DO
116
117  END SUBROUTINE INIT_VARySV
118
119END MODULE VARySV
120
Note: See TracBrowser for help on using the repository browser.