SUEWS API Site
Documentation of SUEWS source code
suews_ctrl_calculations.f95
Go to the documentation of this file.
1!This subroutine does the actual calculations of the SUEWS code (mostly old SUEWS_Temporal).
2!Made by LJ and HW Oct 2014
3!Gives in the grid ID (Gridiv) and number of line in the met forcing data to be analyzed (ir)
4!Last modification
5!
6!Last modification:
7! TS 02 May 2018 - added explict interfaces
8! LJ 15 Jun 2017 - Add parts where NonWaterFraction=0 is used as a divider so that the code won't crash with 100% water
9! TS 20 May 2017 - Add surface-level diagonostics: T2, Q2, U10
10! HCW 09 Dec 2016 - Add zenith and azimuth to output file
11! HCw 24 Aug 2016 - smd and state for each surface set to NAN in output file if surface does not exist.
12! - Added Fc to output file
13! HCW 21 Jul 2016 - Set soil variables to -999 in output when grid is 100% water surface.
14! HCW 29 Jun 2016 - Commented out StateDay and SoilMoistDay as creates jumps and should not be needed.
15! Would not work unless each met block consists of a whole day for each grid.
16! TS 09 Mar 2016 - Added AnOHM subroutine to calculate heat storage
17! HCW 10 Mar 2016 - Calculation of soil moisture deficit of vegetated surfaces added (vsmd)
18! LJ 2 Jan 2016 - Calculation of snow fraction moved from SUEWS_Calculations to SUEWS_Snow
19! HCW 12 Nov 2015 - Added z0m and zdm to output file
20! HCW 05 Nov 2015 - Changed Kawai et al. (2007) z0v calculation so VegFraction(veg+soil) rather than veg_fr(veg+soil+water) is used.
21! HCW 29 Jun 2015 - Added albEveTr and albGrass
22! HCW 25 Jun 2015 - Fixed bug in LAI calculation at year change.
23! - Changed AlbDec to use id value (not (id-1) value) to avoid issues at year change.
24! HCW 27 Apr 2015 - Correction to tot_chang_per_tstep calculation (water balance should now close)
25! HCW 16 Feb 2015 - Updated water balance calculations
26! - Corrected area-averaged calculations (soil moisture, drain, two versions of state with/out water)
27! - Replaced soilmoist_state variable with SoilState (as seems to be duplicate)
28! HCW 15 Jan 2015 - Added switch OHMIncQF to calculate QS with (1, default) or without (0) QF added to QSTAR
29! To do
30! - add iy and imin to output files (may impact LUMPS_RunoffFromGrid)
31! - phase out _per_interval water balance variables
32! - renormalise by NonWaterFraction where necessary
33! - Update Snow subroutines similarly in terms of water balance
34!==================================================================
35
36SUBROUTINE suews_calculations(Gridiv, ir, iMB, irMax)
50 USE time, ONLY: iy, id, it, imin, isec, dectime, dt_since_start
51 USE allocatearray, ONLY: &
52 alb, &
58 tmin_id, &
59 tmax_id, &
60 lenday_id, &
61 sdd_id, &
62 hdd_id, &
71 tair_av, &
72 dqndt, qn_av, &
73 dqnsdt, qn_s_av, &
76 soildepth, &
81 waterdist, &
82 wuday_id, &
96 sfr_surf, &
100 k_roof, k_wall, k_surf, &
103 height, building_frac, veg_frac, building_scale, veg_scale, & !input: SPARTACUS
106 USE sues_data, ONLY: &
108 h_maintain, &
116 USE snowmod, ONLY: &
122 USE gis_data, ONLY: &
125 USE mod_z, ONLY: z, z0m_in, zdm_in
127 USE resist, ONLY: g1, g2, g3, g4, g5, g6, gsmodel, kmax, s1, s2, th, tl
130 USE bluews_module, ONLY: cbl
131 USE moist, ONLY: avcp, avdens, es_hpa, lv_j_kg
132
133 IMPLICIT NONE
134
135 INTEGER :: Gridiv
136 INTEGER :: ir
137 INTEGER :: iMB
138 INTEGER :: irMax
139
140 !==================================================================
141
142 !Translate all data to the variables used in the model calculations
143 IF (diagnose == 1) WRITE (*, *) 'Calling SUEWS_Translate...'
144 CALL suews_translate(gridiv, ir, imb)
145
146 ! ASSOCIATE (v => dz_roof(1, 1:5))
147 ! PRINT *, 'dz_roof before driver', v, SIZE(v)
148 ! END ASSOCIATE
149 !PRINT *,''
150 !PRINT *,'Calling SUEWS_cal_Main'
151 IF (diagnose == 1) print *, 'Calling SUEWS_cal_Main...'
152 CALL suews_cal_main( &
153 aerodynamicresistancemethod, ah_min, ahprof_24hr, ah_slope_cooling, & ! input&inout in alphabetical order
158 basetmethod, &
172 isec, it, evapmethod, &
177 nlayer, &
184 height, building_frac, veg_frac, building_scale, veg_scale, & !input: SPARTACUS
191 qf0_beu, qf_a, qf_b, qf_c, &
203 tmax_id, tmin_id, &
206 sfr_surf, &
210 k_roof, k_wall, k_surf, &
218 dataoutlinebeers, & !output
221 dailystateline) !output
222
223 !============ update and write out SUEWS_cal_DailyState ===============
224 ! only works at the last timestep of a day
226 id, datetimeline, & !input
227 gridiv, numberofgrids, &
229 dataoutdailystate) !inout
230
231 !============ write out results ===============
232 ! works at each timestep
233 ! print *, 'writing out results',dataOutLineESTMExt
234 ! PRINT *, '=====writing out results====='
235 CALL suews_update_output( &
236 snowuse, storageheatmethod, & !input
242
243 ! NB: CBL disabled for the moment for interface improvement
244 ! NB: CBL be decoupled from SUEWS TS 10 Jun 2018
245
246 IF (qh_choice == 1) THEN !use QH and QE from SUEWS
247 qhforcbl(gridiv) = dataoutlinesuews(9)
248 qeforcbl(gridiv) = dataoutlinesuews(10)
249 ELSEIF (qh_choice == 2) THEN !use QH and QE from LUMPS
250 qhforcbl(gridiv) = dataoutlinesuews(11)
251 qeforcbl(gridiv) = dataoutlinesuews(12)
252 ELSEIF (qh_choice == 3) THEN !use QH and QE from OBS
253 qhforcbl(gridiv) = qh_obs
254 qeforcbl(gridiv) = qe_obs
255 IF (qh_obs < -900 .OR. qe_obs < -900) THEN ! observed data has a problem
256
257 CALL errorhint(22, 'Unrealistic observed qh or qe_value for CBL.', qh_obs, qe_obs, qh_choice)
258
259 END IF
260 END IF
261 IF (cbluse >= 1) THEN ! If CBL is used, calculated Temp_C and RH are replaced with the obs.
262 IF (diagnose == 1) WRITE (*, *) 'Calling CBL...'
263
265 !ir=1 indicates first row of each met data block
266 CALL cbl(iy, id, it, imin, ir, gridiv, qh_choice, dectime, &
270 END IF
271
272 ! NB: SOLWEIG can be treated as a separate part:
273 ! NB: SOLWEIG is disabled for v2018a TS 10 Jun 2018
274 ! ===================SOLWEIG OUTPUT ========================================
275 ! IF (SOLWEIGuse==1) THEN
276 ! IF (OutInterval==imin) THEN
277 ! IF (RunForGrid==-999) THEN
278 ! IF(Diagnose==1) WRITE(*,*) 'Calling SOLWEIG_2014a_core...'
279 ! CALL SOLWEIG_2014a_core(iMB)
280 ! SolweigCount=SolweigCount+1
281 ! ELSE
282 ! IF (Gridiv == RunForGrid) THEN
283 ! IF(Diagnose==1) WRITE(*,*) 'Calling SOLWEIG_2014a_core...'
284 ! CALL SOLWEIG_2014a_core(iMB)
285 ! SolweigCount=SolweigCount+1
286 ! ENDIF
287 ! ENDIF
288 ! ENDIF
289 ! ELSE
290 ! SOLWEIGpoi_out=0 ! NB: turn off SOLWEIG for the moment
291 ! ENDIF
292 ! ===================SOLWEIG END================================
293
294 IF (diagnose == 1) WRITE (*, *) 'Calling SUEWS_TranslateBack...'
295 CALL suews_translateback(gridiv, ir, irmax)
296
297END SUBROUTINE suews_calculations
real(kind(1d0)), dimension(:), allocatable statelimit_wall
real(kind(1d0)) pormax_dec
real(kind(1d0)), dimension(nsurf) snowpack
real(kind(1d0)), dimension(:, :, :), allocatable dataoutspartacus
real(kind(1d0)), dimension(:, :), allocatable dz_wall
real(kind(1d0)) tmin_id
real(kind(1d0)), dimension(nsurf) chanohm
real(kind(1d0)), dimension(nsurf) soilstorecap_surf
real(kind(1d0)), dimension(0:23, 2) humactivity_24hr
real(kind(1d0)) albmin_grass
real(kind(1d0)), dimension(ncolumnsdataoutdebug) dataoutlinedebug
real(kind(1d0)), dimension(nsurf) statelimit_surf
real(kind(1d0)), dimension(:, :), allocatable metforcingdata_grid
real(kind(1d0)) ground_albedo_dir_mult_fact
real(kind(1d0)), dimension(nsurf) state_surf
real(kind(1d0)), dimension(nsurf) soildepth
real(kind(1d0)), dimension(nsurf) cpanohm
real(kind(1d0)), dimension(:), allocatable soilstorecap_wall
real(kind(1d0)) porosity_id
real(kind(1d0)), dimension(nvegsurf) basete
real(kind(1d0)), dimension(:), allocatable tsfc_roof
real(kind(1d0)), dimension(:), allocatable alb_roof
real(kind(1d0)), dimension(:), allocatable state_wall
real(kind(1d0)), dimension(0:23, 2) traffprof_24hr
integer n_vegetation_region_urban
real(kind(1d0)), dimension(nsurf+1) ohm_threshwd
real(kind(1d0)), dimension(0:23, 2) ahprof_24hr
real(kind(1d0)), dimension(nvegsurf) beta_enh_bioco2
real(kind(1d0)), dimension(nsurf+1, 4, 3) ohm_coef
real(kind(1d0)), dimension(:, :), allocatable temp_roof
real(kind(1d0)), dimension(nvegsurf) min_res_bioco2
real(kind(1d0)), dimension(:), allocatable wetthresh_wall
real(kind(1d0)), dimension(nvegsurf) lai_id
real(kind(1d0)) veg_ssa_lw
real(kind(1d0)), dimension(:, :, :), allocatable dataoutestmext
real(kind(1d0)), dimension(:), allocatable state_roof
real(kind(1d0)) air_ext_sw
real(kind(1d0)), dimension(ncolumnsdataoutsnow - 5) dataoutlinesnow
real(kind(1d0)), dimension(nvegsurf) baset
real(kind(1d0)), dimension(:, :), allocatable dz_surf
real(kind(1d0)), dimension(nvegsurf) resp_b
real(kind(1d0)), dimension(nsurf) sfr_surf
real(kind(1d0)), dimension(4, nvegsurf) laipower
real(kind(1d0)), dimension(nvegsurf) resp_a
real(kind(1d0)), dimension(ncolumnsdataoutsuews - 5) dataoutlinesuews
real(kind(1d0)), dimension(nvegsurf) laimax
real(kind(1d0)), dimension(0:23, 2) wuprofm_24hr
real(kind(1d0)), dimension(:, :, :), allocatable dataoutsnow
real(kind(1d0)), dimension(nsurf) emis
real(kind(1d0)), dimension(nvegsurf) gdd_id
real(kind(1d0)), dimension(ncolumnsdataoutestm - 5) dataoutlineestm
real(kind(1d0)) capmin_dec
real(kind(1d0)) decidcap_id
real(kind(1d0)), dimension(:), allocatable building_scale
real(kind(1d0)) tair_av
real(kind(1d0)) air_ssa_sw
real(kind(1d0)), dimension(:), allocatable soilstore_roof
real(kind(1d0)), dimension(:, :, :), allocatable dataoutdebug
real(kind(1d0)), dimension(nsurf) soilstore_surf
integer, dimension(nvegsurf) laitype
logical use_sw_direct_albedo
real(kind(1d0)), dimension(:), allocatable height
real(kind(1d0)) veg_ssa_sw
real(kind(1d0)), dimension(:, :), allocatable cp_wall
real(kind(1d0)), dimension(nsurf) sathydraulicconduct
real(kind(1d0)) lenday_id
real(kind(1d0)) qn_s_av
real(kind(1d0)), dimension(6, nsurf) storedrainprm
real(kind(1d0)) air_ssa_lw
real(kind(1d0)), dimension(nsurf) kkanohm
real(kind(1d0)), dimension(nvegsurf) theta_bioco2
real(kind(1d0)), dimension(nsurf+1, nsurf - 1) waterdist
real(kind(1d0)), dimension(:), allocatable emis_roof
real(kind(1d0)), dimension(:), allocatable tsfc_surf
real(kind(1d0)) narp_emis_snow
real(kind(1d0)), dimension(:, :), allocatable cp_surf
real(kind(1d0)), dimension(nsurf+1) ohm_threshsw
real(kind(1d0)), dimension(:, :), allocatable k_roof
real(kind(1d0)), dimension(:), allocatable wetthresh_roof
real(kind(1d0)), dimension(:), allocatable soilstorecap_roof
real(kind(1d0)), dimension(:), allocatable veg_frac
real(kind(1d0)), dimension(:, :), allocatable k_surf
real(kind(1d0)) capmax_dec
real(kind(1d0)), dimension(nvegsurf) beta_bioco2
real(kind(1d0)), dimension(5) datetimeline
real(kind(1d0)), dimension(ncolumnsdataoutrsl - 5+12) dataoutlinersl
real(kind(1d0)), dimension(ncolumnsdataoutestmext - 5) dataoutlineestmext
real(kind(1d0)) albevetr_id
real(kind(1d0)), dimension(nsurf) snowwater
real(kind(1d0)) narp_trans_site
real(kind(1d0)), dimension(0:23, 2) wuprofa_24hr
real(kind(1d0)), dimension(ncolumnsdataoutspartacus) dataoutlinespartacus
real(kind(1d0)), dimension(:, :), allocatable roof_albedo_dir_mult_fact
real(kind(1d0)) albmax_dectr
real(kind(1d0)), dimension(nvegsurf) sdd_id
real(kind(1d0)), dimension(nvegsurf) sddfull
real(kind(1d0)), dimension(nvegsurf) laimin
real(kind(1d0)) albmin_dectr
real(kind(1d0)), dimension(:), allocatable tin_surf
real(kind(1d0)), dimension(:), allocatable tin_wall
real(kind(1d0)), dimension(nvegsurf) alpha_enh_bioco2
real(kind(1d0)), dimension(:), allocatable statelimit_roof
real(kind(1d0)), dimension(:), allocatable alb_wall
real(kind(1d0)), dimension(9) wuday_id
real(kind(1d0)) albdectr_id
real(kind(1d0)), dimension(:, :, :), allocatable dataoutrsl
real(kind(1d0)), dimension(:, :, :), allocatable dataoutsuews
real(kind(1d0)), dimension(:, :), allocatable wall_specular_frac
real(kind(1d0)) qn_av
real(kind(1d0)) veg_fsd_const
real(kind(1d0)), dimension(:, :), allocatable dz_roof
real(kind(1d0)) albmax_grass
real(kind(1d0)) pormin_dec
real(kind(1d0)), dimension(:, :), allocatable temp_surf
real(kind(1d0)), dimension(nsurf) snowdens
real(kind(1d0)), dimension(:), allocatable emis_wall
real(kind(1d0)), dimension(:), allocatable ts5mindata_ir
real(kind(1d0)) dqnsdt
real(kind(1d0)) tmax_id
real(kind(1d0)), dimension(ncolumnsdataoutbeers - 5) dataoutlinebeers
real(kind(1d0)) albmax_evetr
real(kind(1d0)), dimension(:, :, :), allocatable dataoutbl
real(kind(1d0)), dimension(nsurf) wetthresh_surf
real(kind(1d0)), dimension(:, :), allocatable cp_roof
real(kind(1d0)), dimension(nsurf) alb
real(kind(1d0)), dimension(:), allocatable building_frac
real(kind(1d0)), dimension(nvegsurf) maxconductance
real(kind(1d0)), dimension(0:23, 2) popprof_24hr
real(kind(1d0)), dimension(nvegsurf) gddfull
real(kind(1d0)), dimension(:, :), allocatable temp_wall
real(kind(1d0)), dimension(:, :, :), allocatable dataoutdailystate
real(kind(1d0)), dimension(:, :, :), allocatable dataoutestm
real(kind(1d0)) air_ext_lw
real(kind(1d0)), dimension(ncolumnsdataoutdailystate - 5) dailystateline
real(kind(1d0)), dimension(:, :), allocatable k_wall
real(kind(1d0)) dqndt
real(kind(1d0)), dimension(:, :, :), allocatable dataoutbeers
real(kind(1d0)), dimension(nsurf) snowpacklimit
real(kind(1d0)), dimension(:), allocatable veg_scale
real(kind(1d0)), dimension(12) hdd_id
real(kind(1d0)), dimension(:), allocatable tin_roof
real(kind(1d0)), dimension(:), allocatable soilstore_wall
real(kind(1d0)) sw_dn_direct_frac
real(kind(1d0)) albgrass_id
real(kind(1d0)), dimension(nsurf) icefrac
real(kind(1d0)) albmin_evetr
real(kind(1d0)), dimension(:), allocatable tsfc_wall
real(kind(1d0)), dimension(nvegsurf) alpha_bioco2
real(kind(1d0)), dimension(nsurf) snowfrac
real(kind(1d0)) veg_contact_fraction_const
subroutine cbl(iy, id, it, imin, ir, Gridiv, qh_choice, dectime, Temp_C, Press_hPa, avkdn, avu1, avrh, avcp, avdens, es_hPa, lv_J_kg, nsh_real, tstep, UStar, psih, is, NumberOfGrids, qhforCBL, qeforCBL, ReadLinesMetdata, dataOutBL)
subroutine suews_update_dailystate(id, datetimeline, Gridiv, NumberOfGrids, DailyStateLine, dataOutDailyState)
real(kind(1d0)) drainrt
integer netradiationmethod
integer evapmethod
real(kind(1d0)) timezone
real(kind(1d0)) minqfmetab
real(kind(1d0)) trafficunits
real(kind(1d0)) rainmaxres
real(kind(1d0)) frfossilfuel_heat
real(kind(1d0)) co2pointsource
integer ohmincqf
real(kind(1d0)) fcld_obs
integer smdmethod
real(kind(1d0)), dimension(2) fcef_v_kgkm
integer laicalcyes
real(kind(1d0)), dimension(2) qf0_beu
real(kind(1d0)) lng
real(kind(1d0)), dimension(2) qf_a
real(kind(1d0)) ef_umolco2perj
real(kind(1d0)) maxqfmetab
real(kind(1d0)) lai_obs
real(kind(1d0)) ldown_obs
real(kind(1d0)) qf_obs
real(kind(1d0)), dimension(2) ah_slope_heating
real(kind(1d0)) snowfrac_obs
real(kind(1d0)) baset_hc
integer emissionsmethod
real(kind(1d0)) precip
real(kind(1d0)), dimension(2) ah_min
real(kind(1d0)), dimension(2) baset_heating
real(kind(1d0)) lat
real(kind(1d0)) wu_m3
integer diagmethod
integer startdls
real(kind(1d0)), dimension(2) ah_slope_cooling
real(kind(1d0)) press_hpa
integer diagnose
integer basetmethod
real(kind(1d0)), dimension(2) qf_c
real(kind(1d0)) qn1_obs
real(kind(1d0)) qs_obs
real(kind(1d0)) qh_obs
integer waterusemethod
real(kind(1d0)) avrh
real(kind(1d0)) maxfcmetab
integer storageheatmethod
real(kind(1d0)), dimension(2) trafficrate
real(kind(1d0)) enef_v_jkm
real(kind(1d0)), dimension(2) qf_b
real(kind(1d0)) raincover
integer roughlenmommethod
real(kind(1d0)) avu1
real(kind(1d0)) temp_c
real(kind(1d0)), dimension(2) popdensdaytime
real(kind(1d0)) xsmd
real(kind(1d0)) minfcmetab
real(kind(1d0)), dimension(2) baset_cooling
real(kind(1d0)) alt
real(kind(1d0)) avkdn
real(kind(1d0)) popdensnighttime
real(kind(1d0)) qe_obs
real(kind(1d0)) frfossilfuel_nonheat
real(kind(1d0)) faibldg
real(kind(1d0)) bldgh
real(kind(1d0)) evetreeh
real(kind(1d0)) dectreeh
real(kind(1d0)) faievetree
real(kind(1d0)) faidectree
integer readlinesmetdata
integer numberofgrids
real(kind(1d0)) z
real(kind(1d0)) zdm_in
real(kind(1d0)) z0m_in
real(kind(1d0)) avcp
real(kind(1d0)) avdens
real(kind(1d0)) es_hpa
real(kind(1d0)) lv_j_kg
real(kind(1d0)) g1
real(kind(1d0)) th
real(kind(1d0)) tl
real(kind(1d0)) s2
real(kind(1d0)) g3
integer gsmodel
real(kind(1d0)) s1
real(kind(1d0)) g4
real(kind(1d0)) g6
real(kind(1d0)) g2
real(kind(1d0)) kmax
real(kind(1d0)) g5
real(kind(1d0)) snowfallcum
real(kind(1d0)), dimension(0:23, 2) snowprof_24hr
real(kind(1d0)) snowlimpaved
real(kind(1d0)) snowalbmax
real(kind(1d0)) radmeltfact
real(kind(1d0)) crwmax
real(kind(1d0)) tau_r
real(kind(1d0)) tempmeltfact
real(kind(1d0)) snowdensmin
real(kind(1d0)) snowlimbldg
real(kind(1d0)) snowalbmin
real(kind(1d0)) preciplimitalb
real(kind(1d0)) snowdensmax
real(kind(1d0)) snowalb
real(kind(1d0)) tau_a
real(kind(1d0)) tau_f
real(kind(1d0)) preciplimit
real(kind(1d0)) crwmin
real(kind(1d0)), dimension(:), allocatable qeforcbl
real(kind(1d0)) irrfracgrass
real(kind(1d0)) irrfracbldgs
real(kind(1d0)) faut
integer stabilitymethod
integer aerodynamicresistancemethod
real(kind(1d0)) ustar
real(kind(1d0)) runofftowater
real(kind(1d0)), dimension(7) daywat
real(kind(1d0)) flowchange
real(kind(1d0)) irrfracevetr
real(kind(1d0)), dimension(:), allocatable qhforcbl
real(kind(1d0)), dimension(3) ie_a
real(kind(1d0)) psih
real(kind(1d0)) irrfracdectr
real(kind(1d0)) irrfracbsoil
real(kind(1d0)) pipecapacity
integer roughlenheatmethod
real(kind(1d0)) irrfracwater
real(kind(1d0)), dimension(7) daywatper
real(kind(1d0)) nsh_real
real(kind(1d0)) irrfracpaved
real(kind(1d0)) surfacearea
real(kind(1d0)), dimension(3) ie_m
real(kind(1d0)) internalwateruse_h
real(kind(1d0)) h_maintain
subroutine suews_update_output(SnowUse, storageheatmethod, ReadLinesMetdata, NumberOfGrids, ir, gridiv, datetimeLine, dataOutLineSUEWS, dataOutLineSnow, dataOutLineESTM, dataoutLineRSL, dataOutLineBEERS, dataoutlineDebug, dataoutlineSPARTACUS, dataOutLineESTMExt, dataOutSUEWS, dataOutSnow, dataOutESTM, dataOutRSL, dataOutBEERS, dataOutDebug, dataOutSPARTACUS, dataOutESTMExt)
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, kdown, avRh, avU1, BaseT, BaseTe, BaseTMethod, BaseT_HC, beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, DecTreeH, DiagMethod, Diagnose, DiagQN, DiagQS, DRAINRT, dt_since_start, dqndt, qn_av, dqnsdt, qn_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, H_maintain, HDD_id, HumActivity_24hr, IceFrac, id, Ie_a, Ie_end, Ie_m, Ie_start, imin, InternalWaterUse_h, IrrFracPaved, IrrFracBldgs, IrrFracEveTr, IrrFracDecTr, IrrFracGrass, IrrFracBSoil, IrrFracWater, isec, it, 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, nlayer, n_vegetation_region_urban, n_stream_sw_urban, n_stream_lw_urban, sw_dn_direct_frac, air_ext_sw, air_ssa_sw, veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, veg_fsd_const, veg_contact_fraction_const, ground_albedo_dir_mult_fact, use_sw_direct_albedo, height, building_frac, veg_frac, building_scale, veg_scale, alb_roof, emis_roof, alb_wall, emis_wall, roof_albedo_dir_mult_fact, wall_specular_frac, OHM_coef, OHMIncQF, OHM_threshSW, OHM_threshWD, PipeCapacity, PopDensDaytime, PopDensNighttime, PopProf_24hr, PorMax_dec, PorMin_dec, Precip, PrecipLimit, PrecipLimitAlb, Press_hPa, QF0_BEU, Qf_A, Qf_B, Qf_C, qn1_obs, qs_obs, qf_obs, RadMeltFact, RAINCOVER, RainMaxRes, resp_a, resp_b, RoughLenHeatMethod, RoughLenMomMethod, RunoffToWater, S1, S2, SatHydraulicConduct, SDDFull, SDD_id, SMDMethod, SnowAlb, SnowAlbMax, SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, SnowLimBldg, SnowLimPaved, snowFrac_obs, SnowPack, SnowProf_24hr, SnowUse, SoilDepth, StabilityMethod, startDLS, soilstore_surf, SoilStoreCap_surf, state_surf, StateLimit_surf, WetThresh_surf, soilstore_roof, SoilStoreCap_roof, state_roof, StateLimit_roof, WetThresh_roof, soilstore_wall, SoilStoreCap_wall, state_wall, StateLimit_wall, WetThresh_wall, StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, Tmax_id, Tmin_id, BaseT_Cooling, BaseT_Heating, Temp_C, TempMeltFact, TH, theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, sfr_surf, tsfc_roof, tsfc_wall, tsfc_surf, temp_roof, temp_wall, temp_surf, tin_roof, tin_wall, tin_surf, k_roof, k_wall, k_surf, cp_roof, cp_wall, cp_surf, dz_roof, dz_wall, dz_surf, TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, WaterDist, WaterUseMethod, wu_m3, WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, WUProfA_24hr, WUProfM_24hr, xsmd, Z, z0m_in, zdm_in, datetimeLine, dataOutLineSUEWS, dataOutLineSnow, dataOutLineESTM, dataoutLineRSL, dataOutLineBEERS, dataOutLineDebug, dataOutLineSPARTACUS, dataOutLineESTMExt, DailyStateLine)
integer dt_since_start
integer isec
integer iy
integer it
integer imin
integer id
real(kind(1d0)) dectime
subroutine suews_calculations(Gridiv, ir, iMB, irMax)
subroutine errorhint(errh, ProblemFile, VALUE, value2, valueI)
subroutine suews_translateback(Gridiv, ir, irMax)
subroutine suews_translate(Gridiv, ir, iMB)