ii = ind+m
mnum=molnum(i)
iti=itype(i,mnum)
- if (mnum.eq.5) then
- mscab=0.0
- else
+! if (mnum.eq.5) then
+! mscab=0.0
+! else
mscab=msc(iabs(iti),mnum)
- endif
+! endif
massvec(ii)=mscab
if (iti.ne.10 .and. iti.ne.ntyp1_molec(mnum) .and. mnum.ne.5) then