X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Funres.F;h=fab67a1d421edabb4c5911d9792561ac6d1d7e27;hb=589a8f3e38b7ce802f8577e862e54849030a5556;hp=7addf74568ce17d035d87d89192c75c50fe3080b;hpb=a70be6e9e57486949af48628f1bc8810bcbe7ce6;p=unres.git

diff --git a/source/unres/src_MD-M/unres.F b/source/unres/src_MD-M/unres.F
index 7addf74..fab67a1 100644
--- a/source/unres/src_MD-M/unres.F
+++ b/source/unres/src_MD-M/unres.F
@@ -56,7 +56,6 @@ c      call memmon_print_usage()
       if (me.eq.king) call cinfo
 C Read force field parameters and job setup data
       call readrtns
-      call flush(iout)
 C
       if (me.eq.king .or. .not. out1file) then
        write (iout,'(2a/)') 
@@ -137,9 +136,15 @@ c--------------------------------------------------------------------------
       include 'COMMON.SETUP'
       include 'COMMON.CONTROL'
       include 'COMMON.IOUNITS'
-      if (me.eq.king .or. .not. out1file)
-     &   write (iout,*) "Calling chainbuild"
+c      if (me.eq.king .or. .not. out1file) then
+c        write (iout,*) "Calling chainbuild"
+c        call flush(iout)
+c      endif
       call chainbuild
+c      if (me.eq.king .or. .not. out1file) then
+c        write (iout,*) "Calling MD"
+c        call flush(iout)
+c      endif
       call MD
       return
       end
@@ -253,6 +258,7 @@ crc overlap test
           time1=tcpu()
 #endif
           call minim_dc(etot,iretcode,nfun)
+          if(iretcode.eq.8) call check_ecartint
         else
           if (indpdb.ne.0) then 
             call bond_regular
@@ -402,7 +408,7 @@ c---------------------------------------------------------------------------
       include 'COMMON.SBRIDGE'
       double precision varia(maxvar)
       dimension ind(6)
-      double precision energy(0:max_ene)
+      double precision energy(0:n_ene)
       logical eof
       eof=.false.
 #ifdef MPI
@@ -697,7 +703,7 @@ c---------------------------------------------------------------------------
       include 'COMMON.MD'
       include 'COMMON.SBRIDGE'
       common /srutu/ icall
-      double precision energy(0:max_ene)
+      double precision energy(0:n_ene)
 c      do i=2,nres
 c        vbld(i)=vbld(i)+ran_number(-0.1d0,0.1d0)
 c        if (itype(i).ne.10) 
@@ -770,7 +776,7 @@ c---------------------------------------------------------------------------
       include 'DIMENSIONS'
       include 'COMMON.IOUNITS'
       include 'COMMON.CONTROL'
-      double precision energy(0:max_ene)
+      double precision energy(0:n_ene)
       call chainbuild
       call etotal(energy(0))
       call enerprint(energy(0))