source: LMDZ5/branches/testing/libf/phylmd/sisvat/VARxSV.F90 @ 5455

Last change on this file since 5455 was 1910, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1860:1909 into testing branch
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 21.2 KB
RevLine 
[1865]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                 
[1872]10!$OMP THREADPRIVATE(LSmask)
[1865]11      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   isotSV  ! Soil       Type                 
[1872]12!$OMP THREADPRIVATE(isotSV)
[1865]13      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   iWaFSV  ! Soil       Drainage:(1,0)=(y,n)
[1872]14!$OMP THREADPRIVATE(iWaFSV)
[1865]15      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   ivgtSV  ! Vegetation Type                 
[1872]16!$OMP THREADPRIVATE(ivgtSV)
[1865]17                                                                             
18      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   coszSV  ! Cosine of Sun zenithal Angle   
[1872]19!$OMP THREADPRIVATE(coszSV)
[1865]20      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   sol_SV  ! Downward  Solar    Radiation   
[1872]21!$OMP THREADPRIVATE(sol_SV)
[1865]22      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   IRd_SV  ! Downward  Longwave Radiation   
[1872]23!$OMP THREADPRIVATE(IRd_SV)
[1865]24                                                                               
25      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   drr_SV  ! Rain  Intensity       [kg/m2/s]
[1872]26!$OMP THREADPRIVATE(drr_SV)
[1865]27      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dsn_SV  ! Snow  Intensity       [kg/m2/s]
[1872]28!$OMP THREADPRIVATE(dsn_SV)
[1865]29      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dsnbSV  ! Idem, fraction, from Drift  [-]
[1872]30!$OMP THREADPRIVATE(dsnbSV)
[1865]31      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   esnbSV  ! Idem, fraction, from Drift  [-]
[1872]32!$OMP THREADPRIVATE(esnbSV)
[1865]33      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dbs_SV  ! Drift Amount            [kg/m2]
[1872]34!$OMP THREADPRIVATE(dbs_SV)
[1865]35      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BrosSV  ! Buffer Snow Layer Density
[1872]36!$OMP THREADPRIVATE(BrosSV)
[1865]37      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BG1sSV  ! Buffer Snow Layer Dendr/Sphe[-]
[1872]38!$OMP THREADPRIVATE(BG1sSV)
[1865]39      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BG2sSV  ! Buffer Snow Layer Spher/Size[-][0.0001m]
[1872]40!$OMP THREADPRIVATE(BG2sSV)
[1865]41      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dz0_SV  ! dz0(Sastrugi dh)            [m]
[1872]42!$OMP THREADPRIVATE(dz0_SV)
[1865]43                                                                               
44      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   cld_SV  ! Cloudiness (seen from SBL)     
[1872]45!$OMP THREADPRIVATE(cld_SV)
[1865]46      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   za__SV  ! SBL Height                     
[1872]47!$OMP THREADPRIVATE(za__SV)
[1865]48      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   VV__SV  !(SBL Top)  Wind Velocity         
[1872]49!$OMP THREADPRIVATE(VV__SV)
[1865]50      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   VVs_SV  !(Sastr,V)  Relevance
[1872]51!$OMP THREADPRIVATE(VVs_SV)
[1865]52      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   RRs_SV  !(Sastr,V)  Counter
[1872]53!$OMP THREADPRIVATE(RRs_SV)
[1865]54      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   DDs_SV  !(Sastr,V)  Angle
[1872]55!$OMP THREADPRIVATE(DDs_SV)
[1865]56      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   TaT_SV  ! SBL Top   Temperature           
[1872]57!$OMP THREADPRIVATE(TaT_SV)
[1865]58      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   ExnrSV  ! Exner     Potential             
[1872]59!$OMP THREADPRIVATE(ExnrSV)
[1865]60      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dSdTSV  ! Sensible Heat Flux T Derivat.   
[1872]61!$OMP THREADPRIVATE(dSdTSV)
[1865]62      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dLdTSV  ! Latent   Heat Flux T Derivat.   
[1872]63!$OMP THREADPRIVATE(dLdTSV)
[1865]64      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rhT_SV  ! SBL Top   Air  Density         
[1872]65!$OMP THREADPRIVATE(rhT_SV)
[1865]66      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   QaT_SV  ! SBL Top   Specific Humidity     
[1872]67!$OMP THREADPRIVATE(QaT_SV)
[1865]68      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   dQa_SV  ! SBL Flux  Limitation of Qa     
[1872]69!$OMP THREADPRIVATE(dQa_SV)
[1865]70      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   qsnoSV  ! SBL Mean  Snow       Content   
[1872]71!$OMP THREADPRIVATE(qsnoSV)
[1865]72                                                                           
73      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   LAI0SV  ! Nominal Leaf Area Index         
[1872]74!$OMP THREADPRIVATE(LAI0SV)
[1865]75      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   glf0SV  ! Green   Leaf Fraction           
[1872]76!$OMP THREADPRIVATE(glf0SV)
[1865]77                                                                               
78      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   alb0SV  ! Soil    Albedo                 
[1872]79!$OMP THREADPRIVATE(alb0SV)
[1865]80      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   slopSV  ! Snow/Ice/Soil-Water Surf. Slope
[1872]81!$OMP THREADPRIVATE(slopSV)
[1865]82
83                                                             
84      REAL,SAVE                      ::   zSBLSV  ! SBL Height (Initial Value)     
[1872]85!$OMP THREADPRIVATE(zSBLSV)
[1865]86      REAL,SAVE                      ::   dt__SV  ! Time Step                       
[1872]87!$OMP THREADPRIVATE(dt__SV)
[1865]88      CHARACTER (len=18),SAVE        ::   daHost  ! Date Host Model                 
[1872]89!$OMP THREADPRIVATE(daHost)
[1865]90                                                                               
91                                                                               
92! +--SISVAT INPUT/OUTPUT Variables                                             
93! +  -----------------------------                                             
94                                                                               
95      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   isnoSV  ! Nb of Ice/Snow Layers           
[1872]96!$OMP THREADPRIVATE(isnoSV)
[1865]97      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   ispiSV  ! Uppermost superimposed ice     
[1872]98!$OMP THREADPRIVATE(ispiSV)
[1865]99      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   iiceSV  ! Nb of Ice      Layers           
[1872]100!$OMP THREADPRIVATE(iiceSV)
[1865]101      INTEGER ,ALLOCATABLE,SAVE    ::   istoSV(:,:)  ! Snow Layer     History         
[1872]102!$OMP THREADPRIVATE(istoSV)
[1865]103                                                                               
104      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   alb_SV  ! Surface-Canopy Albedo           
[1872]105!$OMP THREADPRIVATE(alb_SV)
[1865]106      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   emi_SV  ! Surface-Canopy Emissivity       
[1872]107!$OMP THREADPRIVATE(emi_SV)
[1865]108      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   IRs_SV  ! Soil           IR Flux         
[1872]109!$OMP THREADPRIVATE(IRs_SV)
[1865]110      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   LMO_SV  ! Monin-Obukhov  Scale           
[1872]111!$OMP THREADPRIVATE(LMO_SV)
[1865]112      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   us__SV  ! Friction       Velocity         
[1872]113!$OMP THREADPRIVATE(us__SV)
[1865]114      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   uts_SV  ! Temperature Turbulent Scale     
[1872]115!$OMP THREADPRIVATE(uts_SV)
[1865]116      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   cutsSV  ! Temperature Turbulent Scale C. 
[1872]117!$OMP THREADPRIVATE(cutsSV)
[1865]118      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   uqs_SV  ! Spec.Humid. Turbulent Scale     
[1872]119!$OMP THREADPRIVATE(uqs_SV)
[1865]120      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   uss_SV  ! Blow.Snow   Turbulent Scale     
[1872]121!$OMP THREADPRIVATE(uss_SV)
[1865]122      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   usthSV  ! Blowing Snow Erosion Thresh.   
[1872]123!$OMP THREADPRIVATE(usthSV)
[1865]124      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rCDmSV  ! Square  Root Contribut. Drag_m 
[1872]125!$OMP THREADPRIVATE(rCDmSV)
[1865]126      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rCDhSV  ! Square  Root Contribut. Drag_h 
[1872]127!$OMP THREADPRIVATE(rCDhSV)
[1865]128      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0m_SV  ! Momentum     Roughness Length   
[1872]129!$OMP THREADPRIVATE(Z0m_SV)
[1865]130      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0mmSV  !  z0(Momentum,    Time Mean) [m]
[1872]131!$OMP THREADPRIVATE(Z0mmSV)
[1865]132      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0mnSV  !  z0(Momentum,    instanta.) [m]
[1872]133!$OMP THREADPRIVATE(Z0mnSV)
[1865]134      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0roSV  ! Subgrid Topo Roughness Length   
[1872]135!$OMP THREADPRIVATE(Z0roSV)
[1865]136      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0SaSV  !  z0(Sastrugi  h)            [m]
[1872]137!$OMP THREADPRIVATE(Z0SaSV)
[1865]138      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0e_SV  !  z0(Snow eroded)            [m]
[1872]139!$OMP THREADPRIVATE(Z0e_SV)
[1865]140      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0emSV  !  z0(Snow eroded, Time Mean) [m]
[1872]141!$OMP THREADPRIVATE(Z0emSV)
[1865]142      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0enSV  !  z0(Snow eroded, instanta.) [m]
[1872]143!$OMP THREADPRIVATE(Z0enSV)
[1865]144      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0h_SV  ! Heat         Roughness Length   
[1872]145!$OMP THREADPRIVATE(Z0h_SV)
[1865]146      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0hmSV  !  z0(Heat,        Time Mean) [m]
[1872]147!$OMP THREADPRIVATE(Z0hmSV)
[1865]148      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   Z0hnSV  !  z0(Heat,        instanta.) [m]
[1872]149!$OMP THREADPRIVATE(Z0hnSV)
[1865]150                                                                         
151      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   snCaSV  ! Canopy  Snow   Thickness       
[1872]152!$OMP THREADPRIVATE(snCaSV)
[1865]153      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rrCaSV  ! Canopy  Water  Content         
[1872]154!$OMP THREADPRIVATE(rrCaSV)
[1865]155      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   psivSV  ! Leaf    Water  Potential       
[1872]156!$OMP THREADPRIVATE(psivSV)
[1865]157      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   TvegSV  ! Vegetation     Temperature     
[1872]158!$OMP THREADPRIVATE(TvegSV)
[1865]159                                                                       
160      REAL ,ALLOCATABLE,SAVE   ::   TsisSV(:,:)  ! Snow/Ice/Soil-Water Temperature
[1872]161!$OMP THREADPRIVATE(TsisSV)
[1865]162      REAL ,ALLOCATABLE,SAVE   ::   ro__SV(:,:)  ! Snow/Ice/Soil-Water VolumicMass
[1872]163!$OMP THREADPRIVATE(ro__SV)
[1865]164      REAL,ALLOCATABLE,SAVE    ::   eta_SV(:,:)  ! Snow/Ice/Soil     Water Content
[1872]165!$OMP THREADPRIVATE(eta_SV)
[1865]166      REAL,ALLOCATABLE,SAVE    ::   G1snSV(:,:)  ! Snow Dendricity/Sphericity     
[1872]167!$OMP THREADPRIVATE(G1snSV)
[1865]168      REAL,ALLOCATABLE,SAVE    ::   G2snSV(:,:)  ! Snow Sphericity/Size           
[1872]169!$OMP THREADPRIVATE(G2snSV)
[1865]170      REAL,ALLOCATABLE,SAVE    ::   dzsnSV(:,:)  ! Snow Layer  Thickness           
[1872]171!$OMP THREADPRIVATE(dzsnSV)
[1865]172      REAL,ALLOCATABLE,SAVE    ::   agsnSV(:,:)  ! Snow Age                       
[1872]173!$OMP THREADPRIVATE(agsnSV)
[1865]174      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   BufsSV  ! Snow Buffer Layer               
[1872]175!$OMP THREADPRIVATE(BufsSV)
[1865]176      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   rusnSV  ! Surficial   Water               
[1872]177!$OMP THREADPRIVATE(rusnSV)
[1865]178      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   SWf_SV  ! Normalized  Decay               
[1872]179!$OMP THREADPRIVATE(SWf_SV)
[1865]180      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   SWS_SV  ! Surficial Water Status         
[1872]181!$OMP THREADPRIVATE(SWS_SV)
[1865]182      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   HFraSV  ! Frazil      Thickness           
[1872]183!$OMP THREADPRIVATE(HFraSV)
[1865]184                                                                           
185      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   zWE_SV  ! Current   Snow Thickness [mmWE]
[1872]186!$OMP THREADPRIVATE(zWE_SV)
[1865]187      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   zWEcSV  ! Compacted Snow Thickness [mmWE]
[1872]188!$OMP THREADPRIVATE(zWEcSV)
[1865]189      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   wem_SV  ! Only Melting             [mmWE]
[1872]190!$OMP THREADPRIVATE(wem_SV)
[1865]191      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   wer_SV  ! Refreezing               [mmWE]
[1872]192!$OMP THREADPRIVATE(wer_SV)
[1865]193      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   wes_SV  ! Sublimation              [mmWE]
[1872]194!$OMP THREADPRIVATE(wes_SV)
[1865]195                                                                             
196                                                                             
197! +--SISVAT OUTPUT       Variables                                           
198! +  -----------------------------                                           
199                                                                               
200      INTEGER,DIMENSION(nb_wri),SAVE ::   no__SV  ! OUTPUT file Unit Number         
[1872]201!$OMP THREADPRIVATE(no__SV)
[1865]202      INTEGER,DIMENSION(nb_wri),SAVE ::   i___SV  ! OUTPUT point   i Coordinate   
[1872]203!$OMP THREADPRIVATE(i___SV)
[1865]204      INTEGER,DIMENSION(nb_wri),SAVE ::   j___SV  ! OUTPUT point   j Coordinate     
[1872]205!$OMP THREADPRIVATE(j___SV)
[1865]206      INTEGER,DIMENSION(nb_wri),SAVE ::   n___SV  ! OUTPUT point   n Coordinate     
[1872]207!$OMP THREADPRIVATE(n___SV)
[1865]208      INTEGER,DIMENSION(nb_wri),SAVE ::   lwriSV  ! OUTPUT point vec Index         
[1872]209!$OMP THREADPRIVATE(lwriSV)
[1865]210!                                                                             
211      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   ii__SV  ! WORK   point   i Coordinate     
[1872]212!$OMP THREADPRIVATE(ii__SV)
[1865]213      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   jj__SV  ! WORK   point   j Coordinate     
[1872]214!$OMP THREADPRIVATE(jj__SV)
[1865]215      INTEGER, DIMENSION(:),ALLOCATABLE,SAVE ::   nn__SV  ! WORK   point   n Coordinate     
[1872]216!$OMP THREADPRIVATE(nn__SV)
[1865]217                                                                               
218      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   IRu_SV  ! UPward    IR Flux (effective)   
[1872]219!$OMP THREADPRIVATE(IRu_SV)
[1865]220      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   hSalSV  ! Saltating Layer Height         
[1872]221!$OMP THREADPRIVATE(hSalSV)
[1865]222      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   qSalSV  ! Saltating Snow  Concentration   
[1872]223!$OMP THREADPRIVATE(qSalSV)
[1865]224      REAL, DIMENSION(:),ALLOCATABLE,SAVE    ::   RnofSV  ! RunOFF    Intensity             
[1872]225!$OMP THREADPRIVATE(RnofSV)
[1865]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.