X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fcluster%2Fwham%2Fsrc%2Fprobabl.F;h=bdca32a5328e5fc234f327ed232065d5f5fa2c7b;hb=6afd4863682c2aea9d9b514beab6404c19164a57;hp=53735e872004bda44f0b7274354d06cc0621029a;hpb=45b3172c2641f2fb789c51059f74ae14e7c82b10;p=unres.git diff --git a/source/cluster/wham/src/probabl.F b/source/cluster/wham/src/probabl.F index 53735e8..bdca32a 100644 --- a/source/cluster/wham/src/probabl.F +++ b/source/cluster/wham/src/probabl.F @@ -176,16 +176,20 @@ cc if (wcorr6.eq.0) ecorr6=0.0d0 #ifdef MPI Fdimless_(i)=beta_h(ib)*etot+entfac(ii) totfree_(i)=etot +#ifdef DEBUG + write (iout,*) "etrop", i,ii,ib, + & 1.0d0/(1.987d-3*beta_h(ib)),totfree(i), + & entfac(ii),Fdimless_(i) +#endif #else Fdimless(i)=beta_h(ib)*etot+entfac(ii) totfree(i)=etot -#endif #ifdef DEBUG - write (iout,*) "etrop", i,ii,ib, & 1.0d0/(1.987d-3*beta_h(ib)),totfree(i), & entfac(ii),Fdimless(i) #endif +#endif enddo ! i #ifdef MPI call MPI_Gatherv(Fdimless_(1),scount(me), @@ -232,14 +236,12 @@ c write (iout,*) "qfree",qfree sumprob=0.0 do i=1,min0(ncon,maxstr_proc)-1 sumprob=sumprob+dexp(dble(-fdimless(i)+fdimless(1)))/qfree -c#define DEBUG #ifdef DEBUG write (iout,*) 'i=',i,ib,beta_h(ib), & 1.0d0/(1.987d-3*beta_h(ib)),list_conf(i), & totfree(list_conf(i)), & -entfac(list_conf(i)),fdimless(i),sumprob #endif -c#undef DEBUG if (sumprob.gt.prob_limit) goto 122 c if (sumprob.gt.1.00d0) goto 122 nlist=nlist+1