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 two sets of records -- one for moisture-content parameters and the other for conductivity (permeability) parameters. Therefore, each set contains NMAT records and the total number of records for this subset is 2xNMAT. The number of lines in each record is determined by the number of parameters used to specify the soil properties per material, NSPPM. This pattern must be repeated NDF times to cover all pore regions. For pore region K, the I-th (I = 1, 2, ..., NMAT) record, in free format, contains the following variables:

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.

For pore region K, the (NMAT + I)-th (I = 1, 2, ..., NMAT) record, in free format, contains the following variables.

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.

Subset B2. Soil Properties in Tabular Form: This subset is needed if and only if KSP .NE. 0. Four sets of records are needed for pressure, water-content, relative conductivity (or relative permeability), and water capacity. Therefore, each set contains NMAT records, and the total number of subsets is 4*NMAT. This pattern must be repeated NDF times to cover all pore regions. The number of lines in each record is determined by the input parameters, NSPPM. For pore region K, the I-th (I = 1, 2, ..., NMAT) record, in free format, contains the following variables:

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 pore region K, the (NMAT + I)-th (I = 1, 2, ..., NMAT) record, in free format, contains the following variables:

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 pore region K, the (2*NMAT + I)-th (I = 1, 2, ..., NMAT) record, in free format, contains the following variables:

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.

For pore region K, the (3*NMAT + I)-th (I = 1, 2, ..., NMAT) record, in FREE FORMAT, contains the following variables:

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.