X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MIN%2Fminimize_p.F;h=2fa0ec895281faf72e8bd7e5d58715155691fc00;hb=ba5b075e907935ee68249583f0a7cf9236cac97d;hp=1721510b3b4ee99cd3d7babca3e5f0c82f71ac1a;hpb=4677231d37cbe20d360cd65e3942b48aa4b2fe26;p=unres.git diff --git a/source/unres/src_MIN/minimize_p.F b/source/unres/src_MIN/minimize_p.F index 1721510..2fa0ec8 100644 --- a/source/unres/src_MIN/minimize_p.F +++ b/source/unres/src_MIN/minimize_p.F @@ -27,7 +27,10 @@ c common /przechowalnia/ v icall = 1 -c print *,"liv",liv," lv",lv +C Following 3 line for diagnostics; comment out if not needed + write(iout,*) "Enter MINIMIZE liv",liv," lv",lv + write (iout,*) "Coordinates before minimization" + call intout call deflt(2,iv,liv,lv,v) * 12 means fresh start, dont call deflt iv(1)=12 @@ -86,8 +89,8 @@ cd print *,'Calling SUMSL' etot=v(10) iretcode=iv(1) cd print *,'Exit SUMSL; return code:',iretcode,' energy:',etot - write (iout,'(/a,i4,a,f15.5/)') 'SUMSL return code:',iv(1), - & " energy",etot +cd write (iout,'(/a,i4,a,f15.5/)') 'SUMSL return code:',iv(1), +cd & " energy",etot call var_to_geom(nvar,x) c call chainbuild c call etotal(energia(0))