P-T-F Parameterization

peridoxia.PTF_functions.PTF_geotherm(TK_pot: float = 1623.15, PGPa_range: list[float] = [4.0, 0.1], P_step: float = -0.01, melting: bool = True)

Calculates a set of PTF values along a geotherm from a given potential temperature. Equations are from Langmuir, Klein, and Plank (1992).

Parameters:

TK_potfloat, optional

Potential temperature, by default 1350+273.15

PGPa_rangelist[float], optional

Pressure range, by default [4.,0.1]

P_stepfloat, optional

Pressure increment, by default -0.01

meltingbool, optional

Whether to include melting, by default True

Returns:

np.ndarray, np.ndarray, np.ndarray

Arrays for P (GPa), T (K), and F

Raises:

ValueError

Raised if solidus pressure is higher than initial pressure