VECFEM3 Reference Manual: vemavs

Type: FORTRAN routine

NAME

vemavs - writes nodal values and the mesh data to an AVS UCD file

SYNOPSIS

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

PURPOSE

vemavs crates an UCD file (unstructured cell data file) that can be read by AVS, see [AVS], and other visualization tools (e.g Explorer). The data file contains the geometrical mesh as well as the result data at the geometrical nodes. In the geometry only point, line, triangle, quadrilateral, tetrahedron, prism and hexahedron elements of order one and two are considered, see mesh, are considered. Order three elements are ignored. The cell-id is given by the first and the mat_id by the second integer vector parameters. If there is no second parameter it is set to one. If there is no first parameter the element group is ignored. The routine writes scalar and vector nodal results. In the last case components missed to form a three dimensional vector are set to zero. The label of the component is given by TEXT1. Moreover data can be write without being interpreted as scalar or vector data. Then the label of the i-th component is TEXT1_<i>. the result data can be created by using the AVS modul ucd_extract. Only values lower than 10.D0**IVEM(4) are considered. To distribute a solution at the global nodes onto the geometrical nodes vemu05 can be used. The model-id is TEXT2 with one model component given by the load case id CASE and two additional components STEP and T if nonsteady data are written.

ARGUMENTS

TEXT1 character*80, scalar, input, local

Component label

TEXT2 character*80, scalar, input, local

model id

T double precision, scalar, input, global

Time.

LU integer, scalar, input, local

Length of the input solution vector, LU>=U1*N.

U double precision, array: U(LU), input, local

The input solution vector at the geometrical nodes. U(U1*(j-1)+i) is the value of the j-th component at the i-th geometrical node on the process MYPROC. If a geometrical node has different values on two processes, one arbitrarily selected value of the two is written.

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.
95	`IVEM` or `U` is too small.
98	write error on unit of UCD file.
99	fatal error.

(4), input, local

Only entries of U lower than 10.D0**IVEM(4) are written.

(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.

(127), input, local

Unit of the UCD file.

(128)=U1, input, local

Leading dimension of the input solution array, U1>=NDEG.

(129)=N, input, global

Number of components of the input solution.

(130)=FTYPE, input, global

Data type

0	unknown (`N` arbitrary)
1	-scalar (`N`=1)
2	-vector (`N`=3)

If N does not fit, the missing components are set to zero and useless values are cut.

(137)=CASE, input, global

Load set counter.

(138)=STEP, input, global

Time step counter. If STEP<0, the steady analysis type is assumed.

(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 vemavs 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. 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.

REFERENCES

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

COPYRIGHTS

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