~ ~ ~ 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