source: LMDZ6/trunk/libf/phylmd/inlandsis/VARxSV.F90 @ 3817

Last change on this file since 3817 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: 20.7 KB
RevLine 
[3792]1MODULE VARxSV                     
2
3USE VAR_SV, only : klonv, nsol, nsno, nb_wri
4                                                         
5IMPLICIT NONE                                                                               
6! +--SISVAT INPUT        Variables                                             
7! +  -----------------------------                                             
8 
9      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   LSmask  ! Land-Sea   Mask                 
10!$OMP THREADPRIVATE(LSmask)
11      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   isotSV  ! Soil       Type                 
12!$OMP THREADPRIVATE(isotSV)
13      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   iWaFSV  ! Soil       Drainage:(1,0)=(y,n)
14!$OMP THREADPRIVATE(iWaFSV)                                                                             
15      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   coszSV  ! Cosine of Sun zenithal Angle   
16!$OMP THREADPRIVATE(coszSV)
17      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   sol_SV  ! Downward  Solar    Radiation   
18!$OMP THREADPRIVATE(sol_SV)
19      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   IRd_SV  ! Downward  Longwave Radiation   
20!$OMP THREADPRIVATE(IRd_SV)
21                                                                               
22      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   drr_SV  ! Rain  Intensity       [kg/m2/s]
23!$OMP THREADPRIVATE(drr_SV)
24      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dsn_SV  ! Snow  Intensity       [kg/m2/s]
25!$OMP THREADPRIVATE(dsn_SV)
26      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dsnbSV  ! Idem, fraction, from Drift  [-]
27!$OMP THREADPRIVATE(dsnbSV)
28      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   esnbSV  ! Idem, fraction, from Drift  [-]
29!$OMP THREADPRIVATE(esnbSV)
30      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dbs_SV  ! Drift Amount            [kg/m2]
31!$OMP THREADPRIVATE(dbs_SV)
32      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dbs_Er  ! BS Erosion              [kg/m2]
33!$OMP THREADPRIVATE(dbs_Er)
34      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dbs_Ac  !  BS Accumulation         [kg/m2]
35!$OMP THREADPRIVATE(dbs_Ac)
36      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BrosSV  ! Buffer Snow Layer Density
37!$OMP THREADPRIVATE(BrosSV)
38      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BG1sSV  ! Buffer Snow Layer Dendr/Sphe[-]
39!$OMP THREADPRIVATE(BG1sSV)
40      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BG2sSV  ! Buffer Snow Layer Spher/Size[-][0.0001m]
41!$OMP THREADPRIVATE(BG2sSV)
42      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dz0_SV  ! dz0(Sastrugi dh)            [m]
43!$OMP THREADPRIVATE(dz0_SV)
44                                                                               
45      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   cld_SV  ! Cloudiness (seen from SBL)     
46!$OMP THREADPRIVATE(cld_SV)
47      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   za__SV  ! SBL Height                     
48!$OMP THREADPRIVATE(za__SV)
49      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   VV__SV  !(SBL Top)  Wind Velocity         
50!$OMP THREADPRIVATE(VV__SV)
51      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   VVs_SV  !(Sastr,V)  Relevance
52!$OMP THREADPRIVATE(VVs_SV)
53      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   RRs_SV  !(Sastr,V)  Counter
54!$OMP THREADPRIVATE(RRs_SV)
55      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   DDs_SV  !(Sastr,V)  Angle
56!$OMP THREADPRIVATE(DDs_SV)
57      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   TaT_SV  ! SBL Top   Temperature           
58!$OMP THREADPRIVATE(TaT_SV)
59      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   ExnrSV  ! Exner     Potential             
60!$OMP THREADPRIVATE(ExnrSV)
61      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dSdTSV  ! Sensible Heat Flux T Derivat.   
62!$OMP THREADPRIVATE(dSdTSV)
63      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dLdTSV  ! Latent   Heat Flux T Derivat.   
64!$OMP THREADPRIVATE(dLdTSV)
65      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rhT_SV  ! SBL Top   Air  Density         
66!$OMP THREADPRIVATE(rhT_SV)
67      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   QaT_SV  ! SBL Top   Specific Humidity     
68!$OMP THREADPRIVATE(QaT_SV)
69      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dQa_SV  ! SBL Flux  Limitation of Qa     
70!$OMP THREADPRIVATE(dQa_SV)
71      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   qsnoSV  ! SBL Mean  Snow       Content   
72!$OMP THREADPRIVATE(qsnoSV)                                                                   
73                                                                               
74      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   alb0SV  ! Soil    Albedo                 
75!$OMP THREADPRIVATE(alb0SV)
76      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   slopSV  ! Snow/Ice/Soil-Water Surf. Slope
77!$OMP THREADPRIVATE(slopSV)
78
79                                                             
80      REAL,SAVE                      ::   zSBLSV  ! SBL Height (Initial Value)     
81!$OMP THREADPRIVATE(zSBLSV)
82      REAL,SAVE                      ::   dt__SV  ! Time Step                       
83!$OMP THREADPRIVATE(dt__SV)
84      CHARACTER (len=18),SAVE        ::   daHost  ! Date Host Model                 
85!$OMP THREADPRIVATE(daHost)
86                                                                               
87                                                                               
88! +--SISVAT INPUT/OUTPUT Variables                                             
89! +  -----------------------------                                             
90                                                                               
91      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   isnoSV  ! Nb of Ice/Snow Layers           
92!$OMP THREADPRIVATE(isnoSV)
93      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   ispiSV  ! Uppermost superimposed ice     
94!$OMP THREADPRIVATE(ispiSV)
95      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   iiceSV  ! Nb of Ice      Layers           
96!$OMP THREADPRIVATE(iiceSV)
97      INTEGER ,ALLOCATABLE,SAVE              ::   istoSV(:,:)  ! Snow Layer     History         
98!$OMP THREADPRIVATE(istoSV)                                                                             
99      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   alb_SV  ! Surface Albedo           
100!$OMP THREADPRIVATE(alb_SV)
101      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   emi_SV  ! Surface Emissivity       
102!$OMP THREADPRIVATE(emi_SV)
103      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   IRs_SV  ! Soil           IR Flux         
104!$OMP THREADPRIVATE(IRs_SV)
105      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   LMO_SV  ! Monin-Obukhov  Scale           
106!$OMP THREADPRIVATE(LMO_SV)
107      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   us__SV  ! Friction       Velocity         
108!$OMP THREADPRIVATE(us__SV)
109      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   uts_SV  ! Temperature Turbulent Scale     
110!$OMP THREADPRIVATE(uts_SV)
111      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   cutsSV  ! Temperature Turbulent Scale C. 
112!$OMP THREADPRIVATE(cutsSV)
113      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   uqs_SV  ! Spec.Humid. Turbulent Scale     
114!$OMP THREADPRIVATE(uqs_SV)
115      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   uss_SV  ! Blow.Snow   Turbulent Scale     
116!$OMP THREADPRIVATE(uss_SV)
117      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   usthSV  ! Blowing Snow Erosion Thresh.   
118!$OMP THREADPRIVATE(usthSV)
119      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rCDmSV  ! Square  Root Contribut. Drag_m 
120!$OMP THREADPRIVATE(rCDmSV)
121      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rCDhSV  ! Square  Root Contribut. Drag_h 
122!$OMP THREADPRIVATE(rCDhSV)
123      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0m_SV  ! Momentum     Roughness Length   
124!$OMP THREADPRIVATE(Z0m_SV)
125      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0mmSV  !  z0(Momentum,    Time Mean) [m]
126!$OMP THREADPRIVATE(Z0mmSV)
127      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0mnSV  !  z0(Momentum,    instanta.) [m]
128!$OMP THREADPRIVATE(Z0mnSV)
129      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0roSV  ! Subgrid Topo Roughness Length   
130!$OMP THREADPRIVATE(Z0roSV)
131      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0SaSV  !  z0(Sastrugi  h)            [m]
132!$OMP THREADPRIVATE(Z0SaSV)
133      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0e_SV  !  z0(Snow eroded)            [m]
134!$OMP THREADPRIVATE(Z0e_SV)
135      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0emSV  !  z0(Snow eroded, Time Mean) [m]
136!$OMP THREADPRIVATE(Z0emSV)
137      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0enSV  !  z0(Snow eroded, instanta.) [m]
138!$OMP THREADPRIVATE(Z0enSV)
139      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0h_SV  ! Heat         Roughness Length   
140!$OMP THREADPRIVATE(Z0h_SV)
141      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0hmSV  !  z0(Heat,        Time Mean) [m]
142!$OMP THREADPRIVATE(Z0hmSV)
143      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0hnSV  !  z0(Heat,        instanta.) [m]
144!$OMP THREADPRIVATE(Z0hnSV)
145                                                                         
146
147                                                                       
148      REAL ,ALLOCATABLE,SAVE   ::   TsisSV(:,:)  ! Snow/Ice/Soil-Water Temperature
149!$OMP THREADPRIVATE(TsisSV)
150      REAL ,ALLOCATABLE,SAVE   ::   ro__SV(:,:)  ! Snow/Ice/Soil-Water VolumicMass
151!$OMP THREADPRIVATE(ro__SV)
152      REAL,ALLOCATABLE,SAVE    ::   eta_SV(:,:)  ! Snow/Ice/Soil     Water Content
153!$OMP THREADPRIVATE(eta_SV)
154      REAL,ALLOCATABLE,SAVE    ::   G1snSV(:,:)  ! Snow Dendricity/Sphericity     
155!$OMP THREADPRIVATE(G1snSV)
156      REAL,ALLOCATABLE,SAVE    ::   G2snSV(:,:)  ! Snow Sphericity/Size           
157!$OMP THREADPRIVATE(G2snSV)
158      REAL,ALLOCATABLE,SAVE    ::   dzsnSV(:,:)  ! Snow Layer  Thickness           
159!$OMP THREADPRIVATE(dzsnSV)
160      REAL,ALLOCATABLE,SAVE    ::   agsnSV(:,:)  ! Snow Age                       
161!$OMP THREADPRIVATE(agsnSV)
162      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BufsSV  ! Snow Buffer Layer               
163!$OMP THREADPRIVATE(BufsSV)
164      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rusnSV  ! Surficial   Water               
165!$OMP THREADPRIVATE(rusnSV)
166      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   SWf_SV  ! Normalized  Decay               
167!$OMP THREADPRIVATE(SWf_SV)
168      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   SWS_SV  ! Surficial Water Status         
169!$OMP THREADPRIVATE(SWS_SV)
170      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   HFraSV  ! Frazil      Thickness           
171!$OMP THREADPRIVATE(HFraSV)
172                                                                           
173      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   zWE_SV  ! Current   Snow Thickness [mmWE]
174!$OMP THREADPRIVATE(zWE_SV)
175      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   zWEcSV  ! Compacted Snow Thickness [mmWE]
176!$OMP THREADPRIVATE(zWEcSV)
177      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   wem_SV  ! Only Melting             [mmWE]
178!$OMP THREADPRIVATE(wem_SV)
179      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   wer_SV  ! Refreezing               [mmWE]
180!$OMP THREADPRIVATE(wer_SV)
181      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   wes_SV  ! Sublimation              [mmWE]
182!$OMP THREADPRIVATE(wes_SV)
183      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   zn4_SV  ! snowheight change         [mm]
184!$OMP THREADPRIVATE(zn4_SV)
185      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   zn5_SV  ! snowheight change         [mm]
186!$OMP THREADPRIVATE(zn5_SV)                                                                           
187                                                                             
188! +--SISVAT OUTPUT       Variables                                           
189! +  -----------------------------                                           
190                                                                               
191      INTEGER,DIMENSION(nb_wri),SAVE ::   no__SV  ! OUTPUT file Unit Number         
192!$OMP THREADPRIVATE(no__SV)
193      INTEGER,DIMENSION(nb_wri),SAVE ::   i___SV  ! OUTPUT point   i Coordinate   
194!$OMP THREADPRIVATE(i___SV)
195      INTEGER,DIMENSION(nb_wri),SAVE ::   j___SV  ! OUTPUT point   j Coordinate     
196!$OMP THREADPRIVATE(j___SV)
197      INTEGER,DIMENSION(nb_wri),SAVE ::   n___SV  ! OUTPUT point   n Coordinate     
198!$OMP THREADPRIVATE(n___SV)
199      INTEGER,DIMENSION(nb_wri),SAVE ::   lwriSV  ! OUTPUT point vec Index         
200!$OMP THREADPRIVATE(lwriSV)
201!                                                                             
202      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   ii__SV  ! WORK   point   i Coordinate     
203!$OMP THREADPRIVATE(ii__SV)
204      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   jj__SV  ! WORK   point   j Coordinate     
205!$OMP THREADPRIVATE(jj__SV)
206      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   nn__SV  ! WORK   point   n Coordinate     
207!$OMP THREADPRIVATE(nn__SV)
208                                                                               
209      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   IRu_SV  ! UPward    IR Flux (effective)   
210!$OMP THREADPRIVATE(IRu_SV)
211      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   hSalSV  ! Saltating Layer Height         
212!$OMP THREADPRIVATE(hSalSV)
213      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   qSalSV  ! Saltating Snow  Concentration   
214!$OMP THREADPRIVATE(qSalSV)
215      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   RnofSV  ! RunOFF    Intensity             
216!$OMP THREADPRIVATE(RnofSV)
217      REAL, DIMENSION(:,:),ALLOCATABLE,SAVE    ::   RuofSV  ! RunOFF    Intensity             
218!$OMP THREADPRIVATE(RuofSV)
219
220
221
222                 
223CONTAINS
224
225
226
227  SUBROUTINE INIT_VARxSV
228  IMPLICIT NONE
229
230      INTEGER                                :: ikl
231
232      ALLOCATE(LSmask(klonv))  ! Land-Sea   Mask                 
233      ALLOCATE(isotSV(klonv))  ! Soil       Type                 
234      ALLOCATE(iWaFSV(klonv))  ! Soil       Drainage:(1,0)=(y,n)                                                                               
235      ALLOCATE(coszSV(klonv))  ! Cosine of Sun zenithal Angle   
236      ALLOCATE(sol_SV(klonv))  ! Downward  Solar    Radiation   
237      ALLOCATE(IRd_SV(klonv))  ! Downward  Longwave Radiation   
238                                                                               
239      ALLOCATE(drr_SV(klonv))  ! Rain  Intensity       [kg/m2/s]
240      ALLOCATE(dsn_SV(klonv))  ! Snow  Intensity       [kg/m2/s]
241      ALLOCATE(dsnbSV(klonv))  ! Idem, from Drift        [kg/m2]
242      ALLOCATE(esnbSV(klonv))  ! Idem, from Drift        [kg/m2]
243      ALLOCATE(dbs_SV(klonv))  ! Drift Amount            [kg/m2]
244      ALLOCATE(dbs_er(klonv))  ! Drift erosion           [kg/m2]
245      ALLOCATE(dbs_ac(klonv))  ! Drift accumulation      [kg/m2]
246      ALLOCATE(BrosSV(klonv))
247      ALLOCATE(BG1sSV(klonv))         
248      ALLOCATE(BG2sSV(klonv))
249      ALLOCATE(dz0_SV(klonv))  ! dz0(Sastrugi dh)            [m]
250                                                                               
251      ALLOCATE(cld_SV(klonv))  ! Cloudiness (seen from SBL)     
252      ALLOCATE(za__SV(klonv))  ! SBL Height                     
253      ALLOCATE(VV__SV(klonv))  !(SBL Top)  Wind Velocity         
254      ALLOCATE(VVs_SV(klonv))
255      ALLOCATE(RRs_SV(klonv))
256      ALLOCATE(DDs_SV(klonv))
257      ALLOCATE(TaT_SV(klonv))  ! SBL Top   Temperature           
258      ALLOCATE(ExnrSV(klonv))  ! Exner     Potential             
259      ALLOCATE(dSdTSV(klonv))  ! Sensible Heat Flux T Derivat.   
260      ALLOCATE(dLdTSV(klonv))  ! Latent   Heat Flux T Derivat.   
261      ALLOCATE(rhT_SV(klonv))  ! SBL Top   Air  Density         
262      ALLOCATE(QaT_SV(klonv))  ! SBL Top   Specific Humidity     
263      ALLOCATE(dQa_SV(klonv))  ! SBL Flux  Limitation of Qa     
264      ALLOCATE(qsnoSV(klonv))  ! SBL Mean  Snow       Content   
265                                                                               
266      ALLOCATE(alb0SV(klonv))  ! Soil    Albedo                 
267      ALLOCATE(slopSV(klonv))  ! Snow/Ice/Soil-Water Surf. Slope
268
269                                                                               
270                                                                               
271! +--SISVAT INPUT/OUTPUT Variables                                             
272! +  -----------------------------                                             
273                                                                               
274      ALLOCATE(isnoSV(klonv))  ! Nb of Ice/Snow Layers           
275      ALLOCATE(ispiSV(klonv))  ! Uppermost superimposed ice     
276      ALLOCATE(iiceSV(klonv))  ! Nb of Ice      Layers           
277      ALLOCATE(istoSV(klonv,0:nsno))  ! Snow Layer     History         
278                                                                               
279      ALLOCATE(alb_SV(klonv))  ! Surface-Canopy Albedo           
280      ALLOCATE(emi_SV(klonv))  ! Surface-Canopy Emissivity       
281      ALLOCATE(IRs_SV(klonv))  ! Soil           IR Flux         
282      ALLOCATE(LMO_SV(klonv))  ! Monin-Obukhov  Scale           
283      ALLOCATE(us__SV(klonv))  ! Friction       Velocity         
284      ALLOCATE(uts_SV(klonv))  ! Temperature Turbulent Scale     
285      ALLOCATE(cutsSV(klonv))  ! Temperature Turbulent Scale C. 
286      ALLOCATE(uqs_SV(klonv))  ! Spec.Humid. Turbulent Scale     
287      ALLOCATE(uss_SV(klonv))  ! Blow.Snow   Turbulent Scale     
288      ALLOCATE(usthSV(klonv))  ! Blowing Snow Erosion Thresh.   
289      ALLOCATE(rCDmSV(klonv))  ! Square  Root Contribut. Drag_m 
290      ALLOCATE(rCDhSV(klonv))  ! Square  Root Contribut. Drag_h 
291      ALLOCATE(Z0m_SV(klonv))  ! Momentum     Roughness Length   
292      ALLOCATE(Z0mmSV(klonv))  !  z0(Momentum,    Time Mean) [m]
293      ALLOCATE(Z0mnSV(klonv))  !  z0(Momentum,    instanta.) [m]
294      ALLOCATE(Z0roSV(klonv))  ! Subgrid Topo Roughness Length   
295      ALLOCATE(Z0SaSV(klonv))  !  z0(Sastrugi  h)            [m]
296      ALLOCATE(Z0e_SV(klonv))  !  z0(Snow eroded)            [m]
297      ALLOCATE(Z0emSV(klonv))  !  z0(Snow eroded, Time Mean) [m]
298      ALLOCATE(Z0enSV(klonv))  !  z0(Snow eroded, instanta.) [m]
299      ALLOCATE(Z0h_SV(klonv))  !  Heat         Roughness Length   
300      ALLOCATE(Z0hmSV(klonv))  !  z0(Heat,        Time Mean) [m]
301      ALLOCATE(Z0hnSV(klonv))  !  z0(Heat,        instanta.) [m]
302                                                                         
303                                                                       
304      ALLOCATE(TsisSV(klonv,-nsol:nsno))  ! Snow/Ice/Soil-Water Temperature
305      ALLOCATE(ro__SV(klonv,-nsol:nsno))  ! Snow/Ice/Soil-Water VolumicMass
306      ALLOCATE(eta_SV(klonv,-nsol:nsno))  ! Snow/Ice/Soil     Water Content
307      ALLOCATE(G1snSV(klonv,    0:nsno))  ! Snow Dendricity/Sphericity     
308      ALLOCATE(G2snSV(klonv,    0:nsno))  ! Snow Sphericity/Size           
309      ALLOCATE(dzsnSV(klonv,    0:nsno))  ! Snow Layer  Thickness           
310      ALLOCATE(agsnSV(klonv,    0:nsno))  ! Snow Age                       
311      ALLOCATE(BufsSV(klonv))  ! Snow Buffer Layer               
312      ALLOCATE(rusnSV(klonv))  ! Surficial   Water               
313      ALLOCATE(SWf_SV(klonv))  ! Normalized  Decay               
314      ALLOCATE(SWS_SV(klonv))  ! Surficial Water Status         
315      ALLOCATE(HFraSV(klonv))  ! Frazil      Thickness           
316                                                                           
317      ALLOCATE(zWE_SV(klonv))  ! Current   Snow Thickness [mmWE]
318      ALLOCATE(zWEcSV(klonv))  ! Compacted Snow Thickness [mmWE]
319      ALLOCATE(wem_SV(klonv)) ! Only Melting             [mmWE]
320      ALLOCATE(wer_SV(klonv)) ! Refreezing               [mmWE]
321      ALLOCATE(wes_SV(klonv))  ! Sublimation              [mmWE]
322      ALLOCATE(zn4_SV(klonv))  ! snow height change       [mm]
323      ALLOCATE(zn5_SV(klonv))  ! snow height change       [mm]                                                                             
324                                                                             
325! +--SISVAT OUTPUT       Variables                                           
326! +  -----------------------------                                           
327                                                                               
328!                                                                             
329      ALLOCATE(ii__SV(klonv))  ! WORK   point   i Coordinate     
330      ALLOCATE(jj__SV(klonv))  ! WORK   point   j Coordinate     
331      ALLOCATE(nn__SV(klonv))  ! WORK   point   n Coordinate     
332                                                                               
333      ALLOCATE(IRu_SV(klonv))  ! UPward    IR Flux (effective)   
334      ALLOCATE(hSalSV(klonv))  ! Saltating Layer Height         
335      ALLOCATE(qSalSV(klonv))  ! Saltating Snow  Concentration   
336      ALLOCATE(RnofSV(klonv))  ! RunOFF    Intensity     
337      ALLOCATE(RuofSV(klonv,4))  ! RunOFF    Intensity     
338
339
340      DO ikl=1,klonv       
341        LSmask(ikl)   = 0
342        isotSV(ikl)   = 0               
343        iWaFSV(ikl)   = 0
344        isnoSV(ikl)   = 0 
345        ispiSV(ikl)   = 0
346        iiceSV(ikl)   = 0         
347        istoSV(ikl,:) = 0
348        ii__SV(ikl)   = 0
349        jj__SV(ikl)   = 0
350        nn__SV(ikl)   = 0
351      END DO
352  END SUBROUTINE INIT_VARxSV
353
354
355                                               
356END MODULE VARxSV
Note: See TracBrowser for help on using the repository browser.