doc/MEDMEM/FIELDcreate.cxx

   1 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   3 //
   4 using namespace std;
   5 #include "MEDMEM_Mesh.hxx"
   6 #include "MEDMEM_Field.hxx"
   7
   8 using namespace MEDMEM;
   9 using namespace MED_EN ;
  10
  11 int main (int argc, char ** argv) {
  12
  13   const string MedFile = "pointe.med" ;
  14   const string MeshName = "maa1" ;
  15
  16   /* read MESH */
  17   MESH * myMesh = new MESH(MED_DRIVER,MedFile,MeshName) ;
  18   //  myMesh->read() ;
  19
  20   // we need a support :
  21   SUPPORT * mySupport = new SUPPORT(myMesh,"Support on all CELLs",MED_CELL);
  22
  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) ;
  28
  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) ;
  35
  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) ;
  41
  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) ;
  47
  48   // Iteration information :
  49   int IterationNumber = 10 ; // set value to MED_NOPDT if undefined (default)
  50   myField.setIterationNumber(IterationNumber) ;
  51
  52   int OrderNumber = 1 ; // set value to MED_NONOR if undefined (default)
  53   myField.setOrderNumber(OrderNumber) ;
  54
  55   double Time = 3.435678 ; // in second
  56   myField.setTime(Time) ;
  57
  58   // Value :
  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);
  64     }
  65
  66   // save this new field
  67   int id = myField.addDriver(MED_DRIVER) ;
  68
  69   return 0 ;
  70 }