Write entropy at every generation in HDF5 state points.

This commit is contained in:
Paul Romano 2012-12-06 21:08:25 -05:00
parent 4afb2baeeb
commit aebce4d86e

View file

@ -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