2 #include "MEDMEM_Mesh.hxx"
3 #include "MEDMEM_Field.hxx"
5 using namespace MEDMEM;
7 int main (int argc, char ** argv) {
9 const string MedFile = "pointe.med" ;
10 const string MeshName = "maa1" ;
11 const string FieldName = "fieldcelldoublevector" ;
14 MESH * myMesh = new MESH(MED_DRIVER,MedFile,MeshName) ;
18 // we need a support :
19 SUPPORT * mySupport = new SUPPORT(myMesh,"Support on all Cells",MED_CELL);
20 FIELD<double> myField(mySupport,MED_DRIVER,MedFile,FieldName) ;
24 // How many components
25 int NumberOfCompoennts = myField.getNumberOfComponents() ;
27 const string * ComponentsNames = myField.getComponentsNames();
28 const string * ComponentsDescriptions = myField.getComponentsDescriptions();
29 const string * ComponentsUnits = myField.getMEDComponentsUnits();
31 for(int i=0;i<NumberOfCompoennts; i++) {
32 cout << "Component " << i << " :" <<endl ;
33 cout << " - name : " << ComponentsNames[i] << endl ;
34 cout << " - description : " << ComponentsDescriptions[i] << endl ;
35 cout << " - unit : " << ComponentsUnits[i] << endl ;
39 int IterationNumber = myField.getIterationNumber() ; // negative mean undefined
40 int OrderNumber = myField.getOrderNumber() ;
41 // internal iteration at this time iteration, negative mean undefined
42 double Time = myField.getTime() ;
44 cout << "Iteration " << IterationNumber << " at time " << Time <<
45 " (and order number " << OrderNumber << ")" << endl ;
48 int NumberOfValue = mySupport->getNumberOfElements(MED_ALL_ELEMENTS);
50 const double * Value = myField.getValue(MED_FULL_INTERLACE);
51 for(int i=0; i<NumberOfValue; i++) {
52 for(int j=0; j<NumberOfCompoennts; j++)
53 cout << Value[i*NumberOfCompoennts+j] << " " ;