rlim Function

function rlim(d1, d2)

Uses

  • proc~~rlim~~UsesGraph proc~rlim advec_kappa.f90::rlim module~modglobal modglobal proc~rlim->module~modglobal

Arguments

Type IntentOptional Attributes Name
real, intent(in) :: d1
real, intent(in) :: d2

Return Value real


Called by

proc~~rlim~~CalledByGraph proc~rlim advec_kappa.f90::rlim proc~advecc_kappa advec_kappa.f90::advecc_kappa proc~advecc_kappa->proc~rlim proc~advection advection.f90::advection proc~advection->proc~advecc_kappa program~dalesurban DALESURBAN program~dalesurban->proc~advection

Contents

Source Code


Source Code

  real function rlim(d1, d2)
     use modglobal, only:eps1
     implicit none
     real, intent(in) :: d1 !< Scalar flux at 1.5 cells upwind
     real, intent(in) :: d2 !< Scalar flux at 0.5 cells upwind

     real ri, phir

     ri = (d2 + eps1)/(d1 + eps1)
     phir = max(0., min(2.*ri, min(1./3.+2./3.*ri, 2.)))
     rlim = 0.5*phir*d1
  end function rlim