1 subroutine mkspl3pb(fun,x,nx,th,nth,ph,nph,fspl,nf4,nf5,
2 > ibcxmin,bcxmin,ibcxmax,bcxmax,nb1,
3 > wk,inwk,ilinx,ilinth,ilinph,ier)
26 real fspl(4,4,4,nf4,nf5,nph)
44 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nf4 nx'')')
48 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nf5 nth'')')
52 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nb1 nth'')')
55 if(inwk.lt.5*max(nx,nth,nph))
then
56 write(6,
'('' ?mkspl3pb -- array dim error, inwk too small'')')
65 fspl(1,1,1,ix,ith,iph)=fun(x(ix),th(ith),ph(iph))
70 call tpsplinb(x,nx,th,nth,ph,nph,fspl,nf4,nf5,
71 > ibcxmin,bcxmin,ibcxmax,bcxmax,nb1,
72 > wk,inwk,ilinx,ilinth,ilinph,ier)