Hi,
If you look at the line 871-872 under commands [newf] and [newf, zero]:
do n = i,nneq-1,ndf
hr(np(28)+i+2*nneq) = 0.0d0 ! f0(n,1) = fixed forces
hr(np(28)+i+3*nneq) = 0.0d0 ! f0(n,2) = fixed displacements
end do ! n
It seems like the loop index i in hr(np(28)+i+2*nneq) and hr(np(28)+i+3*nneq) = 0.0d0 should be replaced with n.
Best,
Ziyu
