1 subroutine r8gridintrp3d(
2 > x_newgrid,nx_new,y_newgrid,ny_new,z_newgrid,nz_new,
4 > nx,xpkg,ny,ypkg,nz,zpkg,
5 > jspline,fspl,icoeff,ixdim,iydim,izdim,
21 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
22 INTEGER ny_new,nz_new,nx_new,iz,iy
24 REAL*8 x_newgrid(nx_new)
25 real*8 y_newgrid(ny_new)
26 real*8 z_newgrid(nz_new)
31 real*8 f_new(if1,if2,nz_new)
46 integer :: ixdim,iydim,izdim
50 real*8 fspl(icoeff,ixdim,iydim,izdim)
65 data ict/1,0,0,0,0,0,0,0,0,0/
73 call r8vecintrp3d(ict,nx_new,x_newgrid,ytmp,ztmp,
74 > nx_new,f_new(1,iy,iz),
75 > nx,xpkg,ny,ypkg,nz,zpkg,
76 > jspline,fspl,icoeff,ixdim,iydim,izdim,