302
Appendix B
13. Subroutine COEFF -forms coefficient matrices of discretized water
flow and mass transport equations.
14. Subroutine FVD-calculates velocity distributions and hydrodynamic
dispersion coefficients from hydraulic head distributions.
15. Subroutine WEIT -generates upstream weighted coefficients for all
nodes.
16. Subroutine SNIM-solves the discretization equations of water flow
and mass transport equations to obtain the distributions of hydraulic
head and solute concentration.
17. Subroutine MCB2D-calls above subroutines to form a complete solution procedure.
18. Subroutine FOB-calculates predicted head and concentration values
for designated observation locations and times.
B.3 Input Data Files
The user is required to prepare the following input data files:
1. DATAN-all control integers.
2. DATAG-geometric data.
3. DATAI-initial conditions.
4. DATAB-boundary conditions.
5. DATAW-pumping and recharge data.
6. DATAS-quantity and quality ofinfiltration and return water.
7. DATAP-model parameters.
8. DAT AO-observation locations and times.
9. DAT AC-computation control numbers.
Control integers in the file DAT AN include:
NP
NE
NB
NZON
IHC
IST
IPR
NW
NTW
NHB
the total number of nodes.
the total number of elements.
the maximum number of neighboring nodes around anode.
the total number of zones for distributed parameters.
setting IHC = 1 when only solving groundwater flow problems;
setting IHC = 2 when solving coupled flow and mass transport
problems.
setting IST = 0 for a transient flow filed; setting IST = 1 for a
steady flow field.
the output control number.
the total number of extraction and injection wells.
the total number of time periods during which the extraction and
injection rates are constant.
the total number of given head boundary nodes.