X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2FMD_A-MTS.F;h=995b558926d07fd68fa8b7171e418326d0b7dfb1;hb=4510e7a28accbbc6f70a49b69e5e01f43b1068a9;hp=d4d6be688b71f0664cb0610d2dab964b8b41a35c;hpb=c50460147146c32657710433c7464783b5347bf5;p=unres.git diff --git a/source/unres/src_MD/MD_A-MTS.F b/source/unres/src_MD/MD_A-MTS.F index d4d6be6..995b558 100644 --- a/source/unres/src_MD/MD_A-MTS.F +++ b/source/unres/src_MD/MD_A-MTS.F @@ -1096,7 +1096,7 @@ c Backup the coordinates, velocities, and accelerations if(tnp .or. tnp1) then #ifndef G77 - write (iout,'(a3,7f)') "TTT",EK,s_np,potE,pi_np,Csplit, + write (iout,'(a3,7f20.10)') "TTT",EK,s_np,potE,pi_np,Csplit, & E_long,energia_short(0) #else write (iout,'(a3,7f20.10)') "TTT",EK,s_np,potE,pi_np,Csplit, @@ -1828,7 +1828,8 @@ c Removing the velocity of the center of mass endif if (.not.rest) then call chainbuild - if(iranconf.ne.0) then + write (iout,*) "PREMINIM ",preminim + if(iranconf.ne.0 .or. preminim) then if (overlapsc) then print *, 'Calling OVERLAP_SC' call overlap_sc(fail) @@ -1850,8 +1851,12 @@ c Removing the velocity of the center of mass call minimize(etot,varia,iretcode,nfun) call var_to_geom(nvar,varia) endif - if(me.eq.king.or..not.out1file) - & write(iout,*) 'SUMSL return code is',iretcode,' eval ',nfun + if(me.eq.king.or..not.out1file) then + write(iout,*) "Minimized energy is",etot + write(iout,*) 'SUMSL return code is',iretcode,' eval ',nfun + call etotal(potEcomp) + call enerprint(potEcomp) + endif endif endif call chainbuild_cart