1 subroutine r8vecintrp2d(ict,ivec,xvec,yvec,ivd,fval,
2 > nx,xpkg,ny,ypkg,jspline,fspl,icoeff,ixdim,iydim,
13 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
71 real*8 fspl(icoeff,ixdim,iydim)
84 integer :: iwarn1,iwarn2
103 write(6,*) .lt.
' ?vecintrp2d: nx2: nx = ',nx
108 write(6,*) .lt.
' ?vecintrp2d: ny2: ny = ',ny
112 call vecin2d_argchk(
'vecintrp2d',jspline,
113 > icoeff,nx,ny,ixdim,iydim,ier)
116 write(6,*) .le.
' ?vecintrp2d: vector dimension 0: ivec = ',
123 >
' ?vecintrp2d: output vector dimension less than input ',
124 >
'vector dimension.'
125 write(6,*)
' ivec=',ivec,
' ivd=',ivd
135 call r8xlookup(ivec,xvec,nx,xpkg,2,ix,dxn,hx,hxi,iwarn1)
136 call r8xlookup(ivec,yvec,ny,ypkg,2,iy,dyn,hy,hyi,iwarn2)
141 call r8fvintrp2d(ict,ivec,ivd,fval,ix,iy,dxn,dyn,
142 > hx,hxi,hy,hyi,jspline,fspl,icoeff,ixdim,iydim)