vegetation_s_curve Subroutine

  1. vegetation.f95
  2. vegetation
  3. vegetation_s_curve

public subroutine vegetation_s_curve(x1, x2, x3, x4)

Arguments

Type IntentOptional AttributesName
real(kind=dp) :: x1
real(kind=dp) :: x2
real(kind=dp) :: x3
real(kind=dp) :: x4

Called by

proc~~vegetation_s_curve~~CalledByGraph proc~vegetation_s_curve vegetation_s_curve proc~crop_growth crop_growth proc~crop_growth->proc~vegetation_s_curve proc~crop_process crop_process proc~crop_process->proc~crop_growth proc~hydrotope_process hydrotope_process proc~hydrotope_process->proc~crop_process proc~runsubbasin runsubbasin proc~runsubbasin->proc~hydrotope_process 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

vegetation_s_curve

Source Code

  subroutine vegetation_s_curve(x1, x2, x3, x4)
    !**** PURPOSE: THIS SUBROUTINE COMPUTES S CURVE PARAMETERS x1 & x2
    !              GIVEN 2 (X, Y) POINTS
    !**** CALLED IN: CRPMD
    !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    real(dp) x1, x2, x3, x4, xx

    xx =log(x3 / x1 - x3)
    x2 = (xx -log(x4 / x2 - x4)) / (x4 - x3)
    x1 = xx + x3 * x2
    return
  end subroutine vegetation_s_curve