add entropy to output H% file

evanoconnor · Jul 13, 2017 · fb4535e · fb4535e
1 parent 7308d47
commit fb4535e
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/restart_H5.F90 b/src/restart_H5.F90
@@ -496,6 +496,14 @@ subroutine restart_output_h5
   call h5sclose_f(dspace_id, error)
   cerror = cerror + error
 
+  call h5screate_simple_f(rank, dims1, dspace_id, error)
+  call h5dcreate_f(file_id, "entropy", H5T_NATIVE_DOUBLE, dspace_id,&
+       & dset_id, error)
+  call h5dwrite_f(dset_id, H5T_NATIVE_DOUBLE, entropy, dims1, error)
+  call h5dclose_f(dset_id, error)
+  call h5sclose_f(dspace_id, error)
+  cerror = cerror + error
+
   call h5screate_simple_f(rank, dims1, dspace_id, error)
   call h5dcreate_f(file_id, "ye", H5T_NATIVE_DOUBLE, dspace_id,&
        & dset_id, error)