psi_eqfile Function

public function psi_eqfile(self, r, theta) result(psi)

Type Bound

abstract_eqfile_geo_type

Arguments

Type IntentOptional Attributes Name
class(abstract_eqfile_geo_type), intent(in) :: self
real, intent(in) :: r
real, intent(in) :: theta

Return Value real


Contents

Source Code


Source Code

  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