DATA SET 6: SOIL PROPERTY PARAMETERS

This data set characterizes water retention curves, water capacity functions, conductivity/permeability - water content/pressure head relations, and mass exchange coefficients. It contains three subsets of data. The subsets B1 and B2 contain two and four sequences of data records, and only one of these two subsets is required depending on the given forms of the above-mentioned relations, namely analytical or tabular.

Subset A. Soil Property Control Integers: This record is in free format and contains the following variables.

1. KSP = Soil property input control, 0 = analytical input, 1 = Tabular data input.
2. NSPPM = Number of points in tabular soil property functions or number of parameters to specify analytical soil functions per material.
3. KCP = Permeability input control; 0 = input saturated hydraulic conductivity, 1 = input saturated permeability.
4. KSMOD = Analytical soil property function used in subroutine SPROP.

Subset B1. Analytical Soil Parameters: This sub set is needed if and only if KSP = 0. For each pore region of each material type, two sets of records -- one for moisture-content parameters and the other for conductivity (permeability) parameters are required. For pore region K and material type I, one needs to enter the following two sets of data: for moisture content:

1. THPROP(1,I,K) = Analytical moisture-content parameter 1 of material I.

2. THPROP(2,I,K) = Analytical moisture-content parameter 2 of material I.

NSPPM. THPROP(NSPPM,I,K) = Analytical moisture-content parameter NSPPM of material I.

and for relative hydraulic conductivity:

1. AKPROP(1,I,K) = Analytical relative conductivity parameter 1 of material I.

2. AKPROP(2,I,K) = Analytical relative conductivity parameter 2 of material I.

NSPPM. AKPROP(NSPPM,I,K) = Analytical relative conductivity parameter NSPPM of material I.

These two sets of data are repeated NMAT times for pore region K. Once the pore region is completely configured, one then advances to pore region K+1 and so on.

Subset B2. Soil Properties in Tabular Form: This subset is needed if and only if KSP .NE. 0. For each material, four sets of records are needed for pressure, water-content, relative conductivity (or relative permeability), and water capacity. The records should be repeated NMAT times for each material type. This pattern must then be repeated NDF times to cover all pore regions. The number of lines in each record is determined by the input parameters, NSPPM. For material type I and pore region K, we there fore have for pressure head:

1. HPROP(1,I,K) = Tabular value of pressure head of the first point in material I.

2. HPROP(2,I,K) = Tabular value of pressure head of the second point in material I.

NSPPM. THPROP(NSPPM,I,K) = Tabular value of pressure head of the NSPPM-th point in material I.

for water content:

1. THPROP(1,I,K) = Tabular value of moisture content of the first point in material I.

2. THPROP(2,I,K) = Tabular value of moisture content of the second point in material I.

NSPPM. THPROP(NSPPM,I,K) = Tabular value of moisture content of the NSPPM-th point in material I.

for relative hydraulic conductivity:

1. AKPROP(1,I,K) = Tabular value of relative conductivity of the first point in material I.

2. AKPROP(2,I,K) = Tabular value of relative conductivity of second point in material I.

NSPPM. KPROP(NSPPM,I,K) = Tabular value of relative conductivity of the NSPPM-th point in material I.

and for moisture-content capacity:

1. CAPROP(1,I,K) = Tabular value of moisture-content capacity of the first point in material I.

2. CAPROP(2,I,K) = Tabular value of moisture-content capacity of the second point in material I.

NSPPM. CAPROP(NSPPM,I,K) = Tabular value of moisture content capacity of the NSPPM-th point in material I.

Subset C. Mass transfer coefficients: For each material type, the mass transfer coefficients between every two pore domains need to be specified.

1. EXCF(I,J,K) = The exchange coefficient between regions J and K in material I (L/T), for J = 1 to NDF and K = 1, to NDF.

Note that theoretically there is no exchange effects within the same pore region, that is, when J=K, the exchange coefficients should be set to zero.

Return to MurfMd home page.
Back to data set 5: Material property.
On to data set 7: Nodal coordinate.