X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fcompare.F90;h=75a9bcea5ba364477580aea32ad59589d626c892;hb=60da3d22fd76779405459cb51d1df9003a87f209;hp=8e6c0afe3cdb3feff10982ed6725a77296eb9bf3;hpb=340832b9074a6903554ccfcf3cac2abf65e09d3a;p=unres4.git

diff --git a/source/unres/compare.F90 b/source/unres/compare.F90
index 8e6c0af..75a9bce 100644
--- a/source/unres/compare.F90
+++ b/source/unres/compare.F90
@@ -46,9 +46,9 @@
       ncont=0
       kkk=3
       do i=nnt+kkk,nct
-        iti=iabs(itype(i,1))
+        iti=iabs(itype(i,molnum(i)))
         do j=nnt,i-kkk
-          itj=iabs(itype(j,1))
+          itj=iabs(itype(j,molnum(i)))
           if (ipot.ne.4) then
 !           rcomp=sigmaii(iti,itj)+1.0D0
             rcomp=facont*sigmaii(iti,itj)
@@ -818,7 +818,7 @@
 !     &             obr,non_conv)
 !        rms=dsqrt(rms)
         call rmsd(rms)
-        print *,"before contact"
+!        print *,"before contact"
 !elte(iout,*) "rms_nacc before contact"
         call contact(.false.,ncont,icont,co)
         frac=contact_fract(ncont,ncont_ref,icont,icont_ref)