enddo
M_SC=0.0d0
do i=nnt,nct
- iti=iabs(itype(i))
+ iti=itype(i)
M_SC=M_SC+msc(iti)
inres=i+nres
do j=1,3
enddo
do i=nnt,nct
- iti=iabs(itype(i))
+ iti=itype(i)
inres=i+nres
do j=1,3
pr(j)=c(j,inres)-cm(j)
do i=nnt,nct
if (itype(i).ne.10) then
- iti=iabs(itype(i))
+ iti=itype(i)
inres=i+nres
Im(1,1)=Im(1,1)+Isc(iti)*(1-dc_norm(1,inres)*
& dc_norm(1,inres))*vbld(inres)*vbld(inres)
incr(j)=d_t(j,0)
enddo
do i=nnt,nct
- iti=iabs(itype(i))
+ iti=itype(i)
inres=i+nres
do j=1,3
pr(j)=c(j,inres)-cm(j)
vcm(j)=vcm(j)+mp*(vv(j)+0.5d0*d_t(j,i))
enddo
endif
- amas=msc(iabs(itype(i)))
+ amas=msc(itype(i))
summas=summas+amas
if (itype(i).ne.10) then
do j=1,3