1 subroutine r8mkspl3zb(fun,x,nx,th,nth,ph,nph,fspl,nf2,nf3,
2 > ibcxmin,bcxmin,ibcxmax,bcxmax,nb1,
3 > wk,inwk,ilinx,ilinth,ilinph,ier)
24 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
25 INTEGER nth,nph,nf2,nf3,nb1,inwk,nx,iph,ith,ix
35 real*8 fspl(8,nf2,nf3,nph)
39 REAL*8 bcxmin(nb1,nph)
41 real*8 bcxmax(nb1,nph)
53 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nf2 nx'')')
57 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nf3 nth'')')
61 write(6,
'('.lt.
' ?mkspl3pb -- array dim error, nb1 nth'')')
70 fspl(1,ix,ith,iph)=fun(x(ix),th(ith),ph(iph))
75 call r8mktricubw(x,nx,th,nth,ph,nph,fspl,nf2,nf3,
76 > ibcxmin,bcxmin,ibcxmax,bcxmax,nb1,
77 > -1,zdummy,-1,zdummy,max(nx,nth,nph),
78 > -1,zdummy,-1,zdummy,max(nx,nth,nph),
79 > wk,inwk,ilinx,ilinth,ilinph,ier)