1 SUBROUTINE setup_int(xknots,smesh,hx,w,w1,u,u1,nk,nots,nmesh)
3 USE vparams,
ONLY: epstan
9 INTEGER,
DIMENSION(nots) :: nk
10 REAL(rprec),
DIMENSION(nots) :: xknots
11 REAL(rprec),
DIMENSION(nmesh) :: smesh
12 REAL(rprec),
DIMENSION(nots) :: hx
13 REAL(rprec),
DIMENSION(nmesh) :: w, w1, u, u1
17 INTEGER :: ksp1, k, i, k1
18 REAL(rprec) :: smesh1, hk6
37 IF (smesh1 .le. xknots(1)) smesh(1) = xknots(1) + epstan
49 IF (smesh(i).gt.xknots(k) .and. smesh(i).le.xknots(k1))
THEN
52 w1(i) = (smesh(i)-xknots(k))/hx(k)
53 IF (w1(i)<(-epstan) .or. w1(i)>1.0+epstan) stop
'w1(i)'
55 u(i) = hk6*w(i)*(w(i)*w(i)-1.0)
56 u1(i) = hk6*w1(i)*(w1(i)*w1(i)-1.0)
59 IF (k .gt. ksp1) stop
'K>KSP1'
66 END SUBROUTINE setup_int