wham debug on

diff --git a/source/wham/src-M/enecalc1.F b/source/wham/src-M/enecalc1.F
index ab0a408..913c2e1 100644 (file)
--- a/source/wham/src-M/enecalc1.F
+++ b/source/wham/src-M/enecalc1.F
@@ -202,6 +202,9 @@ 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))
 c              call intout
+          write (iout,*) "The Cartesian geometry is:"
+          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)
               call pdbout(indstart(me1)+iii,
      & 1.0d0/(1.987D-3*beta_h(ib,ipar)),
      &energia(0),eini,0.0d0,0.0d0)

diff --git a/source/wham/src-M/energy_p_new.F b/source/wham/src-M/energy_p_new.F
index 44e5cba..2418fd3 100644 (file)
--- a/source/wham/src-M/energy_p_new.F
+++ b/source/wham/src-M/energy_p_new.F
@@ -823,8 +823,8 @@ C
             IF (dyn_ss_mask(i).and.dyn_ss_mask(j)) THEN
               call dyn_ssbond_ene(i,j,evdwij)
               evdw=evdw+evdwij
-C            write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
-C     &                        'evdw',i,j,evdwij,' ss',evdw,evdw_t
+            write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
+     &                        'evdw',i,j,evdwij,' ss',evdw,evdw_t
 C triple bond artifac removal
              do k=j+1,iend(i,iint)
 C search over all next residues
@@ -835,8 +835,8 @@ C              write(iout,*) 'k=',k
 C call the energy function that removes the artifical triple disulfide
 C bond the soubroutine is located in ssMD.F
               evdw=evdw+evdwij
-C             write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
-C     &                        'evdw',i,j,evdwij,'tss',evdw,evdw_t
+             write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
+     &                        'evdw',i,j,evdwij,'tss',evdw,evdw_t
               endif!dyn_ss_mask(k)
              enddo! k
             ELSE