1 files changed, 87 insertions, 0 deletions
diff --git a/src/AmrGridReader.cc b/src/AmrGridReader.cc
new file mode 100644
index 0000000..3ea8f90
--- /dev/null
+++ b/src/AmrGridReader.cc
@@ -0,0 +1,87 @@
+// AmrGridReader
+#include "AmrGridReader.hh"
+#include <stdio.h>
+#include <stdlib.h>
+
+// Other stuff
+AmrGrid *AmrGridReader::getGridInfo(AmrGrid &g,int index){
+  g.dataveclen=1;
+  if(file.seek(index)<index)
+    return 0; // fail if index past end
+  IObase::DataType dt;
+  file.readInfo(dt,g.rank,g.dims);
+  g.datatype = dt;
+  g.nbytes = IObase::nBytes(dt,g.rank,g.dims);
+  // find the deepest level (finest time resolution)
+  // Attrib Names?
+  IObase::DataType atype;
+  int length;
+  int attrnum=file.readAttributeInfo("level",atype,length);
+  if(attrnum>=0){
+    int lev; // should be Int level
+    file.readAttribute(attrnum,&lev);
+    if(lev>g.maxlevel) g.maxlevel=lev;
+    g.level=lev;
+  }
+  // check for existence of attribute named "time_refinement"
+  attrnum=file.readAttributeInfo("time_refinement",atype,length);
+  // "time_refinement" exists, so lets read it into the member g.timerefinement
+  if(attrnum>=0){
+    file.readAttribute(attrnum,&(g.timerefinement));
+  }
+  // check for existence of the attribute named "spatial_refinement"
+  attrnum=file.readAttributeInfo("spatial_refinement",atype,length);
+  // it exists, so lets read it into the datastruct member g.spatialrefinement
+  if(attrnum>=0){
+    file.readAttribute(attrnum,&(g.spatialrefinement));
+  }
+
+  attrnum=file.readAttributeInfo("timestep",atype,length);
+  if(attrnum>=0){
+    file.readAttribute(attrnum,&(g.timestep));
+  }
+
+  attrnum=file.readAttributeInfo("origin",atype,length);
+  if(attrnum>=0)
+    file.readAttribute(attrnum,(g.origin));
+
+  attrnum=file.readAttributeInfo("delta",atype,length);
+  if(attrnum>=0)
+    file.readAttribute(attrnum,(g.delta));
+
+  attrnum=file.readAttributeInfo("min_ext",atype,length);
+  if(attrnum>=0)
+    file.readAttribute(attrnum,(g.minext));
+
+  attrnum=file.readAttributeInfo("max_ext",atype,length);
+  if(attrnum>=0)
+    file.readAttribute(attrnum,(g.maxext));
+
+  attrnum=file.readAttributeInfo("persistence",atype,length);
+  if(attrnum>=0){
+    file.readAttribute(attrnum,&(g.persistence));
+    g.maxtime = g.timestep + g.persistence;
+  }
+
+  attrnum=file.readAttributeInfo("iorigin",atype,length);
+  if(attrnum>=0)
+    file.readAttribute(attrnum,&(g.iorigin));
+  
+  attrnum=file.readAttributeInfo("grid_placement_refinement",atype,length);
+  if(attrnum>=0)
+    file.readAttribute(attrnum,&(g.gridplacementrefinement));
+  
+  g.data=0;
+  return &g;
+} // done 
+  
+AmrGrid *AmrGridReader::getGridData(AmrGrid &g,int index){
+  IObase::DataType atype;
+  // if(data) free(data); data=0; // make certain it is empty first
+  g.data = malloc(g.nbytes);
+  file.seek(index);
+  file.readInfo(atype,g.rank,g.dims);
+  g.datatype=atype;
+  file.read(g.data);
+  return &g;
+}