1 files changed, 124 insertions, 0 deletions
diff --git a/CarpetDev/CarpetIOF5_standalone/src/timestep.cc b/CarpetDev/CarpetIOF5_standalone/src/timestep.cc
new file mode 100644
index 000000000..409f47221
--- /dev/null
+++ b/CarpetDev/CarpetIOF5_standalone/src/timestep.cc
@@ -0,0 +1,124 @@
+#include <cassert>
+#include <iomanip>
+#include <limits>
+#include <string>
+#include <sstream>
+
+#include <hdf5.h>
+
+#include "cctk.h"
+#include "cctk_Functions.h"
+
+#include "defs.hh"
+
+#include "timestep.hh"
+#include "utils.hh"
+
+
+
+namespace CarpetIOF5 {
+  
+  namespace F5 {
+    
+    using std::numeric_limits;
+    using std::setprecision;
+    using std::ostringstream;
+    
+    
+    
+    timestep_t::
+    timestep_t (file_t & file,
+                CCTK_REAL const time,
+                char const * const name)
+      : m_file (file),
+        m_time (time)
+    {
+      // Construct a group name from the simulation time
+      ostringstream buf;
+      if (name != 0)
+      {
+        // Use the name argument if it is present
+        buf << name;
+      }
+      else
+      {
+        // Create a string from the time without losing information
+        int const precision = numeric_limits<CCTK_REAL>::digits10 + 2;
+        buf << setprecision (precision) << "T=" << time;
+      }
+      m_name = buf.str();
+      
+      m_hdf5_timestep
+        = open_or_create_group (m_file.get_hdf5_file(), m_name.c_str());
+      assert (m_hdf5_timestep >= 0);
+      write_or_check_attribute (m_hdf5_timestep, "Time", time);
+      
+      assert (invariant());
+    }
+    
+    
+    
+    timestep_t::
+    ~ timestep_t()
+    {
+      check (not H5Gclose (m_hdf5_timestep));
+    }
+    
+    
+      
+    file_t & timestep_t::
+    get_file ()
+      const
+    {
+      return m_file;
+    }
+    
+    
+    
+    CCTK_REAL timestep_t::
+    get_time ()
+      const
+    {
+      return m_time;
+    }
+    
+    
+    
+    hid_t timestep_t::
+    get_hdf5_timestep()
+      const
+    {
+      return m_hdf5_timestep;
+    }
+    
+    
+    
+    void timestep_t::
+    get_link_destination (int const proc,
+                          string & filename,
+                          string & objectname)
+      const
+    {
+      get_file().get_link_destination (proc, filename, objectname);
+      if (objectname.empty())
+      {
+        objectname = m_name;
+      }
+      else
+      {
+        objectname += string ("/") + m_name;
+      }
+    }
+    
+    
+    
+    bool timestep_t::
+    invariant()
+      const
+    {
+      return m_hdf5_timestep >= 0;
+    }
+    
+  } // namespace F5
+
+} // namespace CarpetIOF5