VECFEM3 Reference Manual: vemens

Type: FORTRAN routine

NAME

vemens - creates EnSight geometry file

SYNOPSIS

CALL VEMENS(: NAME, LIVEM, IVEM, LNEK, NEK, LRPARM, RPARM, LIPARM, IPARM, LDNOD, DNOD, LRDPRM, RDPARM, LIDPRM, IDPARM, LNODN, NODNUM, LNOD, NOD, LNOPRM, NOPARM, LBIG, RBIG, IBIG)
INTEGER: LIVEM, LNEK, LRPARM, LIPARM, LDNOD, LRDPRM, LIDPRM, LNODN, LNOPRM, LBIG
CHARACTER*80: NAME
INTEGER: IVEM(LIVEM), NEK(LNEK), IPARM(LIPARM), DNOD(LDNOD), IDPARM(LIDPRM), NODNUM(LNODN), IBIG(*)
DOUBLE PRECISION: RPARM(LRPARM), RDPARM(LRDPRM), NOD(LDNOD), NOPARM(LNOPRM), RBIG(LBIG)

PURPOSE

EnSight is a program package for the postprocessing of the FEM results. The geometrical mesh and solutions can be plotted. The data are handed over via special files. vemens writes the geometrical nodes to EnSight geometry file. The solution is handed over to EnSight by the variable file, which is created by veen97. To read the geometry files and variable files into EnSight a result file has to be created, see [EnSight]. vemdis has to be called before the first call of vemens.

ARGUMENTS

NAME character*80, scalar, input, local

First description line

LIVEM integer, scalar, input, local

Length of the integer information vector, LIVEM>= MESH+ NINFO.

IVEM integer, array: IVEM(LIVEM), input/output, local/global

Integer information vector.

(1)=MESH, input, local

Start address of the mesh informations in IVEM, MESH>203+ NPROC.

(2)=ERR, output, global

Error number.

0	program terminated without error.
90	`LBIG` is too small.
99	fatal error.

(5)=NIVEM, output, local

Used length of IVEM.

(120)=LOUT, input, local

Unit number of the standard output file, normally 6.

(121)=OUTCNT, input, local

Output control flag, normally 1.

0	only error messages are printed.
>0	a protocol is printed.

(125), input, local

Unit of the node coordinate file.

(200)=NPROC, input, global

Number of processes, see combgn.

(201)=MYPROC, input, local

Logical process id number, see combgn.

(202)=NMSG, input/output, global

Message counter. The difference of the input and the output values gives the number of communications during the vemens call.

(204)=TIDS(1), input, global

Begin of the list TIDS which defines the mapping of the logical process ids to the physical process ids. See combgn.

(MESH), input, local

Start of mesh informations, see mesh.

LNEK integer, scalar, input, local

Length of the element array.

NEK integer, array: NEK(LNEK), input, local

Array of the elements, see mesh.

LRPARM integer, scalar, input, local

Length of the real parameter array.

RPARM double precision, array: RPARM(LRPARM), input, local

Real parameter array, see mesh.

LIPARM integer, scalar, input, local

Length of the integer parameter array.

IPARM integer, array: IPARM(LIPARM), input, local

Integer parameter array, see mesh.

LDNOD integer, scalar, input, local

Length of the array of the Dirichlet nodes.

DNOD integer, array: DNOD(LDNOD), input, local

Array of the Dirichlet nodes, see mesh.

LRDPRM integer, scalar, input, local

Length of the real Dirichlet parameter array.

RDPARM double precision, array: RDPARM(LRDPRM), input, local

Array of the real Dirichlet parameters, see mesh.

LIDPRM integer, scalar, input, local

Length of the integer Dirichlet parameter array.

IDPARM integer, array: IDPARM(LIDPRM), input, local

Array of the integer Dirichlet parameters, see mesh.

LNODN integer, scalar, input, local

Length of the array of the id numbers of the geometrical nodes.

NODNUM integer, array: NODNUM(LNODN), input, local

Array of the id numbers of the geometrical nodes, see mesh.

LNOD integer, scalar, input, local

Length of the array of the coordinates of the geometrical nodes.

NOD double precision, array: NOD(LNOD), input, local

Array of the coordinates of the geometrical nodes, see mesh.

LNOPRM integer, scalar, input, local

Length of the array of the node parameters.

NOPARM double precision, array: NOPARM(LNOPRM), input, local

Array of the node parameters, see mesh.

LBIG integer, scalar, input, local

Length of the real work array. The needed length of LBIG depends on the given mesh. A minimal length of LBIG cannot be given. It should be as large as possible.

RBIG double precision, array: RBIG(LBIG), work array, local

Real work array.

IBIG integer, array: IBIG(*), work array, local

Integer work array, RBIG and IBIG have to be defined by the EQUIVALENCE statement.

EXAMPLE

See vemexamples.

METHOD

vemens writes the files on the first process only. vemens writes only the hexahedron, tetrahedron, quadrilateral and triangle elements of order 1 and order 2. Elements, which are not allowed in EnSight are skipped. The EnSight parts are defined by the groups. If DIM=3, vemens writes the x-,y- and z-coordinates from the array NOD (if DIM=2 or DIM=1, missing coordinates are filled by zeroes). Only node coordintes are written which are referred by allowed elements.

REFERENCES

[FAQ], [DATAMAN], [DATAMAN2], [EnSight], [P_MPI].

COPYRIGHTS

by L. Grosz, Auckland , 11. June, 2000.