subroutine pest_enrsb !! ~ ~ ~ PURPOSE ~ ~ ~ !! this subroutine calculates the enrichment ratio for nutrient and !! pesticide transport with runoff !! ~ ~ ~ INCOMING VARIABLES ~ ~ ~ !! name |units |definition !! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ !! surfq(:) |mm H2O |surface runoff generated on day in HRU !! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ !! ~ ~ ~ OUTGOING VARIABLES ~ ~ ~ !! name |units |definition !! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ !! enratio |none |enrichment ratio !! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ !! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ !! ~ ~ ~ ~ ~ ~ END SPECIFICATIONS ~ ~ ~ ~ ~ ~ use hru_module, only : hru, surfq, sedyld, sanyld, silyld, clayld, sagyld, lagyld, ihru, enratio implicit none integer :: j = 0 !none |HRU number real :: cy = 0. ! | j = ihru if (sedyld(j) < 1.e-4) then sedyld(j) = 0.0 sanyld(j) = 0.0 silyld(j) = 0.0 clayld(j) = 0.0 sagyld(j) = 0.0 lagyld(j) = 0.0 endif !! CREAMS method for calculating enrichment ratio cy = .1 * sedyld(j) / (hru(j)%area_ha * surfq(j) + 1.e-6) if (cy > 1.e-6) then enratio = .78 * cy ** (-.2468) else enratio = 0. endif if (enratio > 3.0) enratio = 3.0 return end subroutine pest_enrsb