Interoperable Physics Driver for NGGPS
set_soilveg.f
Go to the documentation of this file.
1  subroutine set_soilveg(me,isot,ivet,nlunit)
2  use namelist_soilveg
3  implicit none
4 
5  integer, intent(in) :: isot,ivet,nlunit
6  integer me
7 !my begin locals
8 !for 20 igbp veg type and 19 stasgo soil type
9  integer i
10  REAL WLTSMC1,REFSMC1
11 ! ----------------------------------------------------------------------
12 ! SET TWO SOIL MOISTURE WILT, SOIL MOISTURE REFERENCE PARAMETERS
13 ! ----------------------------------------------------------------------
14  REAL SMLOW
15  REAL SMLOW_DATA
16  DATA smlow_data /0.5/
17 
18  REAL SMHIGH
19  REAL SMHIGH_DATA
20 ! changed in 2.6 from 3 to 6 on June 2nd 2003
21 ! DATA SMHIGH_DATA /3.0/
22  DATA smhigh_data /6.0/
23  namelist /soil_veg/ slope_data, rsmtbl, rgltbl, hstbl, snupx,
24  & bb, drysmc, f11, maxsmc, refsmc, satpsi, satdk, satdw,
25  & wltsmc, qtz, lparam, zbot_data, salp_data, cfactr_data,
26  & cmcmax_data, sbeta_data, rsmax_data, topt_data,
27  & refdk_data, frzk_data, bare, defined_veg, defined_soil,
28  & defined_slope, fxexp_data, nroot_data, refkdt_data, z0_data,
29  & czil_data, lai_data, csoil_data
30 
31 cmy end locals
32  if(ivet.eq.2) then
33 
34 !using umd veg table
35  slope_data =(/0.1, 0.6, 1.0, 0.35, 0.55, 0.8,
36  & 0.63, 0.0, 0.0, 0.0, 0.0, 0.0,
37  & 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0,
38  & 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0,
39  & 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0/)
40  rsmtbl =(/300.0, 175.0, 175.0, 300.0, 300.0, 70.0,
41  & 20.0, 225.0, 225.0, 225.0, 400.0, 20.0,
42  & 150.0, 0.0, 0.0, 0.0, 0.0, 0.0,
43  & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
44  & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
45 c-----------------------------
46  rgltbl =(/30.0, 30.0, 30.0, 30.0, 30.0, 65.0,
47  & 100.0, 100.0, 100.0, 100.0, 100.0, 100.0,
48  & 100.0, 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  hstbl =(/41.69, 54.53, 51.93, 47.35, 47.35, 54.53,
52  & 36.35, 42.00, 42.00, 42.00, 42.00, 36.35,
53  & 42.00, 0.00, 0.00, 0.00, 0.00, 0.00,
54  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,
55  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
56 ! changed for version 2.6 on june 2nd 2003
57 ! data snupx /0.080, 0.080, 0.080, 0.080, 0.080, 0.080,
58 ! & 0.040, 0.040, 0.040, 0.040, 0.025, 0.040,
59 ! & 0.025, 0.000, 0.000, 0.000, 0.000, 0.000,
60  snupx =(/0.040, 0.040, 0.040, 0.040, 0.040, 0.040,
61  * 0.020, 0.020, 0.020, 0.020, 0.013, 0.020,
62  * 0.013, 0.000, 0.000, 0.000, 0.000, 0.000,
63  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
64  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
65 
66  bare =11
67 
68 c---------------------------------------------------------------------
69 ! number of defined veg used.
70 ! ----------------------------------------------------------------------
71  defined_veg=13
72  nroot_data =(/4,4,4,4,4,4,3,3,3,2,3,3,2,0,0,
73  & 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/)
74 ! ----------------------------------------------------------------------
75 ! vegetation class-related arrays
76 ! ----------------------------------------------------------------------
77  z0_data =(/2.653, 0.826, 0.563, 1.089, 0.854, 0.856,
78  & 0.035, 0.238, 0.065, 0.076, 0.011, 0.125,
79  & 0.011, 0.000, 0.000, 0.000, 0.000, 0.000,
80  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
81  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
82 clu: change to 3 or 2 oct 15, 2004
83  lai_data =(/3.0, 3.0, 3.0, 3.0, 3.0, 3.0,
84  & 3.0, 3.0, 3.0, 3.0, 3.0, 3.0,
85  & 3.0, 0.0, 0.0, 0.0, 0.0, 0.0,
86  & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
87  & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
88 ! use igbp table
89  elseif(ivet.eq.1)then
90 
91  slope_data =(/1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
92  & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,
93  & 1.0 , 1.0, 1.0, 1.0, 1.0, 1.0,
94  & 1.0 , 1.0, 0.0, 0.0, 0.0, 0.0,
95  & 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0/)
96  rsmtbl =(/300.0, 300.0, 70.0, 175.0, 175.0, 70.0,
97  & 70.0, 70.0, 70.0, 20.0, 40.0, 20.0,
98  & 400.0, 35.0, 200.0, 70.0, 100.0, 70.0,
99  & 150.0, 200.0, 0.0, 0.0,0.0, 0.0,
100  & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
101 c-----------------------------
102  rgltbl =(/30.0, 30.0, 30.0, 30.0, 30.0, 100.0,
103  & 100.0, 65.0, 65.0, 100.0, 100.0, 100.0,
104  & 100.0, 100.0, 100.0,100.0,30.0, 100.0,
105  & 100.0, 100.0, 0.0,0.0,0.0,0.0,
106  & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
107  hstbl =(/47.35, 41.69, 47.35, 54.53, 51.93, 42.00,
108  & 42.00, 42.00, 42.00, 36.35, 60.00, 36.25,
109  & 42.00, 36.25, 42.00, 42.00, 51.75, 42.00,
110  & 42.00, 42.00, 0.00, 0.00, 0.00, 0.00,
111  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
112  snupx =(/0.080, 0.080, 0.080, 0.080, 0.080, 0.020,
113  * 0.020, 0.060, 0.040, 0.020, 0.010, 0.020,
114  * 0.020, 0.020, 0.013, 0.013, 0.010, 0.020,
115  & 0.020, 0.020, 0.000, 0.000, 0.000, 0.000,
116  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
117 
118  bare =16
119 
120 !---------------------------------------------------------------------
121 ! number of defined veg used.
122 ! ----------------------------------------------------------------------
123  defined_veg=20
124 
125  nroot_data =(/4,4,4,4,4,3,3,3,3,3,3,3,1,3,2,
126  & 3,0,3,3,2,0,0,0,0,0,0,0,0,0,0/)
127 ! ----------------------------------------------------------------------
128 ! VEGETATION CLASS-RELATED ARRAYS
129 ! ----------------------------------------------------------------------
130  z0_data =(/1.089, 2.653, 0.854, 0.826, 0.80, 0.05,
131  & 0.03, 0.856, 0.856, 0.15, 0.04, 0.13,
132  & 1.00, 0.25, 0.011, 0.011, 0.001, 0.076,
133  & 0.05, 0.03, 0.000, 0.000, 0.000, 0.000,
134  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
135 
136  lai_data =(/3.0, 3.0, 3.0, 3.0, 3.0, 3.0,
137  & 3.0, 3.0, 3.0, 3.0, 3.0, 3.0,
138  & 3.0, 3.0, 3.0, 3.0, 3.0, 3.0,
139  & 3.0, 3.0, 0.0, 0.0, 0.0, 0.0,
140  & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/)
141 
142 ! end if veg table
143  endif
144 
145  if(isot.eq.0) then
146 
147  bb =(/4.26, 8.72, 11.55, 4.74, 10.73, 8.17,
148  & 6.77, 5.25, 4.26, 0.00, 0.00, 0.00,
149  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,
150  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,
151  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
152 ! !!!!!!!!!!!!!! the following values in the table are not used
153 ! !!!!!!!!!!!!!! and are just given for reference
154  drysmc=(/0.029, 0.119, 0.139, 0.047, 0.100, 0.103,
155  & 0.069, 0.066, 0.029, 0.000, 0.000, 0.000,
156  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
157  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
158  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
159 ! !!!!!!!!!!!!!! the following values in the table are not used
160 ! !!!!!!!!!!!!!! and are just given for reference
161  f11 =(/-0.999, -1.116, -2.137, -0.572, -3.201, -1.302,
162  & -1.519, -0.329, -0.999, 0.000, 0.000, 0.000,
163  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
164  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
165  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
166  maxsmc=(/0.421, 0.464, 0.468, 0.434, 0.406, 0.465,
167  & 0.404, 0.439, 0.421, 0.000, 0.000, 0.000,
168  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
169  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
170  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
171 !
172 ! ----------------------------------------------------------------------
173 ! the following 5 parameters are derived later in redprm.f from the soil
174 ! data, and are just given here for reference and to force static
175 ! storage allocation. -dag lohmann, feb. 2001
176 ! ----------------------------------------------------------------------
177 ! data refsmc/0.283, 0.387, 0.412, 0.312, 0.338, 0.382,
178 ! & 0.315, 0.329, 0.283, 0.000, 0.000, 0.000,
179 ! !!!!!!!!!!!!!! the following values in the table are not used
180 ! !!!!!!!!!!!!!! and are just given for reference
181  refsmc=(/0.248, 0.368, 0.398, 0.281, 0.321, 0.361,
182  & 0.293, 0.301, 0.248, 0.000, 0.000, 0.000,
183  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
184  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
185  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
186 ! ----------------------------------------------------------------------
187 ! soil texture-related arrays.
188 ! ----------------------------------------------------------------------
189  satpsi=(/0.04, 0.62, 0.47, 0.14, 0.10, 0.26,
190  & 0.14, 0.36, 0.04, 0.00, 0.00, 0.00,
191  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,
192  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,
193  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
194  satdk =(/1.41e-5, 0.20e-5, 0.10e-5, 0.52e-5, 0.72e-5,
195  & 0.25e-5, 0.45e-5, 0.34e-5, 1.41e-5, 0.00,
196  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
197  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
198  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
199  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
200  qtz =(/0.82, 0.10, 0.25, 0.60, 0.52, 0.35,
201  & 0.60, 0.40, 0.82, 0.00, 0.00, 0.00,
202  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,
203  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,
204  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
205 
206 ! !!!!!!!!!!!!!! the following values in the table are not used
207 ! !!!!!!!!!!!!!! and are just given for reference
208  wltsmc=(/0.029, 0.119, 0.139, 0.047, 0.100, 0.103,
209  & 0.069, 0.066, 0.029, 0.000, 0.000, 0.000,
210  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
211  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
212  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
213 ! !!!!!!!!!!!!!! the following values in the table are not used
214 ! !!!!!!!!!!!!!! and are just given for reference
215  satdw =(/5.71e-6, 2.33e-5, 1.16e-5, 7.95e-6, 1.90e-5,
216  & 1.14e-5, 1.06e-5, 1.46e-5, 5.71e-6, 0.00,
217  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
218  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
219  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
220  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
221 
222 ! ----------------------------------------------------------------------
223 ! number of defined soiltyps used.
224 ! ----------------------------------------------------------------------
225 
226  defined_soil=9
227 
228  else
229 
230 ! using stasgo table
231  bb =(/4.05, 4.26, 4.74, 5.33, 5.33, 5.25,
232  & 6.77, 8.72, 8.17, 10.73, 10.39, 11.55,
233  & 5.25, 0.00, 4.05, 4.26, 11.55, 4.05,
234  & 4.05, 0.00, 0.00, 0.00, 0.00, 0.00,
235  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
236 ! !!!!!!!!!!!!!! The following values in the table are NOT used
237 ! !!!!!!!!!!!!!! and are just given for reference
238  drysmc=(/0.023, 0.028, 0.047, 0.084, 0.084, 0.066,
239  & 0.069, 0.120, 0.103, 0.100, 0.126, 0.135,
240  & 0.069, 0.000, 0.012, 0.028, 0.135, 0.012,
241  & 0.023, 0.000, 0.000, 0.000, 0.000, 0.000,
242  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
243 ! !!!!!!!!!!!!!! The following values in the table are NOT used
244 ! !!!!!!!!!!!!!! and are just given for reference
245  f11 =(/-1.090, -1.041, -0.568, 0.162, 0.162, -0.327,
246  & -1.535, -1.118, -1.297, -3.211, -1.916, -2.258,
247  & -0.201, 0.000, -2.287, -1.041, -2.258, -2.287,
248  & -1.090, 0.000, 0.000, 0.000, 0.000, 0.000,
249  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
250  maxsmc=(/0.395, 0.421, 0.434, 0.476, 0.476, 0.439,
251  & 0.404, 0.464, 0.465, 0.406, 0.468, 0.457,
252  & 0.464, 0.000, 0.200, 0.421, 0.457, 0.200,
253  & 0.395, 0.000, 0.000, 0.000, 0.000, 0.000,
254  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
255 !
256 ! ----------------------------------------------------------------------
257 ! THE FOLLOWING 5 PARAMETERS ARE DERIVED LATER IN REDPRM.F FROM THE SOIL
258 ! DATA, AND ARE JUST GIVEN HERE FOR REFERENCE AND TO FORCE STATIC
259 ! STORAGE ALLOCATION. -DAG LOHMANN, FEB. 2001
260 ! ----------------------------------------------------------------------
261 ! DATA REFSMC/0.283, 0.387, 0.412, 0.312, 0.338, 0.382,
262 ! & 0.315, 0.329, 0.283, 0.000, 0.000, 0.000,
263 ! !!!!!!!!!!!!!! The following values in the table are NOT used
264 ! !!!!!!!!!!!!!! and are just given for reference
265  refsmc=(/0.236, 0.283, 0.312, 0.360, 0.360, 0.329,
266  & 0.315, 0.387, 0.382, 0.338, 0.404, 0.403,
267  & 0.348, 0.000, 0.133, 0.283, 0.403, 0.133,
268  & 0.236, 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 ! SOIL TEXTURE-RELATED ARRAYS.
272 ! ----------------------------------------------------------------------
273  satpsi=(/0.035, 0.0363, 0.1413, 0.7586, 0.7586, 0.3548,
274  & 0.1349, 0.6166, 0.2630, 0.0977, 0.3236, 0.4677,
275  & 0.3548, 0.00, 0.0350, 0.0363, 0.4677, 0.0350,
276  & 0.0350, 0.00, 0.00, 0.00, 0.00, 0.00,
277  & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00/)
278  satdk =(/1.76e-4, 1.4078e-5, 5.2304e-6, 2.8089e-6, 2.8089e-6,
279  & 3.377e-6, 4.4518e-6, 2.0348e-6, 2.4464e-6, 7.2199e-6,
280  & 1.3444e-6, 9.7394e-7, 3.377e-6, 0.00, 1.4087e-05,
281  & 1.4078e-5, 9.7394e-7, 1.4078e-5, 1.760e-4, 0.00,
282  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
283  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
284  qtz =(/0.92, 0.82, 0.60, 0.25, 0.10, 0.40,
285  & 0.60, 0.10, 0.35, 0.52, 0.10, 0.25,
286  & 0.05, 0.00, 0.07, 0.25, 0.60, 0.52,
287  & 0.92, 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 ! !!!!!!!!!!!!!! The following values in the table are NOT used
291 ! !!!!!!!!!!!!!! and are just given for reference
292  wltsmc=(/0.023, 0.028, 0.047, 0.084, 0.084, 0.066,
293  & 0.069, 0.120, 0.103, 0.100, 0.126, 0.135,
294  & 0.069, 0.000, 0.012, 0.028, 0.135, 0.012,
295  & 0.023, 0.000, 0.000, 0.000, 0.000, 0.000,
296  & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000/)
297 ! !!!!!!!!!!!!!! The following values in the table are NOT used
298 ! !!!!!!!!!!!!!! and are just given for reference
299  satdw =(/0.6316e-4, 0.5171e-5, 0.8072e-5, 0.2386e-4, 0.2386e-4,
300  & 0.1433e-4, 0.1006e-4, 0.2358e-4, 0.1130e-4,0.1864e-04,
301  & 0.9658e-05, 0.1151e-04, 0.1356e-04, 0.00, 0.9978e-05,
302  & 0.5171e-05, 0.1151e-04, 0.9978e-05, 0.6316e-04, 0.00,
303  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00,
304  & 0.00 , 0.00 , 0.00 , 0.00 , 0.00/)
305 ! ----------------------------------------------------------------------
306 ! number of defined soiltyps used.
307 ! ----------------------------------------------------------------------
308 
309  defined_soil=19
310 ! end if soil table
311  endif
312 
313 
314 ! the values shared by different veg/soil type data
315 
316 ! PT 5/18/2015 - changed to FALSE to match atm_namelist setting
317 ! PT LPARAM is not used anywhere
318 ! LPARAM =.TRUE.
319  lparam =.false.
320 
321 ! changed for version 2.5.2
322 ! data zbot_data /-3.0/
323  zbot_data =-8.0
324 ! changed for version 2.6 june 2nd 2003
325 ! data salp_data /2.6/
326  salp_data =4.0
327  cfactr_data =0.5
328  cmcmax_data =0.5e-3
329  sbeta_data =-2.0
330  rsmax_data =5000.0
331  topt_data =298.0
332  refdk_data =2.0e-6
333  frzk_data =0.15
334 
335  defined_slope=9
336  fxexp_data =2.0
337  refkdt_data =3.0
338 ! changed in version 2.6 june 2nd 2003
339 ! data czil_data /0.2/
340  czil_data =0.075
341 
342 ! DATA CSOIL_DATA /1.26E+6/
343  csoil_data = 2.00e+6
344 ! ----------------------------------------------------------------------
345 ! READ NAMELIST FILE TO OVERRIDE DEFAULT PARAMETERS ONLY ONCE.
346 ! NAMELIST_NAME must be 50 characters or less.
347 ! ----------------------------------------------------------------------
348 !lu: namelist is set up in run script
349 !PT if (me == 0) write(0,*) 'read namelist cwsoilvegSOIL_VEG'
350 !$$$ READ(5, SOIL_VEG)
351 !PT rewind(nlunit)
352 !PT READ(nlunit, SOIL_VEG)
353 
354 !* WRITE(6, SOIL_VEG)
355 ! OPEN(58, FILE = 'namelist_filename.txt')
356 ! READ(58,'(A)') NAMELIST_NAME
357 ! CLOSE(58)
358 ! WRITE(0,*) 'Namelist Filename is ', NAMELIST_NAME
359 ! OPEN(59, FILE = NAMELIST_NAME)
360 ! 50 CONTINUE
361 ! READ(59, SOIL_VEG, END=100)
362 ! IF (LPARAM) GOTO 50
363 ! 100 CONTINUE
364 ! CLOSE(59)
365 
366  IF (defined_soil .GT. max_soiltyp) THEN
367  WRITE(0,*) 'Warning: DEFINED_SOIL too large in namelist'
368  stop 222
369  ENDIF
370  IF (defined_veg .GT. max_vegtyp) THEN
371  WRITE(0,*) 'Warning: DEFINED_VEG too large in namelist'
372  stop 222
373  ENDIF
374  IF (defined_slope .GT. max_slopetyp) THEN
375  WRITE(0,*) 'Warning: DEFINED_SLOPE too large in namelist'
376  stop 222
377  ENDIF
378 
379  smlow = smlow_data
380  smhigh = smhigh_data
381 
382  DO i = 1,defined_soil
383  satdw(i) = bb(i)*satdk(i)*(satpsi(i)/maxsmc(i))
384  f11(i) = alog10(satpsi(i)) + bb(i)*alog10(maxsmc(i)) + 2.0
385  refsmc1 = maxsmc(i)*(5.79e-9/satdk(i))
386  & **(1.0/(2.0*bb(i)+3.0))
387  refsmc(i) = refsmc1 + (maxsmc(i)-refsmc1) / smhigh
388  wltsmc1 = maxsmc(i) * (200.0/satpsi(i))**(-1.0/bb(i))
389  wltsmc(i) = wltsmc1 - smlow * wltsmc1
390 
391 ! ----------------------------------------------------------------------
392 ! CURRENT VERSION DRYSMC VALUES THAT EQUATE TO WLTSMC.
393 ! FUTURE VERSION COULD LET DRYSMC BE INDEPENDENTLY SET VIA NAMELIST.
394 ! ----------------------------------------------------------------------
395  drysmc(i) = wltsmc(i)
396  END DO
397 
398 ! if (me == 0) write(6,soil_veg)
399  return
400  end
set_soilveg
subroutine set_soilveg(me, isot, ivet, nlunit)
Definition: set_soilveg.f:2