cluster_wham
[unres.git] / source / cluster / wham / src-HCD / energy_p_new.F
index e969ea3..6286514 100644 (file)
@@ -2439,12 +2439,12 @@ C          fac_shield(j)=0.6
           el1=el1*fac_shield(i)**2*fac_shield(j)**2
           el2=el2*fac_shield(i)**2*fac_shield(j)**2
           eesij=(el1+el2)
-          ees=ees+eesij*faclipij2
+          ees=ees+eesij*sss*faclipij2
           else
           fac_shield(i)=1.0
           fac_shield(j)=1.0
           eesij=(el1+el2)
-          ees=ees+eesij*faclipij2
+          ees=ees+eesij*sss*faclipij2
           endif
           evdw1=evdw1+evdwij*sss*faclipij2
 cd          write(iout,'(2(2i3,2x),7(1pd12.4)/2(3(1pd12.4),5x)/)')