19 SnowUse, EmissionsMethod, NetRadiationMethod, RoughLenHeatMethod, & ! model options
20 RoughLenMomMethod, StorageHeatMethod, AerodynamicResistanceMethod, LAIType, SnowPackLimit, &
21 OHM_threshSW, OHM_threshWD, th, tl, Kmax, g1, g2, g3, g4, g5, g6, s1, s2, LaiPower, &
37 OHMIncQF, & ! model options
38 iy, id, it, imin, isec, dt_since_start, tstep, tstep_prev, startDLS, endDLS, & ! time-related input
39 alt, lat, lng, Z, timezone, SurfaceArea, sfr_surf, & ! site-specific geographical settings
40 z0m_in, zdm_in, & ! roughness related settings
41 alb, emis, SnowAlb, OHM_coef, WaterDist, & ! surface properties
42 AHProf_24hr, HumActivity_24hr, PopProf_24hr, TraffProf_24hr, WUProfA_24hr, WUProfM_24hr, snowProf_24hr, & ! hourly profile values
43 qn_av, dqndt, qn_s_av, dqnsdt, & ! OHM related Qn quantities
44 surf_var_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, &
60 AH_SLOPE_Cooling_id, &
61 AH_SLOPE_Heating_id, &
66 T_CRITIC_Cooling_id, &
67 T_CRITIC_Heating_id, &
69 surf_attr_MinStorCap_id, &
70 surf_attr_DrainEquat_id, &
71 surf_attr_DrainCoef1_id, &
72 surf_attr_DrainCoef2_id, &
73 surf_attr_MaxStorCap_id, &
76 SatHydraulicConduct_id, &
102 PopDensDaytime_id, PopDensNighttime_id, &
104 GDD_id, HDD_id, LAI_id, WUDay_id, soilstore_id, state_id, SnowWater, &
105 avkdn, avRh, avU1, Press_hPa, Temp_C, Precip, & ! forcing variables
106 qn, qf, qs, qh, qe, qsfc, tsk, CHKLOWQ, &
111 tsurf_o, smd_o, glw1d, cldfra1d)
114 INTEGER,
INTENT(in) :: SnowUse
115 INTEGER,
INTENT(in) :: EmissionsMethod
116 INTEGER,
INTENT(in) :: NetRadiationMethod
117 INTEGER,
INTENT(IN) :: RoughLenHeatMethod
118 INTEGER,
INTENT(IN) :: RoughLenMomMethod
119 INTEGER,
INTENT(IN) :: StorageHeatMethod
120 INTEGER,
INTENT(IN) :: AerodynamicResistanceMethod
121 INTEGER,
DIMENSION(3),
INTENT(IN) :: LAIType
122 REAL(KIND(1D0)),
DIMENSION(7 + 1),
INTENT(IN) :: SnowPackLimit
123 REAL(KIND(1D0)),
DIMENSION(7 + 1),
INTENT(IN) :: OHM_threshSW
124 REAL(KIND(1D0)),
DIMENSION(7 + 1),
INTENT(IN) :: OHM_threshWD
125 REAL(KIND(1D0)),
INTENT(IN) :: th
126 REAL(KIND(1D0)),
INTENT(IN) :: tl
127 REAL(KIND(1D0)),
INTENT(IN) :: Kmax
128 REAL(KIND(1D0)),
INTENT(IN) :: g1
129 REAL(KIND(1D0)),
INTENT(IN) :: g2
130 REAL(KIND(1D0)),
INTENT(IN) :: g3
131 REAL(KIND(1D0)),
INTENT(IN) :: g4
132 REAL(KIND(1D0)),
INTENT(IN) :: g5
133 REAL(KIND(1D0)),
INTENT(IN) :: g6
134 REAL(KIND(1D0)),
INTENT(IN) :: s1
135 REAL(KIND(1D0)),
INTENT(IN) :: s2
136 REAL(KIND(1D0)),
DIMENSION(4, 3),
INTENT(IN) :: LaiPower
137 REAL(KIND(1D0)),
INTENT(IN) :: CRWmax
138 REAL(KIND(1D0)),
INTENT(IN) :: CRWmin
139 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimit
140 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimitAlb
141 REAL(KIND(1D0)),
INTENT(IN) :: RadMeltFact
142 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMax
143 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMin
144 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMax
145 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMin
146 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimBldg
147 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimPaved
148 REAL(KIND(1D0)),
INTENT(IN) :: tau_a
149 REAL(KIND(1D0)),
INTENT(IN) :: tau_f
150 REAL(KIND(1D0)),
INTENT(IN) :: tau_r
151 REAL(KIND(1D0)),
INTENT(IN) :: TempMeltFact
153 INTEGER,
INTENT(IN) :: OHMIncQF
156 INTEGER,
INTENT(IN) :: iy
157 INTEGER,
INTENT(IN) :: id
158 INTEGER,
INTENT(IN) :: it
159 INTEGER,
INTENT(IN) :: imin
160 INTEGER,
INTENT(in) :: isec
161 INTEGER,
INTENT(in) :: dt_since_start
162 INTEGER,
INTENT(IN) :: tstep
163 INTEGER,
INTENT(IN) :: tstep_prev
164 INTEGER,
INTENT(IN) :: startDLS
165 INTEGER,
INTENT(IN) :: endDLS
168 REAL(KIND(1D0)),
INTENT(IN) :: alt
169 REAL(KIND(1D0)),
INTENT(IN) :: lat
170 REAL(KIND(1D0)),
INTENT(IN) :: lng
171 REAL(KIND(1D0)),
INTENT(IN) :: Z
172 REAL(KIND(1D0)),
INTENT(IN) :: timezone
173 REAL(KIND(1D0)),
INTENT(IN) :: SurfaceArea
174 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: sfr_surf
177 REAL(KIND(1D0)),
INTENT(in) :: z0m_in
178 REAL(KIND(1D0)),
INTENT(in) :: zdm_in
181 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: alb
182 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: emis
183 REAL(KIND(1D0)),
INTENT(INOUT) :: SnowAlb
186 REAL(KIND(1D0)),
DIMENSION(7 + 1, 4, 3),
INTENT(IN) :: OHM_coef
189 REAL(KIND(1D0)),
DIMENSION(7 + 1, 7 - 1),
INTENT(IN) :: WaterDist
192 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(in) :: AHProf_24hr
193 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(in) :: HumActivity_24hr
194 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(in) :: PopProf_24hr
195 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(in) :: TraffProf_24hr
196 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(in) :: WUProfA_24hr
197 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(in) :: WUProfM_24hr
198 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(in) :: snowProf_24hr
201 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_av
202 REAL(KIND(1D0)),
INTENT(INOUT) :: dqndt
203 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_s_av
204 REAL(KIND(1D0)),
INTENT(INOUT) :: dqnsdt
205 REAL(KIND(1D0)),
INTENT(INOUT) :: DecidCap_id
206 REAL(KIND(1D0)),
INTENT(INOUT) :: albDecTr_id
207 REAL(KIND(1D0)),
INTENT(INOUT) :: albEveTr_id
208 REAL(KIND(1D0)),
INTENT(INOUT) :: albGrass_id
209 REAL(KIND(1D0)),
INTENT(INOUT) :: NumCapita_id
211 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: BaseT_id
212 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: BaseTe_id
213 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: GDDFull_id
214 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: SDDFull_id
215 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: LaiMin_id
216 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: LaiMax_id
217 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: MaxConductance_id
218 REAL(KIND(1D0)),
INTENT(INOUT) :: FAIbldg_id
219 REAL(KIND(1D0)),
INTENT(INOUT) :: FAIEveTree_id
220 REAL(KIND(1D0)),
INTENT(INOUT) :: FAIDecTree_id
221 REAL(KIND(1D0)),
INTENT(INOUT) :: bldgH_id
222 REAL(KIND(1D0)),
INTENT(INOUT) :: EveTreeH_id
223 REAL(KIND(1D0)),
INTENT(INOUT) :: DecTreeH_id
224 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: AH_MIN_id
225 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: AH_SLOPE_Cooling_id
226 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: AH_SLOPE_Heating_id
227 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: QF0_BEU_id
228 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: Qf_A_id
229 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: Qf_B_id
230 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: Qf_C_id
231 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: T_CRITIC_Cooling_id
232 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: T_CRITIC_Heating_id
233 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(INOUT) :: TrafficRate_id
234 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: surf_attr_MinStorCap_id
235 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: surf_attr_DrainEquat_id
236 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: surf_attr_DrainCoef1_id
237 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: surf_attr_DrainCoef2_id
238 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: surf_attr_MaxStorCap_id
240 REAL(KIND(1D0)),
INTENT(INOUT) :: PopDensDaytime_id
241 REAL(KIND(1D0)),
INTENT(INOUT) :: PopDensNighttime_id
242 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: SoilStoreCap_id
243 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: SoilDepth_id
244 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: SatHydraulicConduct_id
245 REAL(KIND(1D0)),
INTENT(INOUT) :: AlbMin_DecTr_id
246 REAL(KIND(1D0)),
INTENT(INOUT) :: AlbMax_DecTr_id
247 REAL(KIND(1D0)),
INTENT(INOUT) :: AlbMin_EveTr_id
248 REAL(KIND(1D0)),
INTENT(INOUT) :: AlbMax_EveTr_id
249 REAL(KIND(1D0)),
INTENT(INOUT) :: AlbMin_Grass_id
250 REAL(KIND(1D0)),
INTENT(INOUT) :: AlbMax_Grass_id
251 REAL(KIND(1D0)),
INTENT(INOUT) :: CapMin_dec_id
252 REAL(KIND(1D0)),
INTENT(INOUT) :: CapMax_dec_id
253 REAL(KIND(1D0)),
INTENT(INOUT) :: PorMin_dec_id
254 REAL(KIND(1D0)),
INTENT(INOUT) :: PorMax_dec_id
255 REAL(KIND(1D0)),
INTENT(INOUT) :: DRAINRT_id
256 REAL(KIND(1D0)),
INTENT(INOUT) :: RAINCOVER_id
257 REAL(KIND(1D0)),
INTENT(INOUT) :: RAINMAXRES_id
258 REAL(KIND(1D0)),
INTENT(INOUT) :: FlowChange_id
259 REAL(KIND(1D0)),
INTENT(INOUT) :: PipeCapacity_id
260 REAL(KIND(1D0)),
INTENT(INOUT) :: RunoffToWater_id
261 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: StateLimit_id
262 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: WetThresh_id
263 REAL(KIND(1D0)),
INTENT(INOUT) :: BaseT_HC_id
264 REAL(KIND(1D0)),
INTENT(INOUT) :: Faut_id
265 REAL(KIND(1D0)),
INTENT(INOUT) :: IrrFracConif_id
266 REAL(KIND(1D0)),
INTENT(INOUT) :: IrrFracDecid_id
267 REAL(KIND(1D0)),
INTENT(INOUT) :: IrrFracGrass_id
268 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: DayWat_id
269 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: DayWatPer_id
270 REAL(KIND(1D0)),
INTENT(INOUT) :: porosity_id
271 REAL(KIND(1D0)),
DIMENSION(5),
INTENT(INOUT) :: GDD_id
272 REAL(KIND(1D0)),
DIMENSION(12),
INTENT(INOUT) :: HDD_id
273 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(INOUT) :: LAI_id
274 REAL(KIND(1D0)),
DIMENSION(9),
INTENT(INOUT) :: WUDay_id
275 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: soilstore_id
276 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: state_id
277 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: surf_var_id
278 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(INOUT) :: SnowWater
281 REAL(KIND(1D0)),
INTENT(IN) :: avkdn
282 REAL(KIND(1D0)),
INTENT(IN) :: glw1d
283 REAL(KIND(1D0)),
INTENT(IN) :: avRh
284 REAL(KIND(1D0)),
INTENT(IN) :: avU1
285 REAL(KIND(1D0)),
INTENT(IN) :: Press_hPa
286 REAL(KIND(1D0)),
INTENT(IN) :: Temp_C
287 REAL(KIND(1D0)),
INTENT(IN) :: Precip
288 REAL(KIND(1D0)),
INTENT(IN) :: cldfra1d
290 REAL(KIND(1D0)),
INTENT(out) :: kdown_o
291 REAL(KIND(1D0)),
INTENT(out) :: kup_o
292 REAL(KIND(1D0)),
INTENT(out) :: ldown_o
293 REAL(KIND(1D0)),
INTENT(out) :: lup_o
294 REAL(KIND(1D0)),
INTENT(out) :: tsurf_o
295 REAL(KIND(1D0)),
INTENT(out) :: smd_o
297 REAL(KIND(1D0)),
INTENT(out) :: qn
298 REAL(KIND(1D0)),
INTENT(out) :: qf
299 REAL(KIND(1D0)),
INTENT(out) :: qs
300 REAL(KIND(1D0)),
INTENT(out) :: qh
301 REAL(KIND(1D0)),
INTENT(out) :: qe
302 REAL(KIND(1D0)),
INTENT(out) :: qsfc
303 REAL(KIND(1D0)),
INTENT(out) :: tsk
304 REAL(KIND(1D0)),
INTENT(out) :: CHKLOWQ
307 INTEGER,
PARAMETER :: veg_type = 1
308 INTEGER,
PARAMETER :: gsModel = 2
309 INTEGER,
PARAMETER :: StabilityMethod = 3
310 INTEGER,
PARAMETER :: SMDMethod = 0
311 INTEGER,
PARAMETER :: DiagQN = 0
312 INTEGER,
PARAMETER :: DiagQS = 0
313 INTEGER,
PARAMETER :: Diagnose = 0
314 INTEGER,
PARAMETER :: EvapMethod = 2
315 INTEGER,
PARAMETER :: LAICalcYes = 1
316 INTEGER,
PARAMETER :: WaterUseMethod = 0
318 REAL(KIND(1D0)),
PARAMETER :: LAI_obs = 0
319 REAL(KIND(1D0)) :: ldown_obs
320 REAL(KIND(1D0)) :: fcld_obs
321 REAL(KIND(1D0)),
PARAMETER :: snowFrac_obs = 0
322 REAL(KIND(1D0)),
PARAMETER :: qn1_obs = 0
323 REAL(KIND(1D0)),
PARAMETER :: qh_obs = 0
324 REAL(KIND(1D0)),
PARAMETER :: qf_obs = 0
325 REAL(KIND(1D0)),
PARAMETER :: qs_obs = 0
332 CHARACTER(len=1024) :: message
427 REAL(KIND(1D0)),
PARAMETER :: xsmd = 0.
429 REAL(KIND(1D0)),
PARAMETER :: wu_m3 = 0
431 REAL(KIND(1D0)),
PARAMETER :: InternalWaterUse_h = 0
438 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: Ie_a = [-84.535, 9.959, 3.674]
439 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: Ie_m = [-25.36, 2.988, 1.102]
442 REAL(KIND(1D0)),
PARAMETER :: NARP_EMIS_SNOW = 0.9
443 REAL(KIND(1D0)) :: NARP_TRANS_SITE
481 REAL(KIND(1D0)),
PARAMETER :: EF_umolCO2perJ = 1.159
482 REAL(KIND(1D0)),
PARAMETER :: EnEF_v_Jkm = 4e6
483 REAL(KIND(1D0)),
PARAMETER :: FcEF_v_kgkm = 0.285
484 REAL(KIND(1D0)),
PARAMETER :: FrFossilFuel_Heat = 0.05
485 REAL(KIND(1D0)),
PARAMETER :: FrFossilFuel_NonHeat = 0
486 REAL(KIND(1D0)),
PARAMETER :: TrafficUnits = 1
487 REAL(KIND(1D0)),
PARAMETER :: MaxQFMetab = 175
488 REAL(KIND(1D0)),
PARAMETER :: MinQFMetab = 75
491 REAL(KIND(1D0)),
DIMENSION(7),
PARAMETER :: chAnOHM = 3
492 REAL(KIND(1D0)),
DIMENSION(7),
PARAMETER :: cpAnOHM = 2e6
493 REAL(KIND(1D0)),
DIMENSION(7),
PARAMETER :: kkAnOHM = 1.2
494 REAL(KIND(1D0)),
DIMENSION(:, :),
ALLOCATABLE :: MetForcingData_grid
497 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: alpha_bioCO2 = 0.005
498 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: alpha_enh_bioCO2 = 0.016
499 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: beta_bioCO2 = 8.747
500 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: beta_enh_bioCO2 = 33.454
501 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: min_res_bioCO2 = 0.6
502 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: resp_a = 2.43
503 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: resp_b = 0
504 REAL(KIND(1D0)),
DIMENSION(3),
PARAMETER :: theta_bioCO2 = 0.96
507 REAL(KIND(1D0)),
DIMENSION(24*3600/tstep) :: Tair24HR
508 REAL(KIND(1D0)),
DIMENSION(:),
ALLOCATABLE :: Ts5mindata_ir
510 REAL(KIND(1D0)) :: SnowfallCum = 0
511 REAL(KIND(1D0)),
DIMENSION(7) :: IceFrac = 0.2
512 REAL(KIND(1D0)),
DIMENSION(7) :: SnowDens = 300
513 REAL(KIND(1D0)),
DIMENSION(7) :: SnowFrac = 0
514 REAL(KIND(1D0)),
DIMENSION(7) :: SnowPack = 0
516 REAL(KIND(1D0)),
DIMENSION(5) :: datetimeLine
517 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSUEWS - 5) :: dataOutLineSUEWS
518 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSnow - 5) :: dataOutLineSnow
519 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutESTM - 5) :: dataOutLineESTM
520 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutDailyState - 5) :: DailyStateLine
523 REAL(KIND(1D0)),
DIMENSION(6, 7) :: StoreDrainPrm
524 storedrainprm(1, :) = surf_attr_minstorcap_id
525 storedrainprm(2, :) = surf_attr_drainequat_id
526 storedrainprm(3, :) = surf_attr_draincoef1_id
527 storedrainprm(4, :) = surf_attr_draincoef2_id
528 storedrainprm(5, :) = surf_attr_maxstorcap_id
530 storedrainprm(6, :) = surf_var_id
539 WRITE (message, *)
'in SuMin, before calculation, OHM_coef:', ohm_coef(1, :, :)
540 CALL wrf_debug(100, message)
544 aerodynamicresistancemethod, ah_min_id, ahprof_24hr, ah_slope_cooling_id, &
545 ah_slope_heating_id, &
546 alb, albmax_dectr_id, albmax_evetr_id, albmax_grass_id, &
547 albmin_dectr_id, albmin_evetr_id, albmin_grass_id, &
548 alpha_bioco2, alpha_enh_bioco2, alt, avkdn, avrh, avu1, baset_id, basete_id, &
549 baset_hc_id, beta_bioco2, beta_enh_bioco2, bldgh_id, capmax_dec_id, capmin_dec_id, &
550 chanohm, cpanohm, crwmax, crwmin, daywat_id, daywatper_id, &
551 dectreeh_id, diagnose, diagqn, diagqs, drainrt_id, &
552 dt_since_start, dqndt, qn_av, dqnsdt, qn_s_av, &
553 ef_umolco2perj, emis, emissionsmethod, enef_v_jkm, enddls, evetreeh_id, faibldg_id, &
554 faidectree_id, faievetree_id, faut_id, fcef_v_kgkm, fcld_obs, flowchange_id, &
555 frfossilfuel_heat, frfossilfuel_nonheat, g1, g2, g3, g4, g5, g6, gdd_id, &
556 gddfull_id, gridiv, gsmodel, hdd_id, humactivity_24hr, &
557 icefrac, id, ie_a, ie_end, ie_m, ie_start, imin, &
558 internalwateruse_h, irrfracconif_id, irrfracdecid_id, irrfracgrass_id, isec, it, evapmethod, &
559 iy, kkanohm, kmax, lai_id, laicalcyes, laimax_id, laimin_id, lai_obs, &
560 laipower, laitype, lat, ldown_obs, lng, maxconductance_id, maxqfmetab, &
561 snowwater, metforcingdata_grid, minqfmetab, min_res_bioco2, &
562 narp_emis_snow, narp_trans_site, netradiationmethod, &
563 numcapita_id, ohm_coef, ohmincqf, ohm_threshsw, &
564 ohm_threshwd, pipecapacity_id, popdensdaytime_id, &
565 popdensnighttime_id, popprof_24hr, pormax_dec_id, pormin_dec_id, &
566 precip, preciplimit, preciplimitalb, press_hpa, &
567 qf0_beu_id, qf_a_id, qf_b_id, qf_c_id, &
568 qn1_obs, qh_obs, qs_obs, qf_obs, &
569 radmeltfact, raincover_id, rainmaxres_id, resp_a, resp_b, &
570 roughlenheatmethod, roughlenmommethod, runofftowater_id, s1, s2, &
571 sathydraulicconduct_id, sddfull_id, sfr_surf, smdmethod, snowalb, snowalbmax, &
572 snowalbmin, snowpacklimit, snowdens, snowdensmax, snowdensmin, snowfallcum, snowfrac, &
573 snowlimbldg, snowlimpaved, snowfrac_obs, snowpack, snowprof_24hr, snowuse, soildepth_id, &
574 soilstore_id, soilstorecap_id, stabilitymethod, startdls, state_id, statelimit_id, &
575 storageheatmethod, storedrainprm, surfacearea, tair24hr, tau_a, tau_f, tau_r, &
576 t_critic_cooling_id, t_critic_heating_id, temp_c, tempmeltfact, th, &
577 theta_bioco2, timezone, tl, trafficrate_id, trafficunits, &
578 traffprof_24hr, ts5mindata_ir, tstep, tstep_prev, veg_type, &
579 waterdist, waterusemethod, wetthresh_id, wu_m3, &
580 wuday_id, decidcap_id, albdectr_id, albevetr_id, albgrass_id, porosity_id, &
581 wuprofa_24hr, wuprofm_24hr, xsmd, z, z0m_in, zdm_in, &
582 datetimeline, dataoutlinesuews, dataoutlinesnow, dataoutlineestm, &
585 surf_var_id = storedrainprm(6, :)
587 kdown_o = dataoutlinesuews(1)
588 kup_o = dataoutlinesuews(2)
589 ldown_o = dataoutlinesuews(3)
590 lup_o = dataoutlinesuews(4)
591 tsurf_o = dataoutlinesuews(5)
592 smd_o = dataoutlinesuews(23)
594 qn = dataoutlinesuews(6)
595 qf = dataoutlinesuews(7)
596 qs = dataoutlinesuews(8)
597 qh = dataoutlinesuews(9)
598 qe = dataoutlinesuews(10)
599 qsfc = dataoutlinesuews(16)
600 tsk = dataoutlinesuews(77) + 273.15
614 WRITE (message, *)
' in SuMin, after calculation, OHM_coef:', ohm_coef(1, :, :)
615 CALL wrf_debug(100, message)
617 WRITE (message, *)
' in SuMin, qn,qf,qs,qh,qe:', dataoutlinesuews(6:10)
618 CALL wrf_debug(100, message)
subroutine suews_cal_main(ah_min, ahprof_24hr, ah_slope_cooling, ah_slope_heating, alb, albmax_dectr, albmax_evetr, albmax_grass, albmin_dectr, albmin_evetr, albmin_grass, alpha_bioco2, alpha_enh_bioco2, alt, kdown, avrh, avu1, baset, basete, beta_bioco2, beta_enh_bioco2, bldgh, capmax_dec, capmin_dec, chanohm, co2pointsource, cpanohm, crwmax, crwmin, daywat, daywatper, dectreeh, diagmethod, diagnose, drainrt, dt_since_start, dqndt, qn_av, dqnsdt, qn_s_av, ef_umolco2perj, emis, emissionsmethod, enef_v_jkm, enddls, evetreeh, faibldg, faidectree, faievetree, faimethod, faut, fcef_v_kgkm, fcld_obs, flowchange, frfossilfuel_heat, frfossilfuel_nonheat, g_max, g_k, g_q_base, g_q_shape, g_t, g_sm, gdd_id, gddfull, gridiv, gsmodel, h_maintain, hdd_id, humactivity_24hr, icefrac, id, ie_a, ie_end, ie_m, ie_start, imin, internalwateruse_h, irrfracpaved, irrfracbldgs, irrfracevetr, irrfracdectr, irrfracgrass, irrfracbsoil, irrfracwater, isec, it, iy, kkanohm, kmax, lai_id, laimax, laimin, lai_obs, laipower, laitype, lat, lenday_id, ldown_obs, lng, maxconductance, maxfcmetab, maxqfmetab, snowwater, minfcmetab, minqfmetab, min_res_bioco2, narp_emis_snow, narp_trans_site, netradiationmethod, nlayer, n_vegetation_region_urban, n_stream_sw_urban, n_stream_lw_urban, sw_dn_direct_frac, air_ext_sw, air_ssa_sw, veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, veg_fsd_const, veg_contact_fraction_const, ground_albedo_dir_mult_fact, use_sw_direct_albedo, height, building_frac, veg_frac, building_scale, veg_scale, alb_roof, emis_roof, alb_wall, emis_wall, roof_albedo_dir_mult_fact, wall_specular_frac, ohm_coef, ohmincqf, ohm_threshsw, ohm_threshwd, pipecapacity, popdensdaytime, popdensnighttime, popprof_24hr, pormax_dec, pormin_dec, precip, preciplimit, preciplimitalb, press_hpa, qf0_beu, qf_a, qf_b, qf_c, qn1_obs, qs_obs, qf_obs, radmeltfact, raincover, rainmaxres, resp_a, resp_b, roughlenheatmethod, roughlenmommethod, runofftowater, s1, s2, sathydraulicconduct, sddfull, sdd_id, smdmethod, snowalb, snowalbmax, snowalbmin, snowpacklimit, snowdens, snowdensmax, snowdensmin, snowfallcum, snowfrac, snowlimbldg, snowlimpaved, snowfrac_obs, snowpack, snowprof_24hr, snowuse, soildepth, stabilitymethod, startdls, soilstore_surf, soilstorecap_surf, state_surf, statelimit_surf, wetthresh_surf, soilstore_roof, soilstorecap_roof, state_roof, statelimit_roof, wetthresh_roof, soilstore_wall, soilstorecap_wall, state_wall, statelimit_wall, wetthresh_wall, storageheatmethod, storedrainprm, surfacearea, tair_av, tau_a, tau_f, tau_r, tmax_id, tmin_id, baset_cooling, baset_heating, temp_c, tempmeltfact, th, theta_bioco2, timezone, tl, trafficrate, trafficunits, sfr_surf, tsfc_roof, tsfc_wall, tsfc_surf, temp_roof, temp_wall, temp_surf, tin_roof, tin_wall, tin_surf, k_roof, k_wall, k_surf, cp_roof, cp_wall, cp_surf, dz_roof, dz_wall, dz_surf, traffprof_24hr, ts5mindata_ir, tstep, tstep_prev, veg_type, waterdist, waterusemethod, wu_m3, wuday_id, decidcap_id, albdectr_id, albevetr_id, albgrass_id, porosity_id, wuprofa_24hr, wuprofm_24hr, xsmd, z, z0m_in, zdm_in, output_line_suews)
subroutine sumin(snowuse, emissionsmethod, netradiationmethod, roughlenheatmethod, roughlenmommethod, storageheatmethod, aerodynamicresistancemethod, laitype, snowpacklimit, ohm_threshsw, ohm_threshwd, th, tl, kmax, g1, g2, g3, g4, g5, g6, s1, s2, laipower, crwmax, crwmin, preciplimit, preciplimitalb, radmeltfact, snowalbmax, snowalbmin, snowdensmax, snowdensmin, snowlimbldg, snowlimpaved, tau_a, tau_f, tau_r, tempmeltfact, ohmincqf, iy, id, it, imin, isec, dt_since_start, tstep, tstep_prev, startdls, enddls, alt, lat, lng, z, timezone, surfacearea, sfr_surf, z0m_in, zdm_in, alb, emis, snowalb, ohm_coef, waterdist, ahprof_24hr, humactivity_24hr, popprof_24hr, traffprof_24hr, wuprofa_24hr, wuprofm_24hr, snowprof_24hr, qn_av, dqndt, qn_s_av, dqnsdt, surf_var_id, decidcap_id, albdectr_id, albevetr_id, albgrass_id, numcapita_id, baset_id, basete_id, gddfull_id, sddfull_id, laimin_id, laimax_id, maxconductance_id, faibldg_id, faievetree_id, faidectree_id, bldgh_id, evetreeh_id, dectreeh_id, ah_min_id, ah_slope_cooling_id, ah_slope_heating_id, qf0_beu_id, qf_a_id, qf_b_id, qf_c_id, t_critic_cooling_id, t_critic_heating_id, trafficrate_id, surf_attr_minstorcap_id, surf_attr_drainequat_id, surf_attr_draincoef1_id, surf_attr_draincoef2_id, surf_attr_maxstorcap_id, soilstorecap_id, soildepth_id, sathydraulicconduct_id, albmin_dectr_id, albmax_dectr_id, albmin_evetr_id, albmax_evetr_id, albmin_grass_id, albmax_grass_id, capmin_dec_id, capmax_dec_id, pormin_dec_id, pormax_dec_id, drainrt_id, raincover_id, rainmaxres_id, flowchange_id, pipecapacity_id, runofftowater_id, statelimit_id, wetthresh_id, baset_hc_id, faut_id, irrfracconif_id, irrfracdecid_id, irrfracgrass_id, daywat_id, daywatper_id, popdensdaytime_id, popdensnighttime_id, porosity_id, gdd_id, hdd_id, lai_id, wuday_id, soilstore_id, state_id, snowwater, avkdn, avrh, avu1, press_hpa, temp_c, precip, qn, qf, qs, qh, qe, qsfc, tsk, chklowq, kdown_o, kup_o, ldown_o, lup_o, tsurf_o, smd_o, glw1d, cldfra1d)