X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Funres.F;h=ee926b35fd617bae852dee56e8a73fd5b1224524;hb=bf8b3ab324569ca7ae7bb9045ad9358f6826a150;hp=127056b203527551f02807f700de6b8b0e0b5294;hpb=862af3cfd0ee00d51a4f901d8141b05ba17a822d;p=unres.git diff --git a/source/unres/src_MD-M/unres.F b/source/unres/src_MD-M/unres.F index 127056b..ee926b3 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/)') @@ -132,9 +131,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 @@ -600,11 +605,11 @@ c Broadcast the order to compute internal coordinates to the slaves. write (istat,'(a5,20a12)')"# ", & (wname(print_order(i)),i=1,nprint_ene) write (istat,'("# ",20(1pe12.4))') - (weights(print_order(i)),i=1,nprint_ene) + & (weights(print_order(i)),i=1,nprint_ene) if (refstr) then write (istat,'(a5,20a12)')"# ", & (ename(print_order(i)),i=1,nprint_ene), - & ,"ETOT total","RMSD","nat.contact","nnt.contact" + & "ETOT total","RMSD","nat.contact","nnt.contact" else write (istat,'(a5,14a12)')"# ", & (ename(print_order(i)),i=1,nprint_ene),"ETOT total" @@ -671,6 +676,7 @@ c--------------------------------------------------------------------------- include 'COMMON.SBRIDGE' common /srutu/ icall double precision energy(0:max_ene) + print *,"A TU?" c do i=2,nres c vbld(i)=vbld(i)+ran_number(-0.1d0,0.1d0) c if (itype(i).ne.10) @@ -696,10 +702,14 @@ c enddo totT=1.d0 eq_time=0.0d0 call read_fragments + print *, "AFTER read fragments" call chainbuild_cart + print *,"chainbuild_cart" call cartprint + print *,"After cartprint" call intout icall=1 + print *,"before ETOT" call etotal(energy(0)) etot = energy(0) call enerprint(energy(0)) @@ -707,6 +717,7 @@ c enddo print *,'icheckgrad=',icheckgrad goto (10,20,30) icheckgrad 10 call check_ecartint + call check_ecartint return 20 call check_cartgrad return