source: LMDZ6/branches/LMDZ-tracers/libf/phylmd/inlandsis/VARdSV.F90 @ 3999

Last change on this file since 3999 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: 4.9 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    ::   roBdSV= 450.            ! Max blowing snow Density [kg/m3]
39      REAL,PARAMETER    ::   roSdSV= 450.            ! Max pure    snow Density [kg/m3]
40      REAL,PARAMETER    ::   ru_dSV= 200.            ! Surficial Water Scale Factor [kg/m2]
41
42!C +--Ice
43!C +  ---
44
45      REAL,PARAMETER    ::   CdidSV= 2.1
46
47
48
49!C +--Soil
50!C +  ----
51
52      REAL,PARAMETER    ::   rcwdSV = 4.180e+6    ! Density * Water Specific Heat
53      REAL              ::   dz_dSV(-nsol:0)      ! Vertical  Discretization MARSV:
54                                                  !/0.72,0.20,0.060,0.019,0.001/
55                                                  ! Layer's Thickness
56      REAL              ::   zz_dSV               ! Soil      Thickness
57
58      REAL,PARAMETER    ::   etadSV(0:nsot) = (/ 1.000,0.395,0.410,     &
59     &     0.435,0.485,0.451,0.420,0.477,0.476,0.426,0.492,0.482,0.001 /)     
60                                                  ! Water Content at Saturation  [m3/m3]
61
62      REAL,PARAMETER    ::   psidSV(0:nsot) = (/ 1.000,0.121,0.090,     &
63     &     0.218,0.786,0.478,0.299,0.356,0.630,0.153,0.490,0.405,0.001 /)
64                                                  ! Water Succion at Saturation      [m]
65
66      REAL,PARAMETER    ::   Ks_dSV(0:nsot) = (/ 0.e00, 176.0e-6,       &
67     &          156.3e-6,  34.1e-6,   7.2e-6,   7.0e-6,   6.3e-6,       &
68     &            1.7e-6,   2.5e-6,   2.2e-6,   1.0e-6,   1.3e-6,0.0e0 /)
69                                                  ! Hydraulic Conductivity
70                                                  !               at Saturation    [m/s]
71      REAL,PARAMETER    ::   bCHdSV(0:nsot) = (/ 1.00, 4.05, 4.38,      &
72     &      4.90, 5.30, 5.39, 7.12, 7.75, 8.52,10.40,10.40,11.40, 0.02 /)
73                                                  ! Clapp-Hornberger Coefficient b   [-]
74
75 !     etadSV,   psidSV,   Ks_dSV    bCHdSV 
76 !     /1.000,    1.000,   0.0e00,     1.00,      !  0 WATER
77 !      0.395,    0.121, 176.0e-6,     4.05,      !  1 SAND
78 !      0.410,    0.090, 156.3e-6,     4.38,      !  2 LOAMY      SAND
79 !      0.435,    0.218,  34.1e-6,     4.90,      !  3 SANDY      LOAM
80 !      0.485,    0.786,   7.2e-6,     5.30,      !  4 SILT       LOAM
81 !      0.451,    0.478,   7.0e-6,     5.39,      !  5            LOAM
82 !      0.420,    0.299,   6.3e-6,     7.12,      !  6 SANDY CLAY LOAM
83 !      0.477,    0.356,   1.7e-6,     7.75,      !  7 SILTY CLAY LOAM
84 !      0.476,    0.630,   2.5e-6,     8.52,      !  8       CLAY LOAM
85 !      0.426,    0.153,   2.2e-6,    10.40,      !  9 SANDY CLAY
86 !      0.492,    0.490,   1.0e-6,    10.40,      ! 10 SILTY CLAY
87 !      0.482,    0.405,   1.3e-6,    11.40,      ! 11       CLAY
88 !      0.001,    0.001,   0.0e00,     0.02/      ! 12       ICE
89
90
91!C +--Water Bodies
92!C +  ------------
93
94      REAL,PARAMETER    ::   vK_dSV = 1000.       ! Diffusivity in Water          [m2/s]
95      REAL,PARAMETER    ::   TSIdSV = 0.50        ! Sea-Ice Fraction: SST Scale      [K]
96
97      INTEGER ::  iso1  !rajout hjp for ini - check if ivg,iso is possible
98
99!
100
101END MODULE VARdSV
Note: See TracBrowser for help on using the repository browser.