Logical Unit 11 Preparation Guide

Instructions for preparation of logical unit 11
(binary input for geometry and hydrological data)

The data have to be written by a Fortran program into an unformatted file with this set of WRITE statements:

 CHARACTER TITLEM*70

 INTEGER NPROBM,NNPB,NELB,NBNP,NBES,NTIM

 REAL*8 (A-H,O-Z)

 WRITE(11) TITLEM,NPROBM,NNPB,NELB,NBNP,NBES,NTIM

 WRITE(11) (X(N),N=1,NNP),(Y(N),N=1,NNP),(Z(N),N=1,NNP),

& ((IE(M,I),M=1,NEL),I=1,9)

 WRITE(11) TIMEM,(H,NP=1,NNP),(HT,NP=1,NNP),

& ((TH(IQ,M),IQ=1,8),M=1,NEL),

& (VDCY(NP,1),NP=1,NNP),(VDCY(NP,2),NP=1,NNP),

& (VDCY(NP,3),NP=1,NNP)

For KVI = 2 (transient flow field), the last statement has to be repeated (NTI + 1) times (initial conditions and the flow field for each time step). NTI is defined in data set 2.

The variables are:

1. TITLEM = A description of the problem.

2. NPROBM = Problem number.

3. NNPB = Number of nodal points. (This has to be the same number as NNP in data set 2.)

4. NELB = Number of elements. (This has to be the same number as NEL in data set 2.)

5. NBNP = Number of boundary nodal points.

6. NBES = Number of boundary element sides.

7. NTIM = Number of time steps.

8. X(N) = x-coordinate of node N [L].

9. Y(N) = y-coordinate of node N [L].

10. Z(N) = z-coordinate of node N [L].

11a. IE(MI,I) = Global node number of the I-th node of element MI. (I = 1,..,8)

11b. IE(MI,9) = Material type of element MI.

12. TIMEM = Time of the next simulation time step [T].

13. H = Hydraulic head. This variable is currently not used by HBGC123D but has to be present.

14. HT = Total hydraulic head. This variable is currently not used by HBGC123D but has to be present.

15. TH(IQ,M) = Water content at the IQ-th Gaussian point of element M.

16. VDCY(NP,1) = Darcian velocity in x-direction at nodal point NP.

17. VDCY(NP,2) = Darcian velocity in y-direction at nodal point NP.

18. VDCY(NP,3) = Darcian velocity in z-direction at nodal point NP.

The material type assignment as given by the output of the flow code can be overridden by using the material type correction option in Data Set 9 and setting NCM in Data Set 2 accordingly.

Note: The data format is the same for one-dimensional, two-dimensional, and three-dimensional simulations. One-dimensional simulations only use the values in x-direction (coordinates and velocities). All IE(MI,I) values should be 0 for I = 3,...,8. Two-dimensional simulations use the values in x- and z-directions. All IE(MI,I) values should be 0 for I = 4,...,8 (triangular elements) or for I = 5,..,8 (quadrilateral elements). All unused coordinates, water contents, and velocities for the one-dimensional and two-dimensional cases may be set to 0 and will be ignored.

v1.0 - 1/20/1999

Jin-Ping Gwo, email: jgwo@umbc.edu