FIXME : Add documentation
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(leq_type), | intent(in) | :: | self | |||
real, | intent(in) | :: | r | |||
real, | intent(in) | :: | theta |
pure real function psi (self, r, theta)
class(leq_type), intent(in) :: self
real, intent (in) :: r, theta
if (r==0.) then
psi = self%psi_0
else if (r==1. .and. theta == 0.0) then
psi = self%psi_a
else
psi = r - self%surf%r
end if
end function psi