|
|
◆ WriteGmsh22()
| Sub WriteGmsh22 |
( |
FName As |
String, |
|
|
N As |
Long, |
|
|
X() As |
Double, |
|
|
Y() As |
Double, |
|
|
Z() As |
Double, |
|
|
Ne As |
Long, |
|
|
Kc() As |
Long, |
|
|
Lb() As |
Long, |
|
|
Info As |
Long |
|
) |
| |
Write FEM mesh information to Gmsh file (Version 2.2) [Experimental version]
- Purpose
- Write FEM mesh information to the file in Gmsh format (Version 2.2).
- Parameters
-
| [in] | Fname | Input file name. |
| [in] | N | Number of points. |
| [in] | X() | Array X(LX - 1) (LX >= N)
X coordinates of points. |
| [in] | Y() | Array Y(LY - 1) (LY >= N)
Y coordinates of points. |
| [in] | Z() | Array Z(LZ - 1) (LZ >= N)
Z coordinates of points. |
| [in] | Ne | Ne As Long Number of elements of FEM (including boundary elements). (Ne >= 1) |
| [out] | Kc() | Array Kc(LKc1 - 1, LKc2 - 1) (LKc1 >= Nn + 1, LKc2 >= Ne)
Connection matrix: Kc(0, K-1) is the type of K-th element, Kc(1, K-1), ..., Kc(Nn, K-1) point to the nodes of K-th element (Nn is the number of nodes in the element).
The type of element value is as follows:
= 1: 2-node line.
= 2: 3-node triangle.
= 3: 4-node quadrangle.
= 4: 4-node tetrahedron.
= 5: 8-node hexahedron.
= 6: 6-node prism.
= 7: 5-node pyramid.
= 8: 3-node second order line.
= 10: 9-node second order quadrangle.
= 11: 10-node second order tetrahedron.
= 12: 27-node second order hexahedron.
= 13: 18-node second order prism.
= 14: 14-node second order pyramid.
= 15: 1-node point.
= 16: 8-node second order quadrangle.
= 17: 20-node second order hexahedron.
= 18: 15-node second order prism.
= 19: 13-node second order pyramid.
= 20: 9-node third order incomplete triangle.
= 21: 10-node third order triangle.
= 22: 12-node fourth order incomplete triangle.
= 23: 15-node fourth order triangle.
= 24: 15-node fifth order incomplete triangle.
= 25: 21-node fifth order complete triangle.
= 26: 4-node third order edge.
= 28: 6-node fifth order edge.
= 29: 20-node third order tetrahedron.
= 30: 35-node fourth order tetrahedron.
= 31: 56-node fifth order tetrahedron.
= 92: 64-node third order hexahedron.
= 93: 125-node fourth order hexahedron. |
| [out] | Lb() | Array Lb(LLb1 - 1, LLb2 - 1) (LLb1 >= Nlb + 1, LLb2 >= Ne)
Label numbers of the elements: Lb(0, K-1) is the number of labels for K-th element, Lb(1, K-1), ..., Lb(Nlb, K-1) are the labels for the K-th element (Nlb is the number of labels for the K-th element). |
| [out] | Info | Return code.
= 0: Normal return.
= i < 0: The (-i)-th argument is invalid or has an illegal value.
= 11: File open error.
= 12: Invalid element type. |
|
1.9.7