Saturation Function Card Syntax
~Saturation Function Card
Choose One Saturation Functionvan Genuchten
Choose One Sub-Optionvan Genuchten
Repeat for each Rock/Soil Zonation Domain
,
van Genuchten,
α Parameter,
Units
(1/m),
n Parameter,
Minimum Saturation,
m Parameter(Optional),
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
van Genuchten with Entrapment
Repeat for each Rock/Soil Zonation Domain,
van Genuchten w/ Entrapment,
α Parameter,
Units(1/m),
n Parameter,
Minimum Saturation,
m Parameter,
Effective Gas Residual Saturation,
Critical Trapping Number,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
van Genuchten with Webb Extension
Repeat for each Rock/Soil Zonation Domain,
van Genuchten w/ Webb,
α Parameter,
Units(1/m),
n Parameter,
Minimum Saturation,
m Parameter,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
van Genuchten with Fractured Systems
InfoIf the key words 'dp,' 'dual porosity,' or 'fractured' appear in the rock/soil name specified in the Rock/Soil Zonation Card, this indicates a dual porosity medium and the reading of both matrix and fracture properties (e.g., Fracture and Matrix Porosity) is triggered.
Repeat for each Rock/Soil Zonation Domain,
Fractured van Genuchten,
Matrix α Parameter,
Units(1/m),
Matrix n Parameter,
Matrix Minimum Saturation,
Fracture α Parameter,
Units(1/m),
Fracture n Parameter,
Fracture Minimum Saturation,
(Optional) Matrix m Parameter,
(Optional) Fracture m Parameter,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
Brooks and Corey
Choose One Sub-OptionBrooks and Corey
Repeat for each Rock/Soil Zonation Domain,
Brooks and Corey,
Entry head,
Units(m),
Lambda Parameter,
Minimum Saturation,
Effective Gas Residual Saturation,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
Brooks and Corey with Entrapment
Repeat for each Rock/Soil Zonation Domain,
Brooks and Corey w/ Entrapment,
Entry head,
Units(m),
Lambda Parameter,
Minimum Saturation,
Effective Gas Residual Saturation,
Critical Trapping Number,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
Brooks and Corey with Webb Extension
Repeat for each Rock/Soil Zonation Domain,
Brooks and Corey w/ Webb,
Entry head,
Units(m),
Lambda Parameter,
Minimum Saturation,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
Brooks and Corey with Fractured Systems
InfoIf the key words 'dp,' 'dual porosity,' or 'fractured' appear in the rock/soil name specified in the Rock/Soil Zonation Card, this indicates a dual porosity medium and the reading of both matrix and fracture properties (e.g., Fracture and Matrix Porosity) is triggered.
Repeat for each Rock/Soil Zonation Domain,
Fractured Brooks and Corey,
Matrix Entry head,
Units(m),
Matrix Lambda Parameter,
Matrix Minimum Saturation,
Fracture Entry head,
Units(m),
Fracture Lambda Parameter,
Fracture Minimum Saturation,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
Haverkamp
Repeat for each Rock/Soil Zonation Domain,
Haverkamp,
Entry head,
Units(m),
α Parameter,
β Parameter,
Minimum Saturation,
IJK IndexingIf the Rock/Soil Name is "IJK Indexing," "JKI Indexing," or "KIJ Indexing" the value specified will be applied to all nodes in the domain. Alternately, any parameter value can be replaced with an external file using the following syntax in the place of the parameter value:
,file:,filename for parameter,
or for binary files:
,binary file:,filename for parameter,
Units shown in the input line will be applied to all parameters in the external file. Values and files may be used in combination within an input line.
Tabular
Choose One Tabular OptionCapillary Head/Aqueous Saturation
Choose One Interpolation SchemeLinear
Repeat for each Rock/Soil Zonation Domain,
Tabular,
Number of Table Entries
, Repeat for number of table entries,
Units(m),
Aqueous Saturation, Log-Linear
Repeat for each Rock/Soil Zonation Domain,
Tabular Log-Linear,
Number of Table Entries, Repeat for number of table entries,
Units(m),
Aqueous Saturation, Tabular Cubic Spline
Repeat for each Rock/Soil Zonation Domain,
Tabular Spline,
Number of Table Entries, Repeat for number of table entries,
Units(m),
Aqueous Saturation, Tabular Cubic Spline Log
Repeat for each Rock/Soil Zonation Domain,
Tabular Spline Log,
Number of Table Entries, Repeat for number of table entries,
Units(m),
Aqueous Saturation, Capillary Head/Water Content
Choose One Interpolation SchemeLinear
Repeat for each Rock/Soil Zonation Domain,
Tabular,
Number of Table Entries, Repeat for number of table entries,
Units(m),
Water Content, Log-Linear
Repeat for each Rock/Soil Zonation Domain,
Tabular Log-Linear,
Number of Table Entries, Repeat for number of table entries,
Units(m),
Water Content, Tabular Cubic Spline
Repeat for each Rock/Soil Zonation Domain,
Tabular Spline,
Number of Table Entries, Repeat for number of table entries,
Units(m),
Water Content, Tabular Cubic Spline Log
Repeat for each Rock/Soil Zonation Domain,
Tabular Spline Log,
Number of Table Entries, Repeat for number of table entries,
Units(m),
Water Content, 
TipNote that the Entrapped Gas sub-option can be combined with the Webb Extension sub-option such that the syntax is "van Gnuchten w/Entrapment w/ Webb" or "Brooks and Corey w/Entrapment w/ Webb"
External FilesFile names must be all lower case. Files must contain text data unless specified as binary in which case they must contain binary data created with the same FORTRAN compiler used to compile STOMP. Files contain parameter values for every grid-cell, ordered according to the indexing scheme (i.e., IJK, JKI, or KIJ).