~ ~ ~ PURPOSE ~ ~ ~
this subroutine is the master soil percolation component.
~ ~ ~ INCOMING VARIABLES ~ ~ ~
name |units |definition
~ ~ ~ OUTGOING VARIABLES ~ ~ ~
name |units |definition
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
new water table depth equations 01/2009
c |none |a factor used to convert airvol to wtd
dg |mm |soil layer thickness in HRU
new water table depth equations 01/2009
latq(:) |mm H2O |total lateral flow in soil profile for the
|day in HRU
lyrtile |mm H2O |drainage tile flow in soil layer for day
new water table depth equations 01/2009
ne_p |mm/hr |effective porosity in HRU for all soil profile layers
ne_w |mm/hr |effective porosity in HRU for soil layers above wtd
new water table depth equations 01/2009
qtile |mm H2O |drainage tile flow in soil profile for the day
sepday |mm H2O |micropore percolation from soil layer
sepbtm(:) |mm H2O |percolation from bottom of soil profile for
|the day in HRU
sw_excess |mm H2O |amount of water in excess of field capacity
|stored in soil layer on the current day
new water table depth equations 01/2009
new water table depth equations 01/2009
wt_shall |mm H2O |shallow water table height above bottom of soil profile
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
w2 |mm |
y1 |mm |dummy variable for wat
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
initialize water entering first soil layer ht1%flo is infiltration from overland flow routing calculate crack flow add water moving into soil layer from overlying layer septic tank inflow to biozone layer J.Jeong
determine gravity drained water in layer initialize variables for current layer calculate tile flow (lyrtile), lateral flow (latlyr) and percolation (sepday) summary calculations redistribute soil water if above saturation (high water table) update soil profile water compute shallow water table depth and tile flow drainmod tile equations 08/11/2006 drainmod wt_shall equations 10/23/2006 drainmod wt_shall equations 10/23/2006
Start Daniel"s tile equations modifications 01/2006 drainmod tile equations compute tile flow using drainmod tile equations compute tile flow using existing tile equations End Daniel"s tile equations modifications 01/2006
update soil profile water after tile drainage update soil profile water