~ ~ ~ PURPOSE ~ ~ ~ this subroutine stores and lags sediment and nutrients in surface runoff
~ ~ ~ INCOMING VARIABLES ~ ~ ~ name |units |definition ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ brt(:) |none |fraction of surface runoff that takes |one day or less to reach the subbasin |outlet
ihru |none |HRU number sedminpa(:) |kg P/ha |amount of active mineral phosphorus sorbed to |sediment in surface runoff in HRU for day sedminps(:) |kg P/ha |amount of stable mineral phosphorus sorbed to |sediment in surface runoff in HRU for day surf_bs(2,:) |metric tons |amount of sediment yield lagged over one |day surf_bs(3,:) |kg N/ha |amount of organic nitrogen loading lagged |over one day surf_bs(4,:) |kg P/ha |amount of organic phosphorus loading lagged |over one day surf_bs(5,:) |kg N/ha |amount of nitrate loading in surface runoff |lagged over one day surf_bs(6,:) |kg P/ha |amount of soluble phosphorus loading lagged |over one day surf_bs(7,:) |kg P/ha |amount of active mineral phosphorus loading |lagged over one day surf_bs(8,:) |kg P/ha |amount of stable mineral phosphorus loading |lagged over one day surf_bs(9,:) |# colonies/ha|amount of less persistent bacteria in |solution lagged over one day surf_bs(10,:) |# colonies/ha|amount of persistent bacteria in |solution lagged over one day surf_bs(11,:) |# colonies/ha|amount of less persistent bacteria |sorbed lagged over one day surf_bs(12,:) |# colonies/ha|amount of persistent bacteria |sorbed lagged over one day ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ OUTGOING VARIABLES ~ ~ ~ name |units |definition ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ sedminpa(:) |kg P/ha |amount of active mineral phosphorus sorbed to |sediment in surface runoff in HRU for day sedminps(:) |kg P/ha |amount of stable mineral phosphorus sorbed to |sediment in surface runoff in HRU for day surf_bs(2,:) |metric tons |amount of sediment yield lagged over one |day surf_bs(3,:) |kg N/ha |amount of organic nitrogen loading lagged |over one day surf_bs(4,:) |kg P/ha |amount of organic phosphorus loading lagged |over one day surf_bs(5,:) |kg N/ha |amount of nitrate loading in surface runoff |lagged over one day surf_bs(6,:) |kg P/ha |amount of soluble phosphorus loading lagged |over one day surf_bs(7,:) |kg P/ha |amount of active mineral phosphorus loading |lagged over one day surf_bs(8,:) |kg P/ha |amount of stable mineral phosphorus loading |lagged over one day surf_bs(9,:) |# colonies/ha|amount of less persistent bacteria in |solution lagged over one day surf_bs(10,:) |# colonies/ha|amount of persistent bacteria in |solution lagged over one day surf_bs(11,:) |# colonies/ha|amount of less persistent bacteria |sorbed lagged over one day surf_bs(12,:) |# colonies/ha|amount of persistent bacteria |sorbed lagged over one day ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ SUBROUTINES/FUNCTIONS CALLED ~ ~ ~ Intrinsic: Max
~ ~ ~ ~ ~ ~ END SPECIFICATIONS ~ ~ ~ ~ ~ ~
Left-over (previous timestep) + inflow (current timestep)
new estimation of sediment reaching the main channel
lagged at the end of time step
daily total sediment yield from the HRU
sedyld(j) = surf_bs(2,j) * brt(j) <--line of code in x 2. fixes sedyld low prob