X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-HCD%2Fenergy_p_new.F;fp=source%2Fwham%2Fsrc-HCD%2Fenergy_p_new.F;h=3a83918ae95d0816526f3b427e2b4009afbd7bce;hb=68aa1d135805a14583695611acfcf637748b4466;hp=e72d5584c13047539fce853ee48eb7aa5bbb7722;hpb=809cbbbaf17b645d881fe5c61387d09109acdf19;p=unres.git

diff --git a/source/wham/src-HCD/energy_p_new.F b/source/wham/src-HCD/energy_p_new.F
index e72d558..3a83918 100644
--- a/source/wham/src-HCD/energy_p_new.F
+++ b/source/wham/src-HCD/energy_p_new.F
@@ -56,7 +56,7 @@ C
        call set_shield_fac2
       endif
       call eelec(ees,evdw1,eel_loc,eello_turn3,eello_turn4)
-C            write(iout,*) 'po eelec'
+c            write(iout,*) 'po eelec eello_turn4',eello_turn4
 
 C Calculate excluded-volume interaction energy between peptide groups
 C and side chains.
@@ -2273,6 +2273,7 @@ c        write(iout,*) "JESTEM W PETLI"
         call eelecij(i,i+3,ees,evdw1,eel_loc)
         if (wturn4.gt.0.0d0 .and. itype(i+2).ne.ntyp1) 
      &   call eturn4(i,eello_turn4)
+c        write (iout,*) "i",i," eello_turn4",eello_turn4
 #ifdef FOURBODY
         num_cont_hb(i)=num_conti
 #endif
@@ -3581,6 +3582,8 @@ C Third- and fourth-order contributions from turns
       common /locel/ a_temp,agg,aggi,aggi1,aggj,aggj1,a22,a23,a32,a33,
      &    dxi,dyi,dzi,dx_normi,dy_normi,dz_normi,xmedi,ymedi,zmedi,
      &    num_conti,j1,j2
+      double precision sslipi,sslipj,ssgradlipi,ssgradlipj,faclipij
+      common /lipcalc/ sslipi,sslipj,ssgradlipi,ssgradlipj,faclipij
       j=i+3
 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
 C