nind=2
endif
do i=nnt+1,nct-1
+ mnum=molnum(i)
! if (iabs(itype(i,1)).eq.ntyp1) cycle
DM(ind)=2*ip4+mp/2
- if (iabs(itype(i,1)).eq.10 .or. iabs(itype(i,1)).eq.ntyp1) then
+ if (iabs(itype(i,1)).eq.10 .or. &
+ iabs(itype(i,mnum)).eq.ntyp1_molec(mnum) .or. mnum.eq.5) then
if (iabs(itype(i,1)).eq.10) DM(ind)=DM(ind)+msc(10)
ind=ind+1
else
endif
do i=nnt+1,nct-1
! if (iabs(itype(i,1)).eq.ntyp1) cycle
+ mnum=molnum(i)
DM(ind)=2*ip4+mp(1)/2
- if (iabs(itype(i,1)).eq.10 .or. iabs(itype(i,1)).eq.ntyp1) then
+ if (iabs(itype(i,1)).eq.10 .or. &
+ iabs(itype(i,mnum)).eq.ntyp1_molec(mnum) .or. mnum.eq.5) then
if (iabs(itype(i,1)).eq.10) DM(ind)=DM(ind)+msc(10,1)
ind=ind+1
else