sq_greenampt Subroutine

subroutine sq_greenampt()

Uses

  • proc~~sq_greenampt~~UsesGraph proc~sq_greenampt sq_greenampt module~basin_module basin_module proc~sq_greenampt->module~basin_module module~climate_module climate_module proc~sq_greenampt->module~climate_module module~hru_module hru_module proc~sq_greenampt->module~hru_module module~hydrograph_module hydrograph_module proc~sq_greenampt->module~hydrograph_module module~soil_module soil_module proc~sq_greenampt->module~soil_module module~time_module time_module proc~sq_greenampt->module~time_module module~urban_data_module urban_data_module proc~sq_greenampt->module~urban_data_module module~hydrograph_module->module~basin_module module~hydrograph_module->module~time_module module~carbon_module carbon_module module~soil_module->module~carbon_module 
~ ~ ~ PURPOSE ~ ~ ~
Predicts daily runoff given breakpoint precipitation and snow melt
using the Green & Ampt technique

~ ~ ~ INCOMING VARIABLES ~ ~ ~
name        |units         |definition
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
wst(:)%weat%ts(:)  |mm H2O        |precipitation for the time step during day
swtrg(:)    |none          |rainfall event flag:
                           |  0: no rainfall event over midnight
                           |  1: rainfall event over midnight
wfsh(:)     |mm            |average capillary suction at wetting front
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~

~ ~ ~ OUTGOING VARIABLES ~ ~ ~
name        |units         |definition
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
hhqday(:)   |mm H2O        |surface runoff generated each hour of day
                           |in HRU
surfq(:)    |mm H2O        |surface runoff for the day in HRU
swtrg(:)    |none          |rainfall event flag:
                           |  0: no rainfall event over midnight
                           |  1: rainfall event over midnight
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~

~ ~ ~ SUBROUTINES/FUNCTIONS CALLED ~ ~ ~
Intrinsic: Sum, Exp, Real, Mod

~ ~ ~ ~ ~ ~ END SPECIFICATIONS ~ ~ ~ ~ ~ ~

reset values for day Update effective hydraulic conductivity calculate total amount of rainfall during day for time step and rainfall intensity for time step if rainfall intensity is less than infiltration rate everything will infiltrate if rainfall intensity is greater than infiltration rate find cumulative infiltration for time step by successive substitution Urban Impervious cover effective pcp is > 0.017 mm/min (or 4/100 inches/hr) the potential for initial dabstraction from paved surface is less than the user input initial dabstraction the potential for initial dabstraction from paved surface increases based on evaporation nbs

nbs daily total runoff calculate new rate of infiltration ** tu Wunused-label: 5000 format(//,"Excess rainfall calculation for day ",i3," of year ", &
i4," for sub-basin",i4,".",/) ** tu Wunused-label: 5001 format(t2,"Time",t9,"Incremental",t22,"Cumulative",t35,"Rainfall", & t45,"Infiltration",t59,"Cumulative",t71,"Cumulative",t82, &
"Incremental",/,t2,"Step",t10,"Rainfall",t23,"Rainfall", &
t35,"Intensity",t49,"Rate",t58,"Infiltration",t73,"Runoff", & t84,"Runoff",/,t12,"(mm)",t25,"(mm)",t36,"(mm/h)",t48, &
"(mm/h)",t62,"(mm)",t74,"(mm)",t85,"(mm)",/) ** tu Wunused-label: 5002 format(i5,t12,f5.2,t24,f6.2,t36,f6.2,t47,f7.2,t61,f6.2,t73,f6.2, &
t84,f6.2)

Arguments

None