X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-M%2Fenergy_p_new.F;h=83adc13bdbe65ff70aaea3b25abb5e4a487cc6af;hb=bc56c61e5f80e1552da5ce3332ab99c8351e67cf;hp=beae087fc0e5cb1ad93819a8fee125a1a45fb4b5;hpb=058a610c249b4fe9ee0a273237c88b168bdc59d9;p=unres.git diff --git a/source/wham/src-M/energy_p_new.F b/source/wham/src-M/energy_p_new.F index beae087..83adc13 100644 --- 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