CCPP SciDoc v7.0.0  v7.0.0
Common Community Physics Package Developed at DTC
 
Loading...
Searching...
No Matches
noahmp_tables.f90
1
3
11use machine , only : kind_phys
12 implicit none
13
14 integer, private, parameter :: mvt = 30 ! use 30 instead of 27
15 integer, private, parameter :: mband = 2
16 integer, private, parameter :: msc = 20
17 integer, private, parameter :: max_soiltyp = 30
18 integer, private, parameter :: ncrop = 5
19 integer, private, parameter :: nstage = 8
20 integer, private, parameter :: num_slope = 9
21
22! mptable.tbl vegetation parameters
23
24 integer :: isurban_table
25 integer :: iswater_table
26 integer :: isbarren_table
27 integer :: isice_table
28 integer :: iscrop_table
29 integer :: eblforest_table
30 integer :: natural_table
31 integer :: lcz_1_table
32 integer :: lcz_2_table
33 integer :: lcz_3_table
34 integer :: lcz_4_table
35 integer :: lcz_5_table
36 integer :: lcz_6_table
37 integer :: lcz_7_table
38 integer :: lcz_8_table
39 integer :: lcz_9_table
40 integer :: lcz_10_table
41 integer :: lcz_11_table
42 real (kind=kind_phys) :: ch2op_table(mvt)
43 real (kind=kind_phys) :: dleaf_table(mvt)
44 real (kind=kind_phys) :: z0mvt_table(mvt)
45 real (kind=kind_phys) :: hvt_table(mvt)
46 real (kind=kind_phys) :: hvb_table(mvt)
47 real (kind=kind_phys) :: z0mhvt_table(mvt)
48 real (kind=kind_phys) :: den_table(mvt)
49 real (kind=kind_phys) :: rc_table(mvt)
50 real (kind=kind_phys) :: mfsno_table(mvt)
51 real (kind=kind_phys) :: scffac_table(mvt)
52 real (kind=kind_phys) :: cbiom_table(mvt)
53 real (kind=kind_phys) :: saim_table(mvt,12)
54 real (kind=kind_phys) :: laim_table(mvt,12)
55 real (kind=kind_phys) :: sla_table(mvt)
56 real (kind=kind_phys) :: dilefc_table(mvt)
57 real (kind=kind_phys) :: dilefw_table(mvt)
58 real (kind=kind_phys) :: fragr_table(mvt)
59 real (kind=kind_phys) :: ltovrc_table(mvt)
60
61 real (kind=kind_phys) :: c3psn_table(mvt)
62 real (kind=kind_phys) :: kc25_table(mvt)
63 real (kind=kind_phys) :: akc_table(mvt)
64 real (kind=kind_phys) :: ko25_table(mvt)
65 real (kind=kind_phys) :: ako_table(mvt)
66 real (kind=kind_phys) :: vcmx25_table(mvt)
67 real (kind=kind_phys) :: avcmx_table(mvt)
68 real (kind=kind_phys) :: bp_table(mvt)
69 real (kind=kind_phys) :: mp_table(mvt)
70 real (kind=kind_phys) :: qe25_table(mvt)
71 real (kind=kind_phys) :: aqe_table(mvt)
72 real (kind=kind_phys) :: rmf25_table(mvt)
73 real (kind=kind_phys) :: rms25_table(mvt)
74 real (kind=kind_phys) :: rmr25_table(mvt)
75 real (kind=kind_phys) :: arm_table(mvt)
76 real (kind=kind_phys) :: folnmx_table(mvt)
77 real (kind=kind_phys) :: tmin_table(mvt)
78
79 real (kind=kind_phys) :: xl_table(mvt)
80 real (kind=kind_phys) :: rhol_table(mvt,mband)
81 real (kind=kind_phys) :: rhos_table(mvt,mband)
82 real (kind=kind_phys) :: taul_table(mvt,mband)
83 real (kind=kind_phys) :: taus_table(mvt,mband)
84
85 real (kind=kind_phys) :: mrp_table(mvt)
86 real (kind=kind_phys) :: cwpvt_table(mvt)
87
88 real (kind=kind_phys) :: wrrat_table(mvt)
89 real (kind=kind_phys) :: wdpool_table(mvt)
90 real (kind=kind_phys) :: tdlef_table(mvt)
91
92 real (kind=kind_phys) :: nroot_table(mvt)
93 real (kind=kind_phys) :: rgl_table(mvt)
94 real (kind=kind_phys) :: rs_table(mvt)
95 real (kind=kind_phys) :: hs_table(mvt)
96 real (kind=kind_phys) :: topt_table(mvt)
97 real (kind=kind_phys) :: rsmax_table(mvt)
98
99! soilparm.tbl parameters
100
101 integer :: slcats
102
103 real (kind=kind_phys) :: bexp_table(max_soiltyp)
104 real (kind=kind_phys) :: smcdry_table(max_soiltyp)
105 real (kind=kind_phys) :: f1_table(max_soiltyp)
106 real (kind=kind_phys) :: smcmax_table(max_soiltyp)
107 real (kind=kind_phys) :: smcref_table(max_soiltyp)
108 real (kind=kind_phys) :: psisat_table(max_soiltyp)
109 real (kind=kind_phys) :: dksat_table(max_soiltyp)
110 real (kind=kind_phys) :: dwsat_table(max_soiltyp)
111 real (kind=kind_phys) :: smcwlt_table(max_soiltyp)
112 real (kind=kind_phys) :: quartz_table(max_soiltyp)
113 real (kind=kind_phys) :: bvic_table(max_soiltyp) !vic model infiltration parameter (-) for opt_run=6
114 real (kind=kind_phys) :: axaj_table(max_soiltyp) !Xinanjiang: Tension water distribution inflection parameter [-] for opt_run=7
115 real (kind=kind_phys) :: bxaj_table(max_soiltyp) !Xinanjiang: Tension water distribution shape parameter [-] for opt_run=7
116 real (kind=kind_phys) :: xxaj_table(max_soiltyp) !Xinanjiang: Free water distribution shape parameter [-] for opt_run=7
117 real (kind=kind_phys) :: bdvic_table(max_soiltyp) !VIC model infiltration parameter (-)
118 real (kind=kind_phys) :: gdvic_table(max_soiltyp) !mean capilary drive (m)
119 real (kind=kind_phys) :: bbvic_table(max_soiltyp) !heterogeniety parameter for DVIC infiltration [-]
120
121! genparm.tbl parameters
122
123 real (kind=kind_phys) :: slope_table(num_slope)
124
125 real (kind=kind_phys) :: csoil_table
126 real (kind=kind_phys) :: refdk_table
127 real (kind=kind_phys) :: refkdt_table
128 real (kind=kind_phys) :: frzk_table
129 real (kind=kind_phys) :: zbot_table
130 real (kind=kind_phys) :: czil_table
131
132! mptable.tbl radiation parameters
133
134 real (kind=kind_phys) :: albsat_table(msc,mband)
135 real (kind=kind_phys) :: albdry_table(msc,mband)
136 real (kind=kind_phys) :: albice_table(mband)
137 real (kind=kind_phys) :: alblak_table(mband)
138 real (kind=kind_phys) :: omegas_table(mband)
139 real (kind=kind_phys) :: betads_table
140 real (kind=kind_phys) :: betais_table
141 real (kind=kind_phys) :: eg_table(2)
142
143! mptable.tbl global parameters
144
145 real (kind=kind_phys) :: co2_table
146 real (kind=kind_phys) :: o2_table
147 real (kind=kind_phys) :: timean_table
148 real (kind=kind_phys) :: fsatmx_table
149 real (kind=kind_phys) :: z0sno_table
150 real (kind=kind_phys) :: ssi_table
151 real (kind=kind_phys) :: snow_ret_fac_table
152 real (kind=kind_phys) :: snow_emis_table
153 real (kind=kind_phys) :: swemx_table
154 real (kind=kind_phys) :: tau0_table
155 real (kind=kind_phys) :: grain_growth_table
156 real (kind=kind_phys) :: extra_growth_table
157 real (kind=kind_phys) :: dirt_soot_table
158 real (kind=kind_phys) :: bats_cosz_table
159 real (kind=kind_phys) :: bats_vis_new_table
160 real (kind=kind_phys) :: bats_nir_new_table
161 real (kind=kind_phys) :: bats_vis_age_table
162 real (kind=kind_phys) :: bats_nir_age_table
163 real (kind=kind_phys) :: bats_vis_dir_table
164 real (kind=kind_phys) :: bats_nir_dir_table
165 real (kind=kind_phys) :: rsurf_snow_table
166 real (kind=kind_phys) :: rsurf_exp_table
167
168! mptable.tbl irrigation parameters
169
170 real (kind=kind_phys) :: irr_frac_table ! irrigation Fraction
171 integer :: irr_har_table ! number of days before harvest date to stop irrigation
172 real (kind=kind_phys) :: irr_lai_table ! Minimum lai to trigger irrigation
173 real (kind=kind_phys) :: irr_mad_table ! management allowable deficit (0-1)
174 real (kind=kind_phys) :: filoss_table ! fraction of flood irrigation loss (0-1)
175 real (kind=kind_phys) :: sprir_rate_table ! mm/h, sprinkler irrigation rate
176 real (kind=kind_phys) :: micir_rate_table ! mm/h, micro irrigation rate
177 real (kind=kind_phys) :: firtfac_table ! flood application rate factor
178 real (kind=kind_phys) :: ir_rain_table ! maximum precipitation to stop irrigation trigger
179
180! mptable.tbl crop parameters
181
182 integer :: default_crop_table ! Default crop index
183 integer :: pltday_table(ncrop)
184 integer :: hsday_table(ncrop)
185 real (kind=kind_phys) :: plantpop_table(ncrop)
186 real (kind=kind_phys) :: irri_table(ncrop)
187
188 real (kind=kind_phys) :: gddtbase_table(ncrop)
189 real (kind=kind_phys) :: gddtcut_table(ncrop)
190 real (kind=kind_phys) :: gdds1_table(ncrop)
191 real (kind=kind_phys) :: gdds2_table(ncrop)
192 real (kind=kind_phys) :: gdds3_table(ncrop)
193 real (kind=kind_phys) :: gdds4_table(ncrop)
194 real (kind=kind_phys) :: gdds5_table(ncrop)
195
196 real (kind=kind_phys) :: c3psni_table(ncrop) !photosynthetic pathway: 0. = c4, 1. = c3 ! Zhe Zhang 2020-07-03
197 real (kind=kind_phys) :: kc25i_table(ncrop) !co2 michaelis-menten constant at 25c (pa)
198 real (kind=kind_phys) :: akci_table(ncrop) !q10 for kc25
199 real (kind=kind_phys) :: ko25i_table(ncrop) !o2 michaelis-menten constant at 25c (pa)
200 real (kind=kind_phys) :: akoi_table(ncrop) !q10 for ko25
201 real (kind=kind_phys) :: vcmx25i_table(ncrop) !maximum rate of carboxylation at 25c (umol co2/m**2/s)
202 real (kind=kind_phys) :: avcmxi_table(ncrop) !q10 for vcmx25
203 real (kind=kind_phys) :: bpi_table(ncrop) !minimum leaf conductance (umol/m**2/s)
204 real (kind=kind_phys) :: mpi_table(ncrop) !slope of conductance-to-photosynthesis relationship
205 real (kind=kind_phys) :: qe25i_table(ncrop) !quantum efficiency at 25c (umol co2 / umol photon)
206 real (kind=kind_phys) :: folnmxi_table(ncrop) !foliage nitrogen concentration when
207
208 integer :: c3c4_table(ncrop)
209 real (kind=kind_phys) :: aref_table(ncrop)
210 real (kind=kind_phys) :: psnrf_table(ncrop)
211 real (kind=kind_phys) :: i2par_table(ncrop)
212 real (kind=kind_phys) :: tassim0_table(ncrop)
213 real (kind=kind_phys) :: tassim1_table(ncrop)
214 real (kind=kind_phys) :: tassim2_table(ncrop)
215 real (kind=kind_phys) :: k_table(ncrop)
216 real (kind=kind_phys) :: epsi_table(ncrop)
217
218 real (kind=kind_phys) :: q10mr_table(ncrop)
219 real (kind=kind_phys) :: foln_mx_table(ncrop)
220 real (kind=kind_phys) :: lefreez_table(ncrop)
221
222 real (kind=kind_phys) :: dile_fc_table(ncrop,nstage)
223 real (kind=kind_phys) :: dile_fw_table(ncrop,nstage)
224 real (kind=kind_phys) :: fra_gr_table(ncrop)
225
226 real (kind=kind_phys) :: lf_ovrc_table(ncrop,nstage)
227 real (kind=kind_phys) :: st_ovrc_table(ncrop,nstage)
228 real (kind=kind_phys) :: rt_ovrc_table(ncrop,nstage)
229 real (kind=kind_phys) :: lfmr25_table(ncrop)
230 real (kind=kind_phys) :: stmr25_table(ncrop)
231 real (kind=kind_phys) :: rtmr25_table(ncrop)
232 real (kind=kind_phys) :: grainmr25_table(ncrop)
233
234 real (kind=kind_phys) :: lfpt_table(ncrop,nstage)
235 real (kind=kind_phys) :: stpt_table(ncrop,nstage)
236 real (kind=kind_phys) :: rtpt_table(ncrop,nstage)
237 real (kind=kind_phys) :: grainpt_table(ncrop,nstage)
238 real (kind=kind_phys) :: lfct_table(ncrop,nstage) ! fraction of carbohydrate translocation from leaf to grain ! Zhe Zhang 2020-07-13
239 real (kind=kind_phys) :: stct_table(ncrop,nstage) ! stem to grain
240 real (kind=kind_phys) :: rtct_table(ncrop,nstage) ! root to grain
241 real (kind=kind_phys) :: bio2lai_table(ncrop)
242
243! tile drainage parameters
244 real (kind=kind_phys) :: tdsmc_fac_table(max_soiltyp)
245 real (kind=kind_phys) :: td_dc_table(max_soiltyp)
246 integer :: td_depth_table(max_soiltyp)
247 integer :: drain_layer_opt_table
248 real (kind=kind_phys) :: td_dcoef_table(max_soiltyp)
249 real (kind=kind_phys) :: td_d_table(max_soiltyp)
250 real (kind=kind_phys) :: td_adepth_table(max_soiltyp)
251 real (kind=kind_phys) :: td_radi_table(max_soiltyp)
252 real (kind=kind_phys) :: td_spac_table(max_soiltyp)
253 real (kind=kind_phys) :: td_ddrain_table(max_soiltyp)
254 real (kind=kind_phys) :: klat_fac_table(max_soiltyp)
255
256! mptable.tbl optional parameters
257
258 !------------------------------------------------------------------------------
259 ! Saxton and Rawls 2006 Pedo-transfer function coefficients
260 !------------------------------------------------------------------------------
261
262 real (kind=kind_phys) :: sr2006_theta_1500t_a
263 real (kind=kind_phys) :: sr2006_theta_1500t_b
264 real (kind=kind_phys) :: sr2006_theta_1500t_c
265 real (kind=kind_phys) :: sr2006_theta_1500t_d
266 real (kind=kind_phys) :: sr2006_theta_1500t_e
267 real (kind=kind_phys) :: sr2006_theta_1500t_f
268 real (kind=kind_phys) :: sr2006_theta_1500t_g
269
270 real (kind=kind_phys) :: sr2006_theta_1500_a
271 real (kind=kind_phys) :: sr2006_theta_1500_b
272
273 real (kind=kind_phys) :: sr2006_theta_33t_a
274 real (kind=kind_phys) :: sr2006_theta_33t_b
275 real (kind=kind_phys) :: sr2006_theta_33t_c
276 real (kind=kind_phys) :: sr2006_theta_33t_d
277 real (kind=kind_phys) :: sr2006_theta_33t_e
278 real (kind=kind_phys) :: sr2006_theta_33t_f
279 real (kind=kind_phys) :: sr2006_theta_33t_g
280
281 real (kind=kind_phys) :: sr2006_theta_33_a
282 real (kind=kind_phys) :: sr2006_theta_33_b
283 real (kind=kind_phys) :: sr2006_theta_33_c
284
285 real (kind=kind_phys) :: sr2006_theta_s33t_a
286 real (kind=kind_phys) :: sr2006_theta_s33t_b
287 real (kind=kind_phys) :: sr2006_theta_s33t_c
288 real (kind=kind_phys) :: sr2006_theta_s33t_d
289 real (kind=kind_phys) :: sr2006_theta_s33t_e
290 real (kind=kind_phys) :: sr2006_theta_s33t_f
291 real (kind=kind_phys) :: sr2006_theta_s33t_g
292
293 real (kind=kind_phys) :: sr2006_theta_s33_a
294 real (kind=kind_phys) :: sr2006_theta_s33_b
295
296 real (kind=kind_phys) :: sr2006_psi_et_a
297 real (kind=kind_phys) :: sr2006_psi_et_b
298 real (kind=kind_phys) :: sr2006_psi_et_c
299 real (kind=kind_phys) :: sr2006_psi_et_d
300 real (kind=kind_phys) :: sr2006_psi_et_e
301 real (kind=kind_phys) :: sr2006_psi_et_f
302 real (kind=kind_phys) :: sr2006_psi_et_g
303
304 real (kind=kind_phys) :: sr2006_psi_e_a
305 real (kind=kind_phys) :: sr2006_psi_e_b
306 real (kind=kind_phys) :: sr2006_psi_e_c
307
308 real (kind=kind_phys) :: sr2006_smcmax_a
309 real (kind=kind_phys) :: sr2006_smcmax_b
310
311contains
312
313 subroutine read_mp_table_parameters(errmsg, errflg)
314 implicit none
315
316 character(len=*), intent(out) :: errmsg
317 integer, intent(out) :: errflg
318
319 ! vegetation parameters
320 character(len=256) :: dataset_identifier
321 character(len=256) :: veg_dataset_description
322 logical :: file_named
323 integer :: ierr, ik, im
324 integer :: nveg, isurban, iswater, isbarren, isice, iscrop, eblforest, natural
325 integer :: lcz_1, lcz_2, lcz_3, lcz_4, lcz_5, lcz_6, lcz_7, lcz_8, lcz_9, lcz_10, lcz_11
326 real (kind=kind_phys), dimension(mvt) :: sai_jan, sai_feb, sai_mar, sai_apr, sai_may, sai_jun, sai_jul, sai_aug, &
327 sai_sep, sai_oct, sai_nov, sai_dec, lai_jan, lai_feb, lai_mar, lai_apr, &
328 lai_may, lai_jun, lai_jul, lai_aug, lai_sep, lai_oct, lai_nov, lai_dec, &
329 rhol_vis, rhol_nir, rhos_vis, rhos_nir, taul_vis, taul_nir, taus_vis, taus_nir,&
330 ch2op, dleaf, z0mvt, hvt, hvb, z0mhvt, &
331 den, rc, mfsno, scffac, cbiom, xl, cwpvt, c3psn, kc25, &
332 akc, ko25, ako, avcmx, aqe, ltovrc, dilefc, dilefw, rmf25, sla, fragr, tmin, &
333 vcmx25, tdlef, bp, mp, qe25, rms25, rmr25, arm, folnmx, wdpool, wrrat, mrp, &
334 nroot, rgl, rs, hs, topt, rsmax, rtovrc, rswoodc, bf, wstrc, laimin, &
335 xsamin, eps1, eps2, eps3, eps4, eps5
336 namelist / noahmp_usgs_veg_categories / veg_dataset_description, nveg
337 namelist / noahmp_usgs_parameters / isurban, iswater, isbarren, isice, iscrop, eblforest, natural, &
338 lcz_1, lcz_2, lcz_3, lcz_4, lcz_5, lcz_6, lcz_7, lcz_8, lcz_9, lcz_10, lcz_11, &
339 ch2op, dleaf, z0mvt, hvt, hvb, z0mhvt, &
340 den, rc, mfsno, scffac, cbiom, xl, cwpvt, c3psn, kc25, &
341 akc, ko25, ako, avcmx, aqe, ltovrc, dilefc, dilefw, rmf25, sla, fragr, tmin, &
342 vcmx25, tdlef, bp, mp, qe25, rms25, rmr25, arm, folnmx, wdpool, wrrat, mrp, &
343 nroot, rgl, rs, hs, topt, rsmax, rtovrc, rswoodc, bf, wstrc, laimin, &
344 xsamin, sai_jan, sai_feb, sai_mar, sai_apr, sai_may, &
345 sai_jun, sai_jul, sai_aug, sai_sep, sai_oct, sai_nov, sai_dec, lai_jan, &
346 lai_feb, lai_mar, lai_apr, lai_may, lai_jun, lai_jul, lai_aug, lai_sep, &
347 lai_oct, lai_nov, lai_dec, rhol_vis, rhol_nir, rhos_vis, rhos_nir, taul_vis, &
348 taul_nir, taus_vis, taus_nir, eps1, eps2, eps3, eps4, eps5
349 namelist / noahmp_modis_veg_categories / veg_dataset_description, nveg
350 namelist / noahmp_modis_parameters / isurban, iswater, isbarren, isice, iscrop, eblforest, natural, &
351 lcz_1, lcz_2, lcz_3, lcz_4, lcz_5, lcz_6, lcz_7, lcz_8, lcz_9, lcz_10, lcz_11, &
352 ch2op, dleaf, z0mvt, hvt, hvb, z0mhvt, &
353 den, rc, mfsno, scffac, cbiom, xl, cwpvt, c3psn, kc25, &
354 akc, ko25, ako, avcmx, aqe, ltovrc, dilefc, dilefw, rmf25, sla, fragr, tmin, &
355 vcmx25, tdlef, bp, mp, qe25, rms25, rmr25, arm, folnmx, wdpool, wrrat, mrp, &
356 nroot, rgl, rs, hs, topt, rsmax, rtovrc, rswoodc, bf, wstrc, laimin, &
357 xsamin, sai_jan, sai_feb, sai_mar, sai_apr, sai_may, &
358 sai_jun, sai_jul, sai_aug, sai_sep, sai_oct, sai_nov, sai_dec, lai_jan, &
359 lai_feb, lai_mar, lai_apr, lai_may, lai_jun, lai_jul, lai_aug, lai_sep, &
360 lai_oct, lai_nov, lai_dec, rhol_vis, rhol_nir, rhos_vis, rhos_nir, taul_vis, &
361 taul_nir, taus_vis, taus_nir, eps1, eps2, eps3, eps4, eps5
362 ! soil parameters
363 character(len=256) :: message
364 character(len=10) :: sltype
365 integer :: slcats
366 real (kind=kind_phys), dimension(max_soiltyp) :: bb, drysmc, maxsmc, refsmc, satpsi, satdk, satdw, wltsmc, qtz, &
367 bvic, axaj, bxaj, xxaj, bdvic, bbvic, gdvic, hc
368 namelist / noahmp_stas_soil_categories / sltype, slcats
369 namelist / noahmp_soil_stas_parameters / bb, drysmc, maxsmc, refsmc, satpsi, satdk, satdw, wltsmc, qtz, &
370 bvic, axaj, bxaj, xxaj, bdvic, bbvic, gdvic
371 namelist / noahmp_soil_stas_ruc_parameters / bb, drysmc, hc, maxsmc, refsmc, satpsi, satdk, satdw, wltsmc, qtz, &
372 bvic, axaj, bxaj, xxaj, bdvic, bbvic, gdvic
373
374 ! general parameters
375 real (kind=kind_phys) :: csoil_data, refdk_data, refkdt_data, frzk_data, zbot_data, czil_data
376 real (kind=kind_phys), dimension(num_slope) :: slope_data
377 namelist / noahmp_general_parameters / slope_data, csoil_data, refdk_data, refkdt_data, frzk_data, zbot_data, &
378 czil_data
379
380 ! radiation parameters
381 real (kind=kind_phys) :: betads, betais, eice
382 real (kind=kind_phys), dimension(mband) :: albice, alblak, omegas
383 real (kind=kind_phys), dimension(2) :: eg
384 real (kind=kind_phys), dimension(msc) :: albsat_vis, albsat_nir, albdry_vis, albdry_nir
385 namelist / noahmp_rad_parameters / albsat_vis, albsat_nir, albdry_vis, albdry_nir, albice, alblak, omegas, &
386 betads, betais, eg, eice
387
388 ! global parameters
389 real (kind=kind_phys) :: co2, o2, timean, fsatmx, z0sno, ssi, snow_ret_fac ,snow_emis, swemx, tau0, &
390 grain_growth, extra_growth, dirt_soot, bats_cosz, bats_vis_new, &
391 bats_nir_new, bats_vis_age, bats_nir_age, bats_vis_dir, bats_nir_dir, &
392 rsurf_snow, rsurf_exp, c2_snowcompact, c3_snowcompact, c4_snowcompact, &
393 c5_snowcompact, dm_snowcompact, eta0_snowcompact, snliqmaxfrac, swemaxgla, &
394 wslmax, rous, cmic, snowden_max, class_alb_ref, class_sno_age, class_alb_new,&
395 psiwlt, z0soil, z0lake
396 namelist / noahmp_global_parameters / co2, o2, timean, fsatmx, z0sno, ssi, snow_ret_fac ,snow_emis, swemx, tau0, &
397 grain_growth, extra_growth, dirt_soot, bats_cosz, bats_vis_new, &
398 bats_nir_new, bats_vis_age, bats_nir_age, bats_vis_dir, bats_nir_dir, &
399 rsurf_snow, rsurf_exp, c2_snowcompact, c3_snowcompact, c4_snowcompact, &
400 c5_snowcompact, dm_snowcompact, eta0_snowcompact, snliqmaxfrac, swemaxgla, &
401 wslmax, rous, cmic, snowden_max, class_alb_ref, class_sno_age, class_alb_new,&
402 psiwlt, z0soil, z0lake
403
404 ! irrigation parameters
405 integer :: irr_har
406 real (kind=kind_phys) :: irr_frac, irr_lai, irr_mad, filoss, sprir_rate, micir_rate, firtfac, ir_rain
407 namelist / noahmp_irrigation_parameters / irr_frac, irr_har, irr_lai, irr_mad, filoss, sprir_rate, micir_rate, firtfac,&
408 ir_rain
409
410 ! crop parameters
411 integer :: default_crop
412 integer , dimension(ncrop) :: pltday, hsday
413 real (kind=kind_phys), dimension(ncrop) :: plantpop, irri, gddtbase, gddtcut, gdds1, gdds2, gdds3, gdds4, gdds5, c3psni,&
414 kc25i, akci, ko25i, akoi, avcmxi, vcmx25i, bpi, mpi, folnmxi, qe25i, aref, &
415 psnrf, i2par, tassim0, tassim1, tassim2, k, epsi, q10mr, lefreez, &
416 dile_fc_s1, dile_fc_s2, dile_fc_s3, dile_fc_s4, dile_fc_s5, dile_fc_s6, &
417 dile_fc_s7, dile_fc_s8, dile_fw_s1, dile_fw_s2, dile_fw_s3, dile_fw_s4, &
418 dile_fw_s5, dile_fw_s6, dile_fw_s7, dile_fw_s8, fra_gr, lf_ovrc_s1, &
419 lf_ovrc_s2, lf_ovrc_s3, lf_ovrc_s4, lf_ovrc_s5, lf_ovrc_s6, lf_ovrc_s7, &
420 lf_ovrc_s8, st_ovrc_s1, st_ovrc_s2, st_ovrc_s3, st_ovrc_s4, st_ovrc_s5, &
421 st_ovrc_s6, st_ovrc_s7, st_ovrc_s8, rt_ovrc_s1, rt_ovrc_s2, rt_ovrc_s3, &
422 rt_ovrc_s4, rt_ovrc_s5, rt_ovrc_s6, rt_ovrc_s7, rt_ovrc_s8, lfmr25, stmr25, &
423 rtmr25, grainmr25, lfpt_s1, lfpt_s2, lfpt_s3, lfpt_s4, lfpt_s5, lfpt_s6, &
424 lfpt_s7, lfpt_s8, stpt_s1, stpt_s2, stpt_s3, stpt_s4, stpt_s5, stpt_s6, &
425 stpt_s7, stpt_s8, rtpt_s1, rtpt_s2, rtpt_s3, rtpt_s4, rtpt_s5, rtpt_s6, &
426 rtpt_s7, rtpt_s8, grainpt_s1, grainpt_s2, grainpt_s3, grainpt_s4, grainpt_s5,&
427 grainpt_s6, grainpt_s7, grainpt_s8, lfct_s1, lfct_s2, lfct_s3, lfct_s4, &
428 lfct_s5, lfct_s6, lfct_s7, lfct_s8, stct_s1, stct_s2, stct_s3, stct_s4, &
429 stct_s5, stct_s6, stct_s7, stct_s8, rtct_s1, rtct_s2, rtct_s3, rtct_s4, &
430 rtct_s5, rtct_s6, rtct_s7, rtct_s8, bio2lai
431 namelist / noahmp_crop_parameters / default_crop, pltday, hsday, plantpop, irri, gddtbase, gddtcut, gdds1, gdds2,&
432 gdds3, gdds4, gdds5, c3psni, kc25i, akci, ko25i, akoi, avcmxi, vcmx25i, bpi, &
433 mpi, folnmxi, qe25i, aref, psnrf, i2par, tassim0, tassim1, tassim2, k, &
434 epsi,q10mr, lefreez, dile_fc_s1, dile_fc_s2, dile_fc_s3, dile_fc_s4, &
435 dile_fc_s5, dile_fc_s6, dile_fc_s7, dile_fc_s8, dile_fw_s1, dile_fw_s2, &
436 dile_fw_s3, dile_fw_s4, dile_fw_s5, dile_fw_s6, dile_fw_s7, dile_fw_s8, &
437 fra_gr, lf_ovrc_s1, lf_ovrc_s2, lf_ovrc_s3, lf_ovrc_s4, lf_ovrc_s5, &
438 lf_ovrc_s6, lf_ovrc_s7, lf_ovrc_s8, st_ovrc_s1, st_ovrc_s2, st_ovrc_s3, &
439 st_ovrc_s4, st_ovrc_s5, st_ovrc_s6, st_ovrc_s7, st_ovrc_s8, rt_ovrc_s1, &
440 rt_ovrc_s2, rt_ovrc_s3, rt_ovrc_s4, rt_ovrc_s5, rt_ovrc_s6, rt_ovrc_s7, &
441 rt_ovrc_s8, lfmr25, stmr25, rtmr25, grainmr25, lfpt_s1, lfpt_s2, lfpt_s3, &
442 lfpt_s4, lfpt_s5, lfpt_s6, lfpt_s7, lfpt_s8, stpt_s1, stpt_s2, stpt_s3, &
443 stpt_s4, stpt_s5, stpt_s6, stpt_s7, stpt_s8, rtpt_s1, rtpt_s2, rtpt_s3, &
444 rtpt_s4, rtpt_s5, rtpt_s6, rtpt_s7, rtpt_s8, grainpt_s1, grainpt_s2, &
445 grainpt_s3, grainpt_s4, grainpt_s5, grainpt_s6, grainpt_s7, grainpt_s8, &
446 lfct_s1, lfct_s2, lfct_s3, lfct_s4, lfct_s5, lfct_s6, lfct_s7, lfct_s8, &
447 stct_s1, stct_s2, stct_s3, stct_s4, stct_s5, stct_s6, stct_s7, stct_s8, &
448 rtct_s1, rtct_s2, rtct_s3, rtct_s4, rtct_s5, rtct_s6, rtct_s7, rtct_s8, &
449 bio2lai
450
451 ! tile drainage parameters
452 integer :: nsoiltype, drain_layer_opt
453 integer , dimension(max_soiltyp) :: td_depth
454 real (kind=kind_phys), dimension(max_soiltyp) :: tdsmc_fac, td_dc, td_dcoef, td_d, td_adepth, td_radi, td_spac, &
455 td_ddrain, klat_fac
456 namelist / noahmp_tiledrain_parameters / nsoiltype, drain_layer_opt, tdsmc_fac, td_depth, td_dc, td_dcoef, td_d,&
457 td_adepth, td_radi, td_spac, td_ddrain, klat_fac
458
459 ! optional parameters
460 real (kind=kind_phys) :: sr2006_theta_1500t_a, sr2006_theta_1500t_b, sr2006_theta_1500t_c, &
461 sr2006_theta_1500t_d, sr2006_theta_1500t_e, sr2006_theta_1500t_f, &
462 sr2006_theta_1500t_g, sr2006_theta_1500_a , sr2006_theta_1500_b, &
463 sr2006_theta_33t_a, sr2006_theta_33t_b, sr2006_theta_33t_c, &
464 sr2006_theta_33t_d, sr2006_theta_33t_e, sr2006_theta_33t_f, &
465 sr2006_theta_33t_g, sr2006_theta_33_a, sr2006_theta_33_b, &
466 sr2006_theta_33_c, sr2006_theta_s33t_a, sr2006_theta_s33t_b, &
467 sr2006_theta_s33t_c, sr2006_theta_s33t_d, sr2006_theta_s33t_e, &
468 sr2006_theta_s33t_f, sr2006_theta_s33t_g, sr2006_theta_s33_a, &
469 sr2006_theta_s33_b, sr2006_psi_et_a, sr2006_psi_et_b, sr2006_psi_et_c, &
470 sr2006_psi_et_d, sr2006_psi_et_e, sr2006_psi_et_f, sr2006_psi_et_g, &
471 sr2006_psi_e_a, sr2006_psi_e_b, sr2006_psi_e_c, sr2006_smcmax_a, &
472 sr2006_smcmax_b
473 namelist / noahmp_optional_parameters / sr2006_theta_1500t_a, sr2006_theta_1500t_b, sr2006_theta_1500t_c, &
474 sr2006_theta_1500t_d, sr2006_theta_1500t_e, sr2006_theta_1500t_f, &
475 sr2006_theta_1500t_g, sr2006_theta_1500_a, sr2006_theta_1500_b, &
476 sr2006_theta_33t_a, sr2006_theta_33t_b, sr2006_theta_33t_c, &
477 sr2006_theta_33t_d, sr2006_theta_33t_e, sr2006_theta_33t_f, &
478 sr2006_theta_33t_g, sr2006_theta_33_a, sr2006_theta_33_b, &
479 sr2006_theta_33_c, sr2006_theta_s33t_a, sr2006_theta_s33t_b, &
480 sr2006_theta_s33t_c, sr2006_theta_s33t_d, sr2006_theta_s33t_e, &
481 sr2006_theta_s33t_f, sr2006_theta_s33t_g, sr2006_theta_s33_a, &
482 sr2006_theta_s33_b, sr2006_psi_et_a, sr2006_psi_et_b, sr2006_psi_et_c, &
483 sr2006_psi_et_d, sr2006_psi_et_e, sr2006_psi_et_f, sr2006_psi_et_g, &
484 sr2006_psi_e_a, sr2006_psi_e_b, sr2006_psi_e_c, sr2006_smcmax_a, &
485 sr2006_smcmax_b
486
487 errmsg = ''
488 errflg = 0
489
490 ! initialize our variables to bad values, so that if the namelist read fails, we come to a screeching halt as soon as we try to use anything.
491 ! vegetation parameters
492 isurban_table = -99999
493 iswater_table = -99999
494 isbarren_table = -99999
495 isice_table = -99999
496 iscrop_table = -99999
497 eblforest_table = -99999
498 natural_table = -99999
499 lcz_1_table = -99999
500 lcz_2_table = -99999
501 lcz_3_table = -99999
502 lcz_4_table = -99999
503 lcz_5_table = -99999
504 lcz_6_table = -99999
505 lcz_7_table = -99999
506 lcz_8_table = -99999
507 lcz_9_table = -99999
508 lcz_10_table = -99999
509 lcz_11_table = -99999
510 ch2op_table = -1.0e36
511 dleaf_table = -1.0e36
512 z0mvt_table = -1.0e36
513 hvt_table = -1.0e36
514 hvb_table = -1.0e36
515 z0mhvt_table = -1.0e36
516 den_table = -1.0e36
517 rc_table = -1.0e36
518 mfsno_table = -1.0e36
519 scffac_table = -1.0e36
520 cbiom_table = -1.0e36
521 rhol_table = -1.0e36
522 rhos_table = -1.0e36
523 taul_table = -1.0e36
524 taus_table = -1.0e36
525 xl_table = -1.0e36
526 cwpvt_table = -1.0e36
527 c3psn_table = -1.0e36
528 kc25_table = -1.0e36
529 akc_table = -1.0e36
530 ko25_table = -1.0e36
531 ako_table = -1.0e36
532 avcmx_table = -1.0e36
533 aqe_table = -1.0e36
534 ltovrc_table = -1.0e36
535 dilefc_table = -1.0e36
536 dilefw_table = -1.0e36
537 rmf25_table = -1.0e36
538 sla_table = -1.0e36
539 fragr_table = -1.0e36
540 tmin_table = -1.0e36
541 vcmx25_table = -1.0e36
542 tdlef_table = -1.0e36
543 bp_table = -1.0e36
544 mp_table = -1.0e36
545 qe25_table = -1.0e36
546 rms25_table = -1.0e36
547 rmr25_table = -1.0e36
548 arm_table = -1.0e36
549 folnmx_table = -1.0e36
550 wdpool_table = -1.0e36
551 wrrat_table = -1.0e36
552 mrp_table = -1.0e36
553 saim_table = -1.0e36
554 laim_table = -1.0e36
555 nroot_table = -1.0e36
556 rgl_table = -1.0e36
557 rs_table = -1.0e36
558 hs_table = -1.0e36
559 topt_table = -1.0e36
560 rsmax_table = -1.0e36
561 ! not used in the current ufs version
562! rtovrc_table = -1.0e36
563! rswoodc_table = -1.0e36
564! bf_table = -1.0e36
565! wstrc_table = -1.0e36
566! laimin_table = -1.0e36
567! xsamin_table = -1.0e36
568
569 ! soil parameters
570
571 bexp_table = -1.0e36
572 smcdry_table = -1.0e36
573 f1_table = -1.0e36
574 smcmax_table = -1.0e36
575 smcref_table = -1.0e36
576 psisat_table = -1.0e36
577 dksat_table = -1.0e36
578 dwsat_table = -1.0e36
579 smcwlt_table = -1.0e36
580 quartz_table = -1.0e36
581 bvic_table = -1.0e36
582 axaj_table = -1.0e36
583 bxaj_table = -1.0e36
584 xxaj_table = -1.0e36
585 bdvic_table = -1.0e36
586 gdvic_table = -1.0e36
587 bbvic_table = -1.0e36
588
589 ! general parameters
590 slope_table = -1.0e36
591 csoil_table = -1.0e36
592 refdk_table = -1.0e36
593 refkdt_table = -1.0e36
594 frzk_table = -1.0e36
595 zbot_table = -1.0e36
596 czil_table = -1.0e36
597
598 ! radiation parameters
599 albsat_table = -1.0e36
600 albdry_table = -1.0e36
601 albice_table = -1.0e36
602 alblak_table = -1.0e36
603 omegas_table = -1.0e36
604 betads_table = -1.0e36
605 betais_table = -1.0e36
606 eg_table = -1.0e36
607! eice_table = -1.0e36
608
609 ! global parameters
610 co2_table = -1.0e36
611 o2_table = -1.0e36
612 timean_table = -1.0e36
613 fsatmx_table = -1.0e36
614 z0sno_table = -1.0e36
615 ssi_table = -1.0e36
616snow_ret_fac_table = -1.0e36
617 snow_emis_table = -1.0e36
618 swemx_table = -1.0e36
619 tau0_table = -1.0e36
620grain_growth_table = -1.0e36
621extra_growth_table = -1.0e36
622 dirt_soot_table = -1.0e36
623 bats_cosz_table = -1.0e36
624bats_vis_new_table = -1.0e36
625bats_nir_new_table = -1.0e36
626bats_vis_age_table = -1.0e36
627bats_nir_age_table = -1.0e36
628bats_vis_dir_table = -1.0e36
629bats_nir_dir_table = -1.0e36
630rsurf_snow_table = -1.0e36
631 rsurf_exp_table = -1.0e36
632
633! c2_snowcompact_table = -1.0e36
634! c3_snowcompact_table = -1.0e36
635! c4_snowcompact_table = -1.0e36
636! c5_snowcompact_table = -1.0e36
637! dm_snowcompact_table = -1.0e36
638! eta0_snowcompact_table = -1.0e36
639! snliqmaxfrac_table = -1.0e36
640! swemaxgla_table = -1.0e36
641! wslmax_table = -1.0e36
642! rous_table = -1.0e36
643! cmic_table = -1.0e36
644! snowden_max_table = -1.0e36
645! class_alb_ref_table = -1.0e36
646! class_sno_age_table = -1.0e36
647! class_alb_new_table = -1.0e36
648! psiwlt_table = -1.0e36
649! z0soil_table = -1.0e36
650! z0lake_table = -1.0e36
651
652 ! irrigation parameters
653 irr_har_table = -99999 ! number of days before harvest date to stop irrigation
654 irr_frac_table = -1.0e36 ! irrigation fraction
655 irr_lai_table = -1.0e36 ! minimum lai to trigger irrigation
656 irr_mad_table = -1.0e36 ! management allowable deficit (0-1)
657 filoss_table = -1.0e36 ! fraction of flood irrigation loss (0-1)
658 sprir_rate_table = -1.0e36 ! mm/h, sprinkler irrigation rate
659 micir_rate_table = -1.0e36 ! mm/h, micro irrigation rate
660 firtfac_table = -1.0e36 ! flood application rate factor
661 ir_rain_table = -1.0e36 ! maximum precipitation to stop irrigation trigger
662
663 ! crop parameters
664 default_crop_table = -99999
665 pltday_table = -99999
666 hsday_table = -99999
667 plantpop_table = -1.0e36
668 irri_table = -1.0e36
669 gddtbase_table = -1.0e36
670 gddtcut_table = -1.0e36
671 gdds1_table = -1.0e36
672 gdds2_table = -1.0e36
673 gdds3_table = -1.0e36
674 gdds4_table = -1.0e36
675 gdds5_table = -1.0e36
676 c3psni_table = -1.0e36 ! parameter from psn copied from stomata ! zhe zhang 2020-07-13
677 kc25i_table = -1.0e36
678 akci_table = -1.0e36
679 ko25i_table = -1.0e36
680 akoi_table = -1.0e36
681 avcmxi_table = -1.0e36
682 vcmx25i_table = -1.0e36
683 bpi_table = -1.0e36
684 mpi_table = -1.0e36
685 folnmxi_table = -1.0e36
686 qe25i_table = -1.0e36 ! ends here
687!??? c3c4_table = -99999
688 aref_table = -1.0e36
689 psnrf_table = -1.0e36
690 i2par_table = -1.0e36
691 tassim0_table = -1.0e36
692 tassim1_table = -1.0e36
693 tassim2_table = -1.0e36
694 k_table = -1.0e36
695 epsi_table = -1.0e36
696 q10mr_table = -1.0e36
697 foln_mx_table = -1.0e36
698 lefreez_table = -1.0e36
699 dile_fc_table = -1.0e36
700 dile_fw_table = -1.0e36
701 fra_gr_table = -1.0e36
702 lf_ovrc_table = -1.0e36
703 st_ovrc_table = -1.0e36
704 rt_ovrc_table = -1.0e36
705 lfmr25_table = -1.0e36
706 stmr25_table = -1.0e36
707 rtmr25_table = -1.0e36
708 grainmr25_table = -1.0e36
709 lfpt_table = -1.0e36
710 stpt_table = -1.0e36
711 rtpt_table = -1.0e36
712 grainpt_table = -1.0e36
713 lfct_table = -1.0e36 ! convert start
714 stct_table = -1.0e36
715 rtct_table = -1.0e36 ! convert end
716 bio2lai_table = -1.0e36
717
718 ! tile drainage parameters
719
720 drain_layer_opt_table = -99999
721 td_depth_table = -99999
722 tdsmc_fac_table = -1.0e36
723 td_dc_table = -1.0e36
724 td_dcoef_table = -1.0e36
725 td_d_table = -1.0e36
726 td_adepth_table = -1.0e36
727 td_radi_table = -1.0e36
728 td_spac_table = -1.0e36
729 td_ddrain_table = -1.0e36
730 klat_fac_table = -1.0e36
731
732 ! optional parameters
733! sr2006_theta_1500t_a_table = -1.0e36
734! sr2006_theta_1500t_b_table = -1.0e36
735! sr2006_theta_1500t_c_table = -1.0e36
736! sr2006_theta_1500t_d_table = -1.0e36
737! sr2006_theta_1500t_e_table = -1.0e36
738! sr2006_theta_1500t_f_table = -1.0e36
739! sr2006_theta_1500t_g_table = -1.0e36
740! sr2006_theta_1500_a_table = -1.0e36
741! sr2006_theta_1500_b_table = -1.0e36
742! sr2006_theta_33t_a_table = -1.0e36
743! sr2006_theta_33t_b_table = -1.0e36
744! sr2006_theta_33t_c_table = -1.0e36
745! sr2006_theta_33t_d_table = -1.0e36
746! sr2006_theta_33t_e_table = -1.0e36
747! sr2006_theta_33t_f_table = -1.0e36
748! sr2006_theta_33t_g_table = -1.0e36
749! sr2006_theta_33_a_table = -1.0e36
750! sr2006_theta_33_b_table = -1.0e36
751! sr2006_theta_33_c_table = -1.0e36
752! sr2006_theta_s33t_a_table = -1.0e36
753! sr2006_theta_s33t_b_table = -1.0e36
754! sr2006_theta_s33t_c_table = -1.0e36
755! sr2006_theta_s33t_d_table = -1.0e36
756! sr2006_theta_s33t_e_table = -1.0e36
757! sr2006_theta_s33t_f_table = -1.0e36
758! sr2006_theta_s33t_g_table = -1.0e36
759! sr2006_theta_s33_a_table = -1.0e36
760! sr2006_theta_s33_b_table = -1.0e36
761! sr2006_psi_et_a_table = -1.0e36
762! sr2006_psi_et_b_table = -1.0e36
763! sr2006_psi_et_c_table = -1.0e36
764! sr2006_psi_et_d_table = -1.0e36
765! sr2006_psi_et_e_table = -1.0e36
766! sr2006_psi_et_f_table = -1.0e36
767! sr2006_psi_et_g_table = -1.0e36
768! sr2006_psi_e_a_table = -1.0e36
769! sr2006_psi_e_b_table = -1.0e36
770! sr2006_psi_e_c_table = -1.0e36
771! sr2006_smcmax_a_table = -1.0e36
772! sr2006_smcmax_b_table = -1.0e36
773
774 !---------------------------------------------------------------
775 ! transfer values from table to input variables
776 !---------------------------------------------------------------
777
778 !---------------- noahmptable.tbl vegetation parameters
779
780 dataset_identifier = "modified_igbp_modis_noah"
781
782 inquire( file='noahmptable.tbl', exist=file_named )
783 if ( file_named ) then
784 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
785 else
786 open(15, status='old', form='formatted', action='read', iostat=ierr)
787 end if
788 if ( ierr /= 0 ) then
789 errmsg = 'warning: cannot find file noahmptable.tbl'
790 errflg = 1
791 return
792! write(*,'("warning: cannot find file noahmptable.tbl")')
793 endif
794
795 if ( trim(dataset_identifier) == "usgs" ) then
796 read(15, noahmp_usgs_veg_categories)
797 read(15, noahmp_usgs_parameters)
798 elseif ( trim(dataset_identifier) == "modified_igbp_modis_noah" ) then
799 read(15,noahmp_modis_veg_categories)
800 read(15,noahmp_modis_parameters)
801 else
802 write(*,'("warning: unrecognized dataset_identifier in subroutine readnoahmptable")')
803 write(*,'("warning: dataset_identifier = ''", a, "''")') trim(dataset_identifier)
804 endif
805 close(15)
806
807
808 ! assign values
809 isurban_table = isurban
810 iswater_table = iswater
811 isbarren_table = isbarren
812 isice_table = isice
813 iscrop_table = iscrop
814 eblforest_table = eblforest
815 natural_table = natural
816 lcz_1_table = lcz_1
817 lcz_2_table = lcz_2
818 lcz_3_table = lcz_3
819 lcz_4_table = lcz_4
820 lcz_5_table = lcz_5
821 lcz_6_table = lcz_6
822 lcz_7_table = lcz_7
823 lcz_8_table = lcz_8
824 lcz_9_table = lcz_9
825 lcz_10_table = lcz_10
826 lcz_11_table = lcz_11
827 ch2op_table(1:nveg) = ch2op(1:nveg)
828 dleaf_table(1:nveg) = dleaf(1:nveg)
829 z0mvt_table(1:nveg) = z0mvt(1:nveg)
830 hvt_table(1:nveg) = hvt(1:nveg)
831 hvb_table(1:nveg) = hvb(1:nveg)
832 z0mhvt_table(1:nveg) = z0mhvt(1:nveg)
833 den_table(1:nveg) = den(1:nveg)
834 rc_table(1:nveg) = rc(1:nveg)
835 mfsno_table(1:nveg) = mfsno(1:nveg)
836 scffac_table(1:nveg) = scffac(1:nveg)
837 cbiom_table(1:nveg) = cbiom(1:nveg)
838 xl_table(1:nveg) = xl(1:nveg)
839 cwpvt_table(1:nveg) = cwpvt(1:nveg)
840 c3psn_table(1:nveg) = c3psn(1:nveg)
841 kc25_table(1:nveg) = kc25(1:nveg)
842 akc_table(1:nveg) = akc(1:nveg)
843 ko25_table(1:nveg) = ko25(1:nveg)
844 ako_table(1:nveg) = ako(1:nveg)
845 avcmx_table(1:nveg) = avcmx(1:nveg)
846 aqe_table(1:nveg) = aqe(1:nveg)
847 ltovrc_table(1:nveg) = ltovrc(1:nveg)
848 dilefc_table(1:nveg) = dilefc(1:nveg)
849 dilefw_table(1:nveg) = dilefw(1:nveg)
850 rmf25_table(1:nveg) = rmf25(1:nveg)
851 sla_table(1:nveg) = sla(1:nveg)
852 fragr_table(1:nveg) = fragr(1:nveg)
853 tmin_table(1:nveg) = tmin(1:nveg)
854 vcmx25_table(1:nveg) = vcmx25(1:nveg)
855 tdlef_table(1:nveg) = tdlef(1:nveg)
856 bp_table(1:nveg) = bp(1:nveg)
857 mp_table(1:nveg) = mp(1:nveg)
858 qe25_table(1:nveg) = qe25(1:nveg)
859 rms25_table(1:nveg) = rms25(1:nveg)
860 rmr25_table(1:nveg) = rmr25(1:nveg)
861 arm_table(1:nveg) = arm(1:nveg)
862 folnmx_table(1:nveg) = folnmx(1:nveg)
863 wdpool_table(1:nveg) = wdpool(1:nveg)
864 wrrat_table(1:nveg) = wrrat(1:nveg)
865 mrp_table(1:nveg) = mrp(1:nveg)
866 nroot_table(1:nveg) = nroot(1:nveg)
867 rgl_table(1:nveg) = rgl(1:nveg)
868 rs_table(1:nveg) = rs(1:nveg)
869 hs_table(1:nveg) = hs(1:nveg)
870 topt_table(1:nveg) = topt(1:nveg)
871 rsmax_table(1:nveg) = rsmax(1:nveg)
872! rtovrc_table (1:nveg) = rtovrc (1:nveg)
873! rswoodc_table(1:nveg) = rswoodc(1:nveg)
874! bf_table (1:nveg) = bf (1:nveg)
875! wstrc_table (1:nveg) = wstrc (1:nveg)
876! laimin_table (1:nveg) = laimin (1:nveg)
877! xsamin_table (1:nveg) = xsamin (1:nveg)
878
879 saim_table(1:nveg, 1) = sai_jan(1:nveg)
880 saim_table(1:nveg, 2) = sai_feb(1:nveg)
881 saim_table(1:nveg, 3) = sai_mar(1:nveg)
882 saim_table(1:nveg, 4) = sai_apr(1:nveg)
883 saim_table(1:nveg, 5) = sai_may(1:nveg)
884 saim_table(1:nveg, 6) = sai_jun(1:nveg)
885 saim_table(1:nveg, 7) = sai_jul(1:nveg)
886 saim_table(1:nveg, 8) = sai_aug(1:nveg)
887 saim_table(1:nveg, 9) = sai_sep(1:nveg)
888 saim_table(1:nveg,10) = sai_oct(1:nveg)
889 saim_table(1:nveg,11) = sai_nov(1:nveg)
890 saim_table(1:nveg,12) = sai_dec(1:nveg)
891 laim_table(1:nveg, 1) = lai_jan(1:nveg)
892 laim_table(1:nveg, 2) = lai_feb(1:nveg)
893 laim_table(1:nveg, 3) = lai_mar(1:nveg)
894 laim_table(1:nveg, 4) = lai_apr(1:nveg)
895 laim_table(1:nveg, 5) = lai_may(1:nveg)
896 laim_table(1:nveg, 6) = lai_jun(1:nveg)
897 laim_table(1:nveg, 7) = lai_jul(1:nveg)
898 laim_table(1:nveg, 8) = lai_aug(1:nveg)
899 laim_table(1:nveg, 9) = lai_sep(1:nveg)
900 laim_table(1:nveg,10) = lai_oct(1:nveg)
901 laim_table(1:nveg,11) = lai_nov(1:nveg)
902 laim_table(1:nveg,12) = lai_dec(1:nveg)
903 rhol_table(1:nveg,1) = rhol_vis(1:nveg) !leaf reflectance: 1=vis, 2=nir
904 rhol_table(1:nveg,2) = rhol_nir(1:nveg) !leaf reflectance: 1=vis, 2=nir
905 rhos_table(1:nveg,1) = rhos_vis(1:nveg) !stem reflectance: 1=vis, 2=nir
906 rhos_table(1:nveg,2) = rhos_nir(1:nveg) !stem reflectance: 1=vis, 2=nir
907 taul_table(1:nveg,1) = taul_vis(1:nveg) !leaf transmittance: 1=vis, 2=nir
908 taul_table(1:nveg,2) = taul_nir(1:nveg) !leaf transmittance: 1=vis, 2=nir
909 taus_table(1:nveg,1) = taus_vis(1:nveg) !stem transmittance: 1=vis, 2=nir
910 taus_table(1:nveg,2) = taus_nir(1:nveg) !stem transmittance: 1=vis, 2=nir
911
912 !---------------- noahmptable.tbl soil parameters
913 inquire( file='noahmptable.tbl', exist=file_named )
914 if ( file_named ) then
915 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
916 else
917 open(15, status='old', form='formatted', action='read', iostat=ierr)
918 end if
919 if ( ierr /= 0 ) then
920 errmsg = 'warning: cannot find file noahmptable.tbl'
921 errflg = 1
922 return
923! write(*,'("warning: cannot find file noahmptable.tbl")')
924 endif
925 read(15, noahmp_stas_soil_categories)
926 if ( trim(sltype) == "stas" ) then
927 read(15, noahmp_soil_stas_parameters)
928 elseif ( trim(sltype) == "stas_ruc" ) then
929 read(15, noahmp_soil_stas_ruc_parameters)
930 else
931 write(*,'("warning: unrecognized soiltype in subroutine readnoahmptable")')
932 write(*,'("warning: dataset_identifier = ''", a, "''")') trim(sltype)
933 endif
934 close(15)
935
936 ! assign values
937! slcats_table = slcats
938 bexp_table(1:slcats) = bb(1:slcats)
939 smcdry_table(1:slcats) = drysmc(1:slcats)
940 smcmax_table(1:slcats) = maxsmc(1:slcats)
941 smcref_table(1:slcats) = refsmc(1:slcats)
942 psisat_table(1:slcats) = satpsi(1:slcats)
943 dksat_table(1:slcats) = satdk(1:slcats)
944 dwsat_table(1:slcats) = satdw(1:slcats)
945 smcwlt_table(1:slcats) = wltsmc(1:slcats)
946 quartz_table(1:slcats) = qtz(1:slcats)
947 bvic_table(1:slcats) = bvic(1:slcats)
948 axaj_table(1:slcats) = axaj(1:slcats)
949 bxaj_table(1:slcats) = bxaj(1:slcats)
950 xxaj_table(1:slcats) = xxaj(1:slcats)
951 bdvic_table(1:slcats) = bdvic(1:slcats)
952 gdvic_table(1:slcats) = gdvic(1:slcats)
953 bbvic_table(1:slcats) = bbvic(1:slcats)
954
955 !---------------- noahmptable.tbl general parameters
956 inquire( file='noahmptable.tbl', exist=file_named )
957 if ( file_named ) then
958 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
959 else
960 open(15, status='old', form='formatted', action='read', iostat=ierr)
961 end if
962 if ( ierr /= 0 ) then
963 errmsg = 'warning: cannot find file noahmptable.tbl'
964 errflg = 1
965 return
966! write(*,'("warning: cannot find file noahmptable.tbl")')
967 endif
968 read(15, noahmp_general_parameters)
969 close(15)
970
971 ! assign values
972 slope_table(1:num_slope) = slope_data(1:num_slope)
973 csoil_table = csoil_data
974 refdk_table = refdk_data
975 refkdt_table = refkdt_data
976 frzk_table = frzk_data
977 zbot_table = zbot_data
978 czil_table = czil_data
979
980 !---------------- noahmptable.tbl radiation parameters
981 inquire( file='noahmptable.tbl', exist=file_named )
982 if ( file_named ) then
983 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
984 else
985 open(15, status='old', form='formatted', action='read', iostat=ierr)
986 end if
987 if (ierr /= 0) then
988 errmsg = 'warning: cannot find file noahmptable.tbl'
989 errflg = 1
990 return
991! write(*,'("warning: cannot find file noahmptable.tbl")')
992 endif
993 read(15,noahmp_rad_parameters)
994 close(15)
995
996 ! assign values
997 albsat_table(:,1) = albsat_vis ! saturated soil albedos: 1=vis, 2=nir
998 albsat_table(:,2) = albsat_nir ! saturated soil albedos: 1=vis, 2=nir
999 albdry_table(:,1) = albdry_vis ! dry soil albedos: 1=vis, 2=nir
1000 albdry_table(:,2) = albdry_nir ! dry soil albedos: 1=vis, 2=nir
1001 albice_table = albice
1002 alblak_table = alblak
1003 omegas_table = omegas
1004 betads_table = betads
1005 betais_table = betais
1006 eg_table = eg
1007! eice_table = eice
1008
1009 !---------------- noahmptable.tbl global parameters
1010 inquire( file='noahmptable.tbl', exist=file_named )
1011 if ( file_named ) then
1012 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
1013 else
1014 open(15, status='old', form='formatted', action='read', iostat=ierr)
1015 end if
1016 if (ierr /= 0) then
1017 errmsg = 'warning: cannot find file noahmptable.tbl'
1018 errflg = 1
1019 return
1020! write(*,'("warning: cannot find file noahmptable.tbl")')
1021 endif
1022 read(15,noahmp_global_parameters)
1023 close(15)
1024
1025 ! assign values
1026 co2_table = co2
1027 o2_table = o2
1028 timean_table = timean
1029 fsatmx_table = fsatmx
1030 z0sno_table = z0sno
1031 ssi_table = ssi
1032 snow_ret_fac_table = snow_ret_fac
1033 snow_emis_table = snow_emis
1034 swemx_table = swemx
1035 tau0_table = tau0
1036 grain_growth_table = grain_growth
1037 extra_growth_table = extra_growth
1038 dirt_soot_table = dirt_soot
1039 bats_cosz_table = bats_cosz
1040 bats_vis_new_table = bats_vis_new
1041 bats_nir_new_table = bats_nir_new
1042 bats_vis_age_table = bats_vis_age
1043 bats_nir_age_table = bats_nir_age
1044 bats_vis_dir_table = bats_vis_dir
1045 bats_nir_dir_table = bats_nir_dir
1046 rsurf_snow_table = rsurf_snow
1047 rsurf_exp_table = rsurf_exp
1048! c2_snowcompact_table = c2_snowcompact
1049! c3_snowcompact_table = c3_snowcompact
1050! c4_snowcompact_table = c4_snowcompact
1051! c5_snowcompact_table = c5_snowcompact
1052! dm_snowcompact_table = dm_snowcompact
1053! eta0_snowcompact_table = eta0_snowcompact
1054! snliqmaxfrac_table = snliqmaxfrac
1055! swemaxgla_table = swemaxgla
1056! wslmax_table = wslmax
1057! rous_table = rous
1058! cmic_table = cmic
1059! snowden_max_table = snowden_max
1060! class_alb_ref_table = class_alb_ref
1061! class_sno_age_table = class_sno_age
1062! class_alb_new_table = class_alb_new
1063! psiwlt_table = psiwlt
1064! z0soil_table = z0soil
1065! z0lake_table = z0lake
1066
1067 !---------------- noahmptable.tbl irrigation parameters
1068 inquire( file='noahmptable.tbl', exist=file_named )
1069 if ( file_named ) then
1070 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
1071 else
1072 open(15, status='old', form='formatted', action='read', iostat=ierr)
1073 end if
1074 if (ierr /= 0) then
1075 errmsg = 'warning: cannot find file noahmptable.tbl'
1076 errflg = 1
1077 return
1078! write(*,'("warning: cannot find file noahmptable.tbl")')
1079 endif
1080 read(15,noahmp_irrigation_parameters)
1081 close(15)
1082
1083 ! assign values
1084 irr_frac_table = irr_frac
1085 irr_har_table = irr_har
1086 irr_lai_table = irr_lai
1087 irr_mad_table = irr_mad
1088 filoss_table = filoss
1089 sprir_rate_table = sprir_rate
1090 micir_rate_table = micir_rate
1091 firtfac_table = firtfac
1092 ir_rain_table = ir_rain
1093
1094 !---------------- noahmptable.tbl crop parameters
1095 inquire( file='noahmptable.tbl', exist=file_named )
1096 if ( file_named ) then
1097 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
1098 else
1099 open(15, status='old', form='formatted', action='read', iostat=ierr)
1100 end if
1101 if (ierr /= 0) then
1102 errmsg = 'warning: cannot find file noahmptable.tbl'
1103 errflg = 1
1104 return
1105! write(*,'("warning: cannot find file noahmptable.tbl")')
1106 endif
1107 read(15,noahmp_crop_parameters)
1108 close(15)
1109
1110 ! assign values
1111 default_crop_table = default_crop
1112 pltday_table = pltday
1113 hsday_table = hsday
1114 plantpop_table = plantpop
1115 irri_table = irri
1116 gddtbase_table = gddtbase
1117 gddtcut_table = gddtcut
1118 gdds1_table = gdds1
1119 gdds2_table = gdds2
1120 gdds3_table = gdds3
1121 gdds4_table = gdds4
1122 gdds5_table = gdds5
1123 c3psni_table(1:5) = c3psni(1:5)
1124 kc25i_table(1:5) = kc25i(1:5)
1125 akci_table(1:5) = akci(1:5)
1126 ko25i_table(1:5) = ko25i(1:5)
1127 akoi_table(1:5) = akoi(1:5)
1128 avcmxi_table(1:5) = avcmxi(1:5)
1129 vcmx25i_table(1:5) = vcmx25i(1:5)
1130 bpi_table(1:5) = bpi(1:5)
1131 mpi_table(1:5) = mpi(1:5)
1132 folnmxi_table(1:5) = folnmxi(1:5)
1133 qe25i_table(1:5) = qe25i(1:5)
1134 aref_table = aref
1135 psnrf_table = psnrf
1136 i2par_table = i2par
1137 tassim0_table = tassim0
1138 tassim1_table = tassim1
1139 tassim2_table = tassim2
1140 k_table = k
1141 epsi_table = epsi
1142 q10mr_table = q10mr
1143 lefreez_table = lefreez
1144 fra_gr_table = fra_gr
1145 lfmr25_table = lfmr25
1146 stmr25_table = stmr25
1147 rtmr25_table = rtmr25
1148 grainmr25_table = grainmr25
1149 bio2lai_table = bio2lai
1150 dile_fc_table(:,1) = dile_fc_s1
1151 dile_fc_table(:,2) = dile_fc_s2
1152 dile_fc_table(:,3) = dile_fc_s3
1153 dile_fc_table(:,4) = dile_fc_s4
1154 dile_fc_table(:,5) = dile_fc_s5
1155 dile_fc_table(:,6) = dile_fc_s6
1156 dile_fc_table(:,7) = dile_fc_s7
1157 dile_fc_table(:,8) = dile_fc_s8
1158 dile_fw_table(:,1) = dile_fw_s1
1159 dile_fw_table(:,2) = dile_fw_s2
1160 dile_fw_table(:,3) = dile_fw_s3
1161 dile_fw_table(:,4) = dile_fw_s4
1162 dile_fw_table(:,5) = dile_fw_s5
1163 dile_fw_table(:,6) = dile_fw_s6
1164 dile_fw_table(:,7) = dile_fw_s7
1165 dile_fw_table(:,8) = dile_fw_s8
1166 lf_ovrc_table(:,1) = lf_ovrc_s1
1167 lf_ovrc_table(:,2) = lf_ovrc_s2
1168 lf_ovrc_table(:,3) = lf_ovrc_s3
1169 lf_ovrc_table(:,4) = lf_ovrc_s4
1170 lf_ovrc_table(:,5) = lf_ovrc_s5
1171 lf_ovrc_table(:,6) = lf_ovrc_s6
1172 lf_ovrc_table(:,7) = lf_ovrc_s7
1173 lf_ovrc_table(:,8) = lf_ovrc_s8
1174 st_ovrc_table(:,1) = st_ovrc_s1
1175 st_ovrc_table(:,2) = st_ovrc_s2
1176 st_ovrc_table(:,3) = st_ovrc_s3
1177 st_ovrc_table(:,4) = st_ovrc_s4
1178 st_ovrc_table(:,5) = st_ovrc_s5
1179 st_ovrc_table(:,6) = st_ovrc_s6
1180 st_ovrc_table(:,7) = st_ovrc_s7
1181 st_ovrc_table(:,8) = st_ovrc_s8
1182 rt_ovrc_table(:,1) = rt_ovrc_s1
1183 rt_ovrc_table(:,2) = rt_ovrc_s2
1184 rt_ovrc_table(:,3) = rt_ovrc_s3
1185 rt_ovrc_table(:,4) = rt_ovrc_s4
1186 rt_ovrc_table(:,5) = rt_ovrc_s5
1187 rt_ovrc_table(:,6) = rt_ovrc_s6
1188 rt_ovrc_table(:,7) = rt_ovrc_s7
1189 rt_ovrc_table(:,8) = rt_ovrc_s8
1190 lfpt_table(:,1) = lfpt_s1
1191 lfpt_table(:,2) = lfpt_s2
1192 lfpt_table(:,3) = lfpt_s3
1193 lfpt_table(:,4) = lfpt_s4
1194 lfpt_table(:,5) = lfpt_s5
1195 lfpt_table(:,6) = lfpt_s6
1196 lfpt_table(:,7) = lfpt_s7
1197 lfpt_table(:,8) = lfpt_s8
1198 stpt_table(:,1) = stpt_s1
1199 stpt_table(:,2) = stpt_s2
1200 stpt_table(:,3) = stpt_s3
1201 stpt_table(:,4) = stpt_s4
1202 stpt_table(:,5) = stpt_s5
1203 stpt_table(:,6) = stpt_s6
1204 stpt_table(:,7) = stpt_s7
1205 stpt_table(:,8) = stpt_s8
1206 rtpt_table(:,1) = rtpt_s1
1207 rtpt_table(:,2) = rtpt_s2
1208 rtpt_table(:,3) = rtpt_s3
1209 rtpt_table(:,4) = rtpt_s4
1210 rtpt_table(:,5) = rtpt_s5
1211 rtpt_table(:,6) = rtpt_s6
1212 rtpt_table(:,7) = rtpt_s7
1213 rtpt_table(:,8) = rtpt_s8
1214 grainpt_table(:,1) = grainpt_s1
1215 grainpt_table(:,2) = grainpt_s2
1216 grainpt_table(:,3) = grainpt_s3
1217 grainpt_table(:,4) = grainpt_s4
1218 grainpt_table(:,5) = grainpt_s5
1219 grainpt_table(:,6) = grainpt_s6
1220 grainpt_table(:,7) = grainpt_s7
1221 grainpt_table(:,8) = grainpt_s8
1222 lfct_table(:,1) = lfct_s1
1223 lfct_table(:,2) = lfct_s2
1224 lfct_table(:,3) = lfct_s3
1225 lfct_table(:,4) = lfct_s4
1226 lfct_table(:,5) = lfct_s5
1227 lfct_table(:,6) = lfct_s6
1228 lfct_table(:,7) = lfct_s7
1229 lfct_table(:,8) = lfct_s8
1230 stct_table(:,1) = stct_s1
1231 stct_table(:,2) = stct_s2
1232 stct_table(:,3) = stct_s3
1233 stct_table(:,4) = stct_s4
1234 stct_table(:,5) = stct_s5
1235 stct_table(:,6) = stct_s6
1236 stct_table(:,7) = stct_s7
1237 stct_table(:,8) = stct_s8
1238 rtct_table(:,1) = rtct_s1
1239 rtct_table(:,2) = rtct_s2
1240 rtct_table(:,3) = rtct_s3
1241 rtct_table(:,4) = rtct_s4
1242 rtct_table(:,5) = rtct_s5
1243 rtct_table(:,6) = rtct_s6
1244 rtct_table(:,7) = rtct_s7
1245 rtct_table(:,8) = rtct_s8
1246
1247 !---------------- noahmptable.tbl tile drainage parameters
1248 inquire( file='noahmptable.tbl', exist=file_named )
1249 if ( file_named ) then
1250 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
1251 else
1252 open(15, status='old', form='formatted', action='read', iostat=ierr)
1253 end if
1254 if (ierr /= 0) then
1255 errmsg = 'warning: cannot find file noahmptable.tbl'
1256 errflg = 1
1257 return
1258! write(*,'("warning: cannot find file noahmptable.tbl")')
1259 endif
1260 read(15,noahmp_tiledrain_parameters)
1261 close(15)
1262
1263 ! assign values
1264 drain_layer_opt_table = drain_layer_opt
1265 tdsmc_fac_table(1:nsoiltype) = tdsmc_fac(1:nsoiltype)
1266 td_depth_table(1:nsoiltype) = td_depth(1:nsoiltype)
1267 td_dc_table(1:nsoiltype) = td_dc(1:nsoiltype)
1268 td_dcoef_table(1:nsoiltype) = td_dcoef(1:nsoiltype)
1269 td_d_table(1:nsoiltype) = td_d(1:nsoiltype)
1270 td_adepth_table(1:nsoiltype) = td_adepth(1:nsoiltype)
1271 td_radi_table(1:nsoiltype) = td_radi(1:nsoiltype)
1272 td_spac_table(1:nsoiltype) = td_spac(1:nsoiltype)
1273 td_ddrain_table(1:nsoiltype) = td_ddrain(1:nsoiltype)
1274 klat_fac_table(1:nsoiltype) = klat_fac(1:nsoiltype)
1275
1276 !---------------- noahmptable.tbl optional parameters
1277 inquire( file='noahmptable.tbl', exist=file_named )
1278 if ( file_named ) then
1279 open(15, file="noahmptable.tbl", status='old', form='formatted', action='read', iostat=ierr)
1280 else
1281 open(15, status='old', form='formatted', action='read', iostat=ierr)
1282 end if
1283 if (ierr /= 0) then
1284 errmsg = 'warning: cannot find file noahmptable.tbl'
1285 errflg = 1
1286 return
1287! write(*,'("warning: cannot find file noahmptable.tbl")')
1288 endif
1289 read(15,noahmp_optional_parameters)
1290 close(15)
1291
1292 ! assign values
1293! sr2006_theta_1500t_a_table = sr2006_theta_1500t_a
1294! sr2006_theta_1500t_b_table = sr2006_theta_1500t_b
1295! sr2006_theta_1500t_c_table = sr2006_theta_1500t_c
1296! sr2006_theta_1500t_d_table = sr2006_theta_1500t_d
1297! sr2006_theta_1500t_e_table = sr2006_theta_1500t_e
1298! sr2006_theta_1500t_f_table = sr2006_theta_1500t_f
1299! sr2006_theta_1500t_g_table = sr2006_theta_1500t_g
1300! sr2006_theta_1500_a_table = sr2006_theta_1500_a
1301! sr2006_theta_1500_b_table = sr2006_theta_1500_b
1302! sr2006_theta_33t_a_table = sr2006_theta_33t_a
1303! sr2006_theta_33t_b_table = sr2006_theta_33t_b
1304! sr2006_theta_33t_c_table = sr2006_theta_33t_c
1305! sr2006_theta_33t_d_table = sr2006_theta_33t_d
1306! sr2006_theta_33t_e_table = sr2006_theta_33t_e
1307! sr2006_theta_33t_f_table = sr2006_theta_33t_f
1308! sr2006_theta_33t_g_table = sr2006_theta_33t_g
1309! sr2006_theta_33_a_table = sr2006_theta_33_a
1310! sr2006_theta_33_b_table = sr2006_theta_33_b
1311! sr2006_theta_33_c_table = sr2006_theta_33_c
1312! sr2006_theta_s33t_a_table = sr2006_theta_s33t_a
1313! sr2006_theta_s33t_b_table = sr2006_theta_s33t_b
1314! sr2006_theta_s33t_c_table = sr2006_theta_s33t_c
1315! sr2006_theta_s33t_d_table = sr2006_theta_s33t_d
1316! sr2006_theta_s33t_e_table = sr2006_theta_s33t_e
1317! sr2006_theta_s33t_f_table = sr2006_theta_s33t_f
1318! sr2006_theta_s33t_g_table = sr2006_theta_s33t_g
1319! sr2006_theta_s33_a_table = sr2006_theta_s33_a
1320! sr2006_theta_s33_b_table = sr2006_theta_s33_b
1321! sr2006_psi_et_a_table = sr2006_psi_et_a
1322! sr2006_psi_et_b_table = sr2006_psi_et_b
1323! sr2006_psi_et_c_table = sr2006_psi_et_c
1324! sr2006_psi_et_d_table = sr2006_psi_et_d
1325! sr2006_psi_et_e_table = sr2006_psi_et_e
1326! sr2006_psi_et_f_table = sr2006_psi_et_f
1327! sr2006_psi_et_g_table = sr2006_psi_et_g
1328! sr2006_psi_e_a_table = sr2006_psi_e_a
1329! sr2006_psi_e_b_table = sr2006_psi_e_b
1330! sr2006_psi_e_c_table = sr2006_psi_e_c
1331! sr2006_smcmax_a_table = sr2006_smcmax_a
1332! sr2006_smcmax_b_table = sr2006_smcmax_b
1333
1334 end subroutine read_mp_table_parameters
1335
1336end module noahmp_tables
1337
Data from MPTABLE.TBL, SOILPARM.TBL, GENPARM.TBL for NoahMP.