1 subroutine r8mkspl2zb(fun,x,nx,th,nth,fspl,nf2,
2 > ibcxmin,bcxmin,ibcxmax,bcxmax,
3 > wk,inwk,ilinx,ilinth,ier)
21 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
22 INTEGER nth,nf2,inwk,nx,ix,ith
31 real*8 fspl(4,nf2,nth)
48 write(6,
'('.lt.
' ?mkspl2pb -- array dim error, nf2nx'')')
54 fspl(1,ix,ith)=fun(x(ix),th(ith))
58 call r8mkbicubw(x,nx,th,nth,fspl,nf2,
59 > ibcxmin,bcxmin,ibcxmax,bcxmax,
60 > -1,zdummy,-1,zdummy,
61 > wk,inwk,ilinx,ilinth,ier)