toaxis.f90

subroutine i2mex_toAxis (nt1, ns, the, psi, power, f, ier)

Cubic extrapolation to axis for 1st, psi(1) node.

*NOTE* Need to have at least 5 nodes

integer, intent(in) :: nt1, ns

real(i2mex_r8), intent(in) :: the(nt1), psi(ns)

real(i2mex_r8), intent(in) :: power     f ~ psi**power as psi-> 0

real(i2mex_r8), intent(inout) :: f(nt1, ns)

integer, intent(out) :: ier

subroutine i2mex_toAxisFraction (nt1, ns, the, psi, power, fraction, f, ier)

Cubic extrapolation to axis for first fraction of points

fraction should be < 0.5

integer, intent(in) :: nt1, ns

real(i2mex_r8), intent(in) :: the(nt1), psi(ns)

real(i2mex_r8), intent(in) :: power     f ~ psi**power as psi-> 0

real(i2mex_r8), intent(in) :: fraction     fraction of grid where fct is to

be extrapolated

real(i2mex_r8), intent(inout) :: f(nt1, ns)

integer, intent(out) :: ier

subroutine i2mex_toEdgeFraction (nt1, ns, the, psi, fraction, f, ier)

Cubic extrapolation to edge for last fraction of points

fraction should be < 0.5

integer, intent(in) :: nt1, ns

real(i2mex_r8), intent(in) :: the(nt1), psi(ns)

real(i2mex_r8), intent(in) :: fraction     fraction of grid where fct is to

be extrapolated

real(i2mex_r8), intent(inout) :: f(nt1, ns)

integer, intent(out) :: ier