1 files changed, 67 insertions, 0 deletions
diff --git a/CarpetDev/CarpetIOF5/src/physical_quantity.cc b/CarpetDev/CarpetIOF5/src/physical_quantity.cc
new file mode 100644
index 000000000..0d6640442
--- /dev/null
+++ b/CarpetDev/CarpetIOF5/src/physical_quantity.cc
@@ -0,0 +1,67 @@
+#include <cassert>
+#include <cstdlib>
+
+#include "cctk.h"
+
+#include "physical_quantity.hh"
+#include "utils.hh"
+
+
+
+namespace CarpetIOF5 {
+  
+  namespace F5 {
+    
+    physical_quantity_t::
+    physical_quantity_t (coordinate_system_t & coordinate_system,
+                         int const group)
+      : m_coordinate_system (coordinate_system),
+        m_group (group)
+    {
+      assert (group >= 0 and group < CCTK_NumGroups());
+      
+      char * const name = CCTK_GroupName (group);
+      assert (name != 0);
+      
+      m_hdf5_physical_quantity
+        = open_or_create_group (m_coordinate_system
+                                .get_hdf5_coordinate_system(),
+                                name);
+      assert (m_hdf5_physical_quantity >= 0);
+      
+      free (name);
+      
+      assert (invariant());
+    }
+    
+    
+    
+    physical_quantity_t::
+    ~ physical_quantity_t ()
+    {
+      herr_t const herr = H5Gclose (m_hdf5_physical_quantity);
+      assert (! herr);
+    }
+    
+    
+    
+    hid_t physical_quantity_t::
+    get_hdf5_physical_quantity ()
+      const
+    {
+      return m_hdf5_physical_quantity;
+    }
+    
+    
+    
+    bool physical_quantity_t::
+    invariant ()
+      const
+    {
+      return (m_group >= 0 and m_group < CCTK_NumGroups()
+              and m_hdf5_physical_quantity >= 0);
+    }
+    
+  } // namespace F5
+
+} // namespace CarpetIOF5