PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ ice_temperature_guess_smb()

double pism::energy::ice_temperature_guess_smb ( EnthalpyConverter::Ptr  EC,
double  H,
double  z,
double  T_surface,
double  G,
double  ice_k,
double  K,
double  SMB 
)

A heuristic formula for the temperature distribution within a column of ice. Used during bootstrapping. Uses the surface mass balance to set the vertical velocity at the top surface.

Parameters
[in]ECenthalpy converter
[in]Hice thickness
[in]zheight above the base of the ice, 0 <= z <= H
[in]T_surfacesurface temperature, in Kelvin
[in]Gupward basal heat flux, in W / meter^2
[in]ice_kthermal conductivity of ice
[in]Ktemperature diffusivity in ice, ice_k / (ice_density * ice_c), where ice_c is the ice specific heat capacity
[in]SMBsurface mass balance in m / second

Definition at line 45 of file bootstrapping.cc.

References G, pism::hydrology::K(), pism::array::min(), and pism::diagnostics::SMB.

Referenced by bootstrap_ice_temperature().