CarpetIOHDF5: zero cctk_nghostzones attribute unless ghost zones are written

to disk

1 files changed, 1 insertions, 1 deletions

diff --git a/Carpet/CarpetIOHDF5/src/OutputSlice.cc b/Carpet/CarpetIOHDF5/src/OutputSlice.cc
index a74c865bb..1dfa369f5 100644
--- a/Carpet/CarpetIOHDF5/src/OutputSlice.cc
+++ b/Carpet/CarpetIOHDF5/src/OutputSlice.cc
@@ -1404,7 +1404,7 @@ namespace CarpetIOHDF5 {
         const b2vect obnds       = vhh.at(m)->outer_boundaries(rl,c);
         const i2vect ghost_width = arrdata.at(gi).at(m).dd->ghost_widths.AT(rl);
         for (int d = 0; d < outdim; d++) {
-          nghostzones[d] = ghost_width[0][dirs[d]];
+          nghostzones[d] = output_ghost_points ? ghost_width[0][dirs[d]] : 0;
           assert (all (ghost_width[0] == ghost_width[1]));
 
           bbox[2*d] = obnds[0][dirs[d]];