FIXME : Add documentation
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(in) | :: | vel |
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