minmod Function

public pure function minmod(a, b)

Standard minmod limiter function.

Arguments

Type IntentOptional Attributes Name
real(kind=rk), intent(in) :: a
real(kind=rk), intent(in) :: b

Return Value real(kind=rk)


Source Code

   pure real(rk) function minmod(a, b)
      real(rk), intent(in) :: a, b
      if (a * b <= zero) then
         minmod = zero
      else
         minmod = sign(one, a) * min(abs(a), abs(b))
      end if
   end function minmod