The computer model presented in this report was used to simulate the principal ground-water reservoir in the Beryl-Enterprise area, Escalante Desert, Beaver, Iron, and Washington Counties, Utah (Mower, 1981). The details of the formulation of the model, testing of its validity, and the results of predictions are discussed in the cited report. This report was prepared as part of a cooperative program with the Utah Department of Natural Resources, Division of Water Rights, to investigate the water resources of the State. It is an addendum to the principal interpretive report, and it is presented in order to make the model available to anyone desiring to use it for additional predictions.

The main program used was the finite-difference model for aquifer simulation in two dimensions documented by Trescott, Pinder, and Larson, (1976). Minor modifications were made to adapt the program to the principal ground-water reservoir in the Beryl-Enterprise area. All the modifications are listed at the top of table 1, and were related to parameter input and output, thus none of the computational subroutines were affected. The parameter arrays (table 1) and map of the area with a grid overlay (pi. 1) are given on following pages.

The model simulates an aquifer- under water-table conditions, mostly composed of unconsoliuated basin-fill deposits. The boundaries of the modeled area (pi. 1) generally coincide with the boundaries of the saturated basin fill. However, in the southwest-central part of the model, permeable consolidated rock is included; and that part of the northern boundary between the Black and Wah Wah Mountains is an arbitrary boundary in basin fill between the Beryl-Enterprise area and the Milford area that lies to the northeast. The ignimbrite at Table Butte also was included in the active part of the model. The model includes simulation of discharge by evapotranspiration from phreatophytes.

The areal recharge array was used to simulate recharge entering the modeled area at its boundaries and from stream infiltration in the southern corner near Enterprise. In addition, this array included discharge by wells operated during the period simulated as being under steady-state conditions (virtually 1937), and discharging wells simulating flow of water northeast to the Milford area. These wells also were included in the transient-state simulation (1937-77), although any changes in this discharge were modeled using the pumpage array (Group IV, table 1). The wells simulating outflow to the Milford area are shown on plate 1, but the wells pumping in 1937 are not shown unless they also were pumped during 1937-77.

The pumpage array was used to simulate: (1) Discharge from wells, (2) discharge after 1977 from a mine in the southwest-central part of the model and recharge resulting form the mine discharge (pi. 1), and (3) changes in discharge in wells operated during the steady-state period. Recharge from irrigation was simulated by reducing pumpage from nodes where irrigation occurs. Discharge from all wells was reduced by 5 percent by multiplying all pumpage by 0.95 in the computer program. North of Newcastle, in T. 35 S., R. 15 W., pumpage was reduced by 35 percent because surface materials are very permeable.