source: LMDZ6/trunk/libf/phylmd/inlandsis/VARtSV.F90 @ 5212

Last change on this file since 5212 was 3900, checked in by evignon, 3 years ago

Commission de la nouvelle interface LMDZ-SISVAT
la prochaine commission consistera a supprimer l'ancien repertoire sisvat
et a faire un peu de nettoyage.

File size: 3.2 KB
RevLine 
[3792]1MODULE VARtSV
2
3USE VAR_SV, only : nsol, nsno, klonv
4
5IMPLICIT NONE
6! +
7      INTEGER, PARAMETER      ::   ntaver = 4
8! +
9      REAL,ALLOCATABLE,SAVE                ::   toicSV(:)    ! Snow to ice balance   
10!$OMP THREADPRIVATE(toicSV)
11
12      REAL,ALLOCATABLE,SAVE                ::   dz1_SV(:,:) ! "inverse" layer thicknes
13!$OMP THREADPRIVATE(dz1_SV)
14      REAL,ALLOCATABLE,SAVE                ::   dz2_SV(:,:) ! layer thickness
15!$OMP THREADPRIVATE(dz2_SV)
16
17
18      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   Tsf_SV  ! surface temperature !hj12032010
19!$OMP THREADPRIVATE(Tsf_SV)
20      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   TsfnSV  ! new surface temperature
21!$OMP THREADPRIVATE(TsfnSV)
22      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   AcoHSV, BcoHSV ! coefficients for Enthalpy
23!$OMP THREADPRIVATE(AcoHSV, BcoHSV)                             ! evolution,from atmosphere
24      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   AcoQSV, BcoQSV ! coefficients for Humidity
25!$OMP THREADPRIVATE( AcoQSV, BcoQSV)                            ! evolution,from atmosphere
26      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   ps__SV  ! surface pressure
27!$OMP THREADPRIVATE(ps__SV)
28      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   p1l_SV  ! 1st layer pressure
29!$OMP THREADPRIVATE(p1l_SV)
30      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   cdH_SV  ! drag coeff Energy (?)
31!$OMP THREADPRIVATE(cdH_SV)
32      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   cdM_SV  ! drag coeff momentum (?)
33!$OMP THREADPRIVATE(cdM_SV)
34      REAL, DIMENSION(:),ALLOCATABLE,SAVE  ::   rsolSV  ! Radiation balance surface
35!$OMP THREADPRIVATE(rsolSV)
36      REAL,SAVE                            ::   lambSV  ! coefficient soil layers
37!$OMP THREADPRIVATE(lambSV)
38
39 
40CONTAINS
41
42  SUBROUTINE INIT_VARtSV
[3900]43 
[3792]44  IMPLICIT NONE
[3900]45
46  INTEGER ikl
47
48
49
50
51
52
53
[3792]54      ALLOCATE(toicSV(klonv))
55
56      ALLOCATE(dz1_SV(klonv,-nsol:nsno))     ! "inverse" layer thicknes
57      ALLOCATE(dz2_SV(klonv,-nsol:nsno))     ! layer thickness
58      ALLOCATE(Tsf_SV(klonv))                ! surface temperature !hj12032010
59      ALLOCATE(TsfnSV(klonv))                ! new surface temperature
60      ALLOCATE(AcoHSV(klonv), BcoHSV(klonv)) ! coefficients for Enthalpy
61                                             ! evolution,from atmosphere
62      ALLOCATE(AcoQSV(klonv), BcoQSV(klonv)) ! coefficients for Humidity
63                                             ! evolution,from atmosphere                                                                                   
64      ALLOCATE(ps__SV(klonv))                ! surface pressure
65      ALLOCATE(p1l_SV(klonv))                ! 1st layer pressure
66      ALLOCATE(cdH_SV(klonv))                ! drag coeff Energy (?)
67      ALLOCATE(cdM_SV(klonv))                ! drag coeff momentum (?)
68      ALLOCATE(rsolSV(klonv))                ! Radiation balance surface
69
[3900]70      DO ikl=1,klonv       
71           
72         toicSV(ikl)   = 0.
73         dz1_SV(ikl,:) = 0.
74         dz2_SV(ikl,:) = 0.
75         Tsf_SV(ikl)   = 0.
76         TsfnSV(ikl)   = 0.
77         AcoHSV(ikl)   = 0.
78         BcoHSV(ikl)   = 0.
79         AcoQSV(ikl)   = 0.
80         ps__SV(ikl)   = 0.
81         p1l_SV(ikl)   = 0.
82         cdH_SV(ikl)   = 0.
83         cdM_SV(ikl)   = 0.
84         rsolSV(ikl)   = 0.
85      END DO
86
87
88
[3792]89  END SUBROUTINE INIT_VARtSV
90
91END MODULE VARtSV
92
Note: See TracBrowser for help on using the repository browser.