DATA INPUT ERRORS

In order to facilitate its application, the SWIFT model performs a number of checks on the input data. If an error is detected, the code then converts to an error mode of operation. In this mode, the reading of data continues through the last recurrent data set. An error summary is then written, and the run is aborted. This appendix provides an explanation of the codes which appear in the summary.

Error Explanation Record

1 Increment numbers are incorrectly specified: M-3
NX < 2, NY < 1 or NZ < 1.

2 Number of aquifer influence function blocks (NABLMX) M-3
exceeds the total number of blocks (NB)

3 Compressibilities assume unphysical values: R1-1
CW < 0 or CR < 0;
or heat capacities are given unphysical values:
CPW < 0 or CPR < 0.

4 Thermal conductivities are given unphysical values: R1-2
UKTX < 0 , UKTY < 0 or UKTZ < 0;
dispersivities are incorrectly specified:
ALPHL < 0 or ALPHT < 0;
the diffusion is erroneously defined:
DMEFF < 0;
or the dimensions of ALPHL and ALPHT have
exceeded the dimensions (NRTMAXM):
NRT > NRTMAXM.

5 Fluid or rock densities are incorrectly specified: R1-3
BROCK ≤ 0, BWRN ≤ 0 or BWRI ≤ 0. 

6 Viscosities are given unphysical values: R1-7
VISRR ≤ 0, VISIR ≤ 0, VCC ≤ 0, VISR ≤ 0 or VISI ≤ 0; R1-8
brine concentrations are incorrectly specified: R1-9
SCC < 0 R1-10
the reference temperatures have been used a second time:
TR = TRR or TI = TIR;
or NCV > NCVMX, NTVR > NVRMX,
NTVI > NVIMX, NDT > NDTMX.

7 The heterogeneity/geometry control is outside M-3
its permitted range:
HTG < 1 or HTG > 3. 

8 Initialization output is incorrectly specified: M-3
KOUT ≠ 0, 1 or 3.

9 The print orientation control is outside its permitted range: M-3
PRT < -1 or PRT > 2.

10 For the over/underburden, either the thermal conductivity is
incorrectly stated: R1-13
KOB < 0 with NZOB > 2;
or KUB < 0 with NZUB > 2. 

11 Grid-block sizes are specified incorrectly: R1-17
DELX ≤ 0, DELY ≤ 0 or DELZ ≤ 0. R1-18,R1-19

12 Homogeneous reservoir properties are prescribed R1-20
incorrectly for the
incorrectly for the permeabilities:
KX < 0, KY < 0 or KZ < 0;
for the porosity:
PHI < 0;
or for the reservoir dip angles:
SINX < -1 or SINX > 1,
SINY < -1 or SINY > 1. 

13 Heterogeneous reservoir properties are prescribed R1-21
incorrectly for the grid-block indices:
I1 < 1 or I2 > NX,
J1 < 1 or J2 > NY,
K1 < 1 or K2 > NZ,
I1 > I2, J1 > J2 or K1 > K2;
for the permeabilities:
KX < 0, KY < 0 or KZ < 0;
for the porosity:
PHI < 0 or PHI > 1
or for the heat capacity of the saturated rock:
UCPR < 0. 

14 For the gridding of the radial coordinate, the radii R1-22
are incorrectly specified:
R1 ≥ RWW with R1 > 0,
R1 ≤ RE,
or ln (RE/R1)/(NX-0.5) will cause overflow.

15 For simulation in cylindrical coordinates, the value R1-23
of the vertical increment is incorrect:
DELZ ≤ 0;
the values of the permeabilities are incorrect:
KYY < 0 or KZZ < 0;
or the value of the porosity is incorrect:
POROS < 0 or POROS > 1. 

16 Radial block centers are erroneously prescribed: R1-25
RR(I) > RE for 1 # I # NX.

17 Block edges are incorrectly specified for the reservoir R1-26
reservoir modifications:
I1 < 1 or I1 > NX,
J1 < 1 or J1 > NY,
K1 < 1 or K2 > NZ,
I1 > I2, J1 > J2 or K1 > K2.

18 Aquifer-influence conditions are in error for the R1-27
option parameter: R1-28
IAQ < 0 or IAQ > 4;
or for the block edges:
I1 < 1 or I2 > NX,
J1 < 1 or J2 > NY,
K1 < 1 or K2 > NZ
I1 > I2, J1 > J2 or K1 > K2;
or NPT > NPTMX. 

19 The number of aquifer-influence functions (NABL) R1-28
summed from the input specifications exceeds that
allocated (NABLMX) in READ M-3 or
the number of points in the table of
influence function (NPT) exceeds that allocated (NPTMX).

20 The number of recharge blocks R2-2.5
(NRCH) summed from the input
specifications exceeds that allocated (NRCHMX)
in READ M-3.

21 For initialization of the brine I-2
concentrations, the block indices
are given incorrectly:
I1 < 1 or I2 > NX,
J1 < 1 or J2 > NY,
K1 < 1 or K2 > NZ
I1 > I2, J1 > J2 or K1 > K2;
or for the concentration:
CINIT < 0.

22 For initialization of the I-4
radionuclide concentrations, the
block indices are given incorrectly,
as specified for Error 21.

23 Some grid block, which has a nonzero pore volume R1-21
also has a zero sum of transmissibilities over its surface. R1-26

24 A grid-block pore volume is negative. R1-26

25 The number of wells is specified incorrectly: R2-4
NWT < 1 or NWT > NWMAX
where the latter parameter is defined internally within routine
MAIN. Currently the value of this quantity is:
NWMAX = 100.

26 A well index I is out of range: R2-6
I < 1 or I > NWT. 

27 Well locators IIW and IJW are outside the reservoir: R2-7
IIW < 1 or IIW > NX,
IJW < 1 or IJW > NY. 

28 The well-completion zone lies outside of the system: R2-7
IIC1 < 1 or IIC2 < 1,
IIC1 > NZ or IIC2 > NZ,
IIC1 > IIC2,
or the top block of the completion interval
has zero pore volume. 

30 The well-specification option is erroneously given: R2-7
IINDW1 ≠ 1, ±2, ±3, or ±4
or IINDW1 = 2, 3, 4 with NCALL
greater than 4.

32 Well index is incorrectly defined: R2-7
WI # 0 with IINDW1 … 1.

35 The completion interval contains R2-7
only blocks having zero pore volumes.

37 Layer-allocation factors are incorrectly prescribed: R2-7
KHL < 0. 

38 All layer-allocation factors, R2-7
KHL, are zero. At least one KHL factor must be nonzero.

39 A well identification number I is out of range: R2-7
I < 1 or I > NWT.

40 Differencing and matrix solution are improperly given: R2-2
METHOD < -2 or METHOD > 2
or WTFAC > 1. 

41 The number of nonlinear iterations is incorrect: R2-11
MINITN < 1 or MINITN > MAXITN.

43 Recurrent time is specified incorrectly R2-12
TCHG ≤ TIME
where TIME is the current time. 

44 The time step is: R2-12
DT = 0 for the first time step.

45 The maximum time step is erroneous: R2-12
DTMAX < DTMIN.

46 Mapping is specified incorrectly R2-13
in that one of the four digits of MAP is greater than unity.

47 Maps are prescribed incorrectly: R2-15
I1 < 1 or I2 > NX,
J1 < 1 or J2 > NY,
K1 < 1 or K2 > NZ,
I1 > I2, J1 > J2 or K1 > K2,
or for cylindrical geometry (HTG = 3, RAD M-3) J2 ≠1. 

48 The number of nuclide monitor blocks encountered R2-16
exceeds that allocated (NDSCMX).

50 Sink/source information contains error. The number R2-9
of blocks containing sinks is greater than that allocated: R2-10
NSS > NSMAX;
the sink number is too large:
I > NSS;
or the block indices are incorrect:
IIS < 1 or IIS > NX,
IJS < 1 or IJS > NY,
IKS < 1 or IKS > NZ. 

51 Radioactive-component specification is incorrect. R0-1
For the component number,
I = 0 or I > NCP,
or for the number of parents,
NP < 0. 

52 The half-life is physically unreasonable: R0-1
DEC < 0.

53 Branching, as specified, is incorrect. For the R0-1
parent index,
KP < 1 or KP > NCP,
or for the branching ratio,
AP < 0. 

54 The Freundlich k factor is given incorrectly: R0-2
DIS(I) < 0, I ≤ NCP.

55 For the rock-type modifications, the block numbers R1A-1
are out of range:
I1A < 1 or I1B > NX,
J1A < 1 or J1B > NY,
K1A < 1 or K1B > NZ.
I1A > I1B, J1A > J1B or K1A > K1B.

56 The salt-dissolution coefficient is prescribed R1A-2
erroneously:
ACS < 0.

57 For canistered storage (ILEVEL = 1, READ R1A-3), the R1A-4
spacing parameters are given incorrectly:
SDRIFT < 0, SCNSTR < 0,
DCNSTR < 0, or HCNSTR < 0.

58 The specified location of the repository is unreasonable. R1A-6
A block number is out of range:
I1A < 1 or I1B > NX,
J1A < 1 or J1B > NY,
K1A < 1 or K1B > NZ;
I1A > I1B, J1A > J1B or K1A > K1B.
or the allocated number of repository blocks,
NREPB, is exceeded. 

59 For the concentration of the unleached wastes R1A-7
or the repository heat load, the interpolation
table is given incorrectly:
CTIME ≤ 0,
or the values are not in ascending order. 

60 For canistered storage (ILEVEL = 1 in - R1A-4
READ R1A-3), the storage parameter values are erroneous:
SDRIFT > Max (DELX, DELY),
SCNSTR > Max (DELX, DELY),
or HCNSTR > DELZ. 

61 Solubilities are given incorrectly: R1A-9
CS < 0.

62 The dual-porosity blocks are not ROD-3
specified correctly.
Block indices are in error:
I1A < 1 or I1B > NX,
J1A < 1 or J1B > NY,
K1A < 1 or K1B > NZ;
I1A > I1B, J1A > J1B or K1A > K1B;
the local rock-type indicator is outside its permitted range:
IR > NRTD;
or the position control for the local unit can not be interpreted:
IFD < -3 or IFD > 3. 

63 The properties of the local units have unphysical values. R1D-1
For the rock compressibility: R1D-2
CRD < 0;
for auto mesh generation the length of the local units:
SAD ≤ 0 or DSD ≤ 0;
for double mesh generation:
DSDO < 0;
for user defined mesh:
DS(I) ≤ 0;
for permeability and porosity:
AKSD < 0 or PHID < 0;
for the dispersivity:
ALPD < 0;
or for the thermal conductivity
and heat capacity:
UKTD < 0 or CPRKD < 0.

64 Distribution coefficients are incorrect: R1D-4
D1SD < 0.

65 Mesh generation is prescribed ROD-2
incorrectly. For single mesh-
generation (outer increment
DSDO = 0), thee is an insufficient R1D-2
number of nodes (KGRD = 1):
NSD ≤ 3;
or for double-mesh generation
(DSDO > 0), there is an insufficient
number of nodes (KGRD = 2): NSD ≤ 4.

66 Grid block volume is less than R1-26
pore volume in a dual porosity block.
Avoid pore volume modification
with dual porosity blocks.

67 Modifications to dual porosity blocks are not. R1D-5
specified correctly Block indices are in error:
I1 < 1 or I1 > NX,
J1 < 1 or J1 > NY,
K1 < 1 or K1 > NZ,
I1 > I1, J1 > J2, or K1 > K2.