X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Fenergy_split-sep.F;h=97442a3f630f68c999566a098d15b8b77ff0a976;hb=73ca05fce4b785367bd9600c24a7f026feff14bf;hp=81e4d81135b79701b0cadc23a7637358a9da5d21;hpb=478a9d9a1c99eb3f4bc4ca676ff3162bdd01d633;p=unres.git

diff --git a/source/unres/src_MD/energy_split-sep.F b/source/unres/src_MD/energy_split-sep.F
index 81e4d81..97442a3 100644
--- a/source/unres/src_MD/energy_split-sep.F
+++ b/source/unres/src_MD/energy_split-sep.F
@@ -263,6 +263,7 @@ cMS$ATTRIBUTES C ::  proc_proc
       include 'COMMON.CHAIN'
       include 'COMMON.VAR'
       include 'COMMON.LOCAL'
+      include 'COMMON.CONTROL'
 
 c      write(iout,'(a,i2)')'Calling etotal_short ipot=',ipot
 c      call flush(iout)
@@ -395,6 +396,16 @@ C
 C Calculate electrostatic (H-bonding) energy of the main chain.
 C
   107 continue
+
+C     BARTEK for dfa test!
+      if (wdfa_dist.gt.0) call edfad(edfadis)
+c      print*, 'edfad is finished!', edfadis
+      if (wdfa_tor.gt.0) call edfat(edfator)
+c      print*, 'edfat is finished!', edfator
+      if (wdfa_nei.gt.0) call edfan(edfanei)
+c      print*, 'edfan is finished!', edfanei
+      if (wdfa_beta.gt.0) call edfab(edfabet)
+c      print*, 'edfab is finished!', edfabet
 c
 c Calculate the short-range part of Evdwpp
 c
@@ -426,6 +437,14 @@ C
 C Calculate the virtual-bond torsional energy.
 C
       call etor(etors,edihcnstr)
+c
+c Homology restraints
+c
+      if (constr_homology.ge.1) then
+        call e_modeller(ehomology_constr)
+      else
+        ehomology_constr=0.0d0
+      endif
 C
 C 6/23/01 Calculate double-torsional energy
 C
@@ -467,6 +486,11 @@ C
       energia(21)=esccor
       energia(22)=evdw_p
       energia(23)=evdw_m
+      energia(24)=ehomology_constr
+      energia(25)=edfadis
+      energia(26)=edfator
+      energia(27)=edfanei
+      energia(28)=edfabet
 c      write (iout,*) "ETOTAL_SHORT before SUM_ENERGY"
       call flush(iout)
       call sum_energy(energia,.true.)