ind=ind+1
ii = ind+m
iti=itype(i)
- massvec(ii)=msc(iti)
+ massvec(ii)=msc(iabs(iti))
if (iti.ne.10 .and. iti.ne.ntyp1) then
ind1=ind1+1
ii1= ind1+m1
A(ii,ii1)=1.0d0
- Gmat(ii1,ii1)=ISC(iti)
+ Gmat(ii1,ii1)=ISC(iabs(iti))
endif
enddo
c Off-diagonal elements of the dX part of A