pyMelt.lithologies package

Submodules

pyMelt.lithologies.katz module

Katz et al. (2003)

Implementation of the anhydrous melting model presented by Katz et al. (2003).

class pyMelt.lithologies.katz.lherzolite(CP=1000.0, alphas=40.0, alphaf=68.0, rhos=3.3, rhof=2.9, DeltaS=300.0, parameters={'A1': 1085.7, 'A2': 132.9, 'A3': -5.1, 'B1': 1475.0, 'B2': 80.0, 'B3': -3.2, 'C1': 1780.0, 'C2': 45.0, 'C3': -2.0, 'Mcpx': 0.17, 'beta1': 1.5, 'beta2': 1.5, 'r1': 0.5, 'r2': 0.08})[source]

Bases: lithology

Implementation of the Katz et al. (2003) anhydrous lherzolite melting model.

To use the same format of parameterisation for another lithology, the parameter values may be changed. They are provided as a dictionary during initialisation of the object.

Mcpx: Mass fraction of cpx in the source. Controls the transition to low-productivity harzburgite-type melting.
A1: Parameter used to define solidus.
A2: Parameter used to define solidus.
A3: Parameter used to define solidus.
B1: Parameter used to define lherzolite-liquidus.
B2: Parameter used to define lherzolite-liquidus.
B3: Parameter used to define lherzolite-liquidus.
C1: Parameter used to define liquidus.
C2: Parameter used to define liquidus.
C3: Parameter used to define liquidus.
beta1: Parameter used to calculate melt fraction during cpx-present melting.
beta2: Parameter used to calculate melt fraction during cpx-absent melting.
r1: Parameter used to define cpx reaction coefficient.
r2: Parameter used to define cpx reaction coefficient.

The thermal expansivities, the heat capacity, the densities, and the entropy of fusion may also be changed during object initialisation.

Parameters:

CP (float, default: pyMelt.lithology_class.default_properties[‘CP’]) – The heat capacity (J K-1 kg-1)
alphas (float, default: pyMelt.lithology_class.default_properties[‘alphas’]) – The thermal expansivity of the solid (1e-6 K-1)
alphaf (float, default: pyMelt.lithology_class.default_properties[‘alphaf’]) – The thermal expansivity of the melt (1e-6 K-1)
rhos (float, default: pyMelt.lithology_class.default_properties[‘rhos’]) – The density of the solid (kg m-3)
rhof (float, default: pyMelt.lithology_class.default_properties[‘rhof’]) – The density of the melt (kg m-3)
DeltaS (float, default: pyMelt.lithology_class.default_properties[‘DeltaS’]) – The entropy of fusion J K-1 kg-1
parameters (dict, default: parameters from Katz et al. (2003)) – The model parameters described above

F(P, T, **kwargs)[source]

Wrapper for the melt fraction functions. If T and P are below the solidus, returns 0, if they are above the liquidus, returns 1. If below the temperature of cpx-exhaustion, calls the Fcpx function, otherwise calls the Fopx function.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

Melt fraction.

Return type:

float

TLiquidus(P, **kwargs)[source]

Returns the temperature of the liquidus at any given pressure. Eqn(10).

Parameters:: P (float) – Pressure (GPa).
Returns:: Liquidus temperature (degC).
Return type:: float

TSolidus(P, **kwargs)[source]

Returns the temperature of the solidus at any given pressure. Eqn(4).

Parameters:: P (float) – Pressure (GPa).
Returns:: Solidus temperature (degC).
Return type:: float

dTdF(P, T, **kwargs)[source]

Calculates dT/dF(const. P). First calculates the melt fraction. If F is zero, returns _np.inf. If F is 1, returns _np.inf. Otherwise uses the appropriate expressions for cpx present or absent melting.

Parameters:

P (float) – Pressure (GPa)
T (float) – Temperature (degC)

Returns:

dT/dF(const. P) (K).

Return type:

float

dTdP(P, T, **kwargs)[source]

Calculates dT/dP(const. F). First calculates F, then chooses the appropriate expression for cpx present or absent melting.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

dT/dP(const. F) (K GPa-1).

Return type:

float

pyMelt.lithologies.matthews module

Matthews et al. (2021)

Implentation of the melting models developed by Matthews et al. (2021).

class pyMelt.lithologies.matthews.eclogite(CP=1000.0, alphas=40.0, alphaf=68.0, rhos=3.3, rhof=2.9, DeltaS=300.0, parameters={'C1': 533.842, 'C2': 4.921, 'C3': 20.148, 'C4': 80.879, 'D1': 994.149, 'D2': 8.092, 'D3': -11.778, 'D4': -862.641, 'beta': 2.134})[source]

Bases: lithology

Implementation of the silica-saturated pyroxenite (or eclogite) melting model from Matthews et al. (2021).

To use the same format of parameterisation for another lithology, the parameter values may be changed. They are provided as a dictionary during initialisation of the object.

C1: Parameter used in solidus definition.
C2: Parameter used in solidus definition.
C3: Parameter used in solidus definition.
C4: Parameter used in solidus definition.
D1: Parameter used in liquidus definition.
D2: Parameter used in liquidus definition.
D3: Parameter used in liquidus definition.
D4: Parameter used in liquidus definition.
beta: Parameter used in melt fraction definition.

The thermal expansivities, the heat capacity, the densities, and the entropy of fusion may also be changed during object initialisation.

Parameters:

CP (float, default: pyMelt.default_properties[‘CP’]) – The heat capacity (J K-1 kg-1)
alphas (float, default: pyMelt.default_properties[‘alphas’]) – The thermal expansivity of the solid (1e-6 K-1)
alphaf (float, default: pyMelt.default_properties[‘alphaf’]) – The thermal expansivity of the melt (1e-6 K-1)
rhos (float, default: pyMelt.default_properties[‘rhos’]) – The density of the solid (kg m-3)
rhof (float, default: pyMelt.default_properties[‘rhof’]) – The density of the melt (kg m-3)
DeltaS (float, default: pyMelt.default_properties[‘DeltaS’]) – The entropy of fusion J K-1 kg-1
parameters (dict, default: parameters from Matthews et al. (2021)) – The model parameters described above

F(P, T, **kwargs)[source]

Calculates melt fraction at a given pressure and temperature using T’**beta, where T is the normalised temperature: (T-Tsolidus)/(T-Tliquidus). If P and T are below the the solidus, 0 is returned, if they are above the liquidus, 1 is returned.

Parameters:

P (float) – Pressure (GPa)
T (float) – Temperature (degC)

Returns:

Melt fraction.

Return type:

float

TLiquidus(P, **kwargs)[source]

Calculates the liquidus temperature at a given pressure.

Parameters:: P (float) – Pressure (GPa)
Returns:: Liquidus temperature (degC).
Return type:: float

TSolidus(P, **kwargs)[source]

Calculates the solidus temperature at a given pressure.

Parameters:: P (float) – Pressure (GPa).
Returns:: Solidus temperature (degC).
Return type:: float

dTdF(P, T, **kwargs)[source]

Calculates dT/dF(const. P) at a given pressure and temperature.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

dT/dF(const. P) (K)

Return type:

float

dTdP(P, T, **kwargs)[source]

Calculates dT/dP(const. F) at a given pressure and temperature.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

dTdP(const. F) (K GPa-1)

Return type:

float

class pyMelt.lithologies.matthews.kg1(CP=1000.0, alphas=40.0, alphaf=68.0, rhos=3.3, rhof=2.9, DeltaS=300.0, parameters={'A1': 450.0, 'A2': 2.098, 'A3': 17.0, 'A4': 623.828, 'B1': 174.566, 'B2': 336.833, 'B3': 66.762, 'B4': 503.101, 'C': 0.506, 'Mcpx': 0.342, 'beta1': 1.382, 'beta2': 1.8, 'r1': 0.342, 'r2': 0.191})[source]

Bases: lithology

Implementation of the KG1 melting model from Matthews et al. (2021).

To use the same format of parameterisation for another lithology, the parameter values may be changed. They are provided as a dictionary during initialisation of the object.

Mcpx: Mass fraction of cpx in the source. Controls the transition to low-productivity harzburgite-type melting.
A1: Parameter used to define solidus.
A2: Parameter used to define solidus.
A3: Parameter used to define solidus.
A4: Parameter used to define solidus.
B1: Parameter used to define liquidus.
B2: Parameter used to define liquidus.
B3: Parameter used to define liquidus.
B4: Parameter used to define liquidus.
C: Parameter used to define lherzolite-liquidus.
beta1: Parameter used to calculate melt fraction during cpx-present melting.
beta2: Parameter used to calculate melt fraction during cpx-absent melting.
r1: Parameter used to define cpx reaction coefficient.
r2: Parameter used to define cpx reaction coefficient.

The thermal expansivities, the heat capacity, the densities, and the entropy of fusion may also be changed during object initialisation.

Parameters:

CP (float, default: pyMelt.default_properties[‘CP’]) – The heat capacity (J K-1 kg-1)
alphas (float, default: pyMelt.default_properties[‘alphas’]) – The thermal expansivity of the solid (1e-6 K-1)
alphaf (float, default: pyMelt.default_properties[‘alphaf’]) – The thermal expansivity of the melt (1e-6 K-1)
rhos (float, default: pyMelt.default_properties[‘rhos’]) – The density of the solid (kg m-3)
rhof (float, default: pyMelt.default_properties[‘rhof’]) – The density of the melt (kg m-3)
DeltaS (float, default: pyMelt.default_properties[‘DeltaS’]) – The entropy of fusion J K-1 kg-1
parameters (dict, default: parameters from Matthews et al. (2021)) – The model parameters described above

F(P, T, **kwargs)[source]

Wrapper for the melt fraction functions. If T and P are below the solidus, returns 0, if they are above the liquidus, returns 1. If below the temperature of cpx-exhaustion, calls the Fcpx function, otherwise calls the Fopx function.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

Melt fraction.

Return type:

float

TLiquidus(P, **kwargs)[source]

Returns the temperature of the liquidus at any given pressure.

Parameters:: P (float) – Pressure (GPa).
Returns:: Liquidus temperature (degC).
Return type:: float

TSolidus(P, **kwargs)[source]

Returns the temperature of the solidus at any given pressure.

Parameters:: P (float) – Pressure (GPa).
Returns:: Solidus temperature (degC).
Return type:: float

dTdF(P, T, **kwargs)[source]

Calculates dT/dF(const. P). First calculates the melt fraction. If F is zero, returns _np.inf. If F is 1, returns _np.inf. Otherwise uses the appropriate expressions for cpx present or absent melting.

Parameters:

P (float) – Pressure (GPa)
T (float) – Temperature (degC)

Returns:

dT/dF(const. P) (K).

Return type:

float

dTdP(P, T, **kwargs)[source]

Calculates dT/dP(const. F). First calculates F, then chooses the appropriate expression for cpx present or absent melting.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

dT/dP(const. F) (K GPa-1).

Return type:

float

class pyMelt.lithologies.matthews.klb1(CP=1000.0, alphas=40.0, alphaf=68.0, rhos=3.3, rhof=2.9, DeltaS=300.0, parameters={'A1': 2445.754, 'A2': 9.511, 'A3': -99.782, 'A4': -4378.581, 'B1': 480.403, 'B2': 672.391, 'B3': 12.275, 'B4': -1242.536, 'C': 0.6873, 'Mcpx': 0.15, 'beta1': 1.5, 'beta2': 1.5, 'r1': 0.5, 'r2': 0.08})[source]

Bases: object

Implementation of the KLB1 melting model from Matthews et al. (2021).

To use the same format of parameterisation for another lithology, the parameter values may be changed. They are provided as a dictionary during initialisation of the object.

Mcpx: Mass fraction of cpx in the source. Controls the transition to low-productivity harzburgite-type melting.
A1: Parameter used to define solidus.
A2: Parameter used to define solidus.
A3: Parameter used to define solidus.
A4: Parameter used to define solidus.
B1: Parameter used to define liquidus.
B2: Parameter used to define liquidus.
B3: Parameter used to define liquidus.
B4: Parameter used to define liquidus.
C: Parameter used to define lherzolite-liquidus.
beta1: Parameter used to calculate melt fraction during cpx-present melting.
beta2: Parameter used to calculate melt fraction during cpx-absent melting.
r1: Parameter used to define cpx reaction coefficient.
r2: Parameter used to define cpx reaction coefficient.

The thermal expansivities, the heat capacity, the densities, and the entropy of fusion may also be changed during object initialisation.

Parameters:

CP (float, default: pyMelt.default_properties[‘CP’]) – The heat capacity (J K-1 kg-1)
alphas (float, default: pyMelt.default_properties[‘alphas’]) – The thermal expansivity of the solid (1e-6 K-1)
alphaf (float, default: pyMelt.default_properties[‘alphaf’]) – The thermal expansivity of the melt (1e-6 K-1)
rhos (float, default: pyMelt.default_properties[‘rhos’]) – The density of the solid (kg m-3)
rhof (float, default: pyMelt.default_properties[‘rhof’]) – The density of the melt (kg m-3)
DeltaS (float, default: pyMelt.default_properties[‘DeltaS’]) – The entropy of fusion J K-1 kg-1
parameters (dict, default: parameters from Matthews et al. (2021))

F(P, T, **kwargs)[source]

Wrapper for the melt fraction functions. If T and P are below the solidus, returns 0, if they are above the liquidus, returns 1. If below the temperature of cpx-exhaustion, calls the Fcpx function, otherwise calls the Fopx function.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

Melt fraction.

Return type:

float

TLiquidus(P, **kwargs)[source]

Returns the temperature of the liquidus at any given pressure.

Parameters:: P (float) – Pressure (GPa).
Returns:: Liquidus temperature (degC).
Return type:: float

TSolidus(P, **kwargs)[source]

Returns the temperature of the solidus at any given pressure.

Parameters:: P (float) – Pressure (GPa).
Returns:: Solidus temperature (degC).
Return type:: float

dTdF(P, T, **kwargs)[source]

Calculates dT/dF(const. P). First calculates the melt fraction. If F is zero, returns _np.inf. If F is 1, returns _np.inf. Otherwise uses the appropriate expressions for cpx present or absent melting.

Parameters:

P (float) – Pressure (GPa)
T (float) – Temperature (degC)

Returns:

dT/dF(const. P) (K).

Return type:

float

dTdP(P, T, **kwargs)[source]

Calculates dT/dP(const. F). First calculates F, then chooses the appropriate expression for cpx present or absent melting.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

dT/dP(const. F) (K GPa-1).

Return type:

float

pyMelt.lithologies.pertermann module

Pertermann and Hirschmann (2003, JGR)

The pertermann module implements the G2 model.

class pyMelt.lithologies.pertermann.g2(CP=1000.0, alphas=40.0, alphaf=68.0, rhos=3.3, rhof=2.9, DeltaS=300.0, parameters={'a': 0.7368, 'b': 0.2632, 'c': 1175.0, 'd': 114.0, 'e': 920.0, 'f': 130.0})[source]

Bases: lithology

Implementation of the Pertermann and Hirschmann (2003, JGR) G2 melting model. To use the same format of parameterisation for another lithology, the parameter values may be changed. They are provided as a dictionary during initialisation of the object, with values:

a: Parameter used in calculating melt fraction.
b: Parameter used in calculating melt fraction.
c: Parameter used in liquidus definition.
d: Parameter used in liquidus definition.
e: Parameter used in solidus definition.
f: Parameter used in solidus definition.

The thermal expansivities, the heat capacity, the densities, and the entropy of fusion may also be changed during object initialisation.

Parameters:

CP (float, default: pyMelt.lithology_class.default_properties[‘CP’]) – The heat capacity (J K-1 kg-1)
alphas (float, default: pyMelt.lithology_class.default_properties[‘alphas’]) – The thermal expansivity of the solid (1e-6 K-1)
alphaf (float, default: pyMelt.lithology_class.default_properties[‘alphaf’]) – The thermal expansivity of the melt (1e-6 K-1)
rhos (float, default: pyMelt.lithology_class.default_properties[‘rhos’]) – The density of the solid (kg m-3)
rhof (float, default: pyMelt.lithology_class.default_properties[‘rhof’]) – The density of the melt (kg m-3)
DeltaS (float, default: pyMelt.lithology_class.default_properties[‘DeltaS’]) – The entropy of fusion J K-1 kg-1
parameters (dict, default: parameters from Pertermann and Hirschmann (2003, JGR)) – The model parameters described above

F(P, T, **kwargs)[source]

Calculates melt fraction at a given pressure and temperature using Equation 1:: a*T’**2 + b*T’,
where T is the normalised temperature (Equation 2):: (T-Tsolidus)/(Tliquidus-Tsolidus).

If P and T are below the the solidus, 0 is returned, if they are above the liquidus, 1 is returned.

Parameters:

P (float) – Pressure (GPa)
T (float) – Temperature (degC)

Returns:

Melt fraction.

Return type:

float

TLiquidus(P, **kwargs)[source]

Calculates the liquidus temperature, at a given pressure, using Equation 3:: c + d*P.

Parameters:: P (float) – Pressure (GPa).
Returns:: Liquidus temperature (degC).
Return type:: float

TSolidus(P, **kwargs)[source]

Calculates the solidus temperature, at a given pressure, using Equation 4:: e + f*P.

Parameters:: P (float) – Pressure (GPa)
Returns:: Solidus temperature (degC).
Return type:: float

dTdF(P, T, **kwargs)[source]

Calculates dT/dF(const. P) at a given pressure and temperature.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

dT/dF(const. P) (K)

Return type:

float

dTdP(P, T, **kwargs)[source]

Calculates dT/dP(const. F) at a given pressure and temperature.

Parameters:

P (float) – Pressure (GPa).
T (float) – Temperature (degC).

Returns:

dTdP(const. F) (K GPa-1)

Return type:

float

pyMelt.lithologies.shorttle module

Shorttle et al. (2014)

Implementation of the new lithologies in Shorttle et al. (2014).

class pyMelt.lithologies.shorttle.harzburgite(CP=1000.0, alphas=30.0, alphaf=30.0, rhos=3.25, rhof=2.9, DeltaS=300, parameters={})[source]

Bases: lithology

Material that does not melt, i.e. Harzburgite in Shorttle et al. (2014) and Matthews et al. (2016). Default constants as defined by Shorttle et al. (2014).

The thermal expansivities, the heat capacity, the densities, and the entropy of fusion may be changed during object initialisation.

Parameters:

CP (float, default: 1000.0) – The heat capacity (J K-1 kg-1)
alphas (float, default: 30.0) – The thermal expansivity of the solid (1e-6 K-1)
alphaf (float, default: 30.0) – Melt thermal expansivity, not used, here for consistency.
rhos (float, default: 3.25) – The density of the solid (kg m-3)
rhof (float, default: 2.9) – Melt density, not used, here for consistency.
DeltaS (float, default: 300) – The entropy of fusion, not used, here for consistency.
parameters (dict, default: {}) – This model does not use any parameters, here for consistency.

F(P, T, **kwargs)[source]

Melt Fraction. Returns 0.0.

Parameters:

P (float) – Pressure. There to maintain consistancy within lithology definitions.
T (float) – Temperature. There to maintain consistancy within lithology definitions.

Returns:

Melt fraction will always be 0.0.

Return type:

float

TLiquidus(P, **kwargs)[source]

Liquidus temperature. Returns _np.inf

Parameters:: P (float) – Pressure. There to maintain consistancy within lithology definitions.
Returns:: The value will always be infinite.
Return type:: numpy.inf

TSolidus(P, **kwargs)[source]

Solidus temperature. Returns _np.inf.

Parameters:: P (float) – Pressure. There to maintain consistancy within lithology definitions.
Returns:: The value will always be infinite.
Return type:: numpy.inf

dTdF(P, T, **kwargs)[source]

dTdF(constP). Returns _np.inf.

Parameters:

P (float) – Pressure. There to maintain consistancy within lithology definitions.
T (float) – Temperature. There to maintain consistancy within lithology definitions.

Returns:

The value will always be infinite.

Return type:

numpy.inf

dTdP(P, T, **kwargs)[source]

dTdP(constF). Returns 0.0.

Parameters:

P (float) – Pressure. There to maintain consistancy within lithology definitions.
T (float) – Temperature. There to maintain consistancy within lithology definitions.

Returns:

The value will always be 0.0.

Return type:

float

class pyMelt.lithologies.shorttle.kg1(CP=1140.0, alphas=30.0, alphaf=68.0, rhos=3.3, rhof=2.9, DeltaS=380.0, parameters={'A1': 1095.4, 'A2': 124.1, 'A3': -4.7, 'B1': 1179.6, 'B2': 157.2, 'B3': -11.1, 'C1': 1780.0, 'C2': 45.0, 'C3': -2.0, 'a': 0.3187864, 'alpha': 2.0, 'b': 0.4154, 'beta': 1.5, 'c': 0.7341864, 'd': 0.2658136})[source]

Bases: lithology

Implementation of the KG1 parameterisation by Shorttle et al. (2014).

To use the same format of parameterisation for another lithology, the parameter values may be changed. They are provided as a dictionary during initialisation of the object, with values:

A1: Constant used in solidus expression.
A2: Constant used in solidus expression.
A3: Constant used in solidus expression.
B1: Constant used in cpx-out expression.
B2: Constant used in cpx-out expression.
B3: Constant used in cpx-out expression.
C1: Constant used in liquidus expression.
C2: Constant used in liquidus expression.
C3: Constant used in liquidus expression.
a: Constant used in cpx-present melt fraction expression.
b: Constant used in cpx-present melt fraction expression.
c: Constant used in cpx-absent melt fraction expression.
d: Constant used in cpx-absent melt fraction expression.
alpha: Exponent used in the cpx-present melt fraction expression.
beta: Exponent used in the cpx-absent melt fraction expression.

The thermal expansivities, the heat capacity, the densities, and the entropy of fusion may also be changed during object initialisation.

Parameters:

CP (float, default: 1140.0) – The heat capacity (J K-1 kg-1)
alphas (float, default: 30.0) – The thermal expansivity of the solid (1e-6 K-1)
alphaf (float, default: 68.0) – The thermal expansivity of the melt (1e-6 K-1)
rhos (float, default: 3.3) – The density of the solid (kg m-3)
rhof (float, default: 2.9) – The density of the melt (kg m-3)
DeltaS (float, default: 380.0) – The entropy of fusion J K-1 kg-1
parameters (dict, default: parameters from Shorttle et al. (2014)) – The model parameters described above

F(P, T, **kwargs)[source]

Returns melt fraction at a given pressure and temperature. If below the solidus, returns 0. If above the liquidus, returns 1.

Prior to cpx exhaustion: F = a*(Tr)**alpha _ b*Tr Tr = (T-Tsol)/(Tliq-Tsol)

After cpx exhaustion: F = d*(Tr)**beta + c Tr = (T-Tcpx)/(Tliq-Tcpx)

Parameters:

P (float) – Pressure in GPa.
T (float) – Temperature in degC.

Returns:

Melt fraction between 0 and 1.

Return type:

float

TLiquidus(P, **kwargs)[source]

Returns liquidus temperature at a given pressure. T = C1 + C2*P + C3*P**2.

Parameters:: P (float, or list like) – Pressure in GPa.
Returns:: Liquidus temperature in degC.
Return type:: float, or list like

TSolidus(P, **kwargs)[source]

Returns solidus temperature at a given pressure. T = A1 + A2*P + A3*P**2.

Parameters:: P (float, or list like) – Pressure in GPa.
Returns:: Solidus Temperature in degC.
Return type:: float, or list like

dTdF(P, T, **kwargs)[source]

Returns dT/dF (constant P) at a given pressure and temperature. If below the solidus, or above the liquidus, _np.inf is returned.

Parameters:

P (float) – Pressure in GPa.
T (float) – Temperature in degC.

Returns:

dT/dF (constant P) in K.

Return type:

float

dTdP(P, T, **kwargs)[source]

Returns dT/dP (constant F) at a given pressure and temperature.

Parameters:

P (float) – Pressure in GPa.
T (float) – Temperature in degC.

Returns:

dT/dP (constant F) in K GPa-1.

Return type:

float

Module contents

Lithologies

The lithologies module provides lithology objects representing different bulk compositions, with models published by a number of authors.