src.spyice.coefficients#

src.spyice.coefficients.update_coefficients(argument, x_initial, w, phi, nz, salinity_initial)[source]#

Updates of coefficients required to solve the Advection Reaction Diffusion Equation for each time step for temperature or salinity

Parameters:

  • argument (str) – Either ‘temperature’ or ‘salinity’

  • x_initial (float) – Initial value for salinity or temperature [ppt] or [K]

  • w (float) – Brine velocity [ms-1]

  • phi (float) – Liquid fraction [-]

  • nz (int) – Number of computational nodes

Returns:

A tuple containing the following coefficients:

  • a (numpy.ndarray): ‘temperature’: heat capacity | ‘salinity’: liquid fraction

  • b (numpy.ndarray): Brine velocity [ms-1]

  • c (numpy.ndarray): ‘temperature’: thermal conductivity | ‘salinity’: salt diffusivity

  • d (numpy.ndarray): ‘temperature’: latent heat | ‘salinity’: factor to determine salinity increase due to liquid fraction decrease

Return type:

tuple