1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
#include "AMRgridreaderPlus.h"
#include <stdio.h>
#include <stdlib.h>
#include <iostream.h>
AMRgridPlus *AMRgridreaderPlus::getGridInfo(AMRgridPlus &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;
}
attrnum=file.readAttributeInfo("time",atype,length);
if(attrnum>=0){
file.readAttribute(attrnum,&(g.time));
}attrnum=file.readAttributeInfo("time_refinement",atype,length);
if(attrnum>=0){
file.readAttribute(attrnum,&(g.timerefinement));
}
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("persistence",atype,length);
if(attrnum>=0){
file.readAttribute(attrnum,&(g.persistence));
g.maxtime = g.timestep + g.persistence;
}
attrnum=file.readAttributeInfo("time_refinement",atype,length);
if(attrnum>=0){
file.readAttribute(attrnum,&(g.timeref));
}
attrnum=file.readAttributeInfo("spatial_refinement",atype,length);
if(attrnum>=0){
file.readAttribute(attrnum,&(g.spaceref));
}
attrnum=file.readAttributeInfo("grid_placement_refinement",atype,length);
if(attrnum>=0){
file.readAttribute(attrnum,&(g.placeref));
}
attrnum=file.readAttributeInfo("range",atype,length);
if(attrnum>=0){
double range[2];
file.readAttribute(attrnum,&range);
g.scalarmin=range[0]; g.scalarmax=range[1];
}
else {
cout<<"Please convert this file with AMRPlusConv"<<endl;
exit(0);
}
g.data=0;
return &g;
} // done
AMRgridPlus *AMRgridreaderPlus::getGridData(AMRgridPlus &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;
}
|
