1 files changed, 91 insertions, 0 deletions
diff --git a/CarpetDev/CarpetIOF5_standalone/src/data_region.cc b/CarpetDev/CarpetIOF5_standalone/src/data_region.cc
new file mode 100644
index 000000000..99b188323
--- /dev/null
+++ b/CarpetDev/CarpetIOF5_standalone/src/data_region.cc
@@ -0,0 +1,91 @@
+#include <cassert>
+#include <cstdlib>
+#include <sstream>
+#include <string>
+
+#include <hdf5.h>
+
+#include "cctk.h"
+
+#include "carpet.hh"
+
+#include "data_region.hh"
+#include "utils.hh"
+
+
+
+namespace CarpetIOF5 {
+  
+  namespace F5 {
+    
+    using std::ostringstream;
+    
+    
+    
+    data_region_t::
+    data_region_t (physical_quantity_t & physical_quantity,
+                   bbox<int, dim> const & region)
+      : m_physical_quantity (physical_quantity),
+        m_region (region),
+        m_name (string ("region-") + F5::name_from_ibbox (region))
+    {
+      assert (not region.empty());
+      
+      m_hdf5_data_region
+        = open_or_create_group (m_physical_quantity
+                                .get_hdf5_physical_quantity(),
+                                m_name.c_str());
+      assert (m_hdf5_data_region >= 0);
+      
+      assert (invariant());
+    }
+    
+    
+    
+    data_region_t::
+    ~ data_region_t ()
+    {
+      herr_t const herr = H5Gclose (m_hdf5_data_region);
+      assert (not herr);
+    }
+    
+    
+    
+    physical_quantity_t & data_region_t::
+    get_physical_quantity ()
+      const
+    {
+      return m_physical_quantity;
+    }
+    
+    
+    
+    bbox<int, dim> const & data_region_t::
+    get_region ()
+      const
+    {
+      return m_region;
+    }
+    
+    
+    
+    hid_t data_region_t::
+    get_hdf5_data_region ()
+      const
+    {
+      return m_hdf5_data_region;
+    }
+    
+    
+    
+    bool data_region_t::
+    invariant ()
+      const
+    {
+      return (not m_region.empty()
+              and m_hdf5_data_region >= 0);
+    }
+    
+  } // namespace F5
+
+} // namespace CarpetIOF5