1 subroutine r8mkspl2z(fun,x,nx,th,nth,fspl,nf2,wk,inwk,
12 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
13 INTEGER nth,nf2,inwk,nx,ix,ith
24 real*8 fspl(4,nf2,nth)
36 write(6,
'('.lt.
' ?mkspl2p -- array dim error, nf2nx'')')
42 fspl(1,ix,ith)=fun(x(ix),th(ith))
46 call r8mkbicubw(x,nx,th,nth,fspl,nf2,
47 > 0,zdummy,0,zdummy,-1,zdummy,-1,zdummy,
48 > wk,inwk,ilinx,ilinth,ier)