9 REAL(rprec) :: dpsic,xval(kpsi),chipsipsi(kpsi)
11 REAL(rprec) :: psival2(nxd),dum2(nxd)
13 REAL(rprec) :: dchi,psitest
16 dchi=(psilim-psiaxis)*percenflux
18 CALL initpsi(psic,psiv,xval,chipsi,chipsipsi,dpsic,npsi,
22 psiv(j)=psiv(j)+psiaxis
25 psiv(npsi)=psiaxis+dchi
26 dpsiv=(psilim-psiaxis)/(nx-1)
28 psival2(i)=psiaxis+(i-1.)*dpsiv
33 CALL spline1d(pprime,psiv,npsi,spp,psival2,nx,dum2)
34 CALL spline1d(ffprime,psiv,npsi,sffp,psival2,nx,dum2)
35 CALL spline1d(fval,psiv,npsi,sf,psival2,nx,dum2)
36 CALL spline1d(qsfin,psiv,npsi,qpsi,psival2,nx,dum2)
37 CALL spline1d(press,psiv,npsi,sp,psival2,nx,dum2)
38 CALL spline1d(pw,psiv,npsi,pressw,psival2,nx,dum2)
39 CALL spline1d(pwp,psiv,npsi,pwprim,psival2,nx,dum2)
40 CALL spline1d(rho,psiv,npsi,rho0,psival2,nx,dum2)
41 CALL spline1d(rhop,psiv,npsi,rho0p,psival2,nx,dum2)
44 CALL bcspline(xgrid,zgrid,psixz,nx,nz,nxd,csplpsi)
56 CALL equalarc(psitest,i)
58 END SUBROUTINE mapperb