river_route_phosphorus Subroutine

  1. river.f95
  2. river
  3. river_route_phosphorus

public subroutine river_route_phosphorus(cpp, da9, yph)

Arguments

Type IntentOptional AttributesName
real(kind=dp), intent(in) :: cpp
real(kind=dp), intent(in) :: da9
real(kind=dp), intent(out) :: yph

Called by

proc~~river_route_phosphorus~~CalledByGraph proc~river_route_phosphorus river_route_phosphorus proc~river_route river_route proc~river_route->proc~river_route_phosphorus proc~time_process_day time_process_day proc~time_process_day->proc~river_route proc~time_process_month time_process_month proc~time_process_month->proc~time_process_day proc~time_process_years time_process_years proc~time_process_years->proc~time_process_month program~swim swim program~swim->proc~time_process_years

Contents

Source Code

river_route_phosphorus

Source Code

  subroutine river_route_phosphorus(cpp, da9, yph)
    !**** PURPOSE: THIS SUBROUTINE CALCULATES PHOSPHORUS ROUTING
    !**** CALLED IN: ROUTE
    !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    !      PARAMETERS & VARIABLES
    !
    !       >>>>> COMMON PARAMETERS & VARIABLES
    !       cpp = xporg * er, g/t: P org. in I layer in subbasin, g/t
    !                  corrected for enrichment
    !       da9 = drainage area, ha
    !       ydi = yd, daily soil loss caused by water erosion, t,
    !                 in route before call rtsed
    !       yph = P org. loss with erosion, kg/ha
    !       >>>>>
    !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    !**** Include common parameters

    real(dp), intent(in) :: cpp
    real(dp), intent(in) :: da9
    real(dp), intent(out) :: yph

    yph = .001 * cpp * ydi / da9
    return
  end subroutine river_route_phosphorus