X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;ds=inline;f=source%2Fwham%2Fsrc%2Fenecalc1.F;h=615fe96a9fb2bc655bbe8296f853c7e7ef61a222;hb=20719f6c7872b3989ffa46371f1621ff012dd344;hp=c8901ad610e8d771ca733f30a2117901fdbf748d;hpb=4c8df5bfd461b8509158a9d5e54f705a99c69fed;p=unres.git

diff --git a/source/wham/src/enecalc1.F b/source/wham/src/enecalc1.F
index c8901ad..615fe96 100644
--- a/source/wham/src/enecalc1.F
+++ b/source/wham/src/enecalc1.F
@@ -46,7 +46,7 @@
       double precision tt
       integer snk_p(MaxR,MaxT_h,Max_parm)
       logical lerr
-      character*64 bprotfile_temp
+      character*128 bprotfile_temp
       integer scount_t(0:maxprocs-1)
       call opentmp(islice,ientout,bprotfile_temp)
       iii=0
@@ -165,8 +165,10 @@ c     &   " kfac",kfac,"quot",quot," fT",fT
      &      wtor_d,wsccor,wbond
 #endif
         call etotal(energia(0),fT)
-        if (constr_homology) energia(0)=energia(0)+
+        if (constr_homology.gt.0) energia(0)=energia(0)+
      &   waga_homology(iset)*energia(22)
+c        write (iout,*) "constr_homology",constr_homology," iset",iset,
+c     &   " waga_homology",waga_homology(iset)
 #ifdef DEBUG
         write (iout,*) "Conformation",i
         call enerprint(energia(0),fT)
@@ -215,8 +217,12 @@ c        call pdbout(ii+1,beta_h(ib,ipar),efree,energia(0),0.0d0,rmsdev)
      &         " the value read in: ",energia(0),eini," point",
      &         iii+1,indstart(me1)+iii," T",
      &         1.0d0/(1.987D-3*beta_h(ib,ipar))
+#define DEBUG
 #ifdef DEBUG
              call enerprint(energia(0),fT)
+#endif
+#undef DEBUG
+#ifdef DEBUG
              write (iout,'(4f10.5,2i5)') 0.0,energia(0),0.0,
      &       1.0d0/(beta_h(ib,ipar)*1.987D-3),
      &       0,0
@@ -382,7 +388,7 @@ c------------------------------------------------------------------------------
       include "COMMON.PROT"
       character*64 nazwa
       character*80 bxname,cxname
-      character*64 bprotfile_temp
+      character*128 bprotfile_temp
       character*3 liczba,licz
       character*2 licz2
       integer i,itj,ii,iii,j,k,l
@@ -519,7 +525,7 @@ c        else
 c     &    potE(i,iparm),-entfac(i),rms_nat,iscore
      &    potE(i,nparmset),-entfac(i),rms_nat,iscore
 c        endif
-        write (iout,'(2i5,3e15.5)') i,me,potE(i,iparm),-entfac(i)
+c        write (iout,'(2i5,3e15.5)') i,me,potE(i,iparm),-entfac(i)
 #ifndef MPI
         if (cxfile) call cxwrite(ixdrf,csingle,potE(i,nparmset),
      &    -entfac(i),rms_nat,iscore)