1 subroutine mkspl3zb(fun,x,nx,th,nth,ph,nph,fspl,nf2,nf3,
2 > ibcxmin,bcxmin,ibcxmax,bcxmax,nb1,
3 > wk,inwk,ilinx,ilinth,ilinph,ier)
26 real fspl(8,nf2,nf3,nph)
44 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nf2 nx'')')
48 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nf3 nth'')')
52 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nb1 nth'')')
61 fspl(1,ix,ith,iph)=fun(x(ix),th(ith),ph(iph))
66 call mktricubw(x,nx,th,nth,ph,nph,fspl,nf2,nf3,
67 > ibcxmin,bcxmin,ibcxmax,bcxmax,nb1,
68 > -1,zdummy,-1,zdummy,max(nx,nth,nph),
69 > -1,zdummy,-1,zdummy,max(nx,nth,nph),
70 > wk,inwk,ilinx,ilinth,ilinph,ier)