doc/MEDMEM/FIELDcreate.py

   1 ######################################################################
   2 #                                                                    #
   3 # This Python script should be executed when the shared library is   #
   4 # generated using SWIG 1.3 (or higher) due to the fact that older    #
   5 # version could not handle the wrapping of several class constructor #
   6 #                                                                    #
   7 ######################################################################
   8 from libMEDMEM_Swig import *
   9
  10 MedFile = "pointe.med"
  11 meshName = "maa1"
  12
  13 myMesh = MESH(MED_DRIVER,MedFile,meshName)
  14
  15 mySupport = SUPPORT(myMesh,"Support on all CELLs",MED_CELL)
  16
  17 numberOfComponents = 3
  18 myField = FIELDDOUBLE(mySupport,numberOfComponents)
  19 fieldName = "fieldcelldouble"
  20 myField.setName(fieldName)
  21
  22 for i in range(numberOfComponents):
  23     if (i == 0):
  24         name = "Vx"
  25         desc = "vitesse selon x"
  26     elif (i == 1):
  27         name = "Vy"
  28         desc = "vitesse selon y"
  29     else:
  30         name = "Vz"
  31         desc = "vitesse selon z"
  32     unit = "m. s-1"
  33     ip1 = i+1
  34     myField.setComponentName(ip1,name)
  35     myField.setComponentDescription(ip1,desc)
  36     myField.setMEDComponentUnit(ip1,unit)
  37
  38 iterationNumber = 10
  39 myField.setIterationNumber(iterationNumber)
  40
  41 orderNumber = 1
  42 myField.setOrderNumber(orderNumber)
  43
  44 time = 3.435678
  45 myField.setTime(time)
  46
  47 numberOfValue = mySupport.getNumberOfElements(MED_ALL_ELEMENTS)
  48
  49 for i in range(numberOfValue):
  50     ip1 = i+1
  51     for j in range(numberOfComponents):
  52         jp1 = j+1
  53         value = (ip1+jp1)*0.1
  54         myField.setValueIJ(ip1,jp1,value)
  55
  56 id = myField.addDriver(MED_DRIVER)