doc/MEDMEM/FIELDgeneral.py

   1 # Copyright (C) 2005  OPEN CASCADE, CEA, EDF R&D, LEG
   2 #           PRINCIPIA R&D, EADS CCR, Lip6, BV, CEDRAT
   3 #
   4 ######################################################################
   5 #                                                                    #
   6 # This Python script should be executed when the shared library is   #
   7 # generated using SWIG 1.3 (or higher) due to the fact that older    #
   8 # version could not handle the wrapping of several class constructor #
   9 #                                                                    #
  10 ######################################################################
  11 from libMEDMEM_Swig import *
  12
  13 MedFile = "pointe.med"
  14 meshName = "maa1"
  15 fieldName = "fieldcelldouble"
  16
  17 myMesh = MESH(MED_DRIVER,MedFile,meshName)
  18
  19 mySupport = SUPPORT(myMesh,"Support on CELLs",MED_CELL)
  20
  21 myField = FIELDDOUBLE(mySupport,MED_DRIVER,MedFile,fieldName)
  22
  23 numberOfComponents = myField.getNumberOfComponents()
  24
  25 for i in range(numberOfComponents):
  26     ip1 = i+1
  27     name = myField.getComponentName(ip1)
  28     desc = myField.getComponentDescription(ip1)
  29     unit = myField.getMEDComponentUnit(ip1)
  30
  31     print "Component ",ip1
  32     print "  - name       : ",name
  33     print "  - decription : ",desc
  34     print "  - unit       : ", unit
  35
  36 iterationNumber = myField.getIterationNumber()
  37 orderNumber = myField.getOrderNumber()
  38 time = myField.getTime()
  39 print "Iteration ",iterationNumber,"  at time ",time,\
  40       " (and order number ",orderNumber,")"
  41
  42 numberOfValue = mySupport.getNumberOfElements(MED_ALL_ELEMENTS)
  43 value = myField.getValue(MED_FULL_INTERLACE)
  44
  45 for i in range(numberOfValue):
  46     print "  * ",value[i*numberOfComponents:(i+1)*numberOfComponents]