1 subroutine r8mkspline(x,nx,fspl,ibcxmin,bcxmin,ibcxmax,bcxmax,
14 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
112 real*8,
dimension(:,:),
allocatable :: fspl4
113 real*8,
dimension(:),
allocatable :: wk
117 allocate(fspl4(4,nx),wk(nx))
130 call r8cspline(x,nx,fspl4,ibcxmin,bcxmin,ibcxmax,bcxmax,
138 fspl(2,i)=2.0_r8*fspl4(3,i)
140 fspl(2,nx)=2.0_r8*fspl4(3,nx-1) +
141 > (x(nx)-x(nx-1))*6.0_r8*fspl4(4,nx-1)