hru_control Subroutine

subroutine hru_control()

~ ~ ~ PURPOSE ~ ~ ~
this subroutine controls the simulation of the land phase of the 
hydrologic cycle

adjust precip and temperature for elevation using lapse rates zero pesticide balance variables zero all variables except pest in soil and in/on plant ht1== deposition: write to deposition.out ht2== outflow from inflow: added to hru generated flows check auto operations check day of future fertilizer application if end of year, reset the one time operation per year increment days since last plant and harvest update base zero total heat units zero stresses calculate albedo for day calculate soil temperature for soil layers compute canopy interception compute snow melt route overland flow across hru - add tile flow if not subirrigation or saturated buffer route incoming surface runoff add surface runon to wetland route across hru - infiltrate and deposit sediment add lateral flow soil water Route incoming lateral soil flow add tile flow to tile (subirrigation and saturated buffer) add aquifer flow to bottom soil layer and redistribute upwards Route incoming aquifer flow compute crack volume compute evapotranspiration perform management operations compute surface runoff processes if wetland - no runoff or sediment yield - all constituents transported in surface runoff and sediment will be zero ht2%sed==sediment routed across hru from surface runon check wetland/paddy continuous irrigation Jaehak 2023 manual irrigation on if ponding depth is lower than the threshold depth wetland/paddy processes compute effective rainfall (amount that percs into soil) for wetland use seepage into soil from ponded water no wetland (no ponded water) add irrigation to subdaily effective precip perform soil water routing compute peak rate similar to swat-deg using SCS triangular unit hydrograph

graze only if adequate biomass in HRU if total above ground biomass is available - graze check to set if grazing period is over compute nitrogen and phosphorus mineralization compute biozone processes in septic HRUs if 1) current is septic hru and 2) soil temperature is above zero compute plant community partitions check irrigation demand decision table for water allocation (after adding irrigation) iauto points to pcom(j)%dtbl(iauto) for days between operation check monsoon season for tropical plants compute plant biomass, leaf, root and seed growth reset harvested biomass and number of harvests for yearly yield output compute total parms for all plants in the community compute aoil water content to 300 mm depth compute total surface residue compute actual ET for day in HRU compute pesticide washoff
compute pesticide uptake compute pesticide degradation compute pesticide movement in soil sum total pesticide in soil Add by zhang ==================== Add by zhang ====================

add nitrate in rainfall to soil profile compute nitrate movement leaching compute phosphorus movement compute pathogen transport compute loadings from urban areas compute sediment loading in lateral flow and add to sedyld lag nutrients and sediment in surface runoff lag subsurface flow and nitrate in subsurface flow compute reduction in pollutants due to edge-of-field filter strip compute reduction in pollutants due to in field grass waterway compute reduction in pollutants due to in fixed BMP eff ht2%flo is outflow from wetland or total saturation excess if no wetland mm = m3/ha ha/10000m2 1000mm/m calculate amount of surface runoff during day (qday) and store the remainder compute water yield for HRU compute chl-a, CBOD and dissolved oxygen loadings add evap from impounded water (wetland) to et and esoil set hydrographs for direct routing or landscape unit

Arguments