source: LMDZ5/branches/LF-private/libf/phylmd/sisvat/VARxSV.F90

Last change on this file was 1872, checked in by Laurent Fairhead, 11 years ago

Correction de bugs suite à l'intégration de SISVAT


Bug corrections following integration of SISVAT

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