X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fgeomout.F;h=f8e8f415639d6932fec57a995f0e1bfd00ca1fb6;hb=f038962fbb96e0c2c6f1ccb910373ceddd5b387b;hp=416122be6fea8612909d459b0c12a0eb295847b3;hpb=0ef28a0babbfafd06d3977c622ecbe98a5f41e86;p=unres.git

diff --git a/source/unres/src_MD-M/geomout.F b/source/unres/src_MD-M/geomout.F
index 416122b..f8e8f41 100644
--- a/source/unres/src_MD-M/geomout.F
+++ b/source/unres/src_MD-M/geomout.F
@@ -454,14 +454,33 @@ c-----------------------------------------------------------------
      &          potEcomp(23),me
           format1="a133"
           print *,'A CHUJ',potEcomp(23)
-        else
           write (line1,'(i10,f15.2,7f12.3,i5,$)')
      &           itime,totT,EK,potE,totE,
      &           kinetic_T,t_bath,gyrate(),
      &           potEcomp(23),me
           format1="a114"
         endif
-       else
+       else if (SELFGUIDE.gt.0) then
+       distance=0.0
+       do j=1,3
+       distance=distance+c(j,afmend)-c(j,afmbeg)
+       enddo
+       distance=dsqrt(distance)
+       if (refstr) then
+         call rms_nac_nnc(rms,frac,frac_nn,co,.false.)
+          write (line1,'(i10,f15.2,3f12.3,f7.2,2f6.3,4f12.3,i5,$)')
+     &          itime,totT,EK,potE,totE,
+     &          rms,frac,frac_nn,kinetic_T,t_bath,gyrate(),
+     &          distance,me
+          format1="a133"
+          print *,'A CHUJ',potEcomp(23)
+          write (line1,'(i10,f15.2,7f12.3,i5,$)')
+     &           itime,totT,EK,potE,totE,
+     &           kinetic_T,t_bath,gyrate(),
+     &           distance,me
+          format1="a114"
+        endif
+
        if (refstr) then
          call rms_nac_nnc(rms,frac,frac_nn,co,.false.)
           write (line1,'(i10,f15.2,3f12.3,f7.2,4f6.3,3f12.3,i5,$)')