rlim Function

function rlim(d1, d2)

Uses

  • proc~~rlim~~UsesGraph proc~rlim rlim module~modglobal modglobal proc~rlim->module~modglobal

Arguments

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

Return Value real


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