X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-HCD-5D%2Fenecalc1.F;h=60addc7a4312601272ce7bac3d2312ac2a04cc73;hb=7ac668d0f9ef861042efc51cbcd60c4e5c0239f4;hp=f037ae8a19920e5a3dd3d6559c9c6047d909fb21;hpb=1bcb3c1de31a7a0295279137b22719a799fa617d;p=unres.git

diff --git a/source/wham/src-HCD-5D/enecalc1.F b/source/wham/src-HCD-5D/enecalc1.F
index f037ae8..60addc7 100644
--- a/source/wham/src-HCD-5D/enecalc1.F
+++ b/source/wham/src-HCD-5D/enecalc1.F
@@ -163,8 +163,8 @@ C        write (iout,*) "tuz przed energia"
 C        write (iout,*) "tuz za energia"
 #ifdef DEBUG
         write (iout,*) "Conformation",i
-c          write (iout,'(8f10.5)') ((c(l,k),l=1,3),k=1,nres),
-c     &                            ((c(l,k+nres),l=1,3),k=nnt,nct)
+        write (iout,'(8f10.5)') ((c(l,k),l=1,3),k=1,nres),
+     &                            ((c(l,k+nres),l=1,3),k=nnt,nct)
         call enerprint(energia(0),fT)
 c        write (iout,'(2i5,21f8.2)') i,iparm,(energia(k),k=1,21)
 c        write (iout,*) "ftors(1)",ftors(1)
@@ -200,6 +200,8 @@ c        call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev)
           if (ipar.eq.iparm) write (iout,*) i,iparm,
      &      1.0d0/(beta_h(ib,ipar)*1.987D-3),eini,energia(0)
 #endif
+c          write (iout,*) "eini",eini,"energia(0)",energia(0)," diff",
+c     &       eini-energia(0)
           if (ipar.eq.iparm .and. einicheck.gt.0 .and. 
 !     &      dabs(eini-energia(0)-energia(27)).gt.tole) then
      &      dabs(eini-energia(0)).gt.tole) then
@@ -262,6 +264,7 @@ c          call enerprint(energia(0),fT)
         iii=iii+1
         if (q(1,iii).le.0.0d0 .and. indpdb.gt.0)
      &    q(1,iii)=qwolynes(0,0,ipermin)
+c        write (iout,*) "iii",iii," q",q(1,iii)
         write (ientout,rec=iii) 
      &   ((csingle(l,k),l=1,3),k=1,nres),
      &   ((csingle(l,k+nres),l=1,3),k=nnt,nct),