do t=1,scount(me1)
c ient=-dlog(entfac(t))-entmin
ient=entfac(t)-entmin
- if (ient.le.2000) histent_p(ient)=histent_p(ient)+1
+ write (iout,*) "ient",ient,entfac(t),entmin
+C if (ient.le.2000) histent_p(ient)=histent_p(ient)+1
enddo
- call MPI_Allreduce(histent_p(0),histent(0),ientmax+1,MPI_INTEGER,
- & MPI_SUM,WHAM_COMM,IERROR)
- if (me1.eq.Master) then
- write (iout,*) "Entropy histogram"
- do i=0,ientmax
- write(iout,'(f15.4,i10)') entmin+i,histent(i)
- enddo
- endif
+C call MPI_Allreduce(histent_p(0),histent(0),ientmax+1,MPI_INTEGER,
+C & MPI_SUM,WHAM_COMM,IERROR)
+C if (me1.eq.Master) then
+C write (iout,*) "Entropy histogram"
+C do i=0,ientmax
+C write(iout,'(f15.4,i10)') entmin+i,histent(i)
+C enddo
+C endif
#else
entmin=1.0d10
entmax=-1.0d10
else
ientmax=2000
endif
- do t=1,ntot(islice)
- ient=entfac(t)-entmin
- if (ient.le.2000) histent(ient)=histent(ient)+1
- enddo
- write (iout,*) "Entropy histogram"
- do i=0,ientmax
- write(iout,'(2f15.4)') entmin+i,histent(i)
- enddo
+C do t=1,ntot(islice)
+C ient=entfac(t)-entmin
+C if (ient.le.2000) histent(ient)=histent(ient)+1
+C enddo
+C write (iout,*) "Entropy histogram"
+C do i=0,ientmax
+C write(iout,'(2f15.4)') entmin+i,histent(i)
+C enddo
#endif
#ifdef MPI