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
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
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