::icebin::gpism::IceCoupler_PISM accumulateFluxes_massContExplicitStep allocate_couplers allocate_subglacial_hydrology base compute_enth2 construct_surface_temp cur dumpToFile elevmask_ice elevmask_land energy_step enthalpy_t ib_surface_model IBIceModel ice_top_senth m_ice_density m_meter_per_s_to_kg_per_m2 mass_t massContExplicitStep massContPostHook massContPreHook misc_setup null_hydrology params post_energy_nc post_mass_nc pre_energy_nc pre_mass_nc prepare_initial_outputs prepare_outputs rate reset_rate set_rate super time_setup ~IBIceModel

compute_enth2() void pism::icebin::IBIceModel::compute_enth2 ( pism::array::Scalar &  enth2, pism::array::Scalar &  mass2  ) Sums over columns to compute enthalpy on 2D m_grid. This is used to track conservation of energy. NOTE: Unfortunately so far PISM does not keep track of enthalpy in "partially-filled" cells, so Enth2(i,j) is not valid at locations like this one. We need to address this, but so far, it seems to me, the best thing we can do is estimate Enth2(i,j) at partially-filled cells by computing the average over icy neighbors. I think you can re-use the idea from IceModel::get_threshold_thickness(...) (iMpartm_grid->cc). Definition at line 372 of file IBIceModel.cc. References pism::Geometry::ice_thickness, pism::k, pism::IceModel::m_config, pism::IceModel::m_energy_model, pism::IceModel::m_geometry, and pism::IceModel::m_grid. Referenced by misc_setup(), and set_rate().