enddo
M_SC=0.0d0
do i=nnt,nct
- iti=itype(i)
- M_SC=M_SC+msc(iti)
+ iti=iabs(itype(i))
+ M_SC=M_SC+msc(iabs(iti))
inres=i+nres
do j=1,3
- cm(j)=cm(j)+msc(iti)*c(j,inres)
+ cm(j)=cm(j)+msc(iabs(iti))*c(j,inres)
enddo
enddo
do j=1,3
do j=1,3
pr(j)=c(j,inres)-cm(j)
enddo
- Im(1,1)=Im(1,1)+msc(iti)*(pr(2)*pr(2)+pr(3)*pr(3))
- Im(1,2)=Im(1,2)-msc(iti)*pr(1)*pr(2)
- Im(1,3)=Im(1,3)-msc(iti)*pr(1)*pr(3)
- Im(2,3)=Im(2,3)-msc(iti)*pr(2)*pr(3)
- Im(2,2)=Im(2,2)+msc(iti)*(pr(3)*pr(3)+pr(1)*pr(1))
- Im(3,3)=Im(3,3)+msc(iti)*(pr(1)*pr(1)+pr(2)*pr(2))
+ Im(1,1)=Im(1,1)+msc(iabs(iti))*(pr(2)*pr(2)+pr(3)*pr(3))
+ Im(1,2)=Im(1,2)-msc(iabs(iti))*pr(1)*pr(2)
+ Im(1,3)=Im(1,3)-msc(iabs(iti))*pr(1)*pr(3)
+ Im(2,3)=Im(2,3)-msc(iabs(iti))*pr(2)*pr(3)
+ Im(2,2)=Im(2,2)+msc(iabs(iti))*(pr(3)*pr(3)+pr(1)*pr(1))
+ Im(3,3)=Im(3,3)+msc(iabs(iti))*(pr(1)*pr(1)+pr(2)*pr(2))
enddo
do i=nnt,nct-1
c write (iout,*) "i",i," iti",iti," pr",(pr(j),j=1,3),
c & " v",(v(j),j=1,3)," vp",(vp(j),j=1,3)
do j=1,3
- L(j)=L(j)+msc(iti)*vp(j)
+ L(j)=L(j)+msc(iabs(iti))*vp(j)
enddo
c write (iout,*) "L",(l(j),j=1,3)
if (itype(i).ne.10 .and. itype(i).ne.21) then