B.4. Output Data Files
305
Input data in the file DATAO include:
ELEMENT the element number where an observation weIl is located.
X
the x coordinate of an observation weIl.
Y
the y co ordinate of an observation weIl.
TIME
the observation times.
Note that the values of head and concentration at observation weIls are
calculated from their nodal values by the finite element interpolation method.
The final value of observation TIME should be equal to the total simulation
period.
Input data in the file DATAC include:
DTI the size of the initial time step.
DTF the factor of increasing the size of time step.
DTL the maximal size of time step.
ALF the relaxation coefficient in the iteration solution.
EPS the criterion of convergency.
PT
the upstream weighting coefficient.
Note that DTF, DTL, ALF, and EPS may have different values when
solving the flow subproblem and the mass transport subproblem. Usually, it
is unnecessary to change the values in this file for solving different problems.
When the user finds that there are oscillations around a concentration front,
he or she may try to increase the value of PT. Decrease the value of DTL may
cause an increase of computation time.
In order to help the user to form these input data files, an auxiliary program FINPUT.FOR is attached. After all control integers are designated in
the file DAT AN, formats of aIl other input data files can be generated automaticaIly by running this program. The user then can follow the prompts
given in these fills to fill in all problem related data. The decimallocation is
denoted for each data. The user may incorporate a mesh generation subroutine into the pro gram to form node coordinates and node numbers of
each element.
B.4 Output Data Files
After running the program MCB2D.FOR, three output files will be generated:
1. RESUL-the computational results of the pro gram.
2. DRA WC-data for drawing contours.
3. DRAWT -data for drawing time-head (time-concentration) curves.
The output of the program is controlled by integer IPR in the data file
DAT AN. First, we can set IPR = 0 to check all input data. Then set IPR = 1
305
Input data in the file DATAO include:
ELEMENT the element number where an observation weIl is located.
X
the x coordinate of an observation weIl.
Y
the y co ordinate of an observation weIl.
TIME
the observation times.
Note that the values of head and concentration at observation weIls are
calculated from their nodal values by the finite element interpolation method.
The final value of observation TIME should be equal to the total simulation
period.
Input data in the file DATAC include:
DTI the size of the initial time step.
DTF the factor of increasing the size of time step.
DTL the maximal size of time step.
ALF the relaxation coefficient in the iteration solution.
EPS the criterion of convergency.
PT
the upstream weighting coefficient.
Note that DTF, DTL, ALF, and EPS may have different values when
solving the flow subproblem and the mass transport subproblem. Usually, it
is unnecessary to change the values in this file for solving different problems.
When the user finds that there are oscillations around a concentration front,
he or she may try to increase the value of PT. Decrease the value of DTL may
cause an increase of computation time.
In order to help the user to form these input data files, an auxiliary program FINPUT.FOR is attached. After all control integers are designated in
the file DAT AN, formats of aIl other input data files can be generated automaticaIly by running this program. The user then can follow the prompts
given in these fills to fill in all problem related data. The decimallocation is
denoted for each data. The user may incorporate a mesh generation subroutine into the pro gram to form node coordinates and node numbers of
each element.
B.4 Output Data Files
After running the program MCB2D.FOR, three output files will be generated:
1. RESUL-the computational results of the pro gram.
2. DRA WC-data for drawing contours.
3. DRAWT -data for drawing time-head (time-concentration) curves.
The output of the program is controlled by integer IPR in the data file
DAT AN. First, we can set IPR = 0 to check all input data. Then set IPR = 1