doc/tutorial/atestMEDCouplingRead.rst

   1
   2 .. _python_testMEDCouplingRead_solution:
   3
   4 Read med File
   5 ~~~~~~~~~~~~~~~~~~~
   6
   7 ::
   8
   9     import medcoupling as mc
  10
  11         medFileName = "MEDCoupling_cube3D.med"
  12         MeshName = "3Dcube"
  13         FieldName = "field"
  14         Field2DName = "fieldBottomFace"
  15
  16         # Retrieving meshes
  17         mesh3D = mc.ReadUMeshFromFile(medFileName,MeshName,0)
  18         mesh2D = mc.ReadUMeshFromFile(medFileName,MeshName,-1)
  19
  20         # Retrieving fields
  21         f = mc.ReadFieldCell(medFileName,mesh3D.getName(),0,FieldName,-1,-1)
  22         f2 = mc.ReadFieldCell(medFileName,mesh2D.getName(),-1,Field2DName,-1,-1)
  23
  24         # Retrieving Coords Mesh
  25         Coords3D = mesh3D.getCoords()
  26         Values = Coords3D.getValuesAsTuple()
  27
  28         # Retrieving field value on 0 tuple
  29         pos= Values[0]
  30         res=f.getValueOn(pos)
  31
  32         # Verify if value is OK
  33         bar = mesh3D.computeCellCenterOfMass()
  34         x=bar.getIJ(1,1)
  35         y=bar.getIJ(1,2)
  36         z=bar.getIJ(1,3)
  37
  38         from math import *
  39         d = sqrt(x*x+y*y+z*z)
  40         sinus = sin(d)
  41
  42         if abs(res[0]-sinus)<1.e-5:
  43                 print("OK")
  44         else:
  45                 print("KO")