source: LMDZ5/trunk/libf/phylmd/sisvat/VARdSV.F90 @ 1907

Last change on this file since 1907 was 1907, checked in by lguez, 10 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • 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: 7.2 KB
Line 
1MODULE VARdSV
2
3USE VAR_SV
4
5IMPLICIT NONE
6
7
8! +--SISVAT Global Variables
9! +  =======================
10
11      LOGICAL           ::   INI_SV=.false.          ! Initialisation Switch     
12      REAL,PARAMETER    ::   eps_21=1.e-21           ! Arbitrary  very small value
13
14
15! +--Snow
16! +  ----
17
18      INTEGER,PARAMETER,DIMENSION(5) ::   istdSV(1:5)=(/1,2,3,4,5/)  ! Snow History
19
20      REAL,PARAMETER    ::   Cn_dSV= 2105.           ! Snow Heat Capacity          [J/kg/K]
21      REAL,PARAMETER    ::   SMndSV= 1.00            ! Minimum Thickness of new Layers
22      REAL,PARAMETER    ::   G1_dSV= 99.             ! Conversion 0/99-->0/1
23      REAL,PARAMETER    ::   DDcdSV= 1.,DFcdSV= 4.,DScdSV= 3.   
24                                                     ! Snow Grains Optical Diameter [1e-4m]
25      REAL,PARAMETER    ::   ADSdSV= 4.              ! Snow Grains Actual  Diameter [1e-4m]
26      REAL,PARAMETER    ::   So1dSV= 0.580,So2dSV= 0.320,So3dSV= 0.100
27                                                     ! Total Solar Irradiance Fractions [-]
28                                                     ! Tuning ETH camp 0.3--0.8mim Interval
29                                                     ! Tuning ETH camp 0.8--1.5mim Interval
30                                                     ! Tuning ETH camp 1.5--2.8mim Interval
31                                                     ! So1dSV=0.606,So2dSV=0.301,So3dSV=0.093
32      REAL,PARAMETER    ::   aI1dSV= 0.40,aI2dSV= 0.45,aI3dSV= 0.65
33                                                     ! Bare Ice Albedo                  [-]
34                                                     ! Minimum/Maximum/ICE lense albedo at
35                                                     ! 800 kg/m3 and minimum pure snow albedo
36      REAL,PARAMETER    ::   ws0dSV= 0.07            ! Irreducible Water Saturation in Snow
37      REAL,PARAMETER    ::   roCdSV= 800.            ! Pore Hole Close OFF Density  [kg/m3]
38      REAL,PARAMETER    ::   ru_dSV= 200.            ! Surficial Water Scale Factor [kg/m2]
39
40!C +--Ice
41!C +  ---
42
43      REAL,PARAMETER    ::   CdidSV= 2.1
44
45!C +--Vegetation
46!C +  ----------
47
48      INTEGER,PARAMETER ::   nvgt=12
49      REAL,PARAMETER    ::   DH_dSV(0:nvgt) = (/ 0.00, 0.07, 0.21,      &
50     &       0.70, 0.07, 0.21, 0.70, 1.40, 5.60,14.00, 1.40, 5.60,14.00/)
51                                                  ! Displacement            Height   [m]
52
53      REAL,PARAMETER    ::   Z0mdSV(0:nvgt) = (/ 0.01, 0.01, 0.03,      &
54     &       0.10, 0.01, 0.03, 0.10, 0.20, 0.80, 2.00, 0.20, 0.80, 2.00/)
55                                                  ! Roughness  Length for Momentum   [m]
56
57      REAL,PARAMETER    ::   StodSV(0:nvgt) = (/5000.,  50.,  50.,      &
58     &        50.,  50.,  50.,  50.,  10.,  10.,  10.,  10.,  10.,  10./)
59                                                  ! Minimum    Stomatal Resistance [s/m]
60
61      REAL,PARAMETER    ::   PR_dSV(0:nvgt) = (/  0.0,0.5e9,0.5e9,      &
62     &      0.5e9,0.5e9,0.5e9,0.5e9,1.0e9,1.0e9,1.0e9,1.0e9,1.0e9,1.0e9/)
63                                                  ! Roots Fraction Beta Coefficient  [-]
64
65      REAL,PARAMETER    ::   rbtdSV(0:nvgt) = (/0.000,0.961,0.961,      &
66     &      0.961,0.943,0.964,0.972,0.968,0.962,0.962,0.971,0.976,0.976/)
67
68!        /0.00,  0.01,5000.,   0.0,   0.000,      !  0 NO     VEGETATION
69!         0.07,  0.01,  50.,   0.5e9, 0.961,      !  1 CROPS      LOW
70!         0.21,  0.03,  50.,   0.5e9, 0.961,      !  2 CROPS      MEDIUM
71!         0.70,  0.10,  50.,   0.5e9, 0.961,      !  3 CROPS      HIGH
72!         0.07,  0.01,  50.,   0.5e9, 0.943,      !  4 GRASS      LOW
73!         0.21,  0.03,  50.,   0.5e9, 0.964,      !  5 GRASS      MEDIUM
74!         0.70,  0.10,  50.,   0.5e9, 0.972,      !  6 GRASS      HIGH
75!         1.40,  0.20,  10.,   1.0e9, 0.968,      !  7 BROADLEAF  LOW
76!         5.60,  0.80,  10.,   1.0e9, 0.962,      !  8 BROADLEAF  MEDIUM
77!        14.00,  2.00,  10.,   1.0e9, 0.962,      !  9 BROADLEAF  HIGH
78!         1.40,  0.20,  10.,   1.0e9, 0.971,      ! 10 NEEDLELEAF LOW
79!         5.60,  0.80,  10.,   1.0e9, 0.976,      ! 11 NEEDLELEAF MEDIUM
80!        14.00,  2.00,  10.,   1.0e9, 0.976/      ! 12 NEEDLELEAF HIGH
81
82                                                  ! Internal Plant      Resistance   [s]
83      REAL,PARAMETER    ::   pscdSV = 250.        ! Critical Leaf Water Potential    [m]
84      REAL,PARAMETER    ::   StxdSV = 5000.       ! maXimum  Stomatal   Resistance [s/m]
85      REAL,PARAMETER    ::   LAIdSV = 4.          ! maximum  LAI
86
87
88!C +--Soil
89!C +  ----
90
91      REAL,PARAMETER    ::   rcwdSV = 4.180e+6    ! Density * Water Specific Heat
92      REAL              ::   dz_dSV(-nsol:0)      ! Vertical  Discretization MARSV:
93                                                  !/0.72,0.20,0.060,0.019,0.001/
94                                                  ! Layer's Thickness
95      REAL              ::   zz_dSV               ! Soil      Thickness
96
97      INTEGER,PARAMETER ::   nsot=12
98      REAL,PARAMETER    ::   etadSV(0:nsot) = (/ 1.000,0.395,0.410,     &
99     &     0.435,0.485,0.451,0.420,0.477,0.476,0.426,0.492,0.482,0.001 /)     
100                                                  ! Water Content at Saturation  [m3/m3]
101
102      REAL,PARAMETER    ::   psidSV(0:nsot) = (/ 1.000,0.121,0.090,     &
103     &     0.218,0.786,0.478,0.299,0.356,0.630,0.153,0.490,0.405,0.001 /)
104                                                  ! Water Succion at Saturation      [m]
105
106      REAL,PARAMETER    ::   Ks_dSV(0:nsot) = (/ 0.e00, 176.0e-6,       &
107     &          156.3e-6,  34.1e-6,   7.2e-6,   7.0e-6,   6.3e-6,       &
108     &            1.7e-6,   2.5e-6,   2.2e-6,   1.0e-6,   1.3e-6,0.0e0 /)
109                                                  ! Hydraulic Conductivity
110                                                  !               at Saturation    [m/s]
111      REAL,PARAMETER    ::   bCHdSV(0:nsot) = (/ 1.00, 4.05, 4.38,      &
112     &      4.90, 5.30, 5.39, 7.12, 7.75, 8.52,10.40,10.40,11.40, 0.02 /)
113                                                  ! Clapp-Hornberger Coefficient b   [-]
114
115 !     etadSV,   psidSV,   Ks_dSV    bCHdSV 
116 !     /1.000,    1.000,   0.0e00,     1.00,      !  0 WATER
117 !      0.395,    0.121, 176.0e-6,     4.05,      !  1 SAND
118 !      0.410,    0.090, 156.3e-6,     4.38,      !  2 LOAMY      SAND
119 !      0.435,    0.218,  34.1e-6,     4.90,      !  3 SANDY      LOAM
120 !      0.485,    0.786,   7.2e-6,     5.30,      !  4 SILT       LOAM
121 !      0.451,    0.478,   7.0e-6,     5.39,      !  5            LOAM
122 !      0.420,    0.299,   6.3e-6,     7.12,      !  6 SANDY CLAY LOAM
123 !      0.477,    0.356,   1.7e-6,     7.75,      !  7 SILTY CLAY LOAM
124 !      0.476,    0.630,   2.5e-6,     8.52,      !  8       CLAY LOAM
125 !      0.426,    0.153,   2.2e-6,    10.40,      !  9 SANDY CLAY
126 !      0.492,    0.490,   1.0e-6,    10.40,      ! 10 SILTY CLAY
127 !      0.482,    0.405,   1.3e-6,    11.40,      ! 11       CLAY
128 !      0.001,    0.001,   0.0e00,     0.02/      ! 12       ICE
129
130
131!C +--Water Bodies
132!C +  ------------
133
134      REAL,PARAMETER    ::   vK_dSV = 1000.       ! Diffusivity in Water          [m2/s]
135      REAL,PARAMETER    ::   TSIdSV = 0.50        ! Sea-Ice Fraction: SST Scale      [K]
136
137      INTEGER ::  ivg1,iso1  !rajout hjp for ini - check if ivg,iso is possible
138
139!
140
141END MODULE VARdSV
Note: See TracBrowser for help on using the repository browser.