From aebce4d86e5cad0efb21c880f0a3e59d147e21a5 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Thu, 6 Dec 2012 21:08:25 -0500 Subject: [PATCH] Write entropy at every generation in HDF5 state points. --- src/hdf5_interface.F90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/hdf5_interface.F90 b/src/hdf5_interface.F90 index 24a5f6b92e..273b3cf6f0 100644 --- a/src/hdf5_interface.F90 +++ b/src/hdf5_interface.F90 @@ -878,6 +878,7 @@ contains dims(1) = current_batch call h5ltmake_dataset_double_f(hdf5_state_point, "k_batch", 1, & dims, k_batch, hdf5_err) + dims(1) = current_batch*gen_per_batch call h5ltmake_dataset_double_f(hdf5_state_point, "entropy", 1, & dims, entropy, hdf5_err) end if @@ -1149,7 +1150,7 @@ contains call h5ltread_dataset_double_f(hdf5_state_point, "k_batch", & k_batch(1:restart_batch), dims, hdf5_err) call h5ltread_dataset_double_f(hdf5_state_point, "entropy", & - entropy(1:restart_batch), dims, hdf5_err) + entropy(1:restart_batch*gen_per_batch), dims, hdf5_err) end if ! Read number of realizations for global tallies