VECFEM3 Reference Manual: vegp97

Type: FORTRAN routine

Google


NAME

vegp97 - writes nodal result file for easy postprocessors

SYNOPSIS

CALL VEGP97(
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

vegp97 writes a nodal result file for easy postprocessors like gnuplot. Only one-dimesnsional meshes can be processed. The written file can be easily read from the standart input unit by the following FORTRAN77 statements, where T denotes the current time step, X denotes the coordinate of the node and U the N real values at node X at time T: In a case of static data indicated by IVEM(138)=STEP<0 use
     
     9999 READ(*,*,END=20) X1,X2,X3,(U(I),I=1,N)
	  GOTO 9999

and for time dependend data indicated by IVEM(138)=STEP=>0 use

     9999 READ(*,*,END=20) STEP,T,X1,X2,X3,(U(I),I=1,N)
	  GOTO 9999

Only values on the input vector U lower than 10.D0**IVEM(4) are considered. To distribute a solution at the global nodes onto the geometrical nodes vemu05 can be used.

ARGUMENTS

TEXT1 character*80, scalar, input, local
unused.
TEXT2 character*80, scalar, input, local
unused.
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.
0program terminated without error.
90LBIG is too small.
95IVEM or U is too small.
98write error on unit of the result file.
99fatal 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.
0only error messages are printed.
>0a protocol is printed.
(127), input, local
Unit of the result 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.
(138)=STEP, input, global
Time step counter. If STEP<0, the static data are 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 vegp97 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], [GNUPLOT].

SEE ALSO

VECFEM, vemcompile, vemrun, vemhint, mesh, vemexamples, vemdis, vemu05.

COPYRIGHTS

Program by L. Grosz, 1996. Copyrights by Universitaet Karlsruhe 1989-1996. Copyrights by Lutz Grosz 1996. All rights reserved. More details see VECFEM.


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