X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc%2Fenecalc1.F;h=459f395933879c0b0b14af65a823a4194c66d9f0;hb=bbd4e14f248304c239adb46a64e66a2f6dfeb191;hp=db671ff708a4de95cc34d05542530b5b0668d2a3;hpb=28c34e29544c2c893d370595a854d963826a1531;p=unres.git diff --git a/source/wham/src/enecalc1.F b/source/wham/src/enecalc1.F index db671ff..459f395 100644 --- a/source/wham/src/enecalc1.F +++ b/source/wham/src/enecalc1.F @@ -195,7 +195,6 @@ c call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev) & iii+1,indstart(me1)+iii," T", & 1.0d0/(1.987D-3*beta_h(ib,ipar)) errmsg_count=errmsg_count+1 - call enerprint(energia(0),fT) if (errmsg_count.gt.maxerrmsg_count) & write (iout,*) "Too many warning messages" if (einicheck.gt.1) then @@ -213,12 +212,11 @@ c call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev) do k=1,21 enetb(k,iii+1,iparm)=energia(k) enddo +#define DEBUG #ifdef DEBUG write (iout,'(2i5,f10.1,3e15.5)') i,iii, & 1.0d0/(beta_h(ib,ipar)*1.987D-3),energia(0),eini,efree -c call enerprint(energia(0),fT) -#endif -#ifdef DEBUG + call enerprint(energia(0),fT) write (iout,'(8f10.5)') ((c(l,k),l=1,3),k=1,nres) write (iout,'(8f10.5)') ((c(l,k+nres),l=1,3),k=nnt,nct) write (iout,'(8f10.4)') (vbld(k),k=nnt+1,nct) @@ -231,8 +229,14 @@ c call enerprint(energia(0),fT) write (iout,'(8f10.5)') (q(k,iii+1),k=1,nQ) write (iout,'(f10.5,i10)') rmsdev,iscor call enerprint(energia(0),fT) + write(liczba,'(bz,i3.3)') me + nazwa="test"//liczba//".pdb" + write (iout,*) "pdb file",nazwa + open (ipdb,file=nazwa,position="append") call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev) + close(ipdb) #endif +#undef DEBUG endif enddo ! iparm