1 // Copyright (C) 2005 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
5 #include "MEDMEM_Mesh.hxx"
6 #include "MEDMEM_Field.hxx"
8 using namespace MEDMEM;
9 using namespace MED_EN ;
11 int main (int argc, char ** argv) {
13 const string MedFile = "pointe.med" ;
14 const string MeshName = "maa1" ;
17 MESH * myMesh = new MESH(MED_DRIVER,MedFile,MeshName) ;
20 // we need a support :
21 SUPPORT * mySupport = new SUPPORT(myMesh,"Support on all CELLs",MED_CELL);
23 /* create FIELD on mySupport, with 3 components */
24 int NumberOfCompoennts = 3 ;
25 FIELD<double> myField(mySupport,NumberOfCompoennts) ;
26 const string FieldName = "fieldcelldouble" ;
27 myField.setName(FieldName) ;
29 // Components information
30 string * ComponentsNames = new string[NumberOfCompoennts] ;
31 ComponentsNames[0] = "Vx" ;
32 ComponentsNames[1] = "Vy" ;
33 ComponentsNames[2] = "Vz" ;
34 myField.setComponentsNames(ComponentsNames) ;
36 string * ComponentsDescriptions = new string[NumberOfCompoennts] ;
37 ComponentsDescriptions[0] = "vitesse selon x" ;
38 ComponentsDescriptions[1] = "vitesse selon y" ;
39 ComponentsDescriptions[2] = "vitesse selon z" ;
40 myField.setComponentsDescriptions(ComponentsDescriptions) ;
42 string * ComponentsUnits = new string[NumberOfCompoennts] ;
43 ComponentsUnits[0] = "m.s-1" ;
44 ComponentsUnits[1] = "m.s-1" ;
45 ComponentsUnits[2] = "m.s-1" ;
46 myField.setMEDComponentsUnits(ComponentsUnits) ;
48 // Iteration information :
49 int IterationNumber = 10 ; // set value to MED_NOPDT if undefined (default)
50 myField.setIterationNumber(IterationNumber) ;
52 int OrderNumber = 1 ; // set value to MED_NONOR if undefined (default)
53 myField.setOrderNumber(OrderNumber) ;
55 double Time = 3.435678 ; // in second
56 myField.setTime(Time) ;
59 int NumberOfValue = mySupport->getNumberOfElements(MED_ALL_ELEMENTS);
60 for(int i=1; i<=NumberOfValue; i++) // i^th element
61 for (int j=1; j<=NumberOfCompoennts; j++) { // j^th component
62 double myValue = (i+j) * 0.1 ;
63 myField.setValueIJ(i,j,myValue);
66 // save this new field
67 int id = myField.addDriver(MED_DRIVER) ;