FOR EACH SIMULATION

Data Set 0

[#Text]
Item 0 is a required comment line —“#” must be in column 1.
Text—is a character variable that starts in column 2. Text is generally used to identify input variables and arrays. For example, #Mode for Text labels the underlying integer as the CFP Mode number.


Data Set 1

MODE—is an integer value that activates conduit pipes and (or) layers.
•If MODE=1, only conduit pipes are active.
•If MODE=2, only conduit layers are active.
•If MODE=3, both conduit pipes and layers are active.
•If MODE=1, datasets 0–29 have to be specified.
•If MODE=2, only datasets 0, 1, and 30–39 have to be specified.
•If MODE=3, datasets 0–39 have to be specified 

Data sets 229 are used only if Mode = 1 or 3.

Data Set 2

[#Text]
Required comment line —“#” must be in column 1.

Data Set 3

[#Text]
Required comment line —“#” must be in column 1.

Data Set 4

NNODES NPIPES NLAYERS (read 1 line)
NNODES— is an integer value for the total number of nodes in the conduit pipe network. Each node is located at the center of a model cell in plan view.

NPIPES— is an integer value for the total number of pipes in the conduit network.

NLAYERS— is an integer value for the total number of model layers.



Data Set 5

[#Text]
Required comment line —“#” must be in column 1.

Data Set 6

TEMPERATURE—is a real number in degrees Celsius representing the average temperature of ground water in the conduit pipes.


Data Set 7

[#Text]
Required comment line —“#” must be in column 1.

Repeat data set 8 for each node.

Data Set 8

NO_N MC MR ML NB1 NB2 NB3 NB4 NB5 NB6 PB1 PB2 PB3 PB4 PB5 PB6
NO_N (column 1)—is the node number. Nodes are at the center of a model cell and define the connections of the pipe network.

MC, MR, and ML (columns 2 to 4)—are the MODFLOW cell column, row, and layer numbers within which node NO_N is located, respectively.

NB1 to NB6 (columns 5 to 10)—are neighbor node numbers, which are connected by pipes to node NO_N (specified in the first column). As many as six adjacent or diagonal neighboring nodes in three dimensions can be listed. Neighbor nodes must be in adjacent cells. If there are less than six neighbor nodes, insert zeros for the remaining columns. There will always be at least one node connected to the specified node NO_N.

PB1 to PB6 (columns 11 to 16)—are the pipe numbers connected to the node NO_N. Each node should be connected to at least 1 pipe, but can be connected to as many as 6 pipes. Pipe numbers represent the pipe segment between nodes. If there are less than six neighboring pipes connected to the node (NO_N), insert zeros in the remaining columns. Again, a zero in the rest of the columns means there are no more pipe connections.



Data Set 9

[#Text]
Required comment line —“#” must be in column 1.

Data Set 10

[#Text]
Required comment line —“#” must be in column 1.

Data Set 11

[#Text]
Required comment line —“#” must be in column 1.

Data Set 12

GEOHEIGHT—is the absolute elevation of the pipe nodes. There are two options for entering node elevations. With option 1 the number of lines read is equal to the total number of nodes (NNODES dataset 4). With option 2 only 1 line is read and the first number must equal NNODES in dataset 4.
Option 1—NO_N ELEVATION— Each line will contain an integer node number (NO_N) and the respective node elevation (ELEVATION) with respect to the user’s model datum. A CFPM1 simulation with 10 nodes, for example, would require 10 lines, each containing an integer node number and the respective node elevation. If the first node number is less than the total number of nodes, this indicates that option 1 has been selected and the CFP program will expect to read in a total number of lines of data equal to the total number of nodes as specified in dataset 4.
Option 2—NNODES ELEVATION— only one line is read describing the elevation of each pipe node within finitedifference cells. If the first integer in this section of input equals the total number of nodes as specified in dataset 4, then only one line will be read. This option allows for the CFPM1 to set each node elevation to the vertical center of the encompassing model cell, if ELEVATION is assigned a value of 0.0. Assigning values other than 0.0 for ELEVATION will raise or lower the node elevations a distance equal to ELEVATION above or below the vertical center of the encompassing model cell.


Data Set 13

[#Text]
Required comment line —“#” must be in column 1.

Data Set 14

SA_EXCHANGE—is an integer that equals either 0 or 1
If SA_EXCHANGE = 0, the user assigns the pipe conductance for each node in the CFP Input File.
If SA_EXCHANGE = 1, the user assigns the conduit wall permeability, and the CFP will compute the surface areas of pipes when assembling pipe conductances for groundwater exchange between pipes and porous media.


Data Set 15

[#Text]
Required comment line —“#” must be in column 1.

Data Set 16

EPSILON—is a real number for the convergence criterion of the NewtonRaphson iteration for pipe flow equations. Use a very small number, such as 0.0001.


Data Set 17

[#Text]
Required comment line —“#” must be in column 1.

Data Set 18

NITER—is an integer number for the maximum number of NewtonRaphson iterations. If convergence cannot be achieved, the program will stop and a warning will be printed in the MODFLOW listing file.


Data Set 19

[#Text]
Required comment line —“#” must be in column 1.

Data Set 20

RELAX—is a real number of relaxation that determines the step length of the NewtonRaphson iterations. Changing RELAX to a value slightly less than 1.0 may facilitate convergence of the pipe flow equations.


Data Set 21

[#Text]
Required comment line —“#” must be in column 1.

Data Set 22

P_NR—is an integer print flag for Newton Raphson iterations. If P_NR equals 0, results from these iterations are not printed. If P_NR equals 1, results from these iterations are printed, including the MODFLOW iteration number, Newton Raphson iteration number, and the maximum node head change.


Data Set 23

[#Text]
Required comment line —“#” must be in column 1.

Data Set 24

[#Text]
Required comment line —“#” must be in column 1.

Repeat data set 25 for each pipe.

Data Set 25

NO_P DIAMETER TORTUOSITY RHEIGHT LCRITREY_P TCRITREY_P
NO_P (column 1)—is the integer pipe number.

DIAMETER (column 2)—is a real number for the pipe diameter.

TORTUOSITY (column 3)—is a real number for the pipe tortuosity. A value of 1.0 indicates a straight pipe. If a pipe length greater than the straightline distance between the two nodes of the pipe is required, a tortuosity greater than 1 can be used.

RHEIGHT (column 4)—is a real number for the internal calculation of roughness and is the mean height of the microtopography of the conduit wall.

LCRITREY_P (column 5)—is the lower critical Reynolds number (turbulent to laminar).

TCRITREY_P (column 6)—is the upper critical Reynolds number (laminar to turbulent).



Data Set 26

[#Text]
Required comment line —“#” must be in column 1.

Repeat data set 27 for each node.

Data Set 27

NO_N N_HEAD
NO_N (column 1)—is the integer node number.

N_HEAD (column 2)—are piezometric heads for nodes with fixed head (positive values) or a flag (1 negative values less than or equal to 1) that indicates that heads will be calculated during simulation.



Data Set 28

[#Text]
Required comment line —“#” must be in column 1.

Repeat data set 29 for each node

Data Set 29

NO_N K_EXCHANGE
NO_N (column 1)—is the node number (in increasing order).



Data sets 3039 are used only if Mode = 2 or 3.

Data Set 30

[#Text]
Required comment line —“#” must be in column 1.
Two required text comment lines read only if MODE equals 2 or 3 from Dataset 1. The following variables are required in the CFP Input File for conduit layers. It is important to remember that when MODE equals 2, the following variables are the only variables needed in the CFP Input File. Also, conduit layers must be convertible layers. For example, conduit layers must be specified as LAYCON = 3 in the BCF Package, LAYTYPE not equal 0 in the LPF Package, or LTHUF not equal 0 in the HUF2 Package.

Data Set 31

[#Text]
Required comment line —“#” must be in column 1.

Data Set 32

NCL—is an integer equal to the total number of conduit layers.


Data Set 33

[#Text]
Required comment line —“#” must be in column 1.

Data Set 34

CL—is a onedimensional integer array entered on a single line of the CFP Input File. This array holds the MODFLOW layer numbers that are conduit layers.


Data Set 35

[#Text]
Required comment line —“#” must be in column 1.

Data Set 36

LTEMP—is the mean water temperature in degrees Celsius of all conduit layers.


Data Set 37

[#Text]
Required comment line —“#” must be in column 1.

Repeat lines 38 and 39 for each conduit flow layer.

Data Set 38

[#Text]
Required comment line —“#” must be in column 1.

Data Set 39

VOID LCRITREY_L TCRITREY_L
VOID—is the mean void diameter. This value is used in calculating the critical head difference from the specified lower and upper critical Reynolds numbers.

LCRITREY_L—is the lower critical Reynolds number for switching from turbulent to laminar flow.

TCRITREY_L—is the upper critical Reynolds number for transitioning from laminar to turbulent flow.


