X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-M%2Fenergy_p_new.F;h=7aa46cb63f9c0a8a53c314c95d4acabe2ba79a9c;hb=e72c149dd2f7744dee9025cb8ce275ae85c6bf64;hp=eb9043d8ea58f38235c042a667dc2f2284344876;hpb=7301f2ef9629b2a4a912914b57e34e138199f87c;p=unres.git diff --git a/source/wham/src-M/energy_p_new.F b/source/wham/src-M/energy_p_new.F index eb9043d..7aa46cb 100644 --- a/source/wham/src-M/energy_p_new.F +++ b/source/wham/src-M/energy_p_new.F @@ -818,7 +818,7 @@ C do j=istart(i,iint),iend(i,iint) IF (dyn_ss_mask(i).and.dyn_ss_mask(j)) THEN call dyn_ssbond_ene(i,j,evdwij) - evdw=evdw+evdwij + evdw_t=evdw_t+evdwij write (iout,'(a6,2i5,0pf7.3,a3)') & 'evdw',i,j,evdwij,' ss' C triple bond artifac removal @@ -830,9 +830,9 @@ C write(iout,*) 'k=',k call triple_ssbond_ene(i,j,k,evdwij) C call the energy function that removes the artifical triple disulfide C bond the soubroutine is located in ssMD.F - evdw=evdw+evdwij -C if (energy_dec) write (iout,'(a6,2i5,0pf7.3,a3)') -C & 'evdw',i,j,evdwij,'tss' + evdw_t=evdw_t+evdwij + write (iout,'(a6,2i5,0pf7.3,a3)') + & 'evdw',i,j,evdwij,'tss' endif!dyn_ss_mask(k) enddo! k ELSE @@ -4658,6 +4658,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",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