1 from libMEDMEM_Swig import *
4 #MedFile = "carre_quad4_3.med"
6 #meshName = "CARRE_EN_QUAD4"
8 myMesh = MESH(MED_DRIVER,MedFile,meshName)
11 nameMesh = myMesh.getName()
13 print "Mesh name : ",nameMesh
15 numberOfTypes = myMesh.getNumberOfTypes(MED_CELL)
16 print "Show Connectivity (Nodal) : "
18 # This example use access with a specified medGeometryElement through
21 for i in range(numberOfTypes):
22 cellType = myMesh.getCellType(MED_CELL,i)
23 nameType = cellType.getName()
24 type = cellType.getType()
25 numberOfElements = myMesh.getNumberOfElements(MED_CELL,type)
26 numberOfNodesPerCell = cellType.getNumberOfNodes()
27 connectivity = myMesh.getConnectivity(MED_FULL_INTERLACE,
28 MED_NODAL,MED_CELL,type)
29 print "For Type ",nameType," : "
30 for j in range(numberOfElements):
31 print "Element ",(j+1)," : ",connectivity[j*numberOfNodesPerCell:
32 (j+1)*numberOfNodesPerCell]
34 print "Show Reverse Nodal Connectivity :"
36 # This example use global access with index array
38 numberOfNodes = myMesh.getNumberOfNodes()
40 reverseNodalConnectivity = myMesh.getReverseConnectivity(MED_NODAL)
41 reverseNodalConnectivityIndex = myMesh.getReverseConnectivityIndex(MED_NODAL)
43 for i in range(numberOfNodes):
44 indexBegin = reverseNodalConnectivityIndex[i]
45 indexEnd = reverseNodalConnectivityIndex[i+1]
47 # Index value begin at 1 so (index-1) is in fact used here
49 print "Node ",(i+1)," : ",reverseNodalConnectivity[(indexBegin-1):
52 print "Show Connectivity (Descending) :"
54 # This example use global access with index array
56 numberOfElements = myMesh.getNumberOfElements(MED_CELL,MED_ALL_ELEMENTS)
57 descendingConnectivity = myMesh.getConnectivity(MED_FULL_INTERLACE,
58 MED_DESCENDING,MED_CELL,
60 descendingConnectivityIndex = myMesh.getConnectivityIndex(MED_DESCENDING,
63 for i in range(numberOfElements):
64 indexBegin = descendingConnectivityIndex[i]
65 indexEnd = descendingConnectivityIndex[i+1]
67 # Index value begin at 1 so (index-1) is in fact used here
69 print "Element ",(i+1)," : ",descendingConnectivity[(indexBegin-1):
72 print "Show Reverse Descending Connectivity :"
74 # This example use global access with index array
76 meshDimension = myMesh.getMeshDimension()
78 if (meshDimension == 1):
79 print "ERROR : Mesh Dimension = 1"
80 print "Then the Reverse Descending Connectivity could not be seen"
82 if (meshDimension == 2):
84 constituentEntity = MED_EDGE
86 if (meshDimension == 3):
88 constituentEntity = MED_FACE
90 numberOfConstituents = myMesh.getNumberOfElements(constituentEntity,
92 reverseDescendingConnectivity = myMesh.getReverseConnectivity(
94 reverseDescendingConnectivityIndex = myMesh.getReverseConnectivityIndex(
97 for i in range(numberOfConstituents):
98 indexBegin = reverseDescendingConnectivityIndex[i]
99 indexEnd = reverseDescendingConnectivityIndex[i+1]
101 # Index value begin at 1 so (index-1) is in fact used here
103 print constituent," : ",(i+1)," : ",reverseDescendingConnectivity[
104 (indexBegin-1):(indexEnd-1)]
106 print "Show ",constituent," Connectivity (Nodal) :"
108 constituentConnectivity = myMesh.getConnectivity(MED_FULL_INTERLACE,
112 constituentConnectivityIndex = myMesh.getConnectivityIndex(MED_NODAL,
115 for i in range(numberOfConstituents):
116 indexBegin = constituentConnectivityIndex[i]
117 indexEnd = constituentConnectivityIndex[i+1]
119 # Index value begin at 1 so (index-1) is in fact used here
121 print constituent," : ",(i+1)," : ",constituentConnectivity[
122 (indexBegin-1):(indexEnd-1)]