A modular three-dimensional finite-difference ground-water flow model

Abstract

This report presents a finite-difference model and its associated .modular computer program.The model simulates flow in three dimensions.The report includes detailed explanations of physical and mathematical concepts on which the model is based and an explanation of how those concepts were incorporated in the modular structure of the computer program.The modular structure consists of a Main Program and a series of highly indepenttent subroutines called "modules."The modules are grouped into "packages."Each package deals with a specific feature of the hydrologic system which is to be simulated, such as flow from rivers or flow into drains, or with a specific method of solving linear equations which describe the flow system, such as the Strongly Implicit Procedure or Slice-Successive Overrelaxation.The division of the program into modules permits the user to examine specific hydrologic features of the model independently.This also facilitates development of additional capabilities because new modules or packages can be added to the program without modifying the existing modules or packages.The input and output systems of the computer program are also designed to permit maximum flexibility.Ground-water flow within the aquifer is simulated using a block-centered finite-difference approach.Layers can be simulated as confined, unconfined, or a combination of confined and unconfined.Flow from external stresses, such as flow to wells, areal recharge, evapotranspiration, flow to drains, and flow through riverbeds, can also be simulated.The finite-difference equations can be solved using either the Strongly Implicit Procedure or SIice-Successive Overrelaxation.The program is written in FORTRAN '66 and will run without modification on most computers which have a FORTRAN '66 compiler.It will also run, without modification, with most extended FORTRAN '77 compilers and with minor modifications on standard FORTRAN '77 compilers.Documentation presented in this report includes a narrative description, a flow chart, a list of variables, and a program listing for each module.' ' .(25) "-m -mm In this equation, head at the beginning of the time step hi,j f k and X(LCDELR),X(LCDELC),NROW,NCOL,NLAY,IUNIT(8),IOUT) IF(IUNIT(5).6T.O) CALL EVT1RP(NEVTOP,X(LCIEVT),X(LCEVTR), X(LCEXDP),X(LCSURF),X(LCDELR),X(LCDELC),NCOL,NROW, NLAY,IUNIT(5),IOUT) IF(IUNIT(4).GT.O) CALL RIV1RP(X(LCRIVR),NRIVER,MXRIVR,IUNIT(4), IOUT) IF(IUNIT(7).6T.O) CALL GHB1RP(X(LCBNDS),NBOUND,MXBND,IUNIT( 7), IOUT) C C7C SIMULATE EACH TIME STEP.DO 200 KSTP=1,NSTP C C7C1 CALCULATE TIME STEP LENGTH.SET HOLD=HNEW..