X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Funres.F;h=5731d5a4cf81dc4b86f8e3c411b97d17639fb610;hb=a43ca72434f567c64e9ab54c6c84ba9c214ab3d0;hp=632374b3cd0d7a0cd886190e32bfcaa887f907b1;hpb=b847bacd7f528a68d640b3d58d4f5b9ce3c6ede2;p=unres.git

diff --git a/source/unres/src_MD/unres.F b/source/unres/src_MD/unres.F
index 632374b..5731d5a 100644
--- a/source/unres/src_MD/unres.F
+++ b/source/unres/src_MD/unres.F
@@ -219,7 +219,6 @@ c---------------------------------------------------------------------------
       time_ene=tcpu()-time00
 #endif
       write (iout,*) "Time for energy evaluation",time_ene
-      print *,"after etotal"
       etota = energy(0)
       etot =etota
       call enerprint(energy(0))
@@ -228,18 +227,14 @@ c---------------------------------------------------------------------------
       if (minim) then
 crc overlap test
         if (overlapsc) then 
-          print *, 'Calling OVERLAP_SC'
           call overlap_sc(fail)
         endif 
 
         if (searchsc) then 
           call sc_move(2,nres-1,10,1d10,nft_sc,etot)
-          print *,'SC_move',nft_sc,etot
-          write(iout,*) 'SC_move',nft_sc,etot
         endif 
 
         if (dccart) then
-          print *, 'Calling MINIM_DC'
 #ifdef MPI
           time1=MPI_WTIME()
 #else
@@ -252,7 +247,6 @@ crc overlap test
             call chainbuild
           endif
           call geom_to_var(nvar,varia)
-          print *,'Calling MINIMIZE.'
 #ifdef MPI
           time1=MPI_WTIME()
 #else
@@ -285,8 +279,15 @@ crc overlap test
         if (refstr) call rms_nac_nnc(rms,frac,frac_nn,co,.true.)
         call briefout(0,etot)
       endif
+      potE=etot
+#ifdef NOXDR
+       if (outx) call cartout(0.0d0)
+#else
+       if (outx) call cartoutx(0.0d0)
+#endif
       if (outpdb) call pdbout(etot,titel(:32),ipdb)
       if (outmol2) call mol2out(etot,titel(:32))
+      if (iretcode.ne.4) call check_ecartint
       return
       end
 c---------------------------------------------------------------------------
@@ -323,6 +324,12 @@ c---------------------------------------------------------------------------
       call enerprint(energy(0))
       call intout
       call briefout(0,etot)
+      potE=etot
+#ifdef NOXDR
+      if (outx) call cartout(0.0d0)
+#else
+      if (outx) call cartoutx(0.0d0)
+#endif
       if (outpdb) call pdbout(etot,titel(:32),ipdb)
       if (outmol2) call mol2out(etot,titel(:32))
       if (refstr) call rms_nac_nnc(rms,frac,frac_nn,co,.true.)