10 subroutine r8bpspline(x,inx,th,inth,fspl,inf3,
11 > wk,nwk,ilinx,ilinth,ier)
14 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
15 integer inx,inth,inf3,nwk
16 real*8 x(inx),th(inth),fspl(4,4,inf3,inth),wk(nwk)
60 if(nwk.lt.5*max(inx,inth))
then
61 write(6,
'('' ?bpspline: workspace too small.'')')
65 write(6,
'('' ?bpspline: at least 2 x points required.'')')
69 write(6,
'('' ?bpspline: need at least 2 theta points.'')')
75 call r8splinck(x,inx,ilinx,1.0e-3_r8,ierx)
79 write(6,
'('' ?bpspline: x axis not strict ascending'')')
84 call r8splinck(th,inth,ilinth,1.0e-3_r8,ierth)
88 write(6,
'('' ?bpspline: th axis not strict ascending'')')
99 call r8bcspline(x,inx,th,inth,fspl,inf3,
102 > wk,nwk,ilinx,ilinth,ier)