1 SUBROUTINE setspline(x,weight,y,h,yfit,y2,wten,tens,nots,nb)
3 USE vparams,
ONLY: zero
10 REAL(rprec),
DIMENSION(*) :: x, weight
11 REAL(rprec),
DIMENSION(nots) :: y
12 REAL(rprec),
DIMENSION(*) :: h
13 REAL(rprec),
DIMENSION(nots) :: yfit
14 REAL(rprec),
DIMENSION(*) :: y2, wten
19 REAL(rprec) :: alsq(nots,nots)
40 CALL initspline (alsq, x, h, weight, nots)
46 CALL add_tension (alsq, wten, h, tens, zero, zero, nots, nb,
52 yfit(:nots) = y(:nots)
53 CALL solver (alsq, yfit, nots, 1, info)
57 CALL gety2 (yfit, y2, h, nots, nb)
59 END SUBROUTINE setspline