Context Navigation

← Previous Revision
Next Revision →
Blame
Revision Log

VARdSV.F90

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

Inclusion de la bibliothèque SISVAT/MAR à LMDZ pour le traitement des surfaces
"land ice"

Menegoz

Integration of the SISVAT/MAR library to LMDZ to model the land ice surfaces

Menegoz

File size: 7.2 KB

Line
1	MODULE VARdSV
2
3	USE VAR_SV
4
5	IMPLICIT 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
141	END MODULE VARdSV

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format