doc/MEDMEM/FIELDcreate.cxx

   1 using namespace std;
   2 #include "MEDMEM_Mesh.hxx"
   3 #include "MEDMEM_Field.hxx"
   4
   5 using namespace MEDMEM;
   6
   7 int main (int argc, char ** argv) {
   8
   9   const string MedFile = "pointe.med" ;
  10   const string MeshName = "maa1" ;
  11
  12   /* read MESH */
  13   MESH * myMesh = new MESH(MED_DRIVER,MedFile,MeshName) ;
  14   //  myMesh->read() ;
  15
  16   // we need a support :
  17   SUPPORT * mySupport = new SUPPORT(myMesh,"Support on all CELLs",MED_CELL);
  18
  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) ;
  24
  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) ;
  31
  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) ;
  37
  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) ;
  43
  44   // Iteration information :
  45   int IterationNumber = 10 ; // set value to MED_NOPDT if undefined (default)
  46   myField.setIterationNumber(IterationNumber) ;
  47
  48   int OrderNumber = 1 ; // set value to MED_NONOR if undefined (default)
  49   myField.setOrderNumber(OrderNumber) ;
  50
  51   double Time = 3.435678 ; // in second
  52   myField.setTime(Time) ;
  53
  54   // Value :
  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);
  60     }
  61
  62   // save this new field
  63   int id = myField.addDriver(MED_DRIVER) ;
  64
  65   return 0 ;
  66 }