spyice.utils.initial_userinput#

This module provides functions to calculate the initial melt temperature of seawater or ice based on boundary salinity. It supports multiple calculation methods: one-phase, Frezchem, and two-phase models.

Functions#

calculate_initial_melt_temperature_onephase#

calculate_initial_melt_temperature_onephase(boundary_salinity)#

Calculates the initial melt temperature using a one-phase model.

Parameters:

boundary_salinity (float) – Salinity at the boundary in ppt.

Returns:

Initial melt temperature in Kelvin.

Formula:

\[T_{melt} = 273.15 - 1.853 \frac{S}{28.0}\]

calculate_initial_temperature_frezchem#

calculate_initial_temperature_frezchem(boundary_salinity)#

Calculates the initial melt temperature using the Frezchem model.

Parameters:

boundary_salinity (float) – Salinity at the boundary in ppt.

Returns:

Initial melt temperature in Kelvin.

Formula:

\[T_{melt} = 272.63617665 - 9.1969758 \times 10^{-5} S^2 - 0.03942059 S\]

calculate_initial_melt_temperature_twophase#

calculate_initial_melt_temperature_twophase(boundary_salinity)#

Calculates the initial melt temperature using a two-phase model.

Parameters:

boundary_salinity (float) – Salinity at the boundary in ppt.

Returns:

Initial melt temperature in Kelvin.

Parameters used:

  • Solidus temperature: T_m = 273.15 K

  • Eutectic temperature: T_s = 252.05 K

  • Brine salinity: S_br = 233 ppt

Formula:

\[T_{melt} = 273.15 + \frac{(T_s - T_m) \cdot S}{S_{br}}\]

calculate_initial_melt_temperature#

calculate_initial_melt_temperature(boundary_salinity, method)#

Calculates the initial melt temperature using the specified method.

Parameters:

  • boundary_salinity (float) – Salinity at the boundary in ppt.

  • method (str) –

    Calculation method. Options are:

    • onephase: One phase model

    • Frezchem: Frezchem model

    • twophase: Two-phase model

Returns:

Initial melt temperature in Kelvin.

Raises:

ValueError – If an invalid method is specified.

Description: This function wraps the three specific calculation methods above and selects the appropriate formula based on the method argument.