1 subroutine tcspvec(ict,ivec,xvec,yvec,zvec,ivd,fval,
2 > nx,xpkg,ny,ypkg,nz,zpkg,fspl,inf4,inf5,
51 real fspl(4,4,4,inf4,inf5,nz)
61 integer,
dimension(:),
allocatable :: ix,iy,iz
62 real,
dimension(:),
allocatable :: dxv,dyv,dzv
71 write(6,*) .lt.
' ?tcspvec: nx2: nx = ',nx
76 write(6,*) .lt.
' ?tcspvec: ny2: ny = ',ny
81 write(6,*) .lt.
' ?tcspvec: nz2: nz = ',nz
86 write(6,*) .le.
' ?tcspvec: vector dimension 0: ivec = ',
93 >
' ?tcspvec: output vector dimension less than input ',
95 write(6,*)
' ivec=',ivec,
' ivd=',ivd
101 allocate(ix(ivec), iy(ivec), iz(ivec),
102 > dxv(ivec), dyv(ivec), dzv(ivec), stat=ier)
106 >
' ?tcspvec: memory allocation failure.'
115 call xlookup(ivec,xvec,nx,xpkg,1,ix,dxv,dxv,dxv,iwarn1)
116 call xlookup(ivec,yvec,ny,ypkg,1,iy,dyv,dyv,dyv,iwarn2)
117 call xlookup(ivec,zvec,nz,zpkg,1,iz,dzv,dzv,dzv,iwarn3)
118 iwarn=max(iwarn1,iwarn2,iwarn3)
122 call tcspevfn(ict,ivec,ivd,fval,ix,iy,iz,dxv,dyv,dzv,
125 deallocate(ix,iy,iz,dxv,dyv,dzv)