X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fmoments.f;h=438962a70f5576fea2c112909af26793fcb87d3b;hb=58b37a8c9dff4d0d448f99de2ea09225f105faa0;hp=007c0891a8769e8bcc5e13986ee94448e5b575df;hpb=478a9d9a1c99eb3f4bc4ca676ff3162bdd01d633;p=unres.git

diff --git a/source/unres/src_MD-M/moments.f b/source/unres/src_MD-M/moments.f
index 007c089..438962a 100644
--- a/source/unres/src_MD-M/moments.f
+++ b/source/unres/src_MD-M/moments.f
@@ -42,11 +42,11 @@ c   calculating the center of the mass of the protein
         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
@@ -71,12 +71,12 @@ c   calculating the center of the mass of the protein
            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
@@ -262,7 +262,7 @@ c  Calculate the angular momentum
 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