X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M-SAXS-homology%2FMD_A-MTS.F;h=4346d0079a50065b9010b00829ce30cc40af8594;hb=34d3ad3987785642be58fb2f26557d3314215577;hp=b65a7e9bb0b86d087d298cf0a7a59131388fbe97;hpb=f690e8b70bab14132839afebf080d4a28363b226;p=unres.git

diff --git a/source/unres/src_MD-M-SAXS-homology/MD_A-MTS.F b/source/unres/src_MD-M-SAXS-homology/MD_A-MTS.F
index b65a7e9..4346d00 100644
--- a/source/unres/src_MD-M-SAXS-homology/MD_A-MTS.F
+++ b/source/unres/src_MD-M-SAXS-homology/MD_A-MTS.F
@@ -417,7 +417,7 @@ c Calculate energy and forces
         t_etotal=t_etotal+tcpu()-tt0
 #endif
 #endif
-        potE=potEcomp(0)-potEcomp(20)
+        potE=potEcomp(0)-potEcomp(27)
         call cartgrad
 c Get the new accelerations
         call lagrangian
@@ -965,7 +965,7 @@ c Compute the complete potential energy
       do i=0,n_ene
         potEcomp(i)=energia_short(i)+energia_long(i)
       enddo
-      potE=potEcomp(0)-potEcomp(20)
+      potE=potEcomp(0)-potEcomp(27)
       if (ntwe.ne.0) then
       if (large.and. mod(itime,ntwe).eq.0) then
         call enerprint(potEcomp)
@@ -1826,6 +1826,7 @@ C 8/22/17 AL Minimize initial structure
 #endif
 #endif
       potE=potEcomp(0)
+c      write (iout,*) "PotE-homology",potE-potEcomp(27)
       call cartgrad
       call lagrangian
       call max_accel
@@ -1838,9 +1839,10 @@ C 8/22/17 AL Minimize initial structure
       if(me.eq.king.or..not.out1file)then 
        write(iout,*) "Potential energy and its components"
        call enerprint(potEcomp)
-       write(iout,*) (potEcomp(i),i=0,n_ene)
+c       write(iout,*) (potEcomp(i),i=0,n_ene)
       endif
-      potE=potEcomp(0)-potEcomp(20)
+      potE=potEcomp(0)-potEcomp(27)
+c      write (iout,*) "PotE-homology",potE
       totE=EK+potE
       itime=0
       if (ntwe.ne.0) call statout(itime)