Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(abstract_eqfile_geo_type), | intent(in) | :: | self | |||
real, | intent(in) | :: | r | |||
real, | intent(in) | :: | theta |
real function psi_eqfile (self, r, theta) result(psi) !Can't mark pure as eqfile_cart variant is not pure
implicit none
class(abstract_eqfile_geo_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
end if
end function psi_eqfile