hsg_func Function

private elemental function hsg_func(vel)

Uses

FIXME : Add documentation

Arguments

Type IntentOptional Attributes Name
real, intent(in) :: vel

Return Value real


Contents

Source Code


Source Code

  elemental real function hsg_func (vel)
    use constants, only: sqrt_pi
    implicit none
    real, intent (in) :: vel
    if (abs(vel) <= epsilon(0.0)) then
       hsg_func = 0.0
    else
       hsg_func = 0.5 * erf(vel) / vel**2 - exp(-vel**2) / (sqrt_pi * vel)
    end if
  end function hsg_func