X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Frmsd.F;fp=source%2Funres%2Fsrc_MD-M%2Frmsd.F;h=58c4a1bc768651db4cf6ce41489dfdfb73b64bf8;hb=481eb42509493d4d1e082eb01a7479377c6db79c;hp=82749b4feb95b727bb2fa30848186bc1e9c12bf8;hpb=2095e643aafa2fb192e066e69e5ca3b1a59c2727;p=unres.git diff --git a/source/unres/src_MD-M/rmsd.F b/source/unres/src_MD-M/rmsd.F index 82749b4..58c4a1b 100644 --- a/source/unres/src_MD-M/rmsd.F +++ b/source/unres/src_MD-M/rmsd.F @@ -43,6 +43,7 @@ c--------------------------------------------------------------------------- rminroz=100d2 c print *,"nz_start",nz_start," nz_end",nz_end c if (symetr.le.1) then + print *,nperm,"nperm" do kkk=1,nperm c do i=nz_start,nz_end c iatom=iatom+1 @@ -81,14 +82,14 @@ c enddo c endif c ----- diagnostics -c do kkk=1,nperm -c write (iout,*) 'Ccopy and CREFcopy' -c print '(i5,3f10.5,5x,3f10.5)',(k,(ccopy(j,k),j=1,3), -c & (crefcopy(j,k),j=1,3),k=1,iatom) -c write (iout,'(i5,3f10.5,5x,3f10.5)') (k,(ccopy(j,k),j=1,3), -c & (crefcopy(j,k),j=1,3),k=1,iatom) -c enddo -c ----- end diagnostics +C do kkk=1,nperm + write (iout,*) 'Ccopy and CREFcopy adasko',iatom + print '(i5,3f10.5,5x,3f10.5)',(k,(ccopy(j,k),j=1,3), + & (crefcopy(j,k),j=1,3),k=1,iatom) + write (iout,'(i5,3f10.5,5x,3f10.5)') (k,(ccopy(j,k),j=1,3), + & (crefcopy(j,k),j=1,3),k=1,iatom) +C enddo +Cc ----- end diagnostics c do kkk=1,nperm call fitsq(roznica,ccopy(1,1),crefcopy(1,1),iatom, & przes,obrot,non_conv) @@ -108,7 +109,7 @@ c call mpi_abort(mpi_comm_world,ierror,ierrcode) stop #endif endif -c write (iout,*) "roznica", roznica,kkk + write (iout,*) "roznica", roznica,kkk if (roznica.le.rminroz) rminroz=roznica enddo drms=dsqrt(dabs(rminroz))