erosion_soil_loss Subroutine

  1. erosion.f95
  2. erosion
  3. erosion_soil_loss

public subroutine erosion_soil_loss(aff, pr, xqd, yd)

Arguments

Type IntentOptional AttributesName
real(kind=dp), intent(in) :: aff
real(kind=dp), intent(in) :: pr
real(kind=dp), intent(in) :: xqd
real(kind=dp), intent(inout) :: yd

Called by

proc~~erosion_soil_loss~~CalledByGraph proc~erosion_soil_loss erosion_soil_loss proc~runsubbasin runsubbasin proc~runsubbasin->proc~erosion_soil_loss proc~time_process_day time_process_day proc~time_process_day->proc~runsubbasin 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

erosion_soil_loss

Source Code

  subroutine erosion_soil_loss(aff, pr, xqd, yd)
    !**** PURPOSE: THIS SUBROUTINE ESTIMATES DAILY SOIL LOSS CAUSED BY WATER
    !              EROSION USING THE MODIFIED UNIVERSAL SOIL LOSS EQUATION
    !**** CALLED IN: SUBBASIN
    !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    !     PARAMETERS & VARIABLES
    !
    !      >>>>> COMMON PARAMETERS & VARIABLES
    !      aff = 1000. * da * flu(j), calc in subbasin
    !      ida = current day
    !      ieros = switch code to print from eros()
    !      iersb = number of subbasin to print from eros(), if ieros = 1
    !      pr = peak runoff rate, m3/sec., calc in peakq
    !      precip = precipitation, mm
    !      sl(j) = USLE slope length/slope steepness factor
    !      xcklsp = combined c, k, ls, p factor for subbasin, calc in subbasin
    !      xqd = surface runoff for subbasin, mm, calc in subbasin
    !      yd = daily soil loss from subbasin caused by water erosion, t
    !      >>>>>
    !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    !**** Include common parameters

    real(dp), intent(in) :: aff
    real(dp), intent(in) :: pr
    real(dp), intent(in) :: xqd
    real(dp), intent(inout) :: yd

    yd = 11.8 * (xqd * pr * aff) ** .56 * xcklsp

    return
  end subroutine erosion_soil_loss