33#include "implicit_f.inc"
37 INTEGER , N, INDX(*),NODREAC(*)
38 my_real a(3,*),fthreac(6,*),dt3,fthdtm(6,*)
45#include "vectorize.inc"
48 fthreac(n+1,nodreac(l)) = fthreac(n+1,nodreac(l))
50 & - fthdtm(n+1,nodreac(l))*dt3/two
51 fthreac(n+2,nodreac(l)) = fthreac(n+2,nodreac(l))
53 & - fthdtm(n+2,nodreac(l))*dt3/two
54 fthreac(n+3,nodreac(l)) = fthreac(n+3,nodreac(l))
56 & - fthdtm(n+3,nodreac(l))*dt3/two
57 fthdtm(n+1,nodreac(l)) = a(1,l)
58 fthdtm(n+2,nodreac(l)) = a(2,l)
59 fthdtm(n+3,nodreac(l)) = a(3,l)
72#include "implicit_f.inc"
76 INTEGER NINDX, N, INDX(*)
83#include "vectorize.inc"
86 fanreac(n+1,l) = - a(1,l)
87 fanreac(n+2,l) = - a(2,l)
88 fanreac(n+3,l) = - a(3,l)
subroutine bcs1th_imp(nindx, indx, a, fthreac, nodreac, n, fthdtm, dt3)