GivenTH init_impl m_salinity_ocean m_theta_ocean max_timestep_impl pointwise_update subshelf_salinity subshelf_salinity_diffusion_only subshelf_salinity_freeze_on subshelf_salinity_melt update_impl ~GivenTH

subshelf_salinity_freeze_on() void pism::ocean::GivenTH::subshelf_salinity_freeze_on ( const Constants &  c, double  sea_water_salinity, double  sea_water_potential_temperature, double  thickness, double *  shelf_base_salinity  ) private Compute basal salinity in the basal freeze-on case. In this case we assume that the temperature gradient at the shelf base is zero: \[ T_{\text{grad}} = 0. \] Please see pointwise_update() for details. In this case the coefficients of the quadratic equation for the basal salinity are: \begin{align*} A &= -b_{0}\,\gamma_T\,c_{pW} \\ B &= \gamma_S\,L+\gamma_T\,c_{pW}\,\left(\Theta^W-b_{2}\,h-b_{1}\right) \\ C &= -\gamma_S\,S^W\,L\\ \end{align*} Parameters [in] c model constants [in] sea_water_salinity sea water salinity [in] sea_water_potential_temperature sea water temperature [in] thickness ice shelf thickness [out] shelf_base_salinity resulting basal salinity Returns 0 on success Definition at line 427 of file GivenTH.cc. References pism::ocean::GivenTH::Constants::b, C, pism::ocean::GivenTH::Constants::gamma_S, pism::ocean::GivenTH::Constants::gamma_T, L, pism::ocean::GivenTH::Constants::sea_water_specific_heat_capacity, and pism::ocean::GivenTH::Constants::water_latent_heat_fusion. Referenced by subshelf_salinity().