Add entropy to the energy diagnostic file. Use an MPI_Reduce instead of MPI_Allreduce to get the values