From 30d2aeeeca6bb0c92def1bd16f5e8b7be0f04e3a Mon Sep 17 00:00:00 2001
From: schnetter <>
Date: Wed, 10 Dec 2003 17:08:00 +0000
Subject: Add the content of the cGH structure to each dataset.  This allows to

Add the content of the cGH structure to each dataset.  This allows to
extract all the information independent of FlexIO.  There is also a
version number "carpet_flexio_version" that will increase for every
incompatible change.

darcs-hash:20031210170853-07bb3-31f1ce9c51d6dfb67d265f1ee806522c7391aac5.gz
---
 CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc | 66 ++++++++++++++++++++++++++++++
 1 file changed, 66 insertions(+)

(limited to 'CarpetAttic')

diff --git a/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc b/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc
index 2a00ebe5c..ccd0af66e 100644
--- a/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc
+++ b/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc
@@ -72,6 +72,15 @@ namespace CarpetIOFlexIO {
 				const char* const varlist, const int vindex);
   static void SetFlag (int index, const char* optstring, void* arg);
   
+  static void WriteAttribute (IObase* writer, const char* name,
+                              int value);
+  static void WriteAttribute (IObase* writer, const char* name,
+                              const int* values, int nvalues);
+  static void WriteAttribute (IObase* writer, const char* name,
+                              CCTK_REAL value);
+  static void WriteAttribute (IObase* writer, const char* name,
+                              const CCTK_REAL* values, int nvalues);
+
   
   
   int CarpetIOFlexIOStartup ()
@@ -317,6 +326,25 @@ namespace CarpetIOFlexIO {
 	  dims[d]   = (ext.shape() / ext.stride())[d];
 	}
 	amrwriter->write (origin, dims, (void*)tmp->storage());
+        
+        // Write some additional attributes
+        WriteAttribute (writer, "carpet_flexio_version", 1);
+        WriteAttribute (writer, "cctk_dim", cgh->cctk_dim);
+        WriteAttribute (writer, "cctk_iteration", cgh->cctk_iteration);
+        WriteAttribute (writer, "cctk_gsh", cgh->cctk_gsh, dim);
+        WriteAttribute (writer, "cctk_lsh", cgh->cctk_lsh, dim);
+        WriteAttribute (writer, "cctk_lbnd", cgh->cctk_lbnd, dim);
+        WriteAttribute (writer, "cctk_delta_time", cgh->cctk_delta_time);
+        WriteAttribute (writer, "cctk_delta_space", cgh->cctk_delta_space, dim);
+        WriteAttribute (writer, "cctk_origin_space", cgh->cctk_origin_space, dim);
+        WriteAttribute (writer, "cctk_bbox", cgh->cctk_bbox, 2*dim);
+        WriteAttribute (writer, "cctk_levfac", cgh->cctk_levfac, dim);
+        WriteAttribute (writer, "cctk_levoff", cgh->cctk_levoff, dim);
+        WriteAttribute (writer, "cctk_levoffdenom", cgh->cctk_levoffdenom, dim);
+        WriteAttribute (writer, "cctk_timefac", cgh->cctk_timefac);
+        WriteAttribute (writer, "cctk_convlevel", cgh->cctk_convlevel);
+        WriteAttribute (writer, "cctk_nghostzones", cgh->cctk_nghostzones, dim);
+        WriteAttribute (writer, "cctk_time", cgh->cctk_time);
       }
       
       // Delete temporary copy
@@ -698,4 +726,42 @@ namespace CarpetIOFlexIO {
   
   
   
+  void WriteAttribute (IObase* writer, const char* name, int value)
+  {
+    WriteAttribute (writer, name, &value, 1);
+  }
+  
+  void WriteAttribute (IObase* writer, const char* name,
+                       const int* values, int nvalues)
+  {
+    assert (writer);
+    assert (name);
+    assert (values);
+    vector<CCTK_INT4> values1(nvalues);
+    for (int i=0; i<nvalues; ++i) {
+      values1[i] = values[i];
+    }
+    writer->writeAttribute (name, IObase::Int32, nvalues, &values1[0]);
+  }
+  
+  void WriteAttribute (IObase* writer, const char* name, CCTK_REAL value)
+  {
+    WriteAttribute (writer, name, &value, 1);
+  }
+  
+  void WriteAttribute (IObase* writer, const char* name,
+                       const CCTK_REAL* values, int nvalues)
+  {
+    assert (writer);
+    assert (name);
+    assert (values);
+    vector<CCTK_REAL8> values1(nvalues);
+    for (int i=0; i<nvalues; ++i) {
+      values1[i] = values[i];
+    }
+    writer->writeAttribute (name, IObase::Float64, nvalues, &values1[0]);
+  }
+  
+  
+  
 } // namespace CarpetIOFlexIO
-- 
cgit v1.2.3