CCPP SciDoc v7.0.0  v7.0.0
Common Community Physics Package Developed at DTC
 
Loading...
Searching...
No Matches
set_soilveg_ruc.F90
1
2
5 module set_soilveg_ruc_mod
6
7 use machine , only : kind_phys
8 use namelist_soilveg_ruc
9
10 implicit none
11
12 private
13
14 public :: set_soilveg_ruc
15
16 contains
17
21 subroutine set_soilveg_ruc(me,isot,ivet,nlunit,errmsg,errflg)
22
23 integer, intent(in) :: isot,ivet,nlunit
24 character(len=*), intent(out) :: errmsg
25 integer, intent(out) :: errflg
26 integer me
27
28 integer i
29 real(kind_phys) refsmc1, wltsmc1
30
31 namelist /soil_veg_ruc/ slope_data, albtbl, z0tbl, lemitbl, &
32 & pctbl, shdtbl, &
33 & ifortbl, rstbl, rgltbl, hstbl, snuptbl, laitbl, maxalb, &
34 & mfsno, sncovfac, &
35 & lparam, topt_data, cmcmax_data, cfactr_data, &
36 & rsmax_data, bare, glacier, natural, crop, urban, &
37 & defined_veg, defined_soil, defined_slope, &
38 & bb, drysmc, hc, maxsmc, refsmc, satpsi, satdk, satdw, &
39 & wltsmc, qtz, refsmcnoah, wltsmcnoah, maxsmcnoah
40
41 ! Initialize error-handling
42 errflg = 0
43 errmsg = ''
44
45 if(ivet.eq.2) then
46! Using umd veg classification
47 slope_data =(/0.1, 0.6, 1.0, 0.35, 0.55, 0.8, &
48 & 0.63, 0.0, 0.0, 0.0, 0.0, 0.0, &
49 & 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0, &
50 & 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0, &
51 & 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0/)
52! ----------------------------------------------------------------------
53! vegetation class-related arrays
54! ----------------------------------------------------------------------
55 rstbl =(/300.0, 175.0, 175.0, 300.0, 300.0, 70.0, &
56 & 20.0, 225.0, 225.0, 225.0, 400.0, 20.0, &
57 & 150.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
58 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
59 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
60
61 rgltbl =(/30.0, 30.0, 30.0, 30.0, 30.0, 65.0, &
62 & 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, &
63 & 100.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
64 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
65 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
66
67 hstbl =(/41.69, 54.53, 51.93, 47.35, 47.35, 54.53, &
68 & 36.35, 42.00, 42.00, 42.00, 42.00, 36.35, &
69 & 42.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
70 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
71 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
72
73 snuptbl =(/0.040, 0.040, 0.040, 0.040, 0.040, 0.040, &
74 & 0.020, 0.020, 0.020, 0.020, 0.013, 0.020, &
75 & 0.013, 0.000, 0.000, 0.000, 0.000, 0.000, &
76 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
77 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
78
79 bare =11
80
81!---------------------------------------------------------------------
82! number of defined veg used.
83! ----------------------------------------------------------------------
84 defined_veg=13
85 defined_slope = 9
86 z0tbl =(/2.653, 0.826, 0.563, 1.089, 0.854, 0.856, &
87 & 0.035, 0.238, 0.065, 0.076, 0.011, 0.125, &
88 & 0.011, 0.000, 0.000, 0.000, 0.000, 0.000, &
89 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
90 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
91 laitbl =(/3.0, 3.0, 3.0, 3.0, 3.0, 3.0, &
92 & 3.0, 3.0, 3.0, 3.0, 3.0, 3.0, &
93 & 3.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
94 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
95 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
96
97! Using igbp ivegetation classification - MODI-RUC
98 elseif(ivet.eq.1)then
99!--------------------------------------------------------------------------------------------------------------------------------------
100!21,1, 'ALBEDO Z0 LEMI PC SHDFAC IFOR RS RGL HS SNUP LAI MAXALB'
101!1 .12, .80, .950, .55, .70, 1, 125., 30., 47.35, 0.08, 6.40, 52., 'Evergreen Needleleaf Forest'
102!2, .12, .80, .950, .55, .95, 2, 150., 30., 41.69, 0.08, 6.48, 35., 'Evergreen Broadleaf Forest'
103!3, .14, .80, .940, .55, .70, 4, 150., 30., 47.35, 0.08, 5.16, 54., 'Deciduous Needleleaf Forest'
104!4, .16, .80, .930, .55, .80, 3, 100., 30., 54.53, 0.08, 3.31, 58., 'Deciduous Broadleaf Forest'
105!5, .13, .80, .940, .55, .80, 2, 125., 30., 51.93, 0.08, 5.50, 53., 'Mixed Forests'
106!6, .22, .10, .930, .40, .70, 4, 300., 100., 42.00, 0.03, 3.66, 60., 'Closed Shrublands'
107!7, .20, .10, .880, .40, .70, 4, 170., 100., 39.18, 0.035, 2.60, 65., 'Open Shrublands'
108!8, .20, .30, .930, .40, .70, 5, 300., 100., 42.00, 0.03, 3.66, 60., 'Woody Savannas'
109!9, .20, .15, .920, .40, .50, 5, 70., 65., 54.53, 0.04, 3.66, 50., 'Savannas'
110!10, .19, .075, .920, .40, .80, 5, 40., 100., 36.35, 0.04, 2.90, 70., 'Grasslands'
111!11 .14, .30, .950, .40, .60, 4, 70., 65., 55.97 0.015 5.72, 59., 'Permanent wetlands'
112!12, .18, .20, .935, .40, .80, 7, 40., 100., 36.25, 0.04, 5.68, 66., 'Croplands'
113!13, .18, 0.5, .880, .40, .10, 9, 200., 999., 999.0, 0.04, 1.00, 46., 'Urban and Built-Up'
114!14 .16, .30, .920, .40, .80, 7, 40., 100., 36.25, 0.04, 4.29, 68., 'cropland/natural vegetation mosaic'
115!15, .55, .011, .980, .00, .00, 9, 999., 999., 999.0, 0.02, 0.01, 82., 'Snow and Ice'
116!16, .25, .065, .850, .30, .01, 5, 999., 999., 999.0, 0.02, 0.75, 75., 'Barren or Sparsely Vegetated'
117!17, .08, .0001, .980, .00, .00, 9, 100., 30., 51.75, 0.01, 0.01, 70., 'Water'
118!18, .15, .15, .930, .40, .60, 5, 150., 100., 42.00, 0.025, 3.35, 55., 'Wooded Tundra'
119!19, .15, .075, .920, .40, .60, 5, 150., 100., 42.00, 0.025, 3.35, 60., 'Mixed Tundra'
120!20, .15, .06, .900, .30, .30, 5, 200., 100., 42.00, 0.02, 3.35, 75., 'Barren Tundra'
121!21, .08, .0001, .980, .00, .00, 9, 100., 30., 51.75, 0.01, 0.01, 70., 'Lakes'
122!--------------------------------------------------------------------------------------------------------------------------------------
123
124! number of defined veg used.
125 defined_veg=21
126 defined_slope=20
127
128 slope_data =(/1.0, 1.0, 1.0, 1.0, 1.0, 1.0, &
129 & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, &
130 & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, &
131 & 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, &
132 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
133! ----------------------------------------------------------------------
134! VEGETATION CLASS-RELATED ARRAYS
135! ----------------------------------------------------------------------
136 rstbl = (/125.0, 150.0, 150.0, 100.0, 125.0, 300.0, &
137 & 170.0, 300.0, 70.0, 40.0, 70.0, 40.0, &
138 & 200.0, 40.0, 999.0, 999.0, 100.0, 150.0, &
139 & 150.0, 200.0, 100.0, 0.0, 0.0, 0.0, &
140 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
141
142 rgltbl = (/30.0, 30.0, 30.0, 30.0, 30.0, 100.0, &
143 & 100.0, 100.0, 65.0, 100.0, 65.0, 100.0, &
144 & 999.0, 100.0, 999.0, 999.0, 30.0, 100.0, &
145 & 100.0, 100.0, 30.0, 0.0, 0.0, 0.0, &
146 & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
147
148 hstbl = (/47.35, 41.69, 47.35, 54.53, 51.93, 42.00, &
149 & 39.18, 42.00, 54.53, 36.35, 55.97, 36.25, &
150 & 999.00, 36.25,999.00,999.00, 51.75, 42.00, &
151 & 42.00, 42.00, 51.75, 0.00, 0.00, 0.00, &
152 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
153
154 snuptbl =(/0.080, 0.080, 0.080, 0.080, 0.080, 0.030, &
155 & 0.035, 0.030, 0.040, 0.040, 0.015, 0.040, &
156 & 0.040, 0.040, 0.020, 0.020, 0.010, 0.025, &
157 & 0.025, 0.020, 0.010, 0.000, 0.000, 0.000, &
158 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
159
160 albtbl =(/0.12, 0.12, 0.14, 0.16, 0.13, 0.22, 0.20, &
161 & 0.20, 0.20, 0.19, 0.14, 0.18, 0.18, 0.16, &
162 & 0.55, 0.25, 0.08, 0.15, 0.15, 0.15, 0.08, &
163 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
164 & 0.00, 0.00/)
165
166 z0tbl =(/0.80, 0.80, 0.80, 0.80, 0.80, 0.10, &
167 & 0.10, 0.30, 0.15, 0.075, 0.30, 0.20, &
168 & 0.50, 0.30, 0.011, 0.065, 0.0001,0.15, &
169 & 0.075, 0.06, 0.0001,0.000, 0.000, 0.000, &
170 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
171
172 lemitbl =(/.950, .950, .940, .930, .940, .930, .880, &
173 & .930, .920, .920, .950, .935, .880, .920, &
174 & .980, .850, .980, .930, .920, .900, .980, &
175 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
176 & 0.00, 0.00/)
177
178 pctbl =(/0.55, 0.55, 0.55, 0.55, 0.55, 0.40, 0.40, &
179 & 0.40, 0.40, 0.40, 0.40, 0.40, 0.40, 0.40, &
180 & 0.00, 0.30, 0.00, 0.40, 0.40, 0.30, 0.00, &
181 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
182 & 0.00, 0.00/)
183
184 shdtbl =(/0.70, 0.95, 0.70, 0.80, 0.80, 0.70, 0.70, &
185 & 0.70, 0.50, 0.80, 0.60, 0.80, 0.10, 0.80, &
186 & 0.00, 0.01, 0.00, 0.60, 0.60, 0.30, 0.00, &
187 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
188 & 0.00, 0.00/)
189
190 ifortbl =(/1, 2, 4, 3, 2, 4, 4, 5, 5, 5, 4, 7, 9, 7, &
191 & 9, 9, 9, 5, 5, 5, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0/)
192
193 laitbl =(/2.80, 5.18, 4.16, 4.81, 4.20, 1.16, 0.90, &
194 & 3.00, 3.00, 1.10, 1.72, 2.58, 1.00, 2.29, &
195 & 0.01, 0.75, 0.01, 1.00, 1.00, 1.00, 0.01, &
196 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
197 & 0.00, 0.00/)
198
199 maxalb =(/52., 35., 54., 58., 53., 60., 65., 60., &
200 & 50., 70., 59., 66., 46., 68., 82., 75., &
201 & 70., 55., 60., 75., 70., 0., 0., 0., &
202 & 0., 0., 0., 0., 0., 0./)
203
204 mfsno = &
205 & (/ 1.00, 1.00, 1.00, 1.00, 2.00, 2.00, &
206 & 2.00, 2.00, 2.00, 2.00, 2.00, 2.00, &
207 & 3.00, 3.00, 2.00, 2.00, 2.00, 2.00, &
208 & 2.00, 2.00, 0.00, 0.00, 0.00, 0.00, &
209! & 3.00, 3.00, 2.00, 3.00, 3.00, 3.00, &
210! & 3.00, 3.00, 0.00, 0.00, 0.00, 0.00, &
211 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 /)
212
213!-- Noah MP snowmelt curve values
214! & (/ 1.00, 1.00, 1.00, 1.00, 1.00, 2.00, &
215! & 2.00, 2.00, 2.00, 2.00, 3.00, 3.00, &
216! & 4.00, 4.00, 2.50, 3.00, 3.00, 3.50, &
217! & 3.50, 3.50, 0.00, 0.00, 0.00, 0.00, &
218! & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 /)
219
220 sncovfac = &
221 & (/ 0.030, 0.030, 0.030, 0.030, 0.030, &
222 !& (/ 0.008, 0.008, 0.008, 0.008, 0.008, &
223 & 0.016, 0.016, 0.020, 0.020, 0.020, &
224 & 0.020, 0.014, 0.042, 0.026, 0.030, &
225 & 0.016, 0.030, 0.030, 0.030, 0.030, &
226 & 0.000, 0.000, 0.000, 0.000, 0.000, &
227 & 0.000, 0.000, 0.000, 0.000, 0.000 /)
228
229 natural = 10
230 crop = 12
231 urban = 13
232 glacier = 15
233 bare = 16
234 endif
235! end if veg table
236
237 topt_data =298.0
238 cmcmax_data =0.2e-3
239 cfactr_data =0.5
240 rsmax_data =5000.0
241
242 if(isot.eq.0) then
243
244! ----------------------------------------------------------------------
245! soil texture-related arrays.
246! ----------------------------------------------------------------------
247 bb =(/4.26, 8.72, 11.55, 4.74, 10.73, 8.17, &
248 & 6.77, 5.25, 4.26, 0.00, 0.00, 0.00, &
249 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
250 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
251 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
252
253 drysmc =(/0.029, 0.119, 0.139, 0.047, 0.100, 0.103, &
254 & 0.069, 0.066, 0.029, 0.000, 0.000, 0.000, &
255 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
256 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
257 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
258
259 maxsmc =(/0.421, 0.464, 0.468, 0.434, 0.406, 0.465, &
260 & 0.404, 0.439, 0.421, 0.000, 0.000, 0.000, &
261 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
262 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
263 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
264
265 refsmc =(/0.248, 0.368, 0.398, 0.281, 0.321, 0.361, &
266 & 0.293, 0.301, 0.248, 0.000, 0.000, 0.000, &
267 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
268 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
269 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
270
271 satpsi =(/0.04, 0.62, 0.47, 0.14, 0.10, 0.26, &
272 & 0.14, 0.36, 0.04, 0.00, 0.00, 0.00, &
273 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
274 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
275 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
276
277 satdk =(/1.41e-5, 0.20e-5, 0.10e-5, 0.52e-5, 0.72e-5, &
278 & 0.25e-5, 0.45e-5, 0.34e-5, 1.41e-5, 0.00, &
279 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
280 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
281 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
282 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
283
284 qtz =(/0.82, 0.10, 0.25, 0.60, 0.52, 0.35, &
285 & 0.60, 0.40, 0.82, 0.00, 0.00, 0.00, &
286 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
287 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
288 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
289
290 wltsmc =(/0.029, 0.119, 0.139, 0.047, 0.100, 0.103, &
291 & 0.069, 0.066, 0.029, 0.000, 0.000, 0.000, &
292 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
293 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, &
294 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
295
296 satdw =(/5.71e-6, 2.33e-5, 1.16e-5, 7.95e-6, 1.90e-5, &
297 & 1.14e-5, 1.06e-5, 1.46e-5, 5.71e-6, 0.00, &
298 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
299 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
300 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
301 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
302
303! ----------------------------------------------------------------------
304! number of defined soiltypes used.
305! ----------------------------------------------------------------------
306
307 defined_soil=9
308
309 else
310
311! using STASGO soil classification
312!--------------------------------------------------------------------------------------------------------------------------------------
313!19,1 'BB DRYSMC HC MAXSMC REFSMC SATPSI SATDK SATDW WLTSMC QTZ '
314!1, 4.05, 0.002, 1.47, 0.395, 0.174, 0.121, 1.76E-4, 0.608E-6, 0.033, 0.92, 'SAND'
315!2, 4.38, 0.035, 1.41, 0.410, 0.179, 0.090, 1.56E-4, 0.514E-5, 0.055, 0.82, 'LOAMY SAND'
316!3, 4.90, 0.041, 1.34, 0.435, 0.249, 0.218, 3.47E-5, 0.805E-5, 0.095, 0.60, 'SANDY LOAM'
317!4, 5.30, 0.034, 1.27, 0.485, 0.369, 0.786, 7.20E-6, 0.239E-4, 0.143, 0.25, 'SILT LOAM'
318!5, 5.30, 0.034, 1.27, 0.485, 0.369, 0.786, 7.20E-6, 0.239E-4, 0.143, 0.10, 'SILT'
319!6, 5.39, 0.050, 1.21, 0.451, 0.314, 0.478, 6.95E-6, 0.143E-4, 0.137, 0.40, 'LOAM'
320!7, 7.12, 0.068, 1.18, 0.420, 0.299, 0.299, 6.30E-6, 0.990E-5, 0.148, 0.60, 'SANDY CLAY LOAM'
321!8, 7.75, 0.060, 1.32, 0.477, 0.357, 0.356, 1.70E-6, 0.237E-4, 0.208, 0.10, 'SILTY CLAY LOAM'
322!9, 8.52, 0.085, 1.23, 0.476, 0.391, 0.630, 2.45E-6, 0.113E-4, 0.230, 0.35, 'CLAY LOAM'
323!10, 10.40, 0.100, 1.18, 0.426, 0.316, 0.153, 2.17E-6, 0.187E-4, 0.210, 0.52, 'SANDY CLAY'
324!11, 10.40, 0.070, 1.15, 0.492, 0.409, 0.490, 1.03E-6, 0.964E-5, 0.250, 0.10, 'SILTY CLAY'
325!12, 11.40, 0.068, 1.09, 0.482, 0.400, 0.405, 1.28E-6, 0.112E-4, 0.268, 0.25, 'CLAY'
326!13, 5.39, 0.027, 1.21, 0.451, 0.314, 0.478, 6.95E-6, 0.143E-4, 0.117, 0.05, 'ORGANIC MATERIAL'
327!14, 0.0, 0.0, 4.18, 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.00, 'WATER'
328!15, 4.05, 0.004, 2.03, 0.200, 0.10 , 0.121, 1.41E-4, 0.136E-3, 0.006, 0.60, 'BEDROCK'
329!16, 4.90, 0.065, 2.10, 0.435, 0.249, 0.218, 3.47E-5, 0.514E-5, 0.114, 0.05, 'OTHER(land-ice)'
330!17, 11.40, 0.030, 1.41, 0.468, 0.454, 0.468, 9.74E-7, 0.112E-4, 0.030, 0.60, 'PLAYA'
331!18, 4.05, 0.006, 1.41, 0.200, 0.17, 0.069, 1.41E-4, 0.136E-3, 0.006, 0.52, 'LAVA'
332!19, 4.05, 0.01, 1.47, 0.339, 0.236, 0.069, 1.76E-4, 0.608E-6, 0.060, 0.92, 'WHITE SAND'
333!--------------------------------------------------------------------------------------------------------------------------------------
334
335! number of defined soiltyps used.
336 defined_soil=19
337! ----------------------------------------------------------------------
338! SOIL TEXTURE-RELATED ARRAYS.
339! ----------------------------------------------------------------------
340 bb =(/4.05, 4.38, 4.90, 5.30, 5.30, 5.39, &
341 & 7.12, 7.75, 5.39, 10.40, 10.40, 11.40, &
342 & 5.39, 0.00, 4.05, 4.90, 11.40, 4.05, &
343 & 4.05, 0.00, 0.00, 0.00, 0.00, 0.00, &
344 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
345
346 drysmc =(/0.002, 0.035, 0.041, 0.034, 0.034, 0.050, &
347 & 0.068, 0.060, 0.050, 0.070, 0.070, 0.068, &
348 & 0.027, 0.000, 0.004, 0.065, 0.030, 0.006, &
349 & 0.010, 0.000, 0.000, 0.000, 0.000, 0.000, &
350 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
351
352 hc =(/1.47, 1.41, 1.34, 1.27, 1.27, 1.21, 1.18, &
353 & 1.32, 1.21, 1.18, 1.15, 1.09, 1.21, 4.18, &
354 & 2.03, 2.10, 1.41, 1.41, 1.47, 0.00, 0.00, &
355 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, &
356 & 0.00, 0.00/)
357
358 maxsmc =(/0.395, 0.410, 0.435, 0.485, 0.485, 0.451, &
359 & 0.420, 0.477, 0.451, 0.426, 0.492, 0.482, &
360 & 0.451, 1.000, 0.200, 0.435, 0.468, 0.200, &
361 & 0.339, 0.000, 0.000, 0.000, 0.000, 0.000, &
362 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
363
364 refsmc =(/0.174, 0.179, 0.249, 0.369, 0.369, 0.314, &
365 & 0.299, 0.357, 0.314, 0.316, 0.409, 0.400, &
366 & 0.314, 1.000, 0.100, 0.249, 0.454, 0.170, &
367 & 0.236, 0.000, 0.000, 0.000, 0.000, 0.000, &
368 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
369
370 satpsi =(/0.121, 0.090, 0.218, 0.786, 0.786, 0.478, &
371 & 0.299, 0.356, 0.478, 0.153, 0.490, 0.405, &
372 & 0.478, 0.000, 0.121, 0.218, 0.468, 0.069, &
373 & 0.069, 0.00, 0.00, 0.00, 0.00, 0.00, &
374 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
375
376 satdk =(/1.76e-4, 1.56e-4, 3.47e-5, 7.20e-6, 7.20e-6, &
377 & 6.95e-6, 6.30e-6, 1.70e-6, 6.95e-6, 2.17e-6, &
378 & 1.03e-6, 1.28e-6, 6.95e-6, 0.0, 1.41e-4, &
379 & 3.47e-5, 9.74e-7, 1.41e-4, 1.76e-4, 0.00, &
380 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
381 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
382
383 satdw =(/0.608e-6, 0.514e-5, 0.805e-5, 0.239e-4, 0.239e-4, &
384 & 0.143e-4, 0.990e-5, 0.237e-4, 0.143e-4, 0.187e-4, &
385 & 0.964e-5, 0.112e-4, 0.143e-4, 0.0, 0.136e-03, &
386 & 0.514e-5, 0.112e-4, 0.136e-3, 0.608e-6, 0.00, &
387 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
388 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
389
390 wltsmc =(/0.033, 0.055, 0.095, 0.143, 0.143, 0.137, &
391 & 0.148, 0.170, 0.137, 0.158, 0.190, 0.198, &
392 & 0.117, 0.000, 0.006, 0.114, 0.030, 0.006, &
393 & 0.060, 0.000, 0.000, 0.000, 0.000, 0.000, &
394 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
395
396 qtz =(/0.92, 0.82, 0.60, 0.25, 0.10, 0.40, &
397 & 0.60, 0.10, 0.40, 0.52, 0.10, 0.25, &
398 & 0.05, 0.00, 0.60, 0.05, 0.60, 0.52, &
399 & 0.92, 0.00, 0.00, 0.00, 0.00, 0.00, &
400 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
401
402! Noah parameter to compute SMCREFnoah and SMCWLTnoah
403 bbnoah =(/4.05, 4.26, 4.74, 5.33, 5.33, 5.25, &
404 & 6.77, 8.72, 8.17, 10.73, 10.39, 11.55, &
405 & 5.25, 4.26, 4.05, 4.26, 11.55, 4.05, &
406 & 4.05, 0.00, 0.00, 0.00, 0.00, 0.00, &
407 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
408 maxsmcnoah=(/0.395, 0.421, 0.434, 0.476, 0.476, 0.439, &
409 & 0.404, 0.464, 0.465, 0.406, 0.468, 0.457, &
410 & 0.464, 0.421, 0.200, 0.421, 0.457, 0.200, &
411 & 0.395, 0.000, 0.000, 0.000, 0.000, 0.000, &
412 & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
413!
414 satpsinoah=(/0.035, 0.0363, 0.1413, 0.7586, 0.7586, 0.3548, &
415 & 0.1349, 0.6166, 0.2630, 0.0977, 0.3236, 0.4677, &
416 & 0.3548, 0.0363, 0.0350, 0.0363, 0.4677, 0.0350, &
417 & 0.0350, 0.00, 0.00, 0.00, 0.00, 0.00, &
418 & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
419 satdknoah =(/1.76e-4, 1.4078e-5, 5.2304e-6, 2.8089e-6, 2.8089e-6,&
420 & 3.377e-6, 4.4518e-6, 2.0348e-6, 2.4464e-6, 7.2199e-6,&
421 & 1.3444e-6, 9.7394e-7, 3.377e-6, 1.4078e-5, 1.4087e-05,&
422 & 1.4078e-5, 9.7394e-7, 1.4078e-5, 1.760e-4, 0.00, &
423 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00, &
424 & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
425
426 endif
427! end if soil table
428
429 DO i = 1,defined_soil
430 if (satdknoah(i) /= 0.0 .and. bbnoah(i) > 0.0) then
431 refsmc1 = maxsmcnoah(i)*(5.79e-9/satdknoah(i)) &
432 & **(1.0/(2.0*bbnoah(i)+3.0))
433 refsmcnoah(i) = refsmc1 + (maxsmcnoah(i)-refsmc1) / 6.
434 wltsmc1 = maxsmcnoah(i) * (200.0/satpsinoah(i))**(-1.0/bbnoah(i))
435 wltsmcnoah(i) = wltsmc1 - 0.5 * wltsmc1
436 endif
437 END DO
438
439! PT 5/18/2015 - changed to FALSE to match atm_namelist setting
440! PT LPARAM is not used anywhere
441 lparam =.false.
442
443 IF (defined_soil .GT. max_soiltyp) THEN
444 WRITE(0,*) 'Warning: DEFINED_SOIL too large in namelist'
445 errflg = 1
446 errmsg = 'ERROR(set_soilveg_ruc): DEFINED_SOIL too large in namelist'
447 return
448 ENDIF
449 IF (defined_veg .GT. max_vegtyp) THEN
450 WRITE(0,*) 'Warning: DEFINED_VEG too large in namelist'
451 errflg = 1
452 errmsg = 'ERROR(set_soilveg_ruc): DEFINED_VEG too large in namelist'
453 return
454 ENDIF
455 IF (defined_slope .GT. max_slopetyp) THEN
456 WRITE(0,*) 'Warning: DEFINED_SLOPE too large in namelist'
457 errflg = 1
458 errmsg = 'ERROR(set_soilveg_ruc): DEFINED_SLOPE too large in namelist'
459 return
460 ENDIF
461
462! if (me == 0) write(6,soil_veg_ruc)
463 return
464 end subroutine set_soilveg_ruc
465 end module set_soilveg_ruc_mod
set_soilveg_ruc_mod::set_soilveg_ruc
subroutine, public set_soilveg_ruc(me, isot, ivet, nlunit, errmsg, errflg)
This subroutine specifies vegetation and soil parameters for a given soil and land-use classification...
Definition set_soilveg_ruc.F90:22
namelist_soilveg_ruc
This module contains the namelist options of soil/vegetation in RUC.
Definition namelist_soilveg_ruc.F90:5
set_soilveg_ruc_mod
This module contains subroutine to specify vegetation and soil parameters for a given soild and land-...
Definition set_soilveg_ruc.F90:5