call the module
CALL COLIB2 (M,NFNOPO,NINOPO,NDIMS,NBS,NA,NBGRO,
. NBFR,XYZ,NTRS,NARET,XYZINT,NTFR,NTGRO)
where
- M is the super array,
- NF(NI)NOPO is the support number of the O.D.S. NOPO and its level,
- NDIMS is the dimension of the output space, if NDIMS = 2, side Z is ignored (without control); this allows
us to create segments, triangles or quadrilaterals in the Rē space,
- NBS is the number of vertices in the coarse mesh,
- NA is the number of edges in the coarse mesh,
- NBGRO is the number of blocks (or elements) in the coarse mesh,
- NBFR is the number of faces to describe in order to assign a non-zero reference number to them,
- XYZ(3,.) is the coordinates of the coarse vertices (XYZ(J,I) is coordinate J of vertex I)
- NARET(5,.) is the array of coarse edges where
- NARET(1,I) is the number of the vertex which is the origin of edge I
- NARET(2,I) is the number of its end-point
- NARET(3,I) is the splitting code of the edge: 0 for a straight edge split automatically into segments of
the same length; 1 if we input the intermediate points (see XYZINT)
- NARET(4,I) is the number of control points to create on the edge (end-points excluded)
- NARET(5,I) is the edge reference
- XYZINT(3,.) are the coordinates of the control points of the edges input by the user. For each edge I for which
NARET(3,I) = 1, we enter these values as follows:
XYZINT(J,I)is coordinate J of point I, where I is the global number of the point; we first count the points
(end-points excluded) of the first edge thus described, followed by those of the second, etc.
- NTFR(6,.) is the array of the faces to specify where:
- NTFR(1,I) is the reference of face I
- NTFR(2,I) is its geometric code (3: triangle, 4: quadrilateral)
- NTFR(3,I) to NTFR(6,I) is the list of the vertices (0 for vertex '4' of a triangle)
- NTGRO(11,.) is the array of coarse elements where:
- NTGRO(1,I) is the geometric type of block I (2: segment, 3: triangle,
4: quadrilateral, 5: tetrahedron,
6: pentahedron, 7: hexahedron)
- NTGRO(2,I) to NTGRO(9,I) is the list of its vertices (0 for non-existing vertices)
- NTGRO(10,I) is the geometric type of the finite elements resulting from the splitting of the block
(see remark below)
- NTGRO(11,I) is the sub-domain number of block I (therefore that of all the sub-elements created)