42 AerodynamicResistanceMethod, AH_MIN, AHProf_24hr, AH_SLOPE_Cooling, & ! input&inout in alphabetical order
44 alb, AlbMax_DecTr, AlbMax_EveTr, AlbMax_Grass, &
45 AlbMin_DecTr, AlbMin_EveTr, AlbMin_Grass, &
46 alpha_bioCO2, alpha_enh_bioCO2, alt, avkdn, avRh, avU1, BaseT, BaseTe, &
47 BaseTHDD, beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, &
48 chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, &
49 DecTreeH, Diagnose, DiagQN, DiagQS, DRAINRT, &
50 dt_since_start, dqndt, qn1_av, dqnsdt, qn1_s_av, &
51 EF_umolCO2perJ, emis, EmissionsMethod, EnEF_v_Jkm, endDLS, EveTreeH, FAIBldg, &
52 FAIDecTree, FAIEveTree, Faut, FcEF_v_kgkm, fcld_obs, FlowChange, &
53 FrFossilFuel_Heat, FrFossilFuel_NonHeat, G1, G2, G3, G4, G5, G6, GDD_id, &
54 GDDFull, Gridiv, gsModel, HDD_id, HumActivity_24hr, &
55 IceFrac, id, Ie_a, Ie_end, Ie_m, Ie_start, imin, &
56 InternalWaterUse_h, IrrFracConif, IrrFracDecid, IrrFracGrass, isec, it, EvapMethod, &
57 iy, kkAnOHM, Kmax, LAI_id, LAICalcYes, LAIMax, LAIMin, LAI_obs, &
58 LAIPower, LAIType, lat, lenDay_id, ldown_obs, lng, MaxConductance, MaxFCMetab, MaxQFMetab, &
59 SnowWater, MetForcingData_grid, MinFCMetab, MinQFMetab, min_res_bioCO2, &
60 NARP_EMIS_SNOW, NARP_TRANS_SITE, NetRadiationMethod, &
61 OHM_coef, OHMIncQF, OHM_threshSW, &
62 OHM_threshWD, PipeCapacity, PopDensDaytime, &
63 PopDensNighttime, PopProf_24hr, PorMax_dec, PorMin_dec, &
64 Precip, PrecipLimit, PrecipLimitAlb, Press_hPa, &
65 QF0_BEU, Qf_A, Qf_B, Qf_C, &
66 qn1_obs, qs_obs, qf_obs, &
67 RadMeltFact, RAINCOVER, RainMaxRes, resp_a, resp_b, &
68 RoughLenHeatMethod, RoughLenMomMethod, RunoffToWater, S1, S2, &
69 SatHydraulicConduct, SDDFull, SDD_id, sfr, SMDMethod, SnowAlb, SnowAlbMax, &
70 SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, &
71 SnowLimBldg, SnowLimPaved, snowFrac_obs, SnowPack, SnowProf_24hr, snowUse, SoilDepth, &
72 soilstore_id, SoilStoreCap, StabilityMethod, startDLS, state_id, StateLimit, &
73 StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, &
75 T_CRITIC_Cooling, T_CRITIC_Heating, Temp_C, TempMeltFact, TH, &
76 theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, &
77 TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, &
78 WaterDist, WaterUseMethod, WetThresh, wu_m3, &
79 WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, &
80 WUProfA_24hr, WUProfM_24hr, xsmd, Z, z0m_in, zdm_in, &
81 datetimeLine, dataOutLineSUEWS, dataOutLineSnow, dataOutLineESTM, dataoutLineRSL, &!output
88 INTEGER,
INTENT(IN)::AerodynamicResistanceMethod
89 INTEGER,
INTENT(IN)::Diagnose
90 INTEGER,
INTENT(IN)::DiagQN
91 INTEGER,
INTENT(IN)::DiagQS
92 INTEGER,
INTENT(IN)::startDLS
93 INTEGER,
INTENT(IN)::endDLS
94 INTEGER,
INTENT(IN)::EmissionsMethod
95 INTEGER,
INTENT(IN)::Gridiv
96 INTEGER,
INTENT(IN)::gsModel
97 INTEGER,
INTENT(IN)::id
98 INTEGER,
INTENT(IN)::Ie_end
99 INTEGER,
INTENT(IN)::Ie_start
100 INTEGER,
INTENT(IN)::isec
101 INTEGER,
INTENT(IN)::imin
102 INTEGER,
INTENT(IN)::it
103 INTEGER,
INTENT(IN)::EvapMethod
104 INTEGER,
INTENT(IN)::iy
105 INTEGER,
INTENT(IN)::LAICalcYes
106 INTEGER,
INTENT(IN)::NetRadiationMethod
107 INTEGER,
INTENT(IN)::OHMIncQF
108 INTEGER,
INTENT(IN)::RoughLenHeatMethod
109 INTEGER,
INTENT(IN)::RoughLenMomMethod
110 INTEGER,
INTENT(IN)::SMDMethod
111 INTEGER,
INTENT(IN)::snowUse
112 INTEGER,
INTENT(IN)::StabilityMethod
113 INTEGER,
INTENT(IN)::StorageHeatMethod
114 INTEGER,
INTENT(IN)::tstep
115 INTEGER,
INTENT(IN)::tstep_prev
116 INTEGER,
INTENT(in)::dt_since_start
117 INTEGER,
INTENT(IN)::veg_type
118 INTEGER,
INTENT(IN)::WaterUseMethod
120 REAL(KIND(1D0)),
INTENT(IN)::AlbMax_DecTr
121 REAL(KIND(1D0)),
INTENT(IN)::AlbMax_EveTr
122 REAL(KIND(1D0)),
INTENT(IN)::AlbMax_Grass
123 REAL(KIND(1D0)),
INTENT(IN)::AlbMin_DecTr
124 REAL(KIND(1D0)),
INTENT(IN)::AlbMin_EveTr
125 REAL(KIND(1D0)),
INTENT(IN)::AlbMin_Grass
126 REAL(KIND(1D0)),
INTENT(IN)::alt
127 REAL(KIND(1D0)),
INTENT(IN)::avkdn
128 REAL(KIND(1D0)),
INTENT(IN)::avRh
129 REAL(KIND(1D0)),
INTENT(IN)::avU1
130 REAL(KIND(1D0)),
INTENT(IN)::BaseTHDD
131 REAL(KIND(1D0)),
INTENT(IN)::bldgH
132 REAL(KIND(1D0)),
INTENT(IN)::CapMax_dec
133 REAL(KIND(1D0)),
INTENT(IN)::CapMin_dec
134 REAL(KIND(1D0)),
INTENT(IN)::CO2PointSource
135 REAL(KIND(1D0)),
INTENT(IN)::CRWmax
136 REAL(KIND(1D0)),
INTENT(IN)::CRWmin
137 REAL(KIND(1D0)),
INTENT(IN)::DecTreeH
138 REAL(KIND(1D0)),
INTENT(IN)::DRAINRT
139 REAL(KIND(1D0)),
INTENT(IN)::EF_umolCO2perJ
140 REAL(KIND(1D0)),
INTENT(IN)::EnEF_v_Jkm
141 REAL(KIND(1D0)),
INTENT(IN)::EveTreeH
142 REAL(KIND(1D0)),
INTENT(IN)::FAIBldg
143 REAL(KIND(1D0)),
INTENT(IN)::FAIDecTree
144 REAL(KIND(1D0)),
INTENT(IN)::FAIEveTree
145 REAL(KIND(1D0)),
INTENT(IN)::Faut
146 REAL(KIND(1D0)),
INTENT(IN)::fcld_obs
147 REAL(KIND(1D0)),
INTENT(IN)::FlowChange
148 REAL(KIND(1D0)),
INTENT(IN)::FrFossilFuel_Heat
149 REAL(KIND(1D0)),
INTENT(IN)::FrFossilFuel_NonHeat
150 REAL(KIND(1D0)),
INTENT(IN)::G1
151 REAL(KIND(1D0)),
INTENT(IN)::G2
152 REAL(KIND(1D0)),
INTENT(IN)::G3
153 REAL(KIND(1D0)),
INTENT(IN)::G4
154 REAL(KIND(1D0)),
INTENT(IN)::G5
155 REAL(KIND(1D0)),
INTENT(IN)::G6
156 REAL(KIND(1D0)),
INTENT(IN)::InternalWaterUse_h
157 REAL(KIND(1D0)),
INTENT(IN)::IrrFracConif
158 REAL(KIND(1D0)),
INTENT(IN)::IrrFracDecid
159 REAL(KIND(1D0)),
INTENT(IN)::IrrFracGrass
160 REAL(KIND(1D0)),
INTENT(IN)::Kmax
161 REAL(KIND(1D0)),
INTENT(IN)::LAI_obs
162 REAL(KIND(1D0)),
INTENT(IN)::lat
163 REAL(KIND(1D0)),
INTENT(IN)::ldown_obs
164 REAL(KIND(1D0)),
INTENT(IN)::lng
165 REAL(KIND(1D0)),
INTENT(IN)::MaxFCMetab
166 REAL(KIND(1D0)),
INTENT(IN)::MaxQFMetab
167 REAL(KIND(1D0)),
INTENT(IN)::MinFCMetab
168 REAL(KIND(1D0)),
INTENT(IN)::MinQFMetab
169 REAL(KIND(1D0)),
INTENT(IN)::NARP_EMIS_SNOW
170 REAL(KIND(1D0)),
INTENT(IN)::NARP_TRANS_SITE
171 REAL(KIND(1D0)),
INTENT(IN)::PipeCapacity
172 REAL(KIND(1D0)),
INTENT(IN)::PopDensNighttime
173 REAL(KIND(1D0)),
INTENT(IN)::PorMax_dec
174 REAL(KIND(1D0)),
INTENT(IN)::PorMin_dec
175 REAL(KIND(1D0)),
INTENT(IN)::Precip
176 REAL(KIND(1D0)),
INTENT(IN)::PrecipLimit
177 REAL(KIND(1D0)),
INTENT(IN)::PrecipLimitAlb
178 REAL(KIND(1D0)),
INTENT(IN)::Press_hPa
179 REAL(KIND(1D0)),
INTENT(IN)::qn1_obs
180 REAL(KIND(1D0)),
INTENT(IN)::qs_obs
181 REAL(KIND(1D0)),
INTENT(IN)::qf_obs
182 REAL(KIND(1D0)),
INTENT(IN)::RadMeltFact
183 REAL(KIND(1D0)),
INTENT(IN)::RAINCOVER
184 REAL(KIND(1D0)),
INTENT(IN)::RainMaxRes
185 REAL(KIND(1D0)),
INTENT(IN)::RunoffToWater
186 REAL(KIND(1D0)),
INTENT(IN)::S1
187 REAL(KIND(1D0)),
INTENT(IN)::S2
188 REAL(KIND(1D0)),
INTENT(IN)::SnowAlbMax
189 REAL(KIND(1D0)),
INTENT(IN)::SnowAlbMin
190 REAL(KIND(1D0)),
INTENT(IN)::SnowDensMax
191 REAL(KIND(1D0)),
INTENT(IN)::SnowDensMin
192 REAL(KIND(1D0)),
INTENT(IN)::SnowLimBldg
193 REAL(KIND(1D0)),
INTENT(IN)::SnowLimPaved
194 REAL(KIND(1D0)),
INTENT(IN)::snowFrac_obs
195 REAL(KIND(1D0)),
INTENT(IN)::SurfaceArea
196 REAL(KIND(1D0)),
INTENT(IN)::tau_a
197 REAL(KIND(1D0)),
INTENT(IN)::tau_f
198 REAL(KIND(1D0)),
INTENT(IN)::tau_r
199 REAL(KIND(1D0)),
INTENT(IN)::Temp_C
200 REAL(KIND(1D0)),
INTENT(IN)::TempMeltFact
201 REAL(KIND(1D0)),
INTENT(IN)::TH
202 REAL(KIND(1D0)),
INTENT(IN)::timezone
203 REAL(KIND(1D0)),
INTENT(IN)::TL
204 REAL(KIND(1D0)),
INTENT(IN)::TrafficUnits
205 REAL(KIND(1D0)),
INTENT(IN)::wu_m3
206 REAL(KIND(1D0)),
INTENT(IN)::xsmd
207 REAL(KIND(1D0)),
INTENT(IN)::Z
208 REAL(KIND(1D0)),
INTENT(IN)::z0m_in
209 REAL(KIND(1D0)),
INTENT(IN)::zdm_in
211 INTEGER,
DIMENSION(NVEGSURF),
INTENT(IN)::LAIType
213 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::AH_MIN
214 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::AH_SLOPE_Cooling
215 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::AH_SLOPE_Heating
216 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::FcEF_v_kgkm
217 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::QF0_BEU
218 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::Qf_A
219 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::Qf_B
220 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::Qf_C
221 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::PopDensDaytime
222 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::T_CRITIC_Cooling
223 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::T_CRITIC_Heating
224 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::TrafficRate
225 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) ::Ie_a
226 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) ::Ie_m
227 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) ::MaxConductance
228 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) ::DayWat
229 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) ::DayWatPer
230 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) ::OHM_threshSW
231 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) ::OHM_threshWD
232 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::chAnOHM
233 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::cpAnOHM
234 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::emis
235 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::kkAnOHM
236 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SatHydraulicConduct
237 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::sfr
238 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SnowPackLimit
239 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SoilDepth
240 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SoilStoreCap
241 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::StateLimit
242 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::WetThresh
243 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::alpha_bioCO2
244 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::alpha_enh_bioCO2
245 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::BaseT
246 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::BaseTe
247 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::beta_bioCO2
248 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::beta_enh_bioCO2
249 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::GDDFull
250 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::LAIMax
251 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::LAIMin
252 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::min_res_bioCO2
253 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::resp_a
254 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::resp_b
255 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::SDDFull
256 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::SnowProf_24hr
257 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::theta_bioCO2
258 REAL(KIND(1D0)),
DIMENSION(4, NVEGSURF),
INTENT(IN) ::LAIPower
259 REAL(KIND(1D0)),
DIMENSION(nsurf + 1, 4, 3),
INTENT(IN) ::OHM_coef
260 REAL(KIND(1D0)),
DIMENSION(NSURF + 1, NSURF - 1),
INTENT(IN) ::WaterDist
261 REAL(KIND(1d0)),
DIMENSION(:),
INTENT(IN) ::Ts5mindata_ir
262 REAL(KIND(1D0)),
DIMENSION(:, :),
INTENT(IN) ::MetForcingData_grid
265 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::AHProf_24hr
266 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::HumActivity_24hr
267 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::PopProf_24hr
268 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::TraffProf_24hr
269 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::WUProfA_24hr
270 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::WUProfM_24hr
277 REAL(KIND(1d0)),
INTENT(INOUT) ::qn1_av
278 REAL(KIND(1d0)),
INTENT(INOUT) ::dqndt
279 REAL(KIND(1d0)),
INTENT(INOUT) ::qn1_s_av
280 REAL(KIND(1d0)),
INTENT(INOUT) ::dqnsdt
283 REAL(KIND(1D0)),
INTENT(INOUT) ::SnowfallCum
284 REAL(KIND(1D0)),
INTENT(INOUT) ::SnowAlb
285 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT)::IceFrac
286 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT)::SnowWater
287 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT)::SnowDens
288 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT)::SnowFrac
289 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT)::SnowPack
292 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::soilstore_id
293 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::state_id
294 REAL(KIND(1D0)),
DIMENSION(6, NSURF),
INTENT(INOUT)::StoreDrainPrm
297 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::alb
298 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(INOUT)::GDD_id
299 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(INout)::SDD_id
300 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(INOUT)::LAI_id
301 REAL(KIND(1d0)),
INTENT(INout)::Tmin_id
302 REAL(KIND(1d0)),
INTENT(INout)::Tmax_id
303 REAL(KIND(1d0)),
INTENT(INout)::lenDay_id
304 REAL(KIND(1d0)),
INTENT(INOUT)::DecidCap_id
305 REAL(KIND(1d0)),
INTENT(INOUT)::albDecTr_id
306 REAL(KIND(1d0)),
INTENT(INOUT)::albEveTr_id
307 REAL(KIND(1d0)),
INTENT(INOUT)::albGrass_id
308 REAL(KIND(1d0)),
INTENT(INOUT)::porosity_id
311 REAL(KIND(1d0)),
DIMENSION(12),
INTENT(INOUT) ::HDD_id
314 REAL(KIND(1d0)),
DIMENSION(9),
INTENT(INOUT) ::WUDay_id
317 REAL(KIND(1d0)),
INTENT(INOUT) ::Tair_av
322 REAL(KIND(1D0)),
DIMENSION(5),
INTENT(OUT) ::datetimeLine
323 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSUEWS - 5),
INTENT(OUT) ::dataOutLineSUEWS
324 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSnow - 5),
INTENT(OUT) ::dataOutLineSnow
325 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutESTM - 5),
INTENT(OUT) ::dataOutLineESTM
326 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutRSL - 5),
INTENT(OUT) ::dataoutLineRSL
327 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutDailyState - 5),
INTENT(OUT)::DailyStateLine
335 REAL(KIND(1D0))::AdditionalWater
336 REAL(KIND(1D0))::U10_ms
337 REAL(KIND(1D0))::azimuth
338 REAL(KIND(1D0))::chSnow_per_interval
340 REAL(KIND(1D0))::dens_dry
341 REAL(KIND(1d0))::deltaLAI
342 REAL(KIND(1D0))::drain_per_tstep
343 REAL(KIND(1D0))::Ea_hPa
344 REAL(KIND(1D0))::QE_LUMPS
345 REAL(KIND(1D0))::es_hPa
346 REAL(KIND(1D0))::ev_per_tstep
347 REAL(KIND(1D0))::ext_wu
349 REAL(KIND(1D0))::Fc_anthro
350 REAL(KIND(1D0))::Fc_biogen
351 REAL(KIND(1D0))::Fc_build
352 REAL(KIND(1D0))::fcld
353 REAL(KIND(1D0))::Fc_metab
354 REAL(KIND(1D0))::Fc_photo
355 REAL(KIND(1D0))::Fc_point
356 REAL(KIND(1D0))::Fc_respi
357 REAL(KIND(1D0))::Fc_traff
358 REAL(KIND(1D0))::gfunc
360 REAL(KIND(1D0))::QH_LUMPS
361 REAL(KIND(1D0))::int_wu
362 REAL(KIND(1D0))::kclear
364 REAL(KIND(1D0))::ldown
366 REAL(KIND(1D0))::L_mod
368 REAL(KIND(1D0))::mwstore
369 REAL(KIND(1D0))::NWstate_per_tstep
370 REAL(KIND(1D0))::planF
372 REAL(KIND(1D0))::q2_gkg
375 REAL(KIND(1D0))::QF_SAHP
377 REAL(KIND(1D0))::qh_residual
378 REAL(KIND(1D0))::qh_resist
380 REAL(KIND(1D0))::QmFreez
381 REAL(KIND(1D0))::QmRain
383 REAL(KIND(1D0))::qn1_S
384 REAL(KIND(1D0))::qn1_snowfree
387 REAL(KIND(1D0))::ResistSurf
389 REAL(KIND(1d0))::runoffAGveg
390 REAL(KIND(1d0))::runoffAGimpervious
391 REAL(KIND(1D0))::runoff_per_tstep
392 REAL(KIND(1D0))::runoffPipes
393 REAL(KIND(1D0))::runoffSoil_per_tstep
394 REAL(KIND(1D0))::runoffwaterbody
396 REAL(KIND(1D0))::SoilState
397 REAL(KIND(1D0))::state_per_tstep
398 REAL(KIND(1D0))::surf_chang_per_tstep
400 REAL(KIND(1D0))::t2_C
401 REAL(KIND(1D0))::TSfc_C
402 REAL(KIND(1D0))::TempVeg
403 REAL(KIND(1D0))::tot_chang_per_tstep
404 REAL(KIND(1D0))::TStar
405 REAL(KIND(1D0))::tsurf
406 REAL(KIND(1D0))::UStar
407 REAL(KIND(1D0))::VPD_Pa
408 REAL(KIND(1D0))::wu_DecTr
409 REAL(KIND(1D0))::wu_EveTr
410 REAL(KIND(1D0))::wu_Grass
413 REAL(KIND(1D0))::ZENITH_deg
416 REAL(KIND(1D0)),
DIMENSION(2)::SnowRemoval
417 REAL(KIND(1D0)),
DIMENSION(NSURF)::FreezMelt
418 REAL(KIND(1d0)),
DIMENSION(nsurf)::kup_ind_snow
419 REAL(KIND(1D0)),
DIMENSION(NSURF)::mw_ind
420 REAL(KIND(1D0)),
DIMENSION(NSURF)::Qm_freezState
421 REAL(KIND(1D0)),
DIMENSION(NSURF)::Qm_melt
422 REAL(KIND(1D0)),
DIMENSION(NSURF)::Qm_rain
423 REAL(KIND(1D0)),
DIMENSION(NSURF)::qn1_ind_snow
424 REAL(KIND(1D0)),
DIMENSION(NSURF)::rainOnSnow
425 REAL(KIND(1D0)),
DIMENSION(NSURF)::runoffSoil
426 REAL(KIND(1D0)),
DIMENSION(NSURF)::smd_nsurf
427 REAL(KIND(1D0)),
DIMENSION(NSURF)::snowDepth
429 REAL(KIND(1d0)),
DIMENSION(nsurf)::Tsurf_ind_snow
431 INTEGER,
DIMENSION(NSURF)::snowCalcSwitch
432 INTEGER,
DIMENSION(3) ::dayofWeek_id
439 REAL(KIND(1D0))::lvS_J_kg
440 REAL(KIND(1D0))::psyc_hPa
441 REAL(KIND(1D0))::RAsnow
443 REAL(KIND(1D0))::runoff_per_interval
444 REAL(KIND(1D0))::s_hPa
445 REAL(KIND(1D0))::sIce_hpa
446 REAL(KIND(1D0))::SoilMoistCap
447 REAL(KIND(1D0))::veg_fr
448 REAL(KIND(1D0))::VegPhenLumps
449 REAL(KIND(1D0))::VPd_hpa
450 REAL(KIND(1D0))::vsmd
453 REAL(KIND(1D0)),
DIMENSION(NSURF)::deltaQi
454 REAL(KIND(1D0)),
DIMENSION(NSURF)::drain
455 REAL(KIND(1D0)),
DIMENSION(NSURF)::FreezState
456 REAL(KIND(1D0)),
DIMENSION(NSURF)::FreezStateVol
457 REAL(KIND(1D0)),
DIMENSION(NSURF)::soilstoreOld
458 REAL(KIND(1D0)),
DIMENSION(NSURF)::stateOld
459 REAL(KIND(1D0)),
DIMENSION(NSURF)::tsurf_ind
464 REAL(KIND(1D0))::addImpervious = 0
465 REAL(KIND(1D0))::addPipes = 0
466 REAL(KIND(1D0))::addVeg = 0
467 REAL(KIND(1D0))::addWaterBody = 0
468 REAL(KIND(1D0)),
DIMENSION(NSURF)::AddWater = 0
469 REAL(KIND(1D0)),
DIMENSION(NSURF)::AddWaterRunoff = 0
473 REAL(KIND(1D0))::nsh_real
474 REAL(KIND(1D0))::tstep_real
475 REAL(KIND(1D0))::dectime
478 REAL(KIND(1D0))::VegFraction
479 REAL(KIND(1D0))::ImpervFraction
480 REAL(KIND(1D0))::PervFraction
481 REAL(KIND(1D0))::NonWaterFraction
485 REAL(KIND(1D0))::alb1
495 REAL(KIND(1D0))::qn1_av_prev, qn1_av_next
496 REAL(KIND(1D0))::dqndt_prev, dqndt_next
497 REAL(KIND(1D0))::qn1_s_av_prev, qn1_s_av_next
498 REAL(KIND(1D0))::dqnsdt_prev, dqnsdt_next
501 REAL(KIND(1D0))::SnowfallCum_prev, SnowfallCum_next
502 REAL(KIND(1D0))::SnowAlb_prev, SnowAlb_next
504 REAL(KIND(1D0)),
DIMENSION(NSURF)::IceFrac_prev, IceFrac_next
505 REAL(KIND(1D0)),
DIMENSION(NSURF)::SnowWater_prev, SnowWater_next
506 REAL(KIND(1D0)),
DIMENSION(NSURF)::SnowDens_prev, SnowDens_next
507 REAL(KIND(1D0)),
DIMENSION(NSURF)::SnowFrac_prev, SnowFrac_next
508 REAL(KIND(1D0)),
DIMENSION(NSURF)::SnowPack_prev, SnowPack_next
511 REAL(KIND(1D0)),
DIMENSION(NSURF) ::soilstore_id_prev, soilstore_id_next
512 REAL(KIND(1D0)),
DIMENSION(NSURF) ::state_id_prev, state_id_next
513 REAL(KIND(1D0)),
DIMENSION(6, NSURF)::StoreDrainPrm_prev, StoreDrainPrm_next
516 REAL(KIND(1D0)),
DIMENSION(NSURF) ::alb_prev, alb_next
517 REAL(KIND(1d0)),
DIMENSION(nvegsurf)::GDD_id_prev, GDD_id_next
518 REAL(KIND(1d0)),
DIMENSION(nvegsurf)::LAI_id_prev, LAI_id_next
519 REAL(KIND(1d0)),
DIMENSION(nvegsurf)::SDD_id_prev, SDD_id_next
521 REAL(KIND(1D0))::DecidCap_id_prev, DecidCap_id_next
522 REAL(KIND(1D0))::albDecTr_id_prev, albDecTr_id_next
523 REAL(KIND(1D0))::albEveTr_id_prev, albEveTr_id_next
524 REAL(KIND(1D0))::albGrass_id_prev, albGrass_id_next
525 REAL(KIND(1D0))::porosity_id_prev, porosity_id_next
527 REAL(KIND(1d0))::Tmin_id_prev, Tmin_id_next
528 REAL(KIND(1d0))::Tmax_id_prev, Tmax_id_next
529 REAL(KIND(1d0))::lenDay_id_prev, lenDay_id_next
532 REAL(KIND(1d0)),
DIMENSION(12)::HDD_id_prev, HDD_id_next
535 REAL(KIND(1d0)),
DIMENSION(9)::WUDay_id_prev, WUDay_id_next
537 REAL(KIND(1D0))::Tair_av_prev, Tair_av_next
541 INTEGER,
PARAMETER :: nz = 90
544 logical:: flag_converge
545 REAL(KIND(1D0)):: Ts_iter
546 REAL(KIND(1D0)):: L_mod_iter
547 REAL(KIND(1D0)):: QH_Init
550 REAL(KIND(1d0)),
DIMENSION(30):: psihatm_z
551 REAL(KIND(1d0)),
DIMENSION(30):: psihath_z
557 qn1_s_av_prev = qn1_s_av
559 snowfallcum_prev = snowfallcum
560 snowalb_prev = snowalb
561 icefrac_prev = icefrac
562 snowwater_prev = snowwater
563 snowdens_prev = snowdens
564 snowfrac_prev = snowfrac
565 snowpack_prev = snowpack
566 soilstore_id_prev = soilstore_id
567 state_id_prev = state_id
568 tair_av_prev = tair_av
572 tmin_id_prev = tmin_id
573 tmax_id_prev = tmax_id
574 lenday_id_prev = lenday_id
575 storedrainprm_prev = storedrainprm
576 decidcap_id_prev = decidcap_id
577 porosity_id_prev = porosity_id
579 albdectr_id_prev = albdectr_id
580 albevetr_id_prev = albevetr_id
581 albgrass_id_prev = albgrass_id
583 wuday_id_prev = wuday_id
588 qn1_s_av_next = qn1_s_av
590 snowfallcum_next = snowfallcum
591 snowalb_next = snowalb
592 icefrac_next = icefrac
593 snowwater_next = snowwater
594 snowdens_next = snowdens
595 snowfrac_next = snowfrac
596 snowpack_next = snowpack
597 soilstore_id_next = soilstore_id
598 state_id_next = state_id
599 tair_av_next = tair_av
603 tmin_id_next = tmin_id
604 tmax_id_next = tmax_id
605 lenday_id_next = lenday_id
606 storedrainprm_next = storedrainprm
607 decidcap_id_next = decidcap_id
608 porosity_id_next = porosity_id
610 albdectr_id_next = albdectr_id
611 albevetr_id_next = albevetr_id
612 albgrass_id_next = albgrass_id
614 wuday_id_next = wuday_id
621 flag_converge = .false.
625 do while (.not. flag_converge)
629 id, it, imin, isec, &
635 nsh, nsh_real, tstep_real)
640 vegfraction, impervfraction, pervfraction, nonwaterfraction)
649 id, startdls, enddls, &
653 tair_av_next =
cal_tair_av(tair_av_prev, dt_since_start, tstep, temp_c)
658 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_RoughnessParameters...' 659 IF (diagnose == 1) print *,
'z0m_in =', z0m_in
661 roughlenmommethod, sfr, &
662 bldgh, evetreeh, dectreeh, &
663 porosity_id_prev, faibldg, faievetree, faidectree, &
669 IF (diagnose == 1)
WRITE (*, *)
'Calling NARP_cal_SunPosition...' 671 REAL(iy, KIND(1d0)), &
673 timezone, lat, lng, alt, &
677 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_DailyState...' 678 CALL suews_cal_dailystate( &
679 iy, id, it, imin, isec, tstep, tstep_prev, dt_since_start, dayofweek_id, &
680 tmin_id_prev, tmax_id_prev, lenday_id_prev, &
681 waterusemethod, ie_start, ie_end, &
682 laicalcyes, laitype, &
683 nsh_real, avkdn, temp_c, precip, basethdd, &
684 lat, faut, lai_obs, &
685 albmax_dectr, albmax_evetr, albmax_grass, &
686 albmin_dectr, albmin_evetr, albmin_grass, &
687 capmax_dec, capmin_dec, pormax_dec, pormin_dec, &
688 ie_a, ie_m, daywatper, daywat, &
689 baset, basete, gddfull, sddfull, laimin, laimax, laipower, &
690 decidcap_id_prev, storedrainprm_prev, lai_id_prev, gdd_id_prev, sdd_id_prev, &
691 albdectr_id_prev, albevetr_id_prev, albgrass_id_prev, porosity_id_prev, &
694 tmin_id_next, tmax_id_next, lenday_id_next, &
695 albdectr_id_next, albevetr_id_next, albgrass_id_next, porosity_id_next, &
696 decidcap_id_next, storedrainprm_next, lai_id_next, gdd_id_next, sdd_id_next, deltalai, wuday_id_next)
699 IF (diagnose == 1)
WRITE (*, *)
'Calling LUMPS_cal_AtmMoist...' 701 temp_c, press_hpa, avrh, dectime, &
703 es_hpa, ea_hpa, vpd_hpa, vpd_pa, dq, dens_dry,
avcp,
avdens)
706 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_update_SoilMoist...' 709 soilstorecap, sfr, soilstore_id_prev, &
710 soilmoistcap, soilstate, &
713 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_WaterUse...' 717 wu_m3, surfacearea, sfr, &
718 irrfracconif, irrfracdecid, irrfracgrass, &
719 dayofweek_id, wuprofa_24hr, wuprofm_24hr, &
720 internalwateruse_h, hdd_id_next, wuday_id_next, &
721 waterusemethod, nsh, it, imin, dls, &
722 wu_evetr, wu_dectr, wu_grass, int_wu, ext_wu)
726 ah_min, ahprof_24hr, ah_slope_cooling, ah_slope_heating, co2pointsource, &
727 dayofweek_id, diagnose, dls, ef_umolco2perj, emissionsmethod, enef_v_jkm, &
728 fcef_v_kgkm, frfossilfuel_heat, frfossilfuel_nonheat, hdd_id_next, humactivity_24hr, &
729 id, imin, it, maxfcmetab, maxqfmetab, minfcmetab, minqfmetab, nsh, &
730 popdensdaytime, popdensnighttime, popprof_24hr, qf, qf0_beu, qf_a, qf_b, qf_c, &
731 qf_obs, qf_sahp, sfr, snowfrac, surfacearea, t_critic_cooling, t_critic_heating, &
732 temp_c, trafficrate, trafficunits, traffprof_24hr, &
733 fc_anthro, fc_build, fc_metab, fc_point, fc_traff)
740 netradiationmethod, snowuse, &
741 tstep, snowpack_prev, tau_a, tau_f, snowalbmax, snowalbmin, &
742 diagnose, snowfrac_obs, ldown_obs, fcld_obs, &
743 dectime, zenith_deg, ts_iter, avkdn, temp_c, avrh, ea_hpa, qn1_obs, &
744 snowalb_prev, snowfrac_prev, diagqn, &
745 narp_trans_site, narp_emis_snow, icefrac_prev, sfr, emis, &
746 alb_prev, albdectr_id_next, albevetr_id_next, albgrass_id_next, &
749 qn1, qn1_snowfree, qn1_s, kclear, kup, lup, tsurf, &
750 qn1_ind_snow, kup_ind_snow, tsurf_ind_snow, tsurf_ind, &
751 alb1, snowfrac_next, snowalb_next)
755 storageheatmethod, qs_obs, ohmincqf, gridiv, &
756 id, tstep, dt_since_start, diagnose, sfr, &
757 ohm_coef, ohm_threshsw, ohm_threshwd, &
758 soilstore_id_prev, soilstorecap, state_id_prev, snowuse, snowfrac_next, diagqs, &
759 hdd_id_next, metforcingdata_grid, ts5mindata_ir, qf, qn1, &
760 avkdn, avu1, temp_c, zenith_deg, avrh, press_hpa, ldown, &
761 bldgh, alb_next, emis, cpanohm, kkanohm, chanohm, emissionsmethod, &
762 tair_av_next, qn1_av_prev, dqndt_prev, qn1_s_av_prev, dqnsdt_prev, &
763 storedrainprm_next, &
764 qn1_s, dataoutlineestm, qs, &
765 qn1_av_next, dqndt_next, qn1_s_av_next, dqnsdt_next, &
769 IF (diagnose == 1)
WRITE (*, *)
'Calling MeltHeat' 773 tstep, tau_r, snowdensmax, &
774 lvs_j_kg,
lv_j_kg, tstep_real, radmeltfact, tempmeltfact, snowalbmax, &
775 snowdensmin, temp_c, precip, preciplimit, preciplimitalb, &
776 nsh_real, sfr, tsurf_ind, tsurf_ind_snow, state_id_prev, qn1_ind_snow, &
777 kup_ind_snow, snowwater_prev, deltaqi, alb1, &
778 snowpack_prev, snowfrac_next, snowalb_next, snowdens_prev, snowfallcum_prev, &
779 snowpack_next, snowfrac_next, snowalb_next, snowdens_next, snowfallcum_next, &
780 mwh, qm, qmfreez, qmrain, &
781 veg_fr, snowcalcswitch, qm_melt, qm_freezstate, qm_rain, freezmelt, &
782 freezstate, freezstatevol, rainonsnow, snowdepth, mw_ind, &
786 IF (diagnose == 1)
WRITE (*, *)
'Calling LUMPS_cal_QHQE...' 787 if (i_iter == 1)
then 791 snowuse, qn1, qf, qs, qm, temp_c, veg_fr,
avcp, press_hpa,
lv_j_kg, &
792 tstep_real, drainrt, nsh_real, &
793 precip, rainmaxres, raincover, sfr, lai_id_next, laimax, laimin, &
795 qe_lumps, psyc_hpa, s_hpa, sice_hpa, tempveg, vegphenlumps)
807 snowuse, nonwaterfraction, addpipes, addimpervious, addveg, addwaterbody, &
808 state_id_prev, soilstore_id_prev, sfr, storedrainprm_next, waterdist, nsh_real, &
810 drain, addwaterrunoff, &
811 additionalwater, runoffpipes, runoff_per_interval, &
812 addwater, stateold, soilstoreold)
818 diagnose, aerodynamicresistancemethod, roughlenheatmethod, snowuse, &
819 id, it, gsmodel, smdmethod, &
821 avu1, temp_c, vegfraction, avkdn, &
826 dq, xsmd, vsmd, maxconductance, laimax, lai_id_next, snowfrac_next, sfr, &
827 ustar, tstar, l_mod, &
828 zl, gsc, resistsurf, ra, rasnow, rb)
833 tstep, imin, it, evapmethod, snowcalcswitch, dayofweek_id, crwmin, crwmax, &
834 dectime, lvs_j_kg, avrh, press_hpa, temp_c, &
835 rasnow, psyc_hpa, sice_hpa, &
836 pervfraction, vegfraction, addimpervious, qn1_snowfree, qf, qs, vpd_hpa, s_hpa, &
837 resistsurf, ra, rb, snowdensmin, precip, pipecapacity, runofftowater, &
838 nonwaterfraction, wu_evetr, wu_dectr, wu_grass, addveg, addwaterbody, snowlimpaved, snowlimbldg, &
839 surfacearea, flowchange, drain, wetthresh, stateold, mw_ind, soilstorecap, rainonsnow, &
840 freezmelt, freezstate, freezstatevol, qm_melt, qm_rain, tsurf_ind, sfr, &
841 statelimit, addwater, addwaterrunoff, storedrainprm_next, snowpacklimit, snowprof_24hr, &
842 snowpack_next, snowfrac_next, snowwater_prev, icefrac_prev, snowdens_next, &
843 runoff_per_interval, state_id_prev, soilstore_id_prev, &
844 state_id_next, soilstore_id_next, &
845 snowpack_next, snowfrac_next, snowwater_next, icefrac_next, snowdens_next, &
848 state_per_tstep, nwstate_per_tstep, qe, &
849 swe, chsnow_per_interval, ev_per_tstep, runoff_per_tstep, &
850 surf_chang_per_tstep, runoffpipes, mwstore, runoffwaterbody, &
851 runoffagveg, runoffagimpervious)
855 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_QH...' 858 qn1, qf, qmrain, qe, qs, qmfreez, qm,
avdens,
avcp, tsurf, temp_c, ra, &
859 qh, qh_residual, qh_resist)
867 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_HorizontalSoilWater...' 872 sathydraulicconduct, &
878 runoffsoil_per_tstep &
882 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_SoilState...' 884 smdmethod, xsmd, nonwaterfraction, soilmoistcap, &
885 soilstorecap, surf_chang_per_tstep, &
886 soilstore_id_next, soilstoreold, sfr, &
887 smd, smd_nsurf, tot_chang_per_tstep, soilstate)
928 if (abs(qh - qh_init) > 0.1)
then 929 flag_converge = .false.
931 flag_converge = .true.
935 if (i_iter > 100)
exit 947 IF (diagnose == 1)
WRITE (*, *)
'Calling RSLProfile...' 950 l_mod, sfr, planf, stabilitymethod, &
952 avu1, temp_c, avrh, press_hpa, z, qh, qe, &
953 t2_c, q2_gkg, u10_ms, rh2, &
958 alpha_bioco2, alpha_enh_bioco2, avkdn, avrh, beta_bioco2, beta_enh_bioco2,
bsoilsurf, &
959 conifsurf,
decidsurf, dectime, diagnose, emissionsmethod, fc_anthro, g1, g2, g3, g4, &
960 g5, g6, gfunc,
grasssurf, gsmodel, id, it,
ivconif,
ivdecid,
ivgrass, kmax, lai_id_next, laimin, &
961 laimax, maxconductance, min_res_bioco2,
nsurf,
nvegsurf, press_hpa, resp_a, &
962 resp_b, s1, s2, sfr, smdmethod, snowfrac, t2_c, temp_c, theta_bioco2, th, tl, vsmd, xsmd, &
963 fc, fc_biogen, fc_photo, fc_respi)
972 qn1_s_av = qn1_s_av_next
974 snowfallcum = snowfallcum_next
975 snowalb = snowalb_next
976 icefrac = icefrac_next
977 snowwater = snowwater_next
978 snowdens = snowdens_next
979 snowfrac = snowfrac_next
980 snowpack = snowpack_next
982 soilstore_id = soilstore_id_next
983 state_id = state_id_next
988 decidcap_id = decidcap_id_next
989 albdectr_id = albdectr_id_next
990 albevetr_id = albevetr_id_next
991 albgrass_id = albgrass_id_next
992 porosity_id = porosity_id_next
993 storedrainprm = storedrainprm_next
994 tair_av = tair_av_next
995 tmin_id = tmin_id_next
996 tmax_id = tmax_id_next
997 lenday_id = lenday_id_next
999 wuday_id = wuday_id_next
1003 additionalwater, alb, avkdn, u10_ms, azimuth, &
1004 chsnow_per_interval, dectime, &
1005 drain_per_tstep, qe_lumps, ev_per_tstep, ext_wu, fc, fc_build, fcld, &
1006 fc_metab, fc_photo, fc_respi, fc_point, fc_traff, flowchange, &
1007 qh_lumps, id, imin, int_wu, it, iy, &
1008 kup, lai_id, ldown, l_mod, lup, mwh, &
1010 nsh_real, nwstate_per_tstep, precip, q2_gkg, &
1011 qe, qf, qh, qh_resist, qm, qmfreez, &
1012 qmrain, qn1, qn1_s, qn1_snowfree, qs, ra, &
1013 resistsurf, rh2, runoffagimpervious, runoffagveg, &
1014 runoff_per_tstep, runoffpipes, runoffsoil_per_tstep, &
1015 runoffwaterbody, sfr, smd, smd_nsurf, snowalb, snowremoval, &
1016 state_id_next, state_per_tstep, surf_chang_per_tstep, swe, t2_c, tsfc_c, &
1017 tot_chang_per_tstep, tsurf, ustar, wu_dectr, &
1018 wu_evetr, wu_grass, z0m, zdm, zenith_deg, &
1019 datetimeline, dataoutlinesuews)
1024 CALL update_dailystateline( &
1025 it, imin, nsh_real, &
1026 gdd_id, hdd_id, lai_id, &
1028 tmin_id, tmax_id, lenday_id, &
1035 deltalai, vegphenlumps, &
1036 snowalb, snowdens, &
1047 AH_MIN, AHProf_24hr, AH_SLOPE_Cooling, AH_SLOPE_Heating, CO2PointSource, &! input:
1048 dayofWeek_id, Diagnose, DLS, EF_umolCO2perJ, EmissionsMethod, EnEF_v_Jkm, &
1049 FcEF_v_kgkm, FrFossilFuel_Heat, FrFossilFuel_NonHeat, HDD_id, HumActivity_24hr, &
1050 id, imin, it, MaxFCMetab, MaxQFMetab, MinFCMetab, MinQFMetab, nsh, &
1051 PopDensDaytime, PopDensNighttime, PopProf_24hr, QF, QF0_BEU, Qf_A, Qf_B, Qf_C, &
1052 QF_obs, QF_SAHP, sfr, SnowFrac, SurfaceArea, T_CRITIC_Cooling, T_CRITIC_Heating, &
1053 Temp_C, TrafficRate, TrafficUnits, TraffProf_24hr, &
1054 Fc_anthro, Fc_build, Fc_metab, Fc_point, Fc_traff)
1058 INTEGER,
INTENT(in)::Diagnose
1059 INTEGER,
INTENT(in)::DLS
1060 INTEGER,
INTENT(in)::EmissionsMethod
1061 INTEGER,
INTENT(in)::id
1062 INTEGER,
INTENT(in)::it
1063 INTEGER,
INTENT(in)::imin
1064 INTEGER,
INTENT(in)::nsh
1065 INTEGER,
DIMENSION(3),
INTENT(in)::dayofWeek_id
1067 REAL(KIND(1d0)),
DIMENSION(6, 2),
INTENT(in)::HDD_id
1069 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::AH_MIN
1070 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::AH_SLOPE_Heating
1071 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::AH_SLOPE_Cooling
1072 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(in)::FcEF_v_kgkm
1074 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::PopDensDaytime
1075 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::QF0_BEU
1076 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::Qf_A
1077 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::Qf_B
1078 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::Qf_C
1079 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::T_CRITIC_Heating
1080 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::T_CRITIC_Cooling
1081 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(in)::TrafficRate
1083 REAL(KIND(1d0)),
DIMENSION(0:23, 2),
INTENT(in)::AHProf_24hr
1084 REAL(KIND(1d0)),
DIMENSION(0:23, 2),
INTENT(in)::HumActivity_24hr
1085 REAL(KIND(1d0)),
DIMENSION(0:23, 2),
INTENT(in)::TraffProf_24hr
1086 REAL(KIND(1d0)),
DIMENSION(0:23, 2),
INTENT(in)::PopProf_24hr
1088 REAL(KIND(1D0)),
INTENT(in)::CO2PointSource
1089 REAL(KIND(1D0)),
INTENT(in)::EF_umolCO2perJ
1090 REAL(KIND(1D0)),
INTENT(in)::EnEF_v_Jkm
1091 REAL(KIND(1D0)),
INTENT(in)::FrFossilFuel_Heat
1092 REAL(KIND(1D0)),
INTENT(in)::FrFossilFuel_NonHeat
1093 REAL(KIND(1D0)),
INTENT(in)::MaxFCMetab
1094 REAL(KIND(1D0)),
INTENT(in)::MaxQFMetab
1095 REAL(KIND(1D0)),
INTENT(in)::MinFCMetab
1096 REAL(KIND(1D0)),
INTENT(in)::MinQFMetab
1097 REAL(KIND(1D0)),
INTENT(in)::PopDensNighttime
1098 REAL(KIND(1D0)),
INTENT(in)::QF_obs
1099 REAL(KIND(1D0)),
INTENT(in)::Temp_C
1100 REAL(KIND(1D0)),
INTENT(in)::TrafficUnits
1102 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::sfr
1103 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowFrac
1104 REAL(KIND(1D0)),
INTENT(IN)::SurfaceArea
1106 REAL(KIND(1D0)),
INTENT(out)::Fc_anthro
1107 REAL(KIND(1D0)),
INTENT(out)::Fc_build
1108 REAL(KIND(1D0)),
INTENT(out)::Fc_metab
1109 REAL(KIND(1D0)),
INTENT(out)::Fc_point
1110 REAL(KIND(1D0)),
INTENT(out)::Fc_traff
1111 REAL(KIND(1D0)),
INTENT(out)::QF
1112 REAL(KIND(1D0)),
INTENT(out)::QF_SAHP
1114 INTEGER,
PARAMETER :: notUsedI = -999
1115 REAL(KIND(1D0)),
PARAMETER::notUsed = -999
1117 IF (emissionsmethod == 0)
THEN 1119 ELSEIF ((emissionsmethod > 0 .AND. emissionsmethod <= 6) .OR. emissionsmethod >= 11)
THEN 1121 co2pointsource, emissionsmethod, &
1122 id, it, imin, dls, nsh, dayofweek_id, &
1123 ef_umolco2perj, fcef_v_kgkm, enef_v_jkm, trafficunits, &
1124 frfossilfuel_heat, frfossilfuel_nonheat, &
1125 minfcmetab, maxfcmetab, minqfmetab, maxqfmetab, &
1126 popdensdaytime, popdensnighttime, &
1127 temp_c, hdd_id, qf_a, qf_b, qf_c, &
1128 ah_min, ah_slope_heating, ah_slope_cooling, &
1129 t_critic_heating, t_critic_cooling, &
1132 fc_anthro, fc_metab, fc_traff, fc_build, fc_point, &
1133 ahprof_24hr, humactivity_24hr, traffprof_24hr, popprof_24hr, surfacearea)
1136 CALL errorhint(73,
'RunControl.nml:EmissionsMethod unusable', notused, notused, emissionsmethod)
1139 IF (emissionsmethod >= 1) qf = qf_sahp
1141 IF (emissionsmethod >= 0 .AND. emissionsmethod <= 6)
THEN 1154 alpha_bioCO2, alpha_enh_bioCO2, avkdn, avRh, beta_bioCO2, beta_enh_bioCO2, BSoilSurf, &! input:
1155 ConifSurf, DecidSurf, dectime, Diagnose, EmissionsMethod, Fc_anthro, G1, G2, G3, G4, &
1156 G5, G6, gfunc, GrassSurf, gsmodel, id, it, ivConif, ivDecid, ivGrass, Kmax, LAI_id, LAIMin, &
1157 LAIMax, MaxConductance, min_res_bioCO2, nsurf, NVegSurf, Press_hPa, resp_a, &
1158 resp_b, S1, S2, sfr, SMDMethod, SnowFrac, t2_C, Temp_C, theta_bioCO2, TH, TL, vsmd, xsmd, &
1159 Fc, Fc_biogen, Fc_photo, Fc_respi)
1163 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::alpha_bioCO2
1164 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::alpha_enh_bioCO2
1165 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::beta_bioCO2
1166 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::beta_enh_bioCO2
1167 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::LAI_id
1168 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::LAIMin
1169 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::LAIMax
1170 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::min_res_bioCO2
1171 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::resp_a
1172 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::resp_b
1173 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(in)::theta_bioCO2
1175 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::sfr
1176 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowFrac
1178 REAL(KIND(1d0)),
DIMENSION(3),
INTENT(in)::MaxConductance
1180 INTEGER,
INTENT(in)::BSoilSurf
1181 INTEGER,
INTENT(in)::ConifSurf
1182 INTEGER,
INTENT(in)::DecidSurf
1183 INTEGER,
INTENT(in)::Diagnose
1184 INTEGER,
INTENT(in)::EmissionsMethod
1185 INTEGER,
INTENT(in)::GrassSurf
1186 INTEGER,
INTENT(in)::gsmodel
1187 INTEGER,
INTENT(in)::id
1188 INTEGER,
INTENT(in)::it
1189 INTEGER,
INTENT(in)::ivConif
1190 INTEGER,
INTENT(in)::ivDecid
1191 INTEGER,
INTENT(in)::ivGrass
1192 INTEGER,
INTENT(in)::nsurf
1193 INTEGER,
INTENT(in)::NVegSurf
1194 INTEGER,
INTENT(in)::SMDMethod
1196 REAL(KIND(1d0)),
INTENT(in)::avkdn
1197 REAL(KIND(1d0)),
INTENT(in)::avRh
1198 REAL(KIND(1d0)),
INTENT(in)::dectime
1199 REAL(KIND(1d0)),
INTENT(in)::Fc_anthro
1200 REAL(KIND(1d0)),
INTENT(in)::G1
1201 REAL(KIND(1d0)),
INTENT(in)::G2
1202 REAL(KIND(1d0)),
INTENT(in)::G3
1203 REAL(KIND(1d0)),
INTENT(in)::G4
1204 REAL(KIND(1d0)),
INTENT(in)::G5
1205 REAL(KIND(1d0)),
INTENT(in)::G6
1206 REAL(KIND(1d0)),
INTENT(in)::gfunc
1207 REAL(KIND(1d0)),
INTENT(in)::Kmax
1208 REAL(KIND(1d0)),
INTENT(in)::Press_hPa
1209 REAL(KIND(1d0)),
INTENT(in)::S1
1210 REAL(KIND(1d0)),
INTENT(in)::S2
1211 REAL(KIND(1d0)),
INTENT(in)::t2_C
1212 REAL(KIND(1d0)),
INTENT(in)::Temp_C
1213 REAL(KIND(1d0)),
INTENT(in)::TH
1214 REAL(KIND(1d0)),
INTENT(in)::TL
1215 REAL(KIND(1d0)),
INTENT(in)::vsmd
1216 REAL(KIND(1d0)),
INTENT(in)::xsmd
1218 REAL(KIND(1d0)),
INTENT(out)::Fc_biogen
1219 REAL(KIND(1d0)),
INTENT(out)::Fc_photo
1220 REAL(KIND(1d0)),
INTENT(out)::Fc_respi
1221 REAL(KIND(1d0)),
INTENT(out)::Fc
1223 REAL(KIND(1d0))::gfunc2
1226 REAL(KIND(1d0))::dummy1
1227 REAL(KIND(1d0))::dummy2
1228 REAL(KIND(1d0))::dummy3
1229 REAL(KIND(1d0))::dummy4
1230 REAL(KIND(1d0))::dummy5
1231 REAL(KIND(1d0))::dummy6
1232 REAL(KIND(1d0))::dummy7
1233 REAL(KIND(1d0))::dummy8
1234 REAL(KIND(1d0))::dummy9
1235 REAL(KIND(1d0))::dummy10
1236 REAL(KIND(1d0))::dummy11
1238 IF (emissionsmethod >= 11)
THEN 1240 IF (gsmodel == 3 .OR. gsmodel == 4)
THEN 1244 IF (abs(temp_c - t2_c) > 5)
THEN 1251 t2, press_hpa, avrh, dectime, &
1253 dummy3, dummy4, dummy5, dummy6, dq, dummy7, dummy8, dummy9)
1257 smdmethod, snowfrac, sfr, avkdn, t2, dq, xsmd, vsmd, maxconductance, &
1258 laimax, lai_id, gsmodel, kmax, &
1259 g1, g2, g3, g4, g5, g6, th, tl, s1, s2, &
1260 gfunc2, dummy10, dummy11)
1264 IF (diagnose == 1)
WRITE (*, *)
'Calling CO2_biogen...' 1266 alpha_bioco2, alpha_enh_bioco2, avkdn, beta_bioco2, beta_enh_bioco2, bsoilsurf, &
1267 conifsurf, decidsurf, dectime, emissionsmethod, gfunc, gfunc2, grasssurf, gsmodel, &
1268 id, it, ivconif, ivdecid, ivgrass, lai_id, laimin, laimax, min_res_bioco2, nsurf, &
1269 nvegsurf, resp_a, resp_b, sfr, snowfrac, t2, temp_c, theta_bioco2, &
1270 fc_biogen, fc_photo, fc_respi)
1273 IF (emissionsmethod >= 0 .AND. emissionsmethod <= 6)
THEN 1279 fc = fc_anthro + fc_biogen
1286 NetRadiationMethod, snowUse, &!input
1287 tstep, SnowPack_prev, tau_a, tau_f, SnowAlbMax, SnowAlbMin, &
1288 Diagnose, snowFrac_obs, ldown_obs, fcld_obs, &
1289 dectime, ZENITH_deg, Tsurf_0, avKdn, Temp_C, avRH, ea_hPa, qn1_obs, &
1290 SnowAlb_prev, snowFrac_prev, DiagQN, &
1291 NARP_TRANS_SITE, NARP_EMIS_SNOW, IceFrac, sfr, emis, &
1292 alb_prev, albDecTr_id, albEveTr_id, albGrass_id, &!input
1294 ldown, fcld, &!output
1295 qn1, qn1_snowfree, qn1_S, kclear, kup, lup, tsurf, &
1296 qn1_ind_snow, kup_ind_snow, Tsurf_ind_snow, Tsurf_ind, &
1297 alb1, snowFrac_next, SnowAlb_next)
1306 INTEGER,
INTENT(in)::NetRadiationMethod
1307 INTEGER,
INTENT(in)::snowUse
1308 INTEGER,
INTENT(in)::Diagnose
1309 INTEGER,
INTENT(in)::DiagQN
1310 INTEGER,
INTENT(in)::tstep
1312 REAL(KIND(1d0)),
INTENT(in)::snowFrac_obs
1313 REAL(KIND(1d0)),
INTENT(in)::ldown_obs
1314 REAL(KIND(1d0)),
INTENT(in)::fcld_obs
1315 REAL(KIND(1d0)),
INTENT(in)::dectime
1316 REAL(KIND(1d0)),
INTENT(in)::ZENITH_deg
1317 REAL(KIND(1d0)),
INTENT(in)::Tsurf_0
1318 REAL(KIND(1d0)),
INTENT(in)::avKdn
1319 REAL(KIND(1d0)),
INTENT(in)::Temp_C
1320 REAL(KIND(1d0)),
INTENT(in)::avRH
1321 REAL(KIND(1d0)),
INTENT(in)::ea_hPa
1322 REAL(KIND(1d0)),
INTENT(in)::qn1_obs
1323 REAL(KIND(1d0)),
INTENT(in)::SnowAlb_prev
1324 REAL(KIND(1d0)),
INTENT(in)::NARP_EMIS_SNOW
1325 REAL(KIND(1d0)),
INTENT(in)::NARP_TRANS_SITE
1326 REAL(KIND(1d0)),
INTENT(in)::tau_a, tau_f, SnowAlbMax, SnowAlbMin
1328 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in):: IceFrac
1329 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in):: sfr
1330 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in):: emis
1332 REAL(KIND(1d0)),
DIMENSION(nsurf) ::alb
1333 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in) ::alb_prev
1334 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out) ::alb_next
1335 REAL(KIND(1d0)),
INTENT(in) ::albDecTr_id
1337 REAL(KIND(1d0)),
INTENT(in) ::albEveTr_id
1338 REAL(KIND(1d0)),
INTENT(in) ::albGrass_id
1342 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowPack_prev
1343 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::snowFrac_prev
1344 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::snowFrac_next
1345 REAL(KIND(1d0)),
DIMENSION(nsurf)::SnowFrac
1347 REAL(KIND(1d0)),
INTENT(out)::ldown
1348 REAL(KIND(1d0)),
INTENT(out)::fcld
1349 REAL(KIND(1d0)),
INTENT(out)::qn1
1350 REAL(KIND(1d0)),
INTENT(out)::qn1_snowfree
1351 REAL(KIND(1d0)),
INTENT(out)::qn1_S
1352 REAL(KIND(1d0)),
INTENT(out)::kclear
1353 REAL(KIND(1d0)),
INTENT(out)::kup
1354 REAL(KIND(1d0)),
INTENT(out)::lup
1355 REAL(KIND(1d0)),
INTENT(out)::tsurf
1356 REAL(KIND(1d0)),
INTENT(out)::alb1
1357 REAL(KIND(1d0)),
INTENT(out)::SnowAlb_next
1358 REAL(KIND(1d0))::alb0
1359 REAL(KIND(1d0))::SnowAlb
1361 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out) ::qn1_ind_snow
1362 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out) ::kup_ind_snow
1363 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out) ::Tsurf_ind_snow
1364 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out):: tsurf_ind
1366 REAL(KIND(1d0)),
DIMENSION(nsurf):: lup_ind
1367 REAL(KIND(1d0)),
DIMENSION(nsurf):: kup_ind
1368 REAL(KIND(1d0)),
DIMENSION(nsurf):: qn1_ind
1370 REAL(KIND(1d0)),
PARAMETER::NAN = -999
1371 INTEGER :: NetRadiationMethod_use
1372 INTEGER::AlbedoChoice, ldown_option
1379 tstep, snowpack_prev, snowalb_prev, temp_c, &
1380 tau_a, tau_f, snowalbmax, snowalbmin)
1383 netradiationmethod, &
1385 netradiationmethod_use, albedochoice, ldown_option)
1387 snowfrac = snowfrac_prev
1388 IF (netradiationmethod_use > 0)
THEN 1391 IF (snowuse == 0) snowfrac = 0
1393 IF (ldown_option == 1)
THEN 1399 IF (ldown_option == 2)
THEN 1412 IF (diagnose == 1)
WRITE (*, *)
'Calling NARP...' 1413 IF (diagqn == 1)
WRITE (*, *)
'NetRadiationMethodX:', netradiationmethod_use
1414 IF (diagqn == 1)
WRITE (*, *)
'AlbedoChoice:', albedochoice
1417 nsurf, sfr, snowfrac, alb, emis, icefrac, &
1418 narp_trans_site, narp_emis_snow, &
1419 dectime, zenith_deg, tsurf_0, avkdn, temp_c, avrh, ea_hpa, qn1_obs, &
1421 albedochoice, ldown_option, netradiationmethod_use, diagqn, &
1422 qn1, qn1_snowfree, qn1_s, kclear, kup, ldown, lup, fcld, tsurf, &
1423 qn1_ind_snow, kup_ind_snow, tsurf_ind_snow, tsurf_ind, alb0, alb1)
1426 snowfrac = snowfrac_obs
1428 qn1_snowfree = qn1_obs
1440 snowfrac_next = snowfrac
1442 IF (ldown_option == 1)
THEN 1448 snowalb_next = snowalb
1455 StorageHeatMethod, qs_obs, OHMIncQF, Gridiv, &!input
1456 id, tstep, dt_since_start, Diagnose, sfr, &
1457 OHM_coef, OHM_threshSW, OHM_threshWD, &
1458 soilstore_id, SoilStoreCap, state_id, SnowUse, SnowFrac, DiagQS, &
1459 HDD_id, MetForcingData_grid, Ts5mindata_ir, qf, qn1, &
1460 avkdn, avu1, temp_c, zenith_deg, avrh, press_hpa, ldown, &
1461 bldgh, alb, emis, cpAnOHM, kkAnOHM, chAnOHM, EmissionsMethod, &
1462 Tair_av, qn1_av_prev, dqndt_prev, qn1_s_av_prev, dqnsdt_prev, &
1464 qn1_S, dataOutLineESTM, qs, &!output
1465 qn1_av_next, dqndt_next, qn1_s_av_next, dqnsdt_next, &
1466 deltaQi, a1, a2, a3)
1470 INTEGER,
INTENT(in) ::StorageHeatMethod
1471 INTEGER,
INTENT(in) ::OHMIncQF
1472 INTEGER,
INTENT(in) ::Gridiv
1473 INTEGER,
INTENT(in) ::id
1474 INTEGER,
INTENT(in) ::tstep
1475 INTEGER,
INTENT(in) ::dt_since_start
1476 INTEGER,
INTENT(in) ::Diagnose
1478 INTEGER,
INTENT(in) ::SnowUse
1479 INTEGER,
INTENT(in) ::DiagQS
1480 INTEGER,
INTENT(in) :: EmissionsMethod
1482 REAL(KIND(1d0)),
INTENT(in)::OHM_coef(
nsurf + 1, 4, 3)
1483 REAL(KIND(1d0)),
INTENT(in)::OHM_threshSW(
nsurf + 1)
1484 REAL(KIND(1d0)),
INTENT(in)::OHM_threshWD(
nsurf + 1)
1485 REAL(KIND(1d0)),
INTENT(in)::soilstore_id(
nsurf)
1486 REAL(KIND(1d0)),
INTENT(in)::SoilStoreCap(
nsurf)
1487 REAL(KIND(1d0)),
INTENT(in)::state_id(
nsurf)
1489 REAL(KIND(1d0)),
DIMENSION(12),
INTENT(in)::HDD_id
1490 REAL(KIND(1d0)),
INTENT(in)::qf
1491 REAL(KIND(1d0)),
INTENT(in)::qn1
1492 REAL(KIND(1d0)),
INTENT(in)::qs_obs
1493 REAL(KIND(1d0)),
INTENT(in)::avkdn, avu1, temp_c, zenith_deg, avrh, press_hpa, ldown
1494 REAL(KIND(1d0)),
INTENT(in)::bldgh
1496 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::sfr
1497 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::alb
1498 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::emis
1499 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::cpAnOHM
1500 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::kkAnOHM
1501 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::chAnOHM
1502 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowFrac
1504 REAL(KIND(1d0)),
DIMENSION(:, :),
INTENT(in)::MetForcingData_grid
1506 REAL(KIND(1d0)),
DIMENSION(:),
INTENT(in)::Ts5mindata_ir
1508 REAL(KIND(1d0)),
INTENT(in)::Tair_av
1509 REAL(KIND(1d0)),
INTENT(in)::qn1_av_prev
1510 REAL(KIND(1d0)),
INTENT(out)::qn1_av_next
1511 REAL(KIND(1d0)),
INTENT(in)::dqndt_prev
1512 REAL(KIND(1d0)),
INTENT(out)::dqndt_next
1513 REAL(KIND(1d0)),
INTENT(in)::qn1_s_av_prev
1514 REAL(KIND(1d0)),
INTENT(out)::qn1_s_av_next
1515 REAL(KIND(1d0)),
INTENT(in)::dqnsdt_prev
1516 REAL(KIND(1d0)),
INTENT(out)::dqnsdt_next
1522 REAL(KIND(1d0)),
DIMENSION(6, nsurf),
INTENT(in)::StoreDrainPrm
1524 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::deltaQi
1526 REAL(KIND(1d0)),
DIMENSION(27),
INTENT(out):: dataOutLineESTM
1527 REAL(KIND(1d0)),
INTENT(out)::qn1_S
1528 REAL(KIND(1d0)),
INTENT(out):: qs
1529 REAL(KIND(1d0)),
INTENT(out):: a1
1530 REAL(KIND(1d0)),
INTENT(out):: a2
1531 REAL(KIND(1d0)),
INTENT(out):: a3
1533 REAL(KIND(1d0))::Tair_mav_5d
1534 REAL(KIND(1d0))::qn1_use
1536 REAL(KIND(1d0)):: moist_surf(
nsurf)
1542 dataoutlineestm = -999
1549 IF (ohmincqf == 1)
THEN 1551 ELSEIF (ohmincqf == 0)
THEN 1555 IF (storageheatmethod == 0)
THEN 1558 ELSEIF (storageheatmethod == 1)
THEN 1559 tair_mav_5d = hdd_id(10)
1560 IF (diagnose == 1)
WRITE (*, *)
'Calling OHM...' 1561 CALL ohm(qn1_use, qn1_av_prev, dqndt_prev, qn1_av_next, dqndt_next, &
1562 qn1_s, qn1_s_av_prev, dqnsdt_prev, qn1_s_av_next, dqnsdt_next, &
1563 tstep, dt_since_start, &
1567 ohm_threshsw, ohm_threshwd, &
1568 soilstore_id, soilstorecap, state_id, &
1570 snowuse, snowfrac, &
1572 a1, a2, a3, qs, deltaqi)
1575 ELSEIF (storageheatmethod == 3)
THEN 1576 IF (diagnose == 1)
WRITE (*, *)
'Calling AnOHM...' 1582 moist_surf = state_id/storedrainprm(6, :)
1584 tstep, dt_since_start, &
1585 qn1_use, qn1_av_prev, dqndt_prev, qf, &
1586 metforcingdata_grid, moist_surf, &
1587 alb, emis, cpanohm, kkanohm, chanohm, &
1588 sfr,
nsurf, emissionsmethod, id, gridiv, &
1589 qn1_av_next, dqndt_next, &
1590 a1, a2, a3, qs, deltaqi)
1593 ELSEIF (storageheatmethod == 4 .OR. storageheatmethod == 14)
THEN 1595 IF (diagnose == 1)
WRITE (*, *)
'Calling ESTM...' 1599 avkdn, avu1, temp_c, zenith_deg, avrh, press_hpa, ldown, &
1600 bldgh, ts5mindata_ir, &
1602 dataoutlineestm, qs)
1613 snowUse, NonWaterFraction, addPipes, addImpervious, addVeg, addWaterBody, &
1614 state_id, soilstore_id, sfr, StoreDrainPrm, WaterDist, nsh_real, &
1615 drain_per_tstep, & !output
1616 drain, AddWaterRunoff, &
1617 AdditionalWater, runoffPipes, runoff_per_interval, &
1618 AddWater, stateOld, soilstoreOld)
1623 INTEGER,
INTENT(in) ::Diagnose
1624 INTEGER,
INTENT(in) ::snowUse
1626 REAL(KIND(1d0)),
INTENT(in)::NonWaterFraction
1627 REAL(KIND(1d0)),
INTENT(in)::addPipes
1628 REAL(KIND(1d0)),
INTENT(in)::addImpervious
1629 REAL(KIND(1d0)),
INTENT(in)::addVeg
1630 REAL(KIND(1d0)),
INTENT(in)::addWaterBody
1631 REAL(KIND(1d0)),
INTENT(in)::nsh_real
1633 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in) ::state_id
1634 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in) ::soilstore_id
1635 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in) ::sfr
1636 REAL(KIND(1d0)),
DIMENSION(6, nsurf),
INTENT(in) ::StoreDrainPrm
1637 REAL(KIND(1d0)),
DIMENSION(nsurf + 1, nsurf - 1),
INTENT(in)::WaterDist
1639 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out):: drain
1640 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out):: AddWaterRunoff
1641 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out):: AddWater
1642 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out):: stateOld
1643 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out):: soilstoreOld
1645 REAL(KIND(1d0)),
INTENT(out)::drain_per_tstep
1646 REAL(KIND(1d0)),
INTENT(out)::AdditionalWater
1647 REAL(KIND(1d0)),
INTENT(out)::runoffPipes
1648 REAL(KIND(1d0)),
INTENT(out)::runoff_per_interval
1653 soilstoreold = soilstore_id
1661 additionalwater = addpipes + addimpervious + addveg + addwaterbody
1664 runoffpipes = addpipes
1666 runoff_per_interval = addpipes
1670 IF (diagnose == 1)
WRITE (*, *)
'Calling Drainage...' 1672 IF (nonwaterfraction /= 0)
THEN 1673 DO is = 1,
nsurf - 1
1678 storedrainprm(6, is), &
1679 storedrainprm(2, is), &
1680 storedrainprm(3, is), &
1681 storedrainprm(4, is), &
1688 drain_per_tstep = dot_product(drain(1:
nsurf - 1), sfr(1:
nsurf - 1))/nonwaterfraction
1690 drain(1:
nsurf - 1) = 0
1697 IF (diagnose == 1)
WRITE (*, *)
'Calling ReDistributeWater...' 1701 snowuse, waterdist, sfr, drain, &
1702 addwaterrunoff, addwater)
1709 avdens, avcp, h_mod, qn1, dectime, &!input
1714 REAL(KIND(1d0)),
INTENT(in)::avdens
1715 REAL(KIND(1d0)),
INTENT(in)::avcp
1716 REAL(KIND(1d0)),
INTENT(in)::h_mod
1717 REAL(KIND(1d0)),
INTENT(in)::qn1
1718 REAL(KIND(1d0)),
INTENT(in)::dectime
1719 REAL(KIND(1d0)),
INTENT(out)::H_init
1721 REAL(KIND(1d0)),
PARAMETER::NAN = -999
1722 INTEGER,
PARAMETER::notUsedI = -999
1728 IF (h_mod /= nan)
THEN 1729 h_init = h_mod/(avdens*avcp)
1731 h_init = (qn1*0.2)/(avdens*avcp)
1732 CALL errorhint(38,
'LUMPS unable to calculate realistic value for H_mod.', h_mod, dectime, notusedi)
1742 Diagnose, snowuse, &!input
1743 tstep, imin, it, EvapMethod, snowCalcSwitch, dayofWeek_id, CRWmin, CRWmax, &
1744 dectime, lvS_J_kg, avRh, Press_hPa, Temp_C, &
1745 RAsnow, psyc_hPa, sIce_hPa, &
1746 PervFraction, vegfraction, addimpervious, qn1_snowfree, qf, qs, vpd_hPa, s_hPa, &
1747 ResistSurf, RA, rb, snowdensmin, precip, PipeCapacity, RunoffToWater, &
1748 NonWaterFraction, wu_EveTr, wu_DecTr, wu_Grass, addVeg, addWaterBody, SnowLimPaved, SnowLimBldg, &
1749 SurfaceArea, FlowChange, drain, WetThresh, stateOld, mw_ind, SoilStoreCap, rainonsnow, &
1750 freezmelt, freezstate, freezstatevol, Qm_Melt, Qm_rain, Tsurf_ind, sfr, &
1751 StateLimit, AddWater, addwaterrunoff, StoreDrainPrm, SnowPackLimit, SnowProf_24hr, &
1752 SnowPack_in, SnowFrac_in, SnowWater_in, iceFrac_in, SnowDens_in, &! input:
1753 runoff_per_interval_in, state_id_in, soilstore_id_in, &! input:
1754 state_id_out, soilstore_id_out, &! output:
1755 SnowPack_out, SnowFrac_out, SnowWater_out, iceFrac_out, SnowDens_out, &! output
1756 runoffSoil, &! output:
1758 state_per_tstep, NWstate_per_tstep, qe, &
1759 swe, chSnow_per_interval, ev_per_tstep, runoff_per_tstep, &
1760 surf_chang_per_tstep, runoffPipes, mwstore, runoffwaterbody, &
1761 runoffAGveg, runoffAGimpervious)
1765 INTEGER,
INTENT(in) ::Diagnose
1766 INTEGER,
INTENT(in) ::snowuse
1767 INTEGER,
INTENT(in) ::tstep
1768 INTEGER,
INTENT(in) ::imin
1769 INTEGER,
INTENT(in) ::it
1770 INTEGER,
INTENT(in) ::EvapMethod
1772 INTEGER,
DIMENSION(nsurf),
INTENT(in)::snowCalcSwitch
1773 INTEGER,
DIMENSION(3),
INTENT(in)::dayofWeek_id
1775 REAL(KIND(1d0)),
INTENT(in)::CRWmin
1776 REAL(KIND(1d0)),
INTENT(in)::CRWmax
1777 REAL(KIND(1d0)),
INTENT(in)::dectime
1778 REAL(KIND(1d0)),
INTENT(in)::lvS_J_kg
1781 REAL(KIND(1d0)),
INTENT(in)::avRh
1782 REAL(KIND(1d0)),
INTENT(in)::Press_hPa
1783 REAL(KIND(1d0)),
INTENT(in)::Temp_C
1784 REAL(KIND(1d0)),
INTENT(in)::RAsnow
1785 REAL(KIND(1d0)),
INTENT(in)::psyc_hPa
1787 REAL(KIND(1d0)),
INTENT(in)::sIce_hPa
1788 REAL(KIND(1d0)),
INTENT(in)::PervFraction
1789 REAL(KIND(1d0)),
INTENT(in)::vegfraction
1790 REAL(KIND(1d0)),
INTENT(in)::addimpervious
1791 REAL(KIND(1d0)),
INTENT(in)::qn1_snowfree
1792 REAL(KIND(1d0)),
INTENT(in)::qf
1793 REAL(KIND(1d0)),
INTENT(in)::qs
1794 REAL(KIND(1d0)),
INTENT(in)::vpd_hPa
1795 REAL(KIND(1d0)),
INTENT(in)::s_hPa
1796 REAL(KIND(1d0)),
INTENT(in)::ResistSurf
1797 REAL(KIND(1d0)),
INTENT(in)::RA
1798 REAL(KIND(1d0)),
INTENT(in)::rb
1799 REAL(KIND(1d0)),
INTENT(in)::snowdensmin
1800 REAL(KIND(1d0)),
INTENT(in)::precip
1801 REAL(KIND(1d0)),
INTENT(in)::PipeCapacity
1802 REAL(KIND(1d0)),
INTENT(in)::RunoffToWater
1803 REAL(KIND(1d0)),
INTENT(in)::NonWaterFraction
1804 REAL(KIND(1d0)),
INTENT(in)::wu_EveTr
1805 REAL(KIND(1d0)),
INTENT(in)::wu_DecTr
1806 REAL(KIND(1d0)),
INTENT(in)::wu_Grass
1807 REAL(KIND(1d0)),
INTENT(in)::addVeg
1808 REAL(KIND(1d0)),
INTENT(in)::addWaterBody
1809 REAL(KIND(1d0)),
INTENT(in)::SnowLimPaved
1810 REAL(KIND(1d0)),
INTENT(in)::SnowLimBldg
1811 REAL(KIND(1d0)),
INTENT(in)::SurfaceArea
1812 REAL(KIND(1d0)),
INTENT(in)::FlowChange
1814 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::drain
1815 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::WetThresh
1816 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::stateOld
1817 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::mw_ind
1818 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SoilStoreCap
1819 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::rainonsnow
1820 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::freezmelt
1821 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::freezstate
1822 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::freezstatevol
1823 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::Qm_Melt
1824 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::Qm_rain
1825 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::Tsurf_ind
1826 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::sfr
1827 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowPackLimit
1828 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::StateLimit
1829 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::AddWater
1830 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::addwaterrunoff
1831 REAL(KIND(1d0)),
DIMENSION(6, nsurf),
INTENT(in)::StoreDrainPrm
1832 REAL(KIND(1d0)),
DIMENSION(0:23, 2),
INTENT(in):: SnowProf_24hr
1835 REAL(KIND(1d0)),
INTENT(in)::runoff_per_interval_in
1836 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::state_id_in
1837 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::soilstore_id_in
1838 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowPack_in
1839 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowFrac_in
1840 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowWater_in
1841 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::iceFrac_in
1842 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowDens_in
1845 REAL(KIND(1d0))::runoff_per_interval_out
1846 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::state_id_out
1847 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::soilstore_id_out
1848 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::SnowPack_out
1849 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::SnowFrac_out
1850 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::SnowWater_out
1851 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::iceFrac_out
1852 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::SnowDens_out
1854 REAL(KIND(1d0)),
DIMENSION(nsurf)::runoffSnow
1855 REAL(KIND(1d0)),
DIMENSION(nsurf)::runoff
1856 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(out)::runoffSoil
1857 REAL(KIND(1d0)),
DIMENSION(nsurf)::chang
1858 REAL(KIND(1d0)),
DIMENSION(nsurf)::changSnow
1859 REAL(KIND(1d0)),
DIMENSION(nsurf)::snowDepth
1860 REAL(KIND(1d0)),
DIMENSION(nsurf)::SnowToSurf
1861 REAL(KIND(1d0)),
DIMENSION(nsurf)::ev_snow
1862 REAL(KIND(1d0)),
DIMENSION(2),
INTENT(out)::SnowRemoval
1863 REAL(KIND(1d0)),
DIMENSION(nsurf)::evap
1864 REAL(KIND(1d0)),
DIMENSION(nsurf)::rss_nsurf
1866 REAL(KIND(1d0))::p_mm
1868 REAL(KIND(1d0))::qe_surf
1869 REAL(KIND(1d0)),
INTENT(out)::state_per_tstep
1870 REAL(KIND(1d0)),
INTENT(out)::NWstate_per_tstep
1871 REAL(KIND(1d0)),
INTENT(out)::qe
1872 REAL(KIND(1d0)),
INTENT(out)::swe
1874 REAL(KIND(1d0)),
INTENT(out)::chSnow_per_interval
1875 REAL(KIND(1d0)),
INTENT(out)::ev_per_tstep
1876 REAL(KIND(1d0))::qe_per_tstep
1877 REAL(KIND(1d0)),
INTENT(out)::runoff_per_tstep
1878 REAL(KIND(1d0)),
INTENT(out)::surf_chang_per_tstep
1879 REAL(KIND(1d0)),
INTENT(out)::runoffPipes
1880 REAL(KIND(1d0)),
INTENT(out)::mwstore
1881 REAL(KIND(1d0)),
INTENT(out)::runoffwaterbody
1882 REAL(KIND(1d0))::runoffWaterBody_m3
1883 REAL(KIND(1d0))::runoffPipes_m3
1884 REAL(KIND(1d0)),
INTENT(out)::runoffAGveg
1885 REAL(KIND(1d0)),
INTENT(out)::runoffAGimpervious
1890 REAL(KIND(1d0))::runoff_per_interval
1891 REAL(KIND(1d0)),
DIMENSION(nsurf)::state_id
1892 REAL(KIND(1d0)),
DIMENSION(nsurf)::soilstore_id
1893 REAL(KIND(1d0)),
DIMENSION(nsurf)::SnowPack
1894 REAL(KIND(1d0)),
DIMENSION(nsurf)::SnowFrac
1895 REAL(KIND(1d0)),
DIMENSION(nsurf)::SnowWater
1896 REAL(KIND(1d0)),
DIMENSION(nsurf)::iceFrac
1897 REAL(KIND(1d0)),
DIMENSION(nsurf)::SnowDens
1899 REAL(KIND(1d0)),
DIMENSION(2) ::SurplusEvap
1900 REAL(KIND(1d0))::surplusWaterBody
1901 REAL(KIND(1d0))::pin
1906 REAL(KIND(1d0))::qn_e
1907 REAL(KIND(1d0))::tlv
1908 REAL(KIND(1d0))::runoffAGimpervious_m3
1909 REAL(KIND(1d0))::runoffAGveg_m3
1910 REAL(KIND(1d0))::nsh_real
1911 REAL(KIND(1d0))::tstep_real
1912 REAL(KIND(1d0))::ev_tot
1913 REAL(KIND(1d0))::qe_tot
1914 REAL(KIND(1d0))::surf_chang_tot
1915 REAL(KIND(1d0))::runoff_tot
1916 REAL(KIND(1d0))::chSnow_tot
1918 REAL(KIND(1d0)),
DIMENSION(7)::capStore
1920 runoff_per_interval = runoff_per_interval_in
1921 state_id = state_id_in
1922 soilstore_id = soilstore_id_in
1923 snowpack = snowpack_in
1924 snowfrac = snowfrac_in
1925 snowwater = snowwater_in
1926 icefrac = icefrac_in
1927 snowdens = snowdens_in
1929 tstep_real = tstep*1.d0
1930 nsh_real = 3600/tstep_real
1936 pin = max(0., precip)
1947 surf_chang_per_tstep = 0
1948 runoff_per_tstep = 0
1950 nwstate_per_tstep = 0
1953 chsnow_per_interval = 0
1956 runoffagimpervious = 0
1957 surpluswaterbody = 0
1965 qn_e = qn1_snowfree + qf - qs
1967 IF (diagnose == 1)
WRITE (*, *)
'Calling evap_SUEWS and SoilStore...' 1969 IF (snowuse == 1 .and. snowcalcswitch(is) == 1)
THEN 1970 IF (sfr(is) /= 0)
THEN 1973 tstep, imin, it, dectime, is, &
1974 evapmethod, crwmin, crwmax, nsh_real, lvs_j_kg,
lv_j_kg,
avdens, &
1975 avrh, press_hpa, temp_c, rasnow, psyc_hpa,
avcp, sice_hpa, &
1976 pervfraction, vegfraction, addimpervious, &
1977 vpd_hpa, qn_e, s_hpa, resistsurf, ra, rb, tlv, snowdensmin, snowprof_24hr, precip, &
1978 pipecapacity, runofftowater, &
1979 addveg, snowlimpaved, snowlimbldg, flowchange, drain, &
1980 wetthresh, stateold, mw_ind, soilstorecap, rainonsnow, &
1981 freezmelt, freezstate, freezstatevol, &
1982 qm_melt, qm_rain, tsurf_ind, sfr, dayofweek_id, storedrainprm, snowpacklimit, &
1983 addwater, addwaterrunoff, &
1984 soilstore_id, snowpack, surplusevap, &
1985 snowfrac, snowwater, icefrac, snowdens, &
1986 runoffagimpervious, runoffagveg, surpluswaterbody, &
1987 rss_nsurf, runoffsnow, &
1988 runoff, runoffsoil, chang, changsnow, snowtosurf, state_id, ev_snow, &
1989 snowdepth, snowremoval, swe, ev, chsnow_tot, &
1990 ev_tot, qe_tot, runoff_tot, surf_chang_tot, &
1991 runoffpipes, mwstore, runoffwaterbody)
1994 ev_per_tstep = ev_per_tstep + ev_tot
1995 qe_per_tstep = qe_per_tstep + qe_tot
1996 runoff_per_tstep = runoff_per_tstep + runoff_tot
1997 surf_chang_per_tstep = surf_chang_per_tstep + surf_chang_tot
1998 chsnow_per_interval = chsnow_per_interval + chsnow_tot
2009 capstore(is) = storedrainprm(6, is)
2012 evapmethod, state_id(is), wetthresh(is), capstore(is), &
2013 vpd_hpa,
avdens,
avcp, qn_e, s_hpa, psyc_hpa, resistsurf, ra, rb, tlv, &
2014 rss_nsurf(is), ev, qe_surf)
2018 is, sfr, pipecapacity, runofftowater, pin, &
2019 wu_evetr, wu_dectr, wu_grass, drain, addwater, addimpervious, nsh_real, stateold, addwaterrunoff, &
2020 pervfraction, addveg, soilstorecap, addwaterbody, flowchange, statelimit, runoffagimpervious, surpluswaterbody, &
2021 runoffagveg, runoffpipes, ev, soilstore_id, surplusevap, runoffwaterbody, &
2022 p_mm, chang, runoff, state_id)
2028 ev_per_tstep = ev_per_tstep + evap(is)*sfr(is)
2031 qe_per_tstep = qe_per_tstep + qe_surf*sfr(is)
2034 surf_chang_per_tstep = surf_chang_per_tstep + (state_id(is) - stateold(is))*sfr(is)
2037 runoff_per_tstep = runoff_per_tstep + runoff(is)*sfr(is)
2040 state_per_tstep = state_per_tstep + state_id(is)*sfr(is)
2043 runoff_per_interval = runoff_per_interval + (runoff(is)*sfr(is))
2045 IF (nonwaterfraction /= 0 .AND. is /=
watersurf)
THEN 2046 nwstate_per_tstep = nwstate_per_tstep + (state_id(is)*sfr(is)/nonwaterfraction)
2060 runoffagimpervious_m3 = runoffagimpervious/1000*surfacearea
2061 runoffagveg_m3 = runoffagveg/1000*surfacearea
2062 runoffwaterbody_m3 = runoffwaterbody/1000*surfacearea
2063 runoffpipes_m3 = runoffpipes/1000*surfacearea
2065 runoff_per_interval_out = runoff_per_interval
2066 state_id_out = state_id
2067 soilstore_id_out = soilstore_id
2068 snowpack_out = snowpack
2069 snowfrac_out = snowfrac
2070 snowwater_out = snowwater
2071 icefrac_out = icefrac
2072 snowdens_out = snowdens
2080 qn1, qf, QmRain, qeOut, qs, QmFreez, qm, avdens, avcp, tsurf, Temp_C, RA, &
2081 qh, qh_residual, qh_resist)
2084 INTEGER,
INTENT(in) :: QHMethod
2086 REAL(KIND(1d0)),
INTENT(in)::qn1
2087 REAL(KIND(1d0)),
INTENT(in)::qf
2088 REAL(KIND(1d0)),
INTENT(in)::QmRain
2089 REAL(KIND(1d0)),
INTENT(in)::qeOut
2090 REAL(KIND(1d0)),
INTENT(in)::qs
2091 REAL(KIND(1d0)),
INTENT(in)::QmFreez
2092 REAL(KIND(1d0)),
INTENT(in)::qm
2093 REAL(KIND(1d0)),
INTENT(in)::avdens
2094 REAL(KIND(1d0)),
INTENT(in)::avcp
2095 REAL(KIND(1d0)),
INTENT(in)::tsurf
2096 REAL(KIND(1d0)),
INTENT(in)::Temp_C
2097 REAL(KIND(1d0)),
INTENT(in)::RA
2099 REAL(KIND(1d0)),
INTENT(out)::qh
2100 REAL(KIND(1d0)),
INTENT(out)::qh_resist
2101 REAL(KIND(1d0)),
INTENT(out)::qh_residual
2103 REAL(KIND(1d0)),
PARAMETER::NAN = -999
2106 qh_residual = (qn1 + qf + qmrain) - (qeout + qs + qm + qmfreez)
2111 qh_resist = avdens*avcp*(tsurf - temp_c)/ra
2117 SELECT CASE (qhmethod)
2129 StabilityMethod, &!input:
2130 Diagnose, AerodynamicResistanceMethod, RoughLenHeatMethod, snowUse, &
2131 id, it, gsModel, SMDMethod, &
2132 avdens, avcp, QH_init, zzd, z0m, zdm, &
2133 avU1, Temp_C, VegFraction, &
2134 avkdn, Kmax, G1, G2, G3, G4, G5, G6, S1, S2, TH, TL, dq, &
2135 xsmd, vsmd, MaxConductance, LAIMax, LAI_id, SnowFrac, sfr, &
2136 UStar, TStar, L_mod, &!output
2137 zL, gsc, ResistSurf, RA, RAsnow, rb)
2141 INTEGER,
INTENT(in)::StabilityMethod
2142 INTEGER,
INTENT(in)::Diagnose
2143 INTEGER,
INTENT(in)::AerodynamicResistanceMethod
2144 INTEGER,
INTENT(in)::RoughLenHeatMethod
2145 INTEGER,
INTENT(in)::snowUse
2146 INTEGER,
INTENT(in)::id
2147 INTEGER,
INTENT(in)::it
2148 INTEGER,
INTENT(in)::gsModel
2149 INTEGER,
INTENT(in)::SMDMethod
2152 REAL(KIND(1d0)),
INTENT(in)::avdens
2153 REAL(KIND(1d0)),
INTENT(in)::avcp
2154 REAL(KIND(1d0)),
INTENT(in)::QH_init
2155 REAL(KIND(1d0)),
INTENT(in)::zzd
2156 REAL(KIND(1d0)),
INTENT(in)::z0m
2157 REAL(KIND(1d0)),
INTENT(in)::zdm
2158 REAL(KIND(1d0)),
INTENT(in)::avU1
2159 REAL(KIND(1d0)),
INTENT(in)::Temp_C
2160 REAL(KIND(1d0)),
INTENT(in)::VegFraction
2161 REAL(KIND(1d0)),
INTENT(in)::avkdn
2162 REAL(KIND(1d0)),
INTENT(in)::Kmax
2163 REAL(KIND(1d0)),
INTENT(in)::G1
2164 REAL(KIND(1d0)),
INTENT(in)::G2
2165 REAL(KIND(1d0)),
INTENT(in)::G3
2166 REAL(KIND(1d0)),
INTENT(in)::G4
2167 REAL(KIND(1d0)),
INTENT(in)::G5
2168 REAL(KIND(1d0)),
INTENT(in)::G6
2169 REAL(KIND(1d0)),
INTENT(in)::S1
2170 REAL(KIND(1d0)),
INTENT(in)::S2
2171 REAL(KIND(1d0)),
INTENT(in)::TH
2172 REAL(KIND(1d0)),
INTENT(in)::TL
2173 REAL(KIND(1d0)),
INTENT(in)::dq
2174 REAL(KIND(1d0)),
INTENT(in)::xsmd
2175 REAL(KIND(1d0)),
INTENT(in)::vsmd
2177 REAL(KIND(1d0)),
DIMENSION(3),
INTENT(in) ::MaxConductance
2178 REAL(KIND(1d0)),
DIMENSION(3),
INTENT(in) ::LAIMax
2179 REAL(KIND(1d0)),
DIMENSION(3),
INTENT(in) ::LAI_id
2181 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::SnowFrac
2182 REAL(KIND(1d0)),
DIMENSION(nsurf),
INTENT(in)::sfr
2184 REAL(KIND(1d0)),
INTENT(out)::TStar
2185 REAL(KIND(1d0)),
INTENT(out) ::UStar
2186 REAL(KIND(1d0)),
INTENT(out) ::zL
2187 REAL(KIND(1d0)),
INTENT(out) ::gsc
2188 REAL(KIND(1d0)),
INTENT(out) ::ResistSurf
2189 REAL(KIND(1d0)),
INTENT(out) ::RA
2190 REAL(KIND(1d0)),
INTENT(out) ::RAsnow
2191 REAL(KIND(1d0)),
INTENT(out) ::rb
2192 REAL(KIND(1d0)),
INTENT(out) ::L_mod
2193 REAL(KIND(1d0)) ::gfunc
2201 IF (diagnose == 1)
WRITE (*, *)
'Calling STAB_lumps...' 2218 IF (diagnose == 1)
WRITE (*, *)
'Calling AerodynamicResistance...' 2226 aerodynamicresistancemethod, &
2228 roughlenheatmethod, &
2231 IF (snowuse == 1)
THEN 2232 IF (diagnose == 1)
WRITE (*, *)
'Calling AerodynamicResistance for snow...' 2240 aerodynamicresistancemethod, &
2246 IF (diagnose == 1)
WRITE (*, *)
'Calling SurfaceResistance...' 2250 smdmethod, snowfrac, sfr, avkdn, temp_c, dq, xsmd, vsmd, maxconductance, &
2251 laimax, lai_id, gsmodel, kmax, &
2252 g1, g2, g3, g4, g5, g6, th, tl, s1, s2, &
2253 gfunc, gsc, resistsurf)
2255 IF (diagnose == 1)
WRITE (*, *)
'Calling BoundaryLayerResistance...' 2268 AdditionalWater, alb, avkdn, avU10_ms, azimuth, &!input
2269 chSnow_per_interval, dectime, &
2270 drain_per_tstep, E_mod, ev_per_tstep, ext_wu, Fc, Fc_build, fcld, &
2271 Fc_metab, Fc_photo, Fc_respi, Fc_point, Fc_traff, FlowChange, &
2272 h_mod, id, imin, int_wu, it, iy, &
2273 kup, LAI_id, ldown, l_mod, lup, mwh, &
2275 nsh_real, NWstate_per_tstep, Precip, q2_gkg, &
2276 qeOut, qf, qh, qh_resist, Qm, QmFreez, &
2277 QmRain, qn1, qn1_S, qn1_snowfree, qs, RA, &
2278 resistsurf, RH2, runoffAGimpervious, runoffAGveg, &
2279 runoff_per_tstep, runoffPipes, runoffSoil_per_tstep, &
2280 runoffWaterBody, sfr, smd, smd_nsurf, SnowAlb, SnowRemoval, &
2281 state_id, state_per_tstep, surf_chang_per_tstep, swe, t2_C, tskin_C, &
2282 tot_chang_per_tstep, tsurf, UStar, wu_DecTr, &
2283 wu_EveTr, wu_Grass, z0m, zdm, zenith_deg, &
2284 datetimeLine, dataOutLineSUEWS)
2287 REAL(KIND(1d0)),
PARAMETER :: NAN = -999
2288 INTEGER,
INTENT(in) :: iy
2290 INTEGER,
INTENT(in) :: id
2292 INTEGER,
INTENT(in) :: it
2293 INTEGER,
INTENT(in) :: imin
2295 REAL(KIND(1d0)),
INTENT(in) :: AdditionalWater
2296 REAL(KIND(1d0)),
INTENT(in) :: alb(
nsurf)
2297 REAL(KIND(1d0)),
INTENT(in) :: avkdn
2298 REAL(KIND(1d0)),
INTENT(in) :: avU10_ms
2299 REAL(KIND(1d0)),
INTENT(in) :: azimuth
2300 REAL(KIND(1d0)),
INTENT(in) :: chSnow_per_interval
2301 REAL(KIND(1d0)),
INTENT(in) :: dectime
2302 REAL(KIND(1d0)),
INTENT(in) :: drain_per_tstep
2303 REAL(KIND(1d0)),
INTENT(in) :: E_mod
2304 REAL(KIND(1d0)),
INTENT(in) :: ev_per_tstep
2305 REAL(KIND(1d0)),
INTENT(in) :: ext_wu
2306 REAL(KIND(1d0)),
INTENT(in) :: Fc
2307 REAL(KIND(1d0)),
INTENT(in) :: Fc_build
2308 REAL(KIND(1d0)),
INTENT(in) :: Fc_metab
2309 REAL(KIND(1d0)),
INTENT(in) :: Fc_photo
2310 REAL(KIND(1d0)),
INTENT(in) :: Fc_respi
2311 REAL(KIND(1d0)),
INTENT(in) :: Fc_point
2312 REAL(KIND(1d0)),
INTENT(in) :: Fc_traff
2313 REAL(KIND(1d0)),
INTENT(in) :: fcld
2314 REAL(KIND(1d0)),
INTENT(in) :: FlowChange
2315 REAL(KIND(1d0)),
INTENT(in) :: h_mod
2316 REAL(KIND(1d0)),
INTENT(in) :: int_wu
2317 REAL(KIND(1d0)),
INTENT(in) :: kup
2318 REAL(KIND(1d0)),
INTENT(in) :: l_mod
2319 REAL(KIND(1d0)),
INTENT(in) :: LAI_id(
nvegsurf)
2320 REAL(KIND(1d0)),
INTENT(in) :: ldown
2321 REAL(KIND(1d0)),
INTENT(in) :: lup
2322 REAL(KIND(1d0)),
INTENT(in) :: mwh
2323 REAL(KIND(1d0)),
INTENT(in) :: MwStore
2324 REAL(KIND(1d0)),
INTENT(in) :: nsh_real
2325 REAL(KIND(1d0)),
INTENT(in) :: NWstate_per_tstep
2326 REAL(KIND(1d0)),
INTENT(in) :: Precip
2327 REAL(KIND(1d0)),
INTENT(in) :: q2_gkg
2328 REAL(KIND(1d0)),
INTENT(in) :: qeOut
2329 REAL(KIND(1d0)),
INTENT(in) :: qf
2330 REAL(KIND(1d0)),
INTENT(in) :: qh
2331 REAL(KIND(1d0)),
INTENT(in) :: qh_resist
2332 REAL(KIND(1d0)),
INTENT(in) :: Qm
2333 REAL(KIND(1d0)),
INTENT(in) :: QmFreez
2334 REAL(KIND(1d0)),
INTENT(in) :: QmRain
2335 REAL(KIND(1d0)),
INTENT(in) :: qn1
2336 REAL(KIND(1d0)),
INTENT(in) :: qn1_S
2337 REAL(KIND(1d0)),
INTENT(in) :: qn1_snowfree
2338 REAL(KIND(1d0)),
INTENT(in) :: qs
2339 REAL(KIND(1d0)),
INTENT(in) :: RA
2340 REAL(KIND(1d0)),
INTENT(in) :: resistsurf
2341 REAL(KIND(1d0)),
INTENT(in) :: RH2
2342 REAL(KIND(1d0)),
INTENT(in) :: runoff_per_tstep
2343 REAL(KIND(1d0)),
INTENT(in) :: runoffAGimpervious
2344 REAL(KIND(1d0)),
INTENT(in) :: runoffAGveg
2345 REAL(KIND(1d0)),
INTENT(in) :: runoffPipes
2346 REAL(KIND(1d0)),
INTENT(in) :: runoffSoil_per_tstep
2347 REAL(KIND(1d0)),
INTENT(in) :: runoffWaterBody
2348 REAL(KIND(1d0)),
INTENT(in) :: sfr(
nsurf)
2349 REAL(KIND(1d0)),
INTENT(in) :: smd
2350 REAL(KIND(1d0)),
INTENT(in) :: smd_nsurf(
nsurf)
2351 REAL(KIND(1d0)),
INTENT(in) :: SnowAlb
2352 REAL(KIND(1d0)),
INTENT(in) :: SnowRemoval(2)
2353 REAL(KIND(1d0)),
INTENT(in) :: state_id(
nsurf)
2354 REAL(KIND(1d0)),
INTENT(in) :: state_per_tstep
2355 REAL(KIND(1d0)),
INTENT(in) :: surf_chang_per_tstep
2356 REAL(KIND(1d0)),
INTENT(in) :: swe
2357 REAL(KIND(1d0)),
INTENT(in) :: t2_C
2358 REAL(KIND(1d0)),
INTENT(in) :: tskin_C
2359 REAL(KIND(1d0)),
INTENT(in) :: tot_chang_per_tstep
2360 REAL(KIND(1d0)),
INTENT(in) :: tsurf
2361 REAL(KIND(1d0)),
INTENT(in) :: UStar
2362 REAL(KIND(1d0)),
INTENT(in) :: wu_DecTr
2363 REAL(KIND(1d0)),
INTENT(in) :: wu_EveTr
2364 REAL(KIND(1d0)),
INTENT(in) :: wu_Grass
2365 REAL(KIND(1d0)),
INTENT(in) :: z0m
2366 REAL(KIND(1d0)),
INTENT(in) :: zdm
2367 REAL(KIND(1d0)),
INTENT(in) :: zenith_deg
2369 REAL(KIND(1D0)),
DIMENSION(5),
INTENT(OUT)::datetimeLine
2370 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutSUEWS - 5),
INTENT(out) :: dataOutLineSUEWS
2374 REAL(KIND(1d0)):: LAI_wt
2375 REAL(KIND(1d0)):: RH2_pct
2378 REAL(KIND(1d0))::ResistSurf_x
2379 REAL(KIND(1d0))::l_mod_x
2380 REAL(KIND(1d0))::bulkalbedo
2381 REAL(KIND(1d0))::smd_nsurf_x(
nsurf)
2382 REAL(KIND(1d0))::state_x(
nsurf)
2389 state_x = unpack(spread(nan, dim=1, ncopies=
SIZE(sfr)), mask=(sfr < 0.00001), field=state_id)
2390 smd_nsurf_x = unpack(spread(nan, dim=1, ncopies=
SIZE(sfr)), mask=(sfr < 0.00001), field=smd_nsurf)
2392 resistsurf_x = min(9999., resistsurf)
2394 l_mod_x = max(min(9999., l_mod), -9999.)
2403 lai_wt = dot_product(lai_id(:), sfr(1 + 2:
nvegsurf + 2))
2406 bulkalbedo = dot_product(alb, sfr)
2432 REAL(iy, KIND(1D0)),
REAL(id, KIND(1D0)), &
2433 REAL(it, KIND(1D0)),
REAL(imin, KIND(1D0)), dectime]
2435 dataoutlinesuews = [ &
2436 avkdn, kup, ldown, lup, tsurf, &
2437 qn1, qf, qs, qh, qeout, &
2438 h_mod, e_mod, qh_resist, &
2439 precip, ext_wu, ev_per_tstep, runoff_per_tstep, tot_chang_per_tstep, &
2440 surf_chang_per_tstep, state_per_tstep, nwstate_per_tstep, drain_per_tstep, smd, &
2441 flowchange/nsh_real, additionalwater, &
2442 runoffsoil_per_tstep, runoffpipes, runoffagimpervious, runoffagveg, runoffwaterbody, &
2443 int_wu, wu_evetr, wu_dectr, wu_grass, &
2444 smd_nsurf_x(1:
nsurf - 1), &
2446 zenith_deg, azimuth, bulkalbedo, fcld, &
2448 ustar, l_mod, ra, resistsurf, &
2450 fc_photo, fc_respi, fc_metab, fc_traff, fc_build, fc_point, &
2451 qn1_snowfree, qn1_s, snowalb, &
2452 qm, qmfreez, qmrain, swe, mwh, mwstore, chsnow_per_interval, &
2454 tskin_c, t2_c, q2_gkg, avu10_ms, rh2_pct &
2466 SnowUse, storageheatmethod, &!input
2467 ReadLinesMetdata, NumberOfGrids, &
2469 datetimeLine, dataOutLineSUEWS, dataOutLineSnow, dataOutLineESTM, dataoutLineURSL, &!input
2470 dataOutSUEWS, dataOutSnow, dataOutESTM, dataOutRSL)
2473 INTEGER,
INTENT(in) ::ReadLinesMetdata
2474 INTEGER,
INTENT(in) ::NumberOfGrids
2475 INTEGER,
INTENT(in) ::Gridiv
2476 INTEGER,
INTENT(in) ::SnowUse
2477 INTEGER,
INTENT(in) ::storageheatmethod
2478 INTEGER,
INTENT(in) ::ir
2480 REAL(KIND(1d0)),
DIMENSION(5),
INTENT(in) :: datetimeLine
2481 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutSUEWS - 5),
INTENT(in) :: dataOutLineSUEWS
2482 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutESTM - 5),
INTENT(in) :: dataOutLineESTM
2483 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutSnow - 5),
INTENT(in) :: dataOutLineSnow
2484 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutRSL - 5),
INTENT(in) :: dataoutLineURSL
2486 REAL(KIND(1d0)),
INTENT(inout) :: dataOutSUEWS(readlinesmetdata,
ncolumnsdataoutsuews, numberofgrids)
2487 REAL(KIND(1d0)),
INTENT(inout) :: dataOutSnow(readlinesmetdata,
ncolumnsdataoutsnow, numberofgrids)
2488 REAL(KIND(1d0)),
INTENT(inout) :: dataOutESTM(readlinesmetdata,
ncolumnsdataoutestm, numberofgrids)
2489 REAL(KIND(1d0)),
INTENT(inout) :: dataOutRSL(readlinesmetdata,
ncolumnsdataoutrsl, numberofgrids)
2498 IF (snowuse == 1)
THEN 2502 IF (storageheatmethod == 4)
THEN 2583 vegfraction, ImpervFraction, PervFraction, NonWaterFraction)
2586 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN)::sfr
2587 REAL(KIND(1D0)),
INTENT(OUT)::VegFraction
2588 REAL(KIND(1D0)),
INTENT(OUT)::ImpervFraction
2589 REAL(KIND(1D0)),
INTENT(OUT)::PervFraction
2590 REAL(KIND(1D0)),
INTENT(OUT)::NonWaterFraction
2594 pervfraction = 1 - impervfraction
2736 ELEMENTAL FUNCTION set_nan(x)
RESULT(xx)
2738 REAL(KIND(1d0)),
PARAMETER::pNAN = 30000
2739 REAL(KIND(1d0)),
PARAMETER::NAN = -999
2740 REAL(KIND(1d0)),
INTENT(in)::x
2743 IF (abs(x) > pnan)
THEN 2753 FUNCTION square(x)
RESULT(xx)
2755 REAL(KIND(1d0)),
PARAMETER::pNAN = 9999
2756 REAL(KIND(1d0)),
PARAMETER::NAN = -999
2757 REAL(KIND(1d0)),
INTENT(in)::x
2760 xx = x**2 + nan/pnan
2767 REAL,
PARAMETER::pNAN = 9999
2768 REAL,
PARAMETER::NAN = -999
2772 xx = x**2 + nan/pnan
2782 INTEGER,
INTENT(in) :: i
2783 CHARACTER(len=15),
INTENT(out) :: name, group, aggreg
2787 IF (i < nvar .AND. i > 0)
THEN 2804 INTEGER,
INTENT(out) :: nVar
2811 MetForcingBlock, len_sim, &
2812 AerodynamicResistanceMethod, AH_MIN, AHProf_24hr, AH_SLOPE_Cooling, & ! input&inout in alphabetical order
2814 alb, AlbMax_DecTr, AlbMax_EveTr, AlbMax_Grass, &
2815 AlbMin_DecTr, AlbMin_EveTr, AlbMin_Grass, &
2816 alpha_bioCO2, alpha_enh_bioCO2, alt, BaseT, BaseTe, &
2817 BaseTHDD, beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, &
2818 chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, &
2819 DecTreeH, Diagnose, DiagQN, DiagQS, DRAINRT, &
2820 dt_since_start, dqndt, qn1_av, dqnsdt, qn1_s_av, &
2821 EF_umolCO2perJ, emis, EmissionsMethod, EnEF_v_Jkm, endDLS, EveTreeH, FAIBldg, &
2822 FAIDecTree, FAIEveTree, Faut, FcEF_v_kgkm, FlowChange, &
2823 FrFossilFuel_Heat, FrFossilFuel_NonHeat, G1, G2, G3, G4, G5, G6, GDD_id, &
2824 GDDFull, Gridiv, gsModel, HDD_id, HumActivity_24hr, &
2825 IceFrac, Ie_a, Ie_end, Ie_m, Ie_start, &
2826 InternalWaterUse_h, IrrFracConif, IrrFracDecid, IrrFracGrass, EvapMethod, &
2827 kkAnOHM, Kmax, LAI_id, LAICalcYes, LAIMax, LAIMin, &
2828 LAIPower, LAIType, lat, lng, MaxConductance, MaxFCMetab, MaxQFMetab, &
2829 SnowWater, MinFCMetab, MinQFMetab, min_res_bioCO2, &
2830 NARP_EMIS_SNOW, NARP_TRANS_SITE, NetRadiationMethod, &
2831 OHM_coef, OHMIncQF, OHM_threshSW, &
2832 OHM_threshWD, PipeCapacity, PopDensDaytime, &
2833 PopDensNighttime, PopProf_24hr, PorMax_dec, PorMin_dec, &
2834 PrecipLimit, PrecipLimitAlb, &
2835 QF0_BEU, Qf_A, Qf_B, Qf_C, &
2836 RadMeltFact, RAINCOVER, RainMaxRes, resp_a, resp_b, &
2837 RoughLenHeatMethod, RoughLenMomMethod, RunoffToWater, S1, S2, &
2838 SatHydraulicConduct, SDDFull, SDD_id, sfr, SMDMethod, SnowAlb, SnowAlbMax, &
2839 SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, &
2840 SnowLimBldg, SnowLimPaved, SnowPack, SnowProf_24hr, snowUse, SoilDepth, &
2841 soilstore_id, SoilStoreCap, StabilityMethod, startDLS, state_id, StateLimit, &
2842 StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, &
2843 T_CRITIC_Cooling, T_CRITIC_Heating, TempMeltFact, TH, &
2844 theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, &
2845 Tmin_id, Tmax_id, lenday_id, &
2846 TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, &
2847 WaterDist, WaterUseMethod, WetThresh, &
2848 WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, &
2849 WUProfA_24hr, WUProfM_24hr, Z, z0m_in, zdm_in, &
2850 dataOutBlockSUEWS, dataOutBlockSnow, dataOutBlockESTM, dataOutBlockRSL, &!output
2856 REAL(KIND(1D0)),
DIMENSION(len_sim, 24),
INTENT(IN) ::MetForcingBlock
2857 INTEGER,
INTENT(IN) :: len_sim
2859 INTEGER,
INTENT(IN)::AerodynamicResistanceMethod
2860 INTEGER,
INTENT(IN)::Diagnose
2861 INTEGER,
INTENT(IN)::DiagQN
2862 INTEGER,
INTENT(IN)::DiagQS
2863 INTEGER,
INTENT(IN)::startDLS
2864 INTEGER,
INTENT(IN)::endDLS
2865 INTEGER,
INTENT(IN)::EmissionsMethod
2866 INTEGER,
INTENT(IN)::Gridiv
2867 INTEGER,
INTENT(IN)::gsModel
2868 INTEGER,
INTENT(IN)::Ie_end
2869 INTEGER,
INTENT(IN)::Ie_start
2870 INTEGER,
INTENT(IN)::EvapMethod
2871 INTEGER,
INTENT(IN)::LAICalcYes
2872 INTEGER,
INTENT(IN)::NetRadiationMethod
2873 INTEGER,
INTENT(IN)::OHMIncQF
2874 INTEGER,
INTENT(IN)::RoughLenHeatMethod
2875 INTEGER,
INTENT(IN)::RoughLenMomMethod
2876 INTEGER,
INTENT(IN)::SMDMethod
2877 INTEGER,
INTENT(IN)::snowUse
2878 INTEGER,
INTENT(IN)::StabilityMethod
2879 INTEGER,
INTENT(IN)::StorageHeatMethod
2880 INTEGER,
INTENT(IN)::tstep
2881 INTEGER,
INTENT(IN)::tstep_prev
2883 INTEGER,
INTENT(inout)::dt_since_start
2884 INTEGER,
INTENT(IN)::veg_type
2885 INTEGER,
INTENT(IN)::WaterUseMethod
2887 INTEGER,
DIMENSION(NVEGSURF),
INTENT(IN)::LAIType
2889 REAL(KIND(1D0)),
INTENT(IN)::AlbMax_DecTr
2890 REAL(KIND(1D0)),
INTENT(IN)::AlbMax_EveTr
2891 REAL(KIND(1D0)),
INTENT(IN)::AlbMax_Grass
2892 REAL(KIND(1D0)),
INTENT(IN)::AlbMin_DecTr
2893 REAL(KIND(1D0)),
INTENT(IN)::AlbMin_EveTr
2894 REAL(KIND(1D0)),
INTENT(IN)::AlbMin_Grass
2895 REAL(KIND(1D0)),
INTENT(IN)::alt
2899 REAL(KIND(1D0)),
INTENT(IN)::BaseTHDD
2900 REAL(KIND(1D0)),
INTENT(IN)::bldgH
2901 REAL(KIND(1D0)),
INTENT(IN)::CapMax_dec
2902 REAL(KIND(1D0)),
INTENT(IN)::CapMin_dec
2903 REAL(KIND(1D0)),
INTENT(IN)::CO2PointSource
2904 REAL(KIND(1D0)),
INTENT(IN)::CRWmax
2905 REAL(KIND(1D0)),
INTENT(IN)::CRWmin
2906 REAL(KIND(1D0)),
INTENT(IN)::DecTreeH
2907 REAL(KIND(1D0)),
INTENT(IN)::DRAINRT
2908 REAL(KIND(1D0)),
INTENT(IN)::EF_umolCO2perJ
2909 REAL(KIND(1D0)),
INTENT(IN)::EnEF_v_Jkm
2910 REAL(KIND(1D0)),
INTENT(IN)::EveTreeH
2911 REAL(KIND(1D0)),
INTENT(IN)::FAIBldg
2912 REAL(KIND(1D0)),
INTENT(IN)::FAIDecTree
2913 REAL(KIND(1D0)),
INTENT(IN)::FAIEveTree
2914 REAL(KIND(1D0)),
INTENT(IN)::Faut
2916 REAL(KIND(1D0)),
INTENT(IN)::FlowChange
2917 REAL(KIND(1D0)),
INTENT(IN)::FrFossilFuel_Heat
2918 REAL(KIND(1D0)),
INTENT(IN)::FrFossilFuel_NonHeat
2919 REAL(KIND(1D0)),
INTENT(IN)::G1
2920 REAL(KIND(1D0)),
INTENT(IN)::G2
2921 REAL(KIND(1D0)),
INTENT(IN)::G3
2922 REAL(KIND(1D0)),
INTENT(IN)::G4
2923 REAL(KIND(1D0)),
INTENT(IN)::G5
2924 REAL(KIND(1D0)),
INTENT(IN)::G6
2925 REAL(KIND(1D0)),
INTENT(IN)::InternalWaterUse_h
2926 REAL(KIND(1D0)),
INTENT(IN)::IrrFracConif
2927 REAL(KIND(1D0)),
INTENT(IN)::IrrFracDecid
2928 REAL(KIND(1D0)),
INTENT(IN)::IrrFracGrass
2929 REAL(KIND(1D0)),
INTENT(IN)::Kmax
2931 REAL(KIND(1D0)),
INTENT(IN)::lat
2933 REAL(KIND(1D0)),
INTENT(IN)::lng
2934 REAL(KIND(1D0)),
INTENT(IN)::MaxFCMetab
2935 REAL(KIND(1D0)),
INTENT(IN)::MaxQFMetab
2936 REAL(KIND(1D0)),
INTENT(IN)::MinFCMetab
2937 REAL(KIND(1D0)),
INTENT(IN)::MinQFMetab
2938 REAL(KIND(1D0)),
INTENT(IN)::NARP_EMIS_SNOW
2939 REAL(KIND(1D0)),
INTENT(IN)::NARP_TRANS_SITE
2940 REAL(KIND(1D0)),
INTENT(IN)::PipeCapacity
2941 REAL(KIND(1D0)),
INTENT(IN)::PopDensNighttime
2942 REAL(KIND(1D0)),
INTENT(IN)::PorMax_dec
2943 REAL(KIND(1D0)),
INTENT(IN)::PorMin_dec
2945 REAL(KIND(1D0)),
INTENT(IN)::PrecipLimit
2946 REAL(KIND(1D0)),
INTENT(IN)::PrecipLimitAlb
2952 REAL(KIND(1D0)),
INTENT(IN)::RadMeltFact
2953 REAL(KIND(1D0)),
INTENT(IN)::RAINCOVER
2954 REAL(KIND(1D0)),
INTENT(IN)::RainMaxRes
2955 REAL(KIND(1D0)),
INTENT(IN)::RunoffToWater
2956 REAL(KIND(1D0)),
INTENT(IN)::S1
2957 REAL(KIND(1D0)),
INTENT(IN)::S2
2958 REAL(KIND(1D0)),
INTENT(IN)::SnowAlbMax
2959 REAL(KIND(1D0)),
INTENT(IN)::SnowAlbMin
2960 REAL(KIND(1D0)),
INTENT(IN)::SnowDensMax
2961 REAL(KIND(1D0)),
INTENT(IN)::SnowDensMin
2962 REAL(KIND(1D0)),
INTENT(IN)::SnowLimBldg
2963 REAL(KIND(1D0)),
INTENT(IN)::SnowLimPaved
2965 REAL(KIND(1D0)),
INTENT(IN)::SurfaceArea
2966 REAL(KIND(1D0)),
INTENT(IN)::tau_a
2967 REAL(KIND(1D0)),
INTENT(IN)::tau_f
2968 REAL(KIND(1D0)),
INTENT(IN)::tau_r
2970 REAL(KIND(1D0)),
INTENT(IN)::TempMeltFact
2971 REAL(KIND(1D0)),
INTENT(IN)::TH
2972 REAL(KIND(1D0)),
INTENT(IN)::timezone
2973 REAL(KIND(1D0)),
INTENT(IN)::TL
2974 REAL(KIND(1D0)),
INTENT(IN)::TrafficUnits
2976 REAL(KIND(1D0)),
INTENT(IN)::Z
2977 REAL(KIND(1D0)),
INTENT(IN)::z0m_in
2978 REAL(KIND(1D0)),
INTENT(IN)::zdm_in
2980 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::AH_MIN
2981 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::AH_SLOPE_Cooling
2982 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::AH_SLOPE_Heating
2983 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::FcEF_v_kgkm
2984 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::QF0_BEU
2985 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::Qf_A
2986 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::Qf_B
2987 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::Qf_C
2989 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::PopDensDaytime
2990 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::T_CRITIC_Cooling
2991 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::T_CRITIC_Heating
2992 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) ::TrafficRate
2993 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) ::Ie_a
2994 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) ::Ie_m
2995 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) ::MaxConductance
2996 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) ::DayWat
2997 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) ::DayWatPer
2998 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN)::OHM_threshSW
2999 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN)::OHM_threshWD
3000 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::chAnOHM
3001 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::cpAnOHM
3002 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::emis
3003 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::kkAnOHM
3004 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SatHydraulicConduct
3005 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::sfr
3006 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SnowPackLimit
3007 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SoilDepth
3008 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::SoilStoreCap
3009 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::StateLimit
3010 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) ::WetThresh
3011 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::alpha_bioCO2
3012 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::alpha_enh_bioCO2
3013 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::BaseT
3014 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::BaseTe
3015 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::beta_bioCO2
3016 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::beta_enh_bioCO2
3017 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::GDDFull
3018 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::LAIMax
3019 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::LAIMin
3020 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::min_res_bioCO2
3021 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::resp_a
3022 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::resp_b
3023 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::SDDFull
3024 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::SnowProf_24hr
3025 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) ::theta_bioCO2
3026 REAL(KIND(1D0)),
DIMENSION(4, NVEGSURF),
INTENT(IN) ::LAIPower
3027 REAL(KIND(1D0)),
DIMENSION(nsurf + 1, 4, 3),
INTENT(IN) ::OHM_coef
3028 REAL(KIND(1D0)),
DIMENSION(NSURF + 1, NSURF - 1),
INTENT(IN) ::WaterDist
3029 REAL(KIND(1d0)),
DIMENSION(:),
INTENT(IN) ::Ts5mindata_ir
3032 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::AHProf_24hr
3033 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::HumActivity_24hr
3034 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::PopProf_24hr
3035 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::TraffProf_24hr
3036 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::WUProfA_24hr
3037 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) ::WUProfM_24hr
3043 REAL(KIND(1d0)),
INTENT(INOUT) ::qn1_av
3044 REAL(KIND(1d0)),
INTENT(INOUT) ::dqndt
3045 REAL(KIND(1d0)),
INTENT(INOUT) ::qn1_s_av
3046 REAL(KIND(1d0)),
INTENT(INOUT) ::dqnsdt
3049 REAL(KIND(1D0)),
INTENT(INOUT) ::SnowfallCum
3050 REAL(KIND(1D0)),
INTENT(INOUT) ::SnowAlb
3051 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::IceFrac
3052 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::SnowWater
3053 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::SnowDens
3054 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::SnowFrac
3055 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::SnowPack
3058 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::soilstore_id
3059 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::state_id
3060 REAL(KIND(1D0)),
DIMENSION(6, NSURF),
INTENT(INOUT) ::StoreDrainPrm
3063 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) ::alb
3064 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(INOUT) ::GDD_id
3065 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(INOUT) ::SDD_id
3066 REAL(KIND(1d0)),
DIMENSION(nvegsurf),
INTENT(INOUT)::LAI_id
3067 REAL(KIND(1d0)),
INTENT(INOUT) :: DecidCap_id
3068 REAL(KIND(1d0)),
INTENT(INOUT) :: albDecTr_id
3069 REAL(KIND(1d0)),
INTENT(INOUT) :: albEveTr_id
3070 REAL(KIND(1d0)),
INTENT(INOUT) :: albGrass_id
3071 REAL(KIND(1d0)),
INTENT(INOUT) :: porosity_id
3072 REAL(KIND(1d0)),
INTENT(INOUT) :: Tmin_id
3073 REAL(KIND(1d0)),
INTENT(INOUT) :: Tmax_id
3074 REAL(KIND(1d0)),
INTENT(INOUT) :: lenday_id
3077 REAL(KIND(1d0)),
DIMENSION(12),
INTENT(INOUT) ::HDD_id
3080 REAL(KIND(1d0)),
DIMENSION(9),
INTENT(INOUT) ::WUDay_id
3083 REAL(KIND(1d0)),
INTENT(INOUT) ::Tair_av
3089 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSUEWS),
INTENT(OUT) ::dataOutBlockSUEWS
3090 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSnow),
INTENT(OUT) ::dataOutBlockSnow
3091 REAL(KIND(1d0)),
DIMENSION(len_sim, ncolumnsDataOutESTM),
INTENT(OUT) ::dataOutBlockESTM
3092 REAL(KIND(1d0)),
DIMENSION(len_sim, ncolumnsDataOutRSL),
INTENT(OUT) ::dataOutBlockRSL
3093 REAL(KIND(1d0)),
DIMENSION(len_sim, ncolumnsDataOutDailyState),
INTENT(OUT) ::DailyStateBlock
3110 INTEGER,
PARAMETER :: gridiv_x = 1
3111 REAL(KIND(1D0))::qn1_obs
3112 REAL(KIND(1D0))::qh_obs
3113 REAL(KIND(1D0))::qe_obs
3114 REAL(KIND(1D0))::qs_obs
3115 REAL(KIND(1D0))::qf_obs
3116 REAL(KIND(1D0))::avu1
3117 REAL(KIND(1D0))::avrh
3118 REAL(KIND(1D0))::Temp_C
3119 REAL(KIND(1D0))::Press_hPa
3120 REAL(KIND(1D0))::Precip
3121 REAL(KIND(1D0))::avkdn
3122 REAL(KIND(1D0))::snowFrac_obs
3123 REAL(KIND(1D0))::ldown_obs
3124 REAL(KIND(1D0))::fcld_obs
3125 REAL(KIND(1D0))::wu_m3
3126 REAL(KIND(1D0))::xsmd
3127 REAL(KIND(1D0))::LAI_obs
3128 REAL(KIND(1D0))::kdiff
3129 REAL(KIND(1D0))::kdir
3130 REAL(KIND(1D0))::wdir
3132 REAL(KIND(1D0)),
DIMENSION(5)::datetimeLine
3133 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSUEWS - 5)::dataOutLineSUEWS
3134 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSnow - 5)::dataOutLineSnow
3135 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutESTM - 5)::dataOutLineESTM
3136 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutRSL - 5)::dataOutLineRSL
3137 REAL(KIND(1d0)),
DIMENSION(ncolumnsDataOutDailyState - 5)::DailyStateLine
3139 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSUEWS, 1) ::dataOutBlockSUEWS_X
3140 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSnow, 1) ::dataOutBlockSnow_X
3141 REAL(KIND(1d0)),
DIMENSION(len_sim, ncolumnsDataOutESTM, 1) ::dataOutBlockESTM_X
3142 REAL(KIND(1d0)),
DIMENSION(len_sim, ncolumnsDataOutRSL, 1) ::dataOutBlockRSL_X
3145 REAL(KIND(1D0)),
DIMENSION(10, 10) ::MetForcingData_grid
3147 CHARACTER(len=150):: FileStateInit
3148 CHARACTER(len=4):: year_txt
3149 CHARACTER(len=3):: id_text
3150 CHARACTER(len=2):: it_text, imin_text
3155 DO ir = 1, len_sim, 1
3159 iy = int(metforcingblock(ir, 1))
3160 id = int(metforcingblock(ir, 2))
3161 it = int(metforcingblock(ir, 3))
3162 imin = int(metforcingblock(ir, 4))
3164 qn1_obs = metforcingblock(ir, 5)
3165 qh_obs = metforcingblock(ir, 6)
3166 qe_obs = metforcingblock(ir, 7)
3167 qs_obs = metforcingblock(ir, 8)
3168 qf_obs = metforcingblock(ir, 9)
3169 avu1 = metforcingblock(ir, 10)
3170 avrh = metforcingblock(ir, 11)
3171 temp_c = metforcingblock(ir, 12)
3172 press_hpa = metforcingblock(ir, 13)
3173 precip = metforcingblock(ir, 14)
3174 avkdn = metforcingblock(ir, 15)
3175 snowfrac_obs = metforcingblock(ir, 16)
3176 ldown_obs = metforcingblock(ir, 17)
3177 fcld_obs = metforcingblock(ir, 18)
3178 wu_m3 = metforcingblock(ir, 19)
3179 xsmd = metforcingblock(ir, 20)
3180 lai_obs = metforcingblock(ir, 21)
3181 kdiff = metforcingblock(ir, 22)
3182 kdir = metforcingblock(ir, 23)
3183 wdir = metforcingblock(ir, 24)
3417 aerodynamicresistancemethod, ah_min, ahprof_24hr, ah_slope_cooling, &
3419 alb, albmax_dectr, albmax_evetr, albmax_grass, &
3420 albmin_dectr, albmin_evetr, albmin_grass, &
3421 alpha_bioco2, alpha_enh_bioco2, alt, avkdn, avrh, avu1, baset, basete, &
3422 basethdd, beta_bioco2, beta_enh_bioco2, bldgh, capmax_dec, capmin_dec, &
3423 chanohm, co2pointsource, cpanohm, crwmax, crwmin, daywat, daywatper, &
3424 dectreeh, diagnose, diagqn, diagqs, drainrt, &
3425 dt_since_start, dqndt, qn1_av, dqnsdt, qn1_s_av, &
3426 ef_umolco2perj, emis, emissionsmethod, enef_v_jkm, enddls, evetreeh, faibldg, &
3427 faidectree, faievetree, faut, fcef_v_kgkm, fcld_obs, flowchange, &
3428 frfossilfuel_heat, frfossilfuel_nonheat, g1, g2, g3, g4, g5, g6, gdd_id, &
3429 gddfull, gridiv, gsmodel, hdd_id, humactivity_24hr, &
3430 icefrac, id, ie_a, ie_end, ie_m, ie_start, imin, &
3431 internalwateruse_h, irrfracconif, irrfracdecid, irrfracgrass, isec, it, evapmethod, &
3432 iy, kkanohm, kmax, lai_id, laicalcyes, laimax, laimin, lai_obs, &
3433 laipower, laitype, lat, lenday_id, ldown_obs, lng, maxconductance, maxfcmetab, maxqfmetab, &
3434 snowwater, metforcingdata_grid, minfcmetab, minqfmetab, min_res_bioco2, &
3435 narp_emis_snow, narp_trans_site, netradiationmethod, &
3436 ohm_coef, ohmincqf, ohm_threshsw, &
3437 ohm_threshwd, pipecapacity, popdensdaytime, &
3438 popdensnighttime, popprof_24hr, pormax_dec, pormin_dec, &
3439 precip, preciplimit, preciplimitalb, press_hpa, &
3440 qf0_beu, qf_a, qf_b, qf_c, &
3441 qn1_obs, qs_obs, qf_obs, &
3442 radmeltfact, raincover, rainmaxres, resp_a, resp_b, &
3443 roughlenheatmethod, roughlenmommethod, runofftowater, s1, s2, &
3444 sathydraulicconduct, sddfull, sdd_id, sfr, smdmethod, snowalb, snowalbmax, &
3445 snowalbmin, snowpacklimit, snowdens, snowdensmax, snowdensmin, snowfallcum, snowfrac, &
3446 snowlimbldg, snowlimpaved, snowfrac_obs, snowpack, snowprof_24hr, snowuse, soildepth, &
3447 soilstore_id, soilstorecap, stabilitymethod, startdls, state_id, statelimit, &
3448 storageheatmethod, storedrainprm, surfacearea, tair_av, tau_a, tau_f, tau_r, &
3450 t_critic_cooling, t_critic_heating, temp_c, tempmeltfact, th, &
3451 theta_bioco2, timezone, tl, trafficrate, trafficunits, &
3452 traffprof_24hr, ts5mindata_ir, tstep, tstep_prev, veg_type, &
3453 waterdist, waterusemethod, wetthresh, wu_m3, &
3454 wuday_id, decidcap_id, albdectr_id, albevetr_id, albgrass_id, porosity_id, &
3455 wuprofa_24hr, wuprofm_24hr, xsmd, z, z0m_in, zdm_in, &
3456 datetimeline, dataoutlinesuews, dataoutlinesnow, dataoutlineestm, dataoutlinersl, &
3460 dt_since_start = dt_since_start + tstep
3463 dailystateblock(ir, :) = [datetimeline, dailystateline]
3468 snowuse, storageheatmethod, &
3470 ir, gridiv_x, datetimeline, dataoutlinesuews, dataoutlinesnow, dataoutlineestm, dataoutlinersl, &
3471 dataoutblocksuews_x, dataoutblocksnow_x, dataoutblockestm_x, dataoutblockrsl_x)
3475 dataoutblocksuews = dataoutblocksuews_x(:, :, 1)
3476 dataoutblocksnow = dataoutblocksnow_x(:, :, 1)
3477 dataoutblockestm = dataoutblockestm_x(:, :, 1)
3478 dataoutblockrsl = dataoutblockrsl_x(:, :, 1)
3485 year, idectime, UTC, locationlatitude, locationlongitude, locationaltitude, & !input
3486 sunazimuth, sunzenith)
3489 REAL(KIND(1D0)),
INTENT(in) :: year, idectime, UTC, &
3490 locationlatitude, locationlongitude, locationaltitude
3491 REAL(KIND(1D0)),
INTENT(out) ::sunazimuth, sunzenith
3494 year, idectime, utc, locationlatitude, locationlongitude, locationaltitude, &
3495 sunazimuth, sunzenith)
3506 function cal_tair_av(tair_av_prev, dt_since_start, tstep, temp_c)
result(tair_av_next)
3510 real(KIND(1D0)),
intent(in) :: tair_av_prev
3511 real(KIND(1D0)),
intent(in) :: temp_c
3512 integer,
intent(in) :: dt_since_start
3513 integer,
intent(in) :: tstep
3515 real(KIND(1D0)) ::tair_av_next
3517 real(KIND(1D0)),
parameter:: len_day_s = 24*3600
3518 real(KIND(1D0)):: len_cal_s
3519 real(KIND(1D0)):: temp_k
3522 if (dt_since_start > len_day_s)
then 3524 len_cal_s = len_day_s
3527 len_cal_s = dt_since_start + tstep
3529 temp_k = temp_c + 273.15
3530 tair_av_next = tair_av_prev*(len_cal_s - tstep*1.)/len_cal_s + temp_k*tstep/len_cal_s
3534 function cal_tsfc(qh, avdens, avcp, RA, temp_c)
result(tsfc_C)
3538 real(KIND(1D0)),
intent(in) :: qh
3539 real(KIND(1D0)),
intent(in) :: avdens
3540 real(KIND(1D0)),
intent(in) :: avcp
3541 real(KIND(1D0)),
intent(in) :: RA
3542 real(KIND(1D0)),
intent(in) :: temp_C
3544 real(KIND(1D0)) ::tsfc_C
3546 tsfc_c = qh/(avdens*avcp)*ra + temp_c
subroutine suews_update_soilmoist(NonWaterFraction, SoilStoreCap, sfr, soilstore_id, SoilMoistCap, SoilState, vsmd, smd)
subroutine anohm(tstep, dt_since_start, qn1, qn1_av_prev, dqndt_prev, qf, MetForcingData_grid, moist_surf, alb, emis, cpAnOHM, kkAnOHM, chAnOHM, sfr, nsurf, EmissionsMethod, id, Gridiv, qn1_av_next, dqndt_next, a1, a2, a3, qs, deltaQi)
High level wrapper for AnOHM calculation.
subroutine suews_cal_soilstate(SMDMethod, xsmd, NonWaterFraction, SoilMoistCap, SoilStoreCap, surf_chang_per_tstep, soilstore_id, soilstoreOld, sfr, smd, smd_nsurf, tot_chang_per_tstep, SoilState)
subroutine co2_biogen(alpha_bioCO2, alpha_enh_bioCO2, avkdn, beta_bioCO2, beta_enh_bioCO2, BSoilSurf, ConifSurf, DecidSurf, dectime, EmissionsMethod, gfunc, gfunc2, GrassSurf, gsmodel, id, it, ivConif, ivDecid, ivGrass, LAI_id, LAIMin, LAIMax, min_res_bioCO2, nsurf, NVegSurf, resp_a, resp_b, sfr, SnowFrac, t2, Temp_C, theta_bioCO2, Fc_biogen, Fc_photo, Fc_respi)
subroutine snowcalc(tstep, imin, it, dectime, is, EvapMethod, CRWmin, CRWmax, nsh_real, lvS_J_kg, lv_j_kg, avdens, avRh, Press_hPa, Temp_C, RAsnow, psyc_hPa, avcp, sIce_hPa, PervFraction, vegfraction, addimpervious, vpd_hPa, qn_e, s_hPa, ResistSurf, RA, rb, tlv, snowdensmin, SnowProf_24hr, precip, PipeCapacity, RunoffToWater, addVeg, SnowLimPaved, SnowLimBldg, FlowChange, drain, WetThresh, stateOld, mw_ind, SoilStoreCap, rainonsnow, freezmelt, freezstate, freezstatevol, Qm_Melt, Qm_rain, Tsurf_ind, sfr, dayofWeek_id, StoreDrainPrm, SnowPackLimit, AddWater, addwaterrunoff, soilstore_id, SnowPack, SurplusEvap, SnowFrac, SnowWater, iceFrac, SnowDens, runoffAGimpervious, runoffAGveg, surplusWaterBody, rss_nsurf, runoffSnow, runoff, runoffSoil, chang, changSnow, SnowToSurf, state_id, ev_snow, SnowDepth, SnowRemoval, swe, ev, chSnow_tot, ev_tot, qe_tot, runoff_tot, surf_chang_tot, runoffPipes, mwstore, runoffwaterbody)
subroutine redistributewater(snowUse, WaterDist, sfr, Drain, AddWaterRunoff, AddWater)
subroutine suews_cal_biogenco2(alpha_bioCO2, alpha_enh_bioCO2, avkdn, avRh, beta_bioCO2, beta_enh_bioCO2, BSoilSurf, ConifSurf, DecidSurf, dectime, Diagnose, EmissionsMethod, Fc_anthro, G1, G2, G3, G4, G5, G6, gfunc, GrassSurf, gsmodel, id, it, ivConif, ivDecid, ivGrass, Kmax, LAI_id, LAIMin, LAIMax, MaxConductance, min_res_bioCO2, nsurf, NVegSurf, Press_hPa, resp_a, resp_b, S1, S2, sfr, SMDMethod, SnowFrac, t2_C, Temp_C, theta_bioCO2, TH, TL, vsmd, xsmd, Fc, Fc_biogen, Fc_photo, Fc_respi)
integer, parameter ncolumnsdataoutestm
subroutine suews_init_qh(avdens, avcp, h_mod, qn1, dectime, H_init)
integer, parameter ivgrass
subroutine suews_cal_dailystate(iy, id, it, imin, isec, tstep, tstep_prev, dt_since_start, DayofWeek_id, Tmin_id_prev, Tmax_id_prev, lenDay_id_prev, WaterUseMethod, Ie_start, Ie_end, LAICalcYes, LAIType, nsh_real, avkdn, Temp_C, Precip, BaseTHDD, lat, Faut, LAI_obs, AlbMax_DecTr, AlbMax_EveTr, AlbMax_Grass, AlbMin_DecTr, AlbMin_EveTr, AlbMin_Grass, CapMax_dec, CapMin_dec, PorMax_dec, PorMin_dec, Ie_a, Ie_m, DayWatPer, DayWat, BaseT, BaseTe, GDDFull, SDDFull, LAIMin, LAIMax, LAIPower, DecidCap_id_prev, StoreDrainPrm_prev, LAI_id_prev, GDD_id_prev, SDD_id_prev, albDecTr_id_prev, albEveTr_id_prev, albGrass_id_prev, porosity_id_prev, HDD_id_prev, HDD_id_next, Tmin_id_next, Tmax_id_next, lenDay_id_next, albDecTr_id_next, albEveTr_id_next, albGrass_id_next, porosity_id_next, DecidCap_id_next, StoreDrainPrm_next, LAI_id_next, GDD_id_next, SDD_id_next, deltaLAI, WUDay_id)
subroutine narp_cal_sunposition(year, idectime, UTC, locationlatitude, locationlongitude, locationaltitude, sunazimuth, sunzenith)
integer, parameter ncolumnsdataoutdailystate
integer, parameter bsoilsurf
subroutine suews_cal_water(Diagnose, snowUse, NonWaterFraction, addPipes, addImpervious, addVeg, addWaterBody, state_id, soilstore_id, sfr, StoreDrainPrm, WaterDist, nsh_real, drain_per_tstep, drain, AddWaterRunoff, AdditionalWater, runoffPipes, runoff_per_interval, AddWater, stateOld, soilstoreOld)
subroutine estm(Gridiv, tstep, avkdn, avu1, temp_c, zenith_deg, avrh, press_hpa, ldown, bldgh, Ts5mindata_ir, Tair_av, dataOutLineESTM, QS)
real(kind(1d0)) function qa2rh(qa_gkg, pres_hPa, Ta_degC)
subroutine snow_cal_meltheat(snowUse, tstep, tau_r, SnowDensMax, lvS_J_kg, lv_J_kg, tstep_real, RadMeltFact, TempMeltFact, SnowAlbMax, SnowDensMin, Temp_C, Precip, PrecipLimit, PrecipLimitAlb, nsh_real, sfr, Tsurf_ind, Tsurf_ind_snow, state_id, qn1_ind_snow, kup_ind_snow, SnowWater, deltaQi, alb1, SnowPack_in, SnowFrac_in, SnowAlb_in, SnowDens_in, SnowfallCum_in, SnowPack_out, SnowFrac_out, SnowAlb_out, SnowDens_out, SnowfallCum_out, mwh, Qm, QmFreez, QmRain, veg_fr, snowCalcSwitch, Qm_melt, Qm_freezState, Qm_rain, FreezMelt, FreezState, FreezStateVol, rainOnSnow, SnowDepth, mw_ind, dataOutLineSnow)
subroutine radmethod(NetRadiationMethod, snowUse, NetRadiationMethod_use, AlbedoChoice, ldown_option)
subroutine suews_cal_sunposition(year, idectime, UTC, locationlatitude, locationlongitude, locationaltitude, sunazimuth, sunzenith)
elemental real(kind(1d0)) function set_nan(x)
subroutine boundarylayerresistance(zzd, z0m, avU1, UStar, rb)
subroutine suews_update_dailystate(id, datetimeline, Gridiv, NumberOfGrids, DailyStateLine, dataOutDailyState)
subroutine cal_atmmoist(Temp_C, Press_hPa, avRh, dectime, lv_J_kg, lvS_J_kg, es_hPa, Ea_hPa, VPd_hpa, VPD_Pa, dq, dens_dry, avcp, air_dens)
subroutine lumps_cal_qhqe(veg_type, snowUse, qn1, qf, qs, Qm, Temp_C, Veg_Fr, avcp, Press_hPa, lv_J_kg, tstep_real, DRAINRT, nsh_real, Precip, RainMaxRes, RAINCOVER, sfr, LAI_id_prev, LAImax, LAImin, QH_LUMPS, QE_LUMPS, psyc_hPa, s_hPa, sIce_hpa, Veg_Fr_temp, VegPhenLumps)
subroutine suews_update_output(SnowUse, storageheatmethod, ReadLinesMetdata, NumberOfGrids, ir, gridiv, datetimeLine, dataOutLineSUEWS, dataOutLineSnow, dataOutLineESTM, dataoutLineURSL, dataOutSUEWS, dataOutSnow, dataOutESTM, dataOutRSL)
subroutine suews_cal_anthropogenicemission(AH_MIN, AHProf_24hr, AH_SLOPE_Cooling, AH_SLOPE_Heating, CO2PointSource, dayofWeek_id, Diagnose, DLS, EF_umolCO2perJ, EmissionsMethod, EnEF_v_Jkm, FcEF_v_kgkm, FrFossilFuel_Heat, FrFossilFuel_NonHeat, HDD_id, HumActivity_24hr, id, imin, it, MaxFCMetab, MaxQFMetab, MinFCMetab, MinQFMetab, nsh, PopDensDaytime, PopDensNighttime, PopProf_24hr, QF, QF0_BEU, Qf_A, Qf_B, Qf_C, QF_obs, QF_SAHP, sfr, SnowFrac, SurfaceArea, T_CRITIC_Cooling, T_CRITIC_Heating, Temp_C, TrafficRate, TrafficUnits, TraffProf_24hr, Fc_anthro, Fc_build, Fc_metab, Fc_point, Fc_traff)
subroutine snowupdate(tstep, Temp_C, tau_a, tau_f, tau_r, SnowDensMax, SnowDensMin, SnowAlbMax, SnowAlbMin, SnowPack_prev, SnowAlb_prev, SnowDens_prev, SnowAlb_next, SnowDens_next)
subroutine suews_cal_qn(NetRadiationMethod, snowUse, tstep, SnowPack_prev, tau_a, tau_f, SnowAlbMax, SnowAlbMin, Diagnose, snowFrac_obs, ldown_obs, fcld_obs, dectime, ZENITH_deg, Tsurf_0, avKdn, Temp_C, avRH, ea_hPa, qn1_obs, SnowAlb_prev, snowFrac_prev, DiagQN, NARP_TRANS_SITE, NARP_EMIS_SNOW, IceFrac, sfr, emis, alb_prev, albDecTr_id, albEveTr_id, albGrass_id, alb_next, ldown, fcld, qn1, qn1_snowfree, qn1_S, kclear, kup, lup, tsurf, qn1_ind_snow, kup_ind_snow, Tsurf_ind_snow, Tsurf_ind, alb1, snowFrac_next, SnowAlb_next)
subroutine ohm(qn1, qn1_av_prev, dqndt_prev, qn1_av_next, dqndt_next, qn1_S, qn1_s_av_prev, dqnsdt_prev, qn1_s_av_next, dqnsdt_next, tstep, dt_since_start, sfr, nsurf, Tair_mav_5d, OHM_coef, OHM_threshSW, OHM_threshWD, soilstore_id, SoilStoreCap, state_id, BldgSurf, WaterSurf, SnowUse, SnowFrac, DiagQS, a1, a2, a3, qs, deltaQi)
subroutine cal_stab(StabilityMethod, zzd, z0m, zdm, avU1, Temp_C, QH_init, avdens, avcp, L_MOD, TStar, UStar, zL)
integer, parameter conifsurf
subroutine suews_cal_tstep(tstep, nsh, nsh_real, tstep_real)
subroutine suews_cal_qh(QHMethod, qn1, qf, QmRain, qeOut, qs, QmFreez, qm, avdens, avcp, tsurf, Temp_C, RA, qh, qh_residual, qh_resist)
subroutine suews_cal_resistance(StabilityMethod, Diagnose, AerodynamicResistanceMethod, RoughLenHeatMethod, snowUse, id, it, gsModel, SMDMethod, avdens, avcp, QH_init, zzd, z0m, zdm, avU1, Temp_C, VegFraction, avkdn, Kmax, G1, G2, G3, G4, G5, G6, S1, S2, TH, TL, dq, xsmd, vsmd, MaxConductance, LAIMax, LAI_id, SnowFrac, sfr, UStar, TStar, L_mod, zL, gsc, ResistSurf, RA, RAsnow, rb)
subroutine surfaceresistance(id, it, SMDMethod, SnowFrac, sfr, avkdn, Temp_C, dq, xsmd, vsmd, MaxConductance, LAIMax, LAI_id, gsModel, Kmax, G1, G2, G3, G4, G5, G6, TH, TL, S1, S2, gfunc, gsc, ResistSurf)
real(kind(1d0)) function update_snow_albedo(tstep, SnowPack_prev, SnowAlb_prev, Temp_C, tau_a, tau_f, SnowAlbMax, SnowAlbMin)
subroutine suews_update_outputline(AdditionalWater, alb, avkdn, avU10_ms, azimuth, chSnow_per_interval, dectime, drain_per_tstep, E_mod, ev_per_tstep, ext_wu, Fc, Fc_build, fcld, Fc_metab, Fc_photo, Fc_respi, Fc_point, Fc_traff, FlowChange, h_mod, id, imin, int_wu, it, iy, kup, LAI_id, ldown, l_mod, lup, mwh, MwStore, nsh_real, NWstate_per_tstep, Precip, q2_gkg, qeOut, qf, qh, qh_resist, Qm, QmFreez, QmRain, qn1, qn1_S, qn1_snowfree, qs, RA, resistsurf, RH2, runoffAGimpervious, runoffAGveg, runoff_per_tstep, runoffPipes, runoffSoil_per_tstep, runoffWaterBody, sfr, smd, smd_nsurf, SnowAlb, SnowRemoval, state_id, state_per_tstep, surf_chang_per_tstep, swe, t2_C, tskin_C, tot_chang_per_tstep, tsurf, UStar, wu_DecTr, wu_EveTr, wu_Grass, z0m, zdm, zenith_deg, datetimeLine, dataOutLineSUEWS)
AnOHM: Analytical Objective Hysteresis Model.
real(kind(1d0)) function square(x)
subroutine update_dailystateline(it, imin, nsh_real, GDD_id, HDD_id, LAI_id, SDD_id, Tmin_id, Tmax_id, lenday_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, WUDay_id, deltaLAI, VegPhenLumps, SnowAlb, SnowDens, a1, a2, a3, DailyStateLine)
subroutine soilstore(is, sfr, PipeCapacity, RunoffToWater, pin, wu_EveTr, wu_DecTr, wu_Grass, drain, AddWater, addImpervious, nsh_real, stateOld, AddWaterRunoff, PervFraction, addVeg, SoilStoreCap, addWaterBody, FlowChange, StateLimit, runoffAGimpervious, surplusWaterBody, runoffAGveg, runoffPipes, ev, soilstore_id, SurplusEvap, runoffWaterBody, p_mm, chang, runoff, state_id)
real(kind(1d0)) function cal_z0v(RoughLenHeatMethod, z0m, VegFraction, UStar)
real(kind(1d0)) function stab_psi_heat(StabilityMethod, ZL)
integer, parameter ncolumnsdataoutrsl
integer, parameter ncolumnsdataoutsnow
subroutine rslprofile(Zh, z0m, zdm, L_MOD, sfr, planF, StabilityMethod, avcp, lv_J_kg, avdens, avU1, Temp_C, avRH, Press_hPa, zMeas, qh, qe, T2_C, q2_gkg, U10_ms, RH2, dataoutLineRSL)
integer, parameter grasssurf
subroutine suews_cal_main(AerodynamicResistanceMethod, 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, avkdn, avRh, avU1, BaseT, BaseTe, BaseTHDD, beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, DecTreeH, Diagnose, DiagQN, DiagQS, DRAINRT, dt_since_start, dqndt, qn1_av, dqnsdt, qn1_s_av, EF_umolCO2perJ, emis, EmissionsMethod, EnEF_v_Jkm, endDLS, EveTreeH, FAIBldg, FAIDecTree, FAIEveTree, Faut, FcEF_v_kgkm, fcld_obs, FlowChange, FrFossilFuel_Heat, FrFossilFuel_NonHeat, G1, G2, G3, G4, G5, G6, GDD_id, GDDFull, Gridiv, gsModel, HDD_id, HumActivity_24hr, IceFrac, id, Ie_a, Ie_end, Ie_m, Ie_start, imin, InternalWaterUse_h, IrrFracConif, IrrFracDecid, IrrFracGrass, isec, it, EvapMethod, iy, kkAnOHM, Kmax, LAI_id, LAICalcYes, LAIMax, LAIMin, LAI_obs, LAIPower, LAIType, lat, lenDay_id, ldown_obs, lng, MaxConductance, MaxFCMetab, MaxQFMetab, SnowWater, MetForcingData_grid, MinFCMetab, MinQFMetab, min_res_bioCO2, NARP_EMIS_SNOW, NARP_TRANS_SITE, NetRadiationMethod, 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, sfr, SMDMethod, SnowAlb, SnowAlbMax, SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, SnowLimBldg, SnowLimPaved, snowFrac_obs, SnowPack, SnowProf_24hr, snowUse, SoilDepth, soilstore_id, SoilStoreCap, StabilityMethod, startDLS, state_id, StateLimit, StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, Tmax_id, Tmin_id, T_CRITIC_Cooling, T_CRITIC_Heating, Temp_C, TempMeltFact, TH, theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, WaterDist, WaterUseMethod, WetThresh, wu_m3, WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, WUProfA_24hr, WUProfM_24hr, xsmd, Z, z0m_in, zdm_in, datetimeLine, dataOutLineSUEWS, dataOutLineSnow, dataOutLineESTM, dataoutLineRSL, DailyStateLine)
real(kind(1d0)) function, dimension(nsurf) update_snow_dens(tstep, SnowPack_prev, SnowDens_prev, tau_r, SnowDensMax, SnowDensMin)
integer, parameter ivconif
subroutine suews_cal_multitsteps(MetForcingBlock, len_sim, AerodynamicResistanceMethod, 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, BaseT, BaseTe, BaseTHDD, beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, DecTreeH, Diagnose, DiagQN, DiagQS, DRAINRT, dt_since_start, dqndt, qn1_av, dqnsdt, qn1_s_av, EF_umolCO2perJ, emis, EmissionsMethod, EnEF_v_Jkm, endDLS, EveTreeH, FAIBldg, FAIDecTree, FAIEveTree, Faut, FcEF_v_kgkm, FlowChange, FrFossilFuel_Heat, FrFossilFuel_NonHeat, G1, G2, G3, G4, G5, G6, GDD_id, GDDFull, Gridiv, gsModel, HDD_id, HumActivity_24hr, IceFrac, Ie_a, Ie_end, Ie_m, Ie_start, InternalWaterUse_h, IrrFracConif, IrrFracDecid, IrrFracGrass, EvapMethod, kkAnOHM, Kmax, LAI_id, LAICalcYes, LAIMax, LAIMin, LAIPower, LAIType, lat, lng, MaxConductance, MaxFCMetab, MaxQFMetab, SnowWater, MinFCMetab, MinQFMetab, min_res_bioCO2, NARP_EMIS_SNOW, NARP_TRANS_SITE, NetRadiationMethod, OHM_coef, OHMIncQF, OHM_threshSW, OHM_threshWD, PipeCapacity, PopDensDaytime, PopDensNighttime, PopProf_24hr, PorMax_dec, PorMin_dec, PrecipLimit, PrecipLimitAlb, QF0_BEU, Qf_A, Qf_B, Qf_C, RadMeltFact, RAINCOVER, RainMaxRes, resp_a, resp_b, RoughLenHeatMethod, RoughLenMomMethod, RunoffToWater, S1, S2, SatHydraulicConduct, SDDFull, SDD_id, sfr, SMDMethod, SnowAlb, SnowAlbMax, SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, SnowLimBldg, SnowLimPaved, SnowPack, SnowProf_24hr, snowUse, SoilDepth, soilstore_id, SoilStoreCap, StabilityMethod, startDLS, state_id, StateLimit, StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, T_CRITIC_Cooling, T_CRITIC_Heating, TempMeltFact, TH, theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, Tmin_id, Tmax_id, lenday_id, TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, WaterDist, WaterUseMethod, WetThresh, WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, WUProfA_24hr, WUProfM_24hr, Z, z0m_in, zdm_in, dataOutBlockSUEWS, dataOutBlockSnow, dataOutBlockESTM, dataOutBlockRSL, DailyStateBlock)
subroutine output_name_n(i, name, group, aggreg)
type(varattr), dimension(500) varlistall
real(kind(1d0)) function rh2qa(RH_dec, pres_hPa, Ta_degC)
subroutine suews_cal_horizontalsoilwater(sfr, SoilStoreCap, SoilDepth, SatHydraulicConduct, SurfaceArea, NonWaterFraction, tstep_real, soilstore_id, runoffSoil, runoffSoil_per_tstep)
subroutine suews_cal_dectime(id, it, imin, isec, dectime)
integer, parameter nvegsurf
subroutine output_size(nVar)
subroutine evap_suews(EvapMethod, state_is, WetThresh_is, capStore_is, vpd_hPa, avdens, avcp, qn_e, s_hPa, psyc_hPa, ResistSurf, RA, rb, tlv, rss, ev, qe)
subroutine suews_cal_weekday(iy, id, lat, dayofWeek_id)
subroutine narp(nsurf, sfr, SnowFrac, alb, emis, IceFrac, NARP_TRANS_SITE, NARP_EMIS_SNOW, DTIME, ZENITH_deg, tsurf_0, kdown, Temp_C, RH, Press_hPa, qn1_obs, SnowAlb, AlbedoChoice, ldown_option, NetRadiationMethod_use, DiagQN, QSTARall, QSTAR_SF, QSTAR_S, kclear, KUPall, LDOWN, LUPall, fcld, TSURFall, qn1_ind_snow, kup_ind_snow, Tsurf_ind_snow, Tsurf_ind, alb0, alb1)
subroutine suews_cal_wateruse(nsh_real, wu_m3, SurfaceArea, sfr, IrrFracConif, IrrFracDecid, IrrFracGrass, DayofWeek_id, WUProfA_24hr, WUProfM_24hr, InternalWaterUse_h, HDD_id, WUDay_id, WaterUseMethod, NSH, it, imin, DLS, wu_EveTr, wu_DecTr, wu_Grass, int_wu, ext_wu)
subroutine suews_cal_qe(Diagnose, snowuse, tstep, imin, it, EvapMethod, snowCalcSwitch, dayofWeek_id, CRWmin, CRWmax, dectime, lvS_J_kg, avRh, Press_hPa, Temp_C, RAsnow, psyc_hPa, sIce_hPa, PervFraction, vegfraction, addimpervious, qn1_snowfree, qf, qs, vpd_hPa, s_hPa, ResistSurf, RA, rb, snowdensmin, precip, PipeCapacity, RunoffToWater, NonWaterFraction, wu_EveTr, wu_DecTr, wu_Grass, addVeg, addWaterBody, SnowLimPaved, SnowLimBldg, SurfaceArea, FlowChange, drain, WetThresh, stateOld, mw_ind, SoilStoreCap, rainonsnow, freezmelt, freezstate, freezstatevol, Qm_Melt, Qm_rain, Tsurf_ind, sfr, StateLimit, AddWater, addwaterrunoff, StoreDrainPrm, SnowPackLimit, SnowProf_24hr, SnowPack_in, SnowFrac_in, SnowWater_in, iceFrac_in, SnowDens_in, runoff_per_interval_in, state_id_in, soilstore_id_in, state_id_out, soilstore_id_out, SnowPack_out, SnowFrac_out, SnowWater_out, iceFrac_out, SnowDens_out, runoffSoil, SnowRemoval, state_per_tstep, NWstate_per_tstep, qe, swe, chSnow_per_interval, ev_per_tstep, runoff_per_tstep, surf_chang_per_tstep, runoffPipes, mwstore, runoffwaterbody, runoffAGveg, runoffAGimpervious)
integer, parameter ncolumnsdataoutsuews
subroutine suews_cal_qs(StorageHeatMethod, qs_obs, OHMIncQF, Gridiv, id, tstep, dt_since_start, Diagnose, sfr, OHM_coef, OHM_threshSW, OHM_threshWD, soilstore_id, SoilStoreCap, state_id, SnowUse, SnowFrac, DiagQS, HDD_id, MetForcingData_grid, Ts5mindata_ir, qf, qn1, avkdn, avu1, temp_c, zenith_deg, avrh, press_hpa, ldown, bldgh, alb, emis, cpAnOHM, kkAnOHM, chAnOHM, EmissionsMethod, Tair_av, qn1_av_prev, dqndt_prev, qn1_s_av_prev, dqnsdt_prev, StoreDrainPrm, qn1_S, dataOutLineESTM, qs, qn1_av_next, dqndt_next, qn1_s_av_next, dqnsdt_next, deltaQi, a1, a2, a3)
real(kind(1d0)) function cal_tsfc(qh, avdens, avcp, RA, temp_c)
subroutine suews_cal_roughnessparameters(RoughLenMomMethod, sfr, bldgH, EveTreeH, DecTreeH, porosity_id, FAIBldg, FAIEveTree, FAIDecTree, z0m_in, zdm_in, Z, planF, Zh, z0m, zdm, ZZD)
real(kind(1d0)) function stab_psi_mom(StabilityMethod, ZL)
subroutine aerodynamicresistance(ZZD, z0m, AVU1, L_mod, UStar, VegFraction, AerodynamicResistanceMethod, StabilityMethod, RoughLenHeatMethod, RA)
integer, parameter decidsurf
integer, parameter pavsurf
real(kind(1d0)) function cal_tair_av(tair_av_prev, dt_since_start, tstep, temp_c)
subroutine errorhint(errh, ProblemFile, VALUE, value2, valueI)
subroutine suews_cal_surf(sfr, vegfraction, ImpervFraction, PervFraction, NonWaterFraction)
subroutine suews_cal_dls(id, startDLS, endDLS, DLS)
integer, parameter bldgsurf
integer, parameter watersurf
real(kind(1d0)) function qsatf(T, PMB)
subroutine drainage(is, state_is, StorCap, DrainEq, DrainCoef1, DrainCoef2, nsh_real, drain_is)
real function square_real(x)
integer, parameter ivdecid
subroutine anthropogenicemissions(CO2PointSource, EmissionsMethod, id, it, imin, DLS, nsh, DayofWeek_id, EF_umolCO2perJ, FcEF_v_kgkm, EnEF_v_Jkm, TrafficUnits, FrFossilFuel_Heat, FrFossilFuel_NonHeat, MinFCMetab, MaxFCMetab, MinQFMetab, MaxQFMetab, PopDensDaytime, PopDensNighttime, Temp_C, HDD_id, Qf_A, Qf_B, Qf_C, AH_MIN, AH_SLOPE_Heating, AH_SLOPE_Cooling, T_CRITIC_Heating, T_CRITIC_Cooling, TrafficRate, QF0_BEU, QF_SAHP, Fc_anthro, Fc_metab, Fc_traff, Fc_build, Fc_point, AHProf_24hr, HumActivity_24hr, TraffProf_24hr, PopProf_24hr, SurfaceArea)