Calculates the normalised radial coordinate
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(in) | :: | rp |
pure real function psifun (rp)
implicit none
real, intent (in) :: rp
psifun = min(1.0, max(0.0, (rp - rpmin) / (rpmax - rpmin)))
end function psifun