Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

Mod_SISVAT_dat.f90 @ 2188

Last change on this file since 2188 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics

File size: 14.9 KB

Line
1	module Mod_SISVAT_dat
2
3
4	!--------------------------------------------------------------------------+
5	! Sun 30-Jun-2013 MAR \|
6	! module Mod_SISVAT_dat contains the constants of the \|
7	! Soil/Ice Snow Vegetation Atmosphere Transfer Scheme \|
8	! \|
9	! \|
10	! CAUTION: Soil Hydraulic Parameters: Please VERIFY ICE DATA \|
11	! ^^^^^^^ \|
12	! \|
13	! \|
14	! Preprocessing Option: \|
15	! ^^^^^^^^^^^^^^^^^^^^^ \|
16	! #fb: Fractions of total Solar Irradiance (Feagle and Businger 1981) \|
17	! \|
18	! \|
19	! version 3.p.4.1 created by H. Gallee, Wed 30-Jan-2013 \|
20	! modified by H. Gallee, 30-Jun-May-2013 \|
21	! \|
22	!--------------------------------------------------------------------------+
23
24
25	! Global Variables
26	! ==================
27
28	use Mod_Real
29	use Mod_SISVAT_dim
30
31
32	IMPLICIT NONE
33
34
35
36	! Internal Variables
37	! ==================
38
39	integer, SAVE :: ivg,iso
40
41
42
43	! SISVAT Global Variables
44	! =======================
45
46	real(kind=real8), SAVE :: eps_21 = 1.e-21 ! Arbitrary very small value
47
48	real(kind=real8), SAVE :: slop1d = 0.00 ! Local Slope (1-d simulations)
49	real(kind=real8), SAVE :: Adz0dt ! Decay of Angle(Wind,Sastrugi) Influence on z0
50
51	real(kind=real8), SAVE :: WatIsv = 273.16 ! used in Evaluation of Saturation Specific Humidity over Ice (see Dudhia)
52	real(kind=real8), SAVE :: ExpIsv = 6150.00 ! used in Evaluation of Saturation Specific Humidity over Ice (see Dudhia)
53
54
55
56	! Snow
57	! ----
58
59	integer, SAVE, dimension(5) :: istdSV ! Snow History
60	data (istdSV(iso),iso=1,5) /1,2,3,4,5/
61	! 1: faceted cristal
62	! 2: liq.watr/no faceted cristal befor
63	! 3: liq.watr/ faceted cristal befor
64
65	real(kind=real8), SAVE :: Cn_dSV = 2105. ! Snow Heat Capacity [J/kg/K]
66	! Loth et al. 1993, JGR 98 D6
67	real(kind=real8), SAVE :: SMndSV = 1.00 ! New Snow Layer Min.Thickn. [mm w.e.]
68	real(kind=real8), SAVE :: G1_dSV = 99.00 ! Conversion 0/99-->0/1
69	! Sphericity/Dendricity
70
71	! Optical Diameter of:
72	real(kind=real8), SAVE :: DDcdSV = 1.00 ! Dendritic Crystals [0.0001 m]
73	real(kind=real8), SAVE :: DFcdSV = 4.00 ! Young Faceted Crystals [0.0001 m]
74	real(kind=real8), SAVE :: DScdSV = 3.00 ! Small Crystals [0.0001 m]
75
76	! Actual Diameter of:
77	real(kind=real8), SAVE :: ADSdSV = 4.00 ! Small Crystals [0.0001 m]
78
79
80	! Fractions of total Solar Irradiance in 3 spectral Intervals
81	! (see Feagle and Businger 1981, Int.Geoph.Ser. 25, p.215-222)
82	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
83	! #fb real(kind=real8), SAVE :: So1dSV = 0.606 ! 0.3--0.8mim Interval
84	! #fb real(kind=real8), SAVE :: So2dSV = 0.301 ! 0.8--1.5mim Interval
85	! #fb real(kind=real8), SAVE :: So3dSV = 0.093 ! 1.5--2.8mim Interval
86
87	real(kind=real8), SAVE :: So1dSV = 0.580 ! Tuning ETH camp 0.3--0.8mim Interval
88	real(kind=real8), SAVE :: So2dSV = 0.320 ! Tuning ETH camp 0.8--1.5mim Interval
89	real(kind=real8), SAVE :: So3dSV = 0.100 ! Tuning ETH camp 1.5--2.8mim Interval
90
91	real(kind=real8), SAVE :: aI1dSV = 0.40 ! Minimum bare ICE albedo [-]
92	real(kind=real8), SAVE :: aI2dSV = 0.45 ! Maximum bare ICE albedo [-]
93	real(kind=real8), SAVE :: aI3dSV = 0.65 ! ICE lense albedo at 800 kg/m3 [-]
94	! and minimum pure snow albedo
95
96	! Water in Snow
97	! -------------
98	real(kind=real8), SAVE :: ws0dSV = 0.07 ! Irreducible Water Saturation in Snow
99	! Coleou et al., 1998, A.Gla.26, 64-68
100	real(kind=real8), SAVE :: roCdSV = 800.00 ! Pore Hole Close OFF Density [kg/m3]
101	! Greuell & Konzelmann (1994)
102	! Glob.Plan.Change 9, 4.5 p.100
103	real(kind=real8), SAVE :: ru_dSV = 200.00 ! Surficial Water Scale Factor [kg/m2]
104
105
106	! Ice
107	! ---
108
109	real(kind=real8), SAVE :: CdidSV = 2.10 ! Conductivity of pure Ice [W/m/K]
110
111
112
113	! Vegetation (SVAT Classification)
114	! -----------------------------------------------
115
116	real(kind=real8), SAVE :: DH_dSV(0:nvgt) ! Displacement Height [m]
117	real(kind=real8), SAVE :: Z0mdSV(0:nvgt) ! Roughness Length for Momentum [m]
118	real(kind=real8), SAVE :: StodSV(0:nvgt) ! Minimum Stomatal Resistance [s/m]
119	real(kind=real8), SAVE :: PR_dSV(0:nvgt) ! Internal Plant Resistance [s]
120	real(kind=real8), SAVE :: rbtdSV(0:nvgt) ! Root Fraction Beta Coefficient [-]
121
122	data (DH_dSV(ivg), &! Displacement Height [m]
123	& Z0mdSV(ivg), &! Roughness Length for Momentum [m]
124	& StodSV(ivg), &! Minimum Stomatal Resistance [s/m]
125	& PR_dSV(ivg), &! Internal Plant Resistance [s]
126	& rbtdSV(ivg), &! Root beta coeffient [-]
127	& ivg=0,nvgt) &
128	& /0.00, 0.01,5000., 0.0, 0.000, &! 0 NO VEGETATION
129	& 0.07, 0.01, 50., 0.5e9, 0.961, &! 1 CROPS LOW
130	& 0.21, 0.03, 50., 0.5e9, 0.961, &! 2 CROPS MEDIUM
131	& 0.70, 0.10, 50., 0.5e9, 0.961, &! 3 CROPS HIGH
132	& 0.07, 0.01, 50., 0.5e9, 0.943, &! 4 GRASS LOW
133	& 0.21, 0.03, 50., 0.5e9, 0.964, &! 5 GRASS MEDIUM
134	& 0.70, 0.10, 50., 0.5e9, 0.972, &! 6 GRASS HIGH
135	& 1.40, 0.20, 10., 1.0e9, 0.968, &! 7 BROADLEAF LOW
136	& 5.60, 0.80, 10., 1.0e9, 0.962, &! 8 BROADLEAF MEDIUM
137	& 14.00, 2.00, 10., 1.0e9, 0.962, &! 9 BROADLEAF HIGH
138	& 1.40, 0.20, 10., 1.0e9, 0.971, &! 10 NEEDLELEAF LOW
139	& 5.60, 0.80, 10., 1.0e9, 0.976, &! 11 NEEDLELEAF MEDIUM
140	& 14.00, 2.00, 10., 1.0e9, 0.976/ ! 12 NEEDLELEAF HIGH
141
142
143	real(kind=real8), SAVE :: pscdSV = 250. ! Critical Leaf Water Potential [m]
144	real(kind=real8), SAVE :: StxdSV = 5000. ! maXimum Stomatal Resistance [s/m]
145	real(kind=real8), SAVE :: LAIdSV = 4. ! maximum LAI
146
147
148	real(kind=real8), SAVE :: f__ust(0:nvgt) !
149
150	data (f__ust(ivg), ivg=0,nvgt) &!
151	& /1.00, &! 0 NO VEGETATION
152	& 1.20, &! 1 CROPS LOW
153	& 5.00, &! 2 CROPS MEDIUM
154	& 10.00, &! 3 CROPS HIGH
155	& 1.20, &! 4 GRASS LOW
156	& 5.00, &! 5 GRASS MEDIUM
157	& 10.00, &! 6 GRASS HIGH
158	& 5.00, &! 7 BROADLEAF LOW
159	& 10.00, &! 8 BROADLEAF MEDIUM
160	& 12.00, &! 9 BROADLEAF HIGH
161	& 10.00, &! 10 NEEDLELEAF LOW
162	& 12.00, &! 11 NEEDLELEAF MEDIUM
163	& 50.00 / ! 12 NEEDLELEAF HIGH
164
165
166	! Soil
167	! ----
168
169	real(kind=real8), SAVE :: rcwdSV = 4.180e+6 ! Density * Water Specific Heat
170
171	! Soil Vertical Discretization
172	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
173	real(kind=real8), SAVE :: dz_dSV(-nsol:0) ! Vertical Discretization
174	data (dz_dSV(iso),iso=-4,0) &!
175	& /0.72,0.20,0.060,0.019,0.001/ ! Layer's Thickness
176
177	real(kind=real8), SAVE :: zz_dSV ! Soil Thickness
178
179	! Soil Hydraulic Parameters (USDA Classification)
180	! ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
181	integer, parameter :: nsot = 12
182
183	real(kind=real8), SAVE :: etadSV(0:nsot) ! Water Content at Saturation [kg/kg]
184	real(kind=real8), SAVE :: psidSV(0:nsot) ! Water Succion at Saturation [m]
185	real(kind=real8), SAVE :: Ks_dSV(0:nsot) ! Hydraulic Conductivity
186	! at Saturation [m/s]
187	real(kind=real8), SAVE :: bCHdSV(0:nsot) ! Clapp-Hornberger Coefficient b [-]
188
189	data (etadSV(iso), &
190	& psidSV(iso), &
191	& Ks_dSV(iso), &
192	& bCHdSV(iso), &
193	& iso=0,nsot) &
194	& / 1.000, 1.000, 0.0e00, 1.00, &! 0 WATER
195	& 0.395, 0.121, 176.0e-6, 4.05, &! 1 SAND
196	& 0.410, 0.090, 156.3e-6, 4.38, &! 2 LOAMY SAND
197	& 0.435, 0.218, 34.1e-6, 4.90, &! 3 SANDY LOAM
198	& 0.485, 0.786, 7.2e-6, 5.30, &! 4 SILT LOAM
199	& 0.451, 0.478, 7.0e-6, 5.39, &! 5 LOAM
200	& 0.420, 0.299, 6.3e-6, 7.12, &! 6 SANDY CLAY LOAM
201	& 0.477, 0.356, 1.7e-6, 7.75, &! 7 SILTY CLAY LOAM
202	& 0.476, 0.630, 2.5e-6, 8.52, &! 8 CLAY LOAM
203	& 0.426, 0.153, 2.2e-6, 10.40, &! 9 SANDY CLAY
204	& 0.492, 0.490, 1.0e-6, 10.40, &! 10 SILTY CLAY
205	& 0.482, 0.405, 1.3e-6, 11.40, &! 11 CLAY
206	& 0.001, 0.001, 0.0e00, 0.02/ ! 12 ICE
207
208
209	real(kind=real8), SAVE :: ustdmn(0:nsot)
210	real(kind=real8), SAVE :: claypc(0:nsot)
211
212	data (ustdmn(iso), &
213	& claypc(iso), &
214	& iso=0,nsot) &
215	& /10.000, 0.0000, &! 0 WATER !
216	& 0.300, 0.0000, &! 1 SAND !
217	& 0.300, 0.0920, &! 2 LOAMY SAND ! Fal99, Table 2
218	& 0.300, 0.1420, &! 3 SANDY LOAM ! Fal99, Table 2
219	& 0.300, 0.1630, &! 4 SILT LOAM ! Guess (Interpol.)
220	& 0.300, 0.1840, &! 5 LOAM ! Fal99, Table 2
221	& 0.300, 0.2280, &! 6 SANDY CLAY LOAM ! Guess (Interpol.)
222	& 0.300, 0.2720, &! 7 SILTY CLAY LOAM ! Guess (Interpol.)
223	& 0.300, 0.3160, &! 8 CLAY LOAM ! Fal99, Table 2
224	& 0.300, 0.3750, &! 9 SANDY CLAY ! Guess (Interpol.)
225	& 0.300, 0.4340, &! 10 SILTY CLAY ! Guess (Interpol.)
226	& 0.300, 0.4920, &! 11 CLAY ! Fal99, Table 2
227	& 10.000, 0.0000/ ! 12 ICE !
228
229
230
231	! Water Bodies
232	! ------------
233
234	real(kind=real8), SAVE :: vK_dSV = 1000. ! Diffusivity in Water [m2/s]
235	real(kind=real8), SAVE :: SIcMIN = 0.1 ! Sea-Ice Layer Min Thickness [m]
236	real(kind=real8), SAVE :: dzSIce(4) ! Sea-Ice Vertical Discretisation [m]
237	data dzSIce /0.5,0.05,0.001,0.0/
238	real(kind=real8), SAVE :: SIc_OK(2) ! Sea-Ice Switch
239	data SIc_OK /1.0,0.00/
240	real(kind=real8), SAVE :: TSIdSV = 0.50 ! Sea-Ice Fraction: SST Scale [K]
241	real(kind=real8), SAVE :: OcnMin = 0.05 ! Open Water Min Fraction (S.Hem.) [-]
242	real(kind=real8), SAVE :: TocnSI = 270.70 ! Ocn Temp. for Full Sea-Ice Cover [K]
243	real(kind=real8), SAVE :: TOF_SV ! Ocn Grid Cell Freez.Temperature [K]
244	real(kind=real8), SAVE :: VarSST ! Variable (0.) / Fixed (1.)
245	real(kind=real8), SAVE :: FixSST ! Fixed (1.) / Variable (0.) SST
246	real(kind=real8), SAVE :: SSTnud ! SST Nudging Parameter
247
248
249	! Auxiliary Variables
250	! -------------------
251
252	integer, parameter :: nkhy=50
253
254	real(kind=real8), SAVE :: rocsSV( 0:nsot) ! Soil Contribution to (ro c)_s
255	real(kind=real8), SAVE :: etamSV( 0:nsot) ! Soil Minimum Humidity
256	real(kind=real8), SAVE :: s1__SV( 0:nsot) ! ... X eta**( b+2), DR97(3.36)
257	real(kind=real8), SAVE :: s2__SV( 0:nsot) ! ... X eta**(2b+3), DR97(3.35)
258	real(kind=real8), SAVE :: aKdtSV( 0:nsot, 0:nkhy) ! Khyd=aeta+b: a dt
259	real(kind=real8), SAVE :: bKdtSV( 0:nsot, 0:nkhy) ! Khyd=aeta+b: b dt
260
261
262
263	end module Mod_SISVAT_dat

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

Download in other formats:

Original Format