hidden part of debug

[unres.git] / source / wham / src-M / energy_p_new.F

diff --git a/source/wham/src-M/energy_p_new.F b/source/wham/src-M/energy_p_new.F
index beae087..83adc13 100644 (file)
--- a/source/wham/src-M/energy_p_new.F
+++ b/source/wham/src-M/energy_p_new.F
@@ -819,8 +819,8 @@ C
             IF (dyn_ss_mask(i).and.dyn_ss_mask(j)) THEN
               call dyn_ssbond_ene(i,j,evdwij)
               evdw=evdw+evdwij
-            write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
-     &                        'evdw',i,j,evdwij,' ss',evdw,evdw_t
+C            write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
+C     &                        '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
@@ -831,8 +831,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
-             write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
-     &                        'evdw',i,j,evdwij,'tss',evdw,evdw_t
+C             write (iout,'(a6,2i5,0pf7.3,a3,2f10.3)')
+C     &                        'evdw',i,j,evdwij,'tss',evdw,evdw_t
               endif!dyn_ss_mask(k)
              enddo! k
             ELSE
@@ -929,7 +929,7 @@ C Calculate the radial part of the gradient
 C Calculate angular part of the gradient.
             call sc_grad
             endif
-            write(iout,*)  "partial sum", evdw, evdw_t
+C            write(iout,*)  "partial sum", evdw, evdw_t
             ENDIF    ! dyn_ss            
           enddo      ! j
         enddo        ! iint
@@ -4659,7 +4659,7 @@ c   3 = SC...Ca...Ca...SCi
            esccor=esccor+v1ij*cosphi+v2ij*sinphi
            gloci=gloci+j*(v2ij*cosphi-v1ij*sinphi)
          enddo
-      write (iout,*)"EBACK_SC_COR",esccor,i
+C      write (iout,*)"EBACK_SC_COR",esccor,i
 c      write (iout,*) "EBACK_SC_COR",i,v1ij*cosphi+v2ij*sinphi,intertyp,
 c     & nterm_sccor(isccori,isccori1),isccori,isccori1
 c        gloc_sc(intertyp,i-3,icg)=gloc_sc(intertyp,i-3,icg)+wsccor*gloci