erosion_organic_nitrate_loss Subroutine

  1. erosion.f95
  2. erosion
  3. erosion_organic_nitrate_loss

public subroutine erosion_organic_nitrate_loss(j, da9, yd, yon)

! if (yon.lt.xnorgp) yon = 0.

Arguments

Type IntentOptional AttributesName
integer :: j
real(kind=dp), intent(in) :: da9
real(kind=dp), intent(in) :: yd
real(kind=dp), intent(inout) :: yon

Called by

proc~~erosion_organic_nitrate_loss~~CalledByGraph proc~erosion_organic_nitrate_loss erosion_organic_nitrate_loss proc~runsubbasin runsubbasin proc~runsubbasin->proc~erosion_organic_nitrate_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_organic_nitrate_loss

Source Code

  subroutine erosion_organic_nitrate_loss(j, da9, yd, yon)
    !**** PURPOSE: COMPUTES organic N loss with erosion
    !**** CALLED IN: SUBBASIN
    !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    !     PARAMETERS & VARIABLES
    !
    !      >>>>> COMMON PARAMETERS & VARIABLES
    !      conn = xnorg * er, g/t
    !      da9 = 100. * da = basin area in ha, from readbas
    !      er = enrichment ratio, from enrsb
    !      xnorg = N org. in I layer for subbasin, g/t
    !      xnorgp = N org. in I layer for subbasin, kg/ha
    !      yd = daily soil loss (erosion), in t, calc in ysed
    !      yon = org N loss with erosion, kg/ha
    !      yone(j) = org N loss with erosion, kg/ha
    !      >>>>>
    !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    !**** Include common parameters

    real(dp), intent(in) :: da9
    real(dp), intent(in) :: yd
    real(dp), intent(inout) :: yon
    integer j

    !**** CALC org N loss with erosion
    conn = xnorg * er
    yon = .001 * conn * yd / da9

    !**** Correction: AnjaH
    !!!      if (yon.lt.xnorgp) yon = 0.
    if (yon .gt. xnorgp) yon = xnorgp
    if (yon .le. 0) yon = 0.

    yone(j) = yon
    return
  end subroutine erosion_organic_nitrate_loss