X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-M%2Frmscalc.f;h=6b03ac504494b1dbb9a492084d3da536df7ccf5e;hb=68aa1d135805a14583695611acfcf637748b4466;hp=826167fce33a2292a2646a0ae5de26b9197c1ec0;hpb=478a9d9a1c99eb3f4bc4ca676ff3162bdd01d633;p=unres.git

diff --git a/source/wham/src-M/rmscalc.f b/source/wham/src-M/rmscalc.f
index 826167f..6b03ac5 100644
--- a/source/wham/src-M/rmscalc.f
+++ b/source/wham/src-M/rmscalc.f
@@ -155,20 +155,28 @@ c-------------------------------------------------------------------------
       rminrms=10.0d10
       rmsminsing=10d10
       nperm=1
+C      write (iout,*) "tu2", nres,nsup
+      noverlap=nres
+      if (nres.gt.nsup+nnt-1) noverlap=nsup+nnt-1
+c      write (iout,*) "tu3,",noverlap
       do i=1,symetr
        nperm=nperm*i
       enddo
       do kkk=1,nperm
+       write (iout,*) "kkk",kkk
        nnsup=0
-       do i=1,nres
-        if (itype(i).ne.21) then
+       do i=1,noverlap
+        if (itype(i).ne.ntyp1) then
           nnsup=nnsup+1
           do j=1,3
             cc(j,nnsup)=c(j,i)
             ccref(j,nnsup)=cref(j,i,kkk)
           enddo
+              write (iout,'(2i5,3f10.5,5x,3f10.5)')
+     &        i,nnsup,(cc(j,nnsup),j=1,3),(ccref(j,nnsup),j=1,3)
         endif
        enddo
+
        call fitsq(rms,cc(1,1),ccref(1,1),nnsup,przes,obrot,non_conv)
        if (non_conv) then
         print *,'Error: FITSQ non-convergent, jcon',jcon,i