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" ;
13 MESH * myMesh = new MESH(MED_DRIVER,MedFile,MeshName) ;
16 // we need a support :
17 SUPPORT * mySupport = new SUPPORT(myMesh,"Support on all CELLs",MED_CELL);
19 /* create FIELD on mySupport, with 3 components */
20 int NumberOfCompoennts = 3 ;
21 FIELD<double> myField(mySupport,NumberOfCompoennts) ;
22 const string FieldName = "fieldcelldouble" ;
23 myField.setName(FieldName) ;
25 // Components information
26 string * ComponentsNames = new string[NumberOfCompoennts] ;
27 ComponentsNames[0] = "Vx" ;
28 ComponentsNames[1] = "Vy" ;
29 ComponentsNames[2] = "Vz" ;
30 myField.setComponentsNames(ComponentsNames) ;
32 string * ComponentsDescriptions = new string[NumberOfCompoennts] ;
33 ComponentsDescriptions[0] = "vitesse selon x" ;
34 ComponentsDescriptions[1] = "vitesse selon y" ;
35 ComponentsDescriptions[2] = "vitesse selon z" ;
36 myField.setComponentsDescriptions(ComponentsDescriptions) ;
38 string * ComponentsUnits = new string[NumberOfCompoennts] ;
39 ComponentsUnits[0] = "m.s-1" ;
40 ComponentsUnits[1] = "m.s-1" ;
41 ComponentsUnits[2] = "m.s-1" ;
42 myField.setMEDComponentsUnits(ComponentsUnits) ;
44 // Iteration information :
45 int IterationNumber = 10 ; // set value to MED_NOPDT if undefined (default)
46 myField.setIterationNumber(IterationNumber) ;
48 int OrderNumber = 1 ; // set value to MED_NONOR if undefined (default)
49 myField.setOrderNumber(OrderNumber) ;
51 double Time = 3.435678 ; // in second
52 myField.setTime(Time) ;
55 int NumberOfValue = mySupport->getNumberOfElements(MED_ALL_ELEMENTS);
56 for(int i=1; i<=NumberOfValue; i++) // i^th element
57 for (int j=1; j<=NumberOfCompoennts; j++) { // j^th component
58 double myValue = (i+j) * 0.1 ;
59 myField.setValueIJ(i,j,myValue);
62 // save this new field
63 int id = myField.addDriver(MED_DRIVER) ;