EnthalpyModel_Regional.cc

Go to the documentation of this file.

/* Copyright (C) 2016, 2017, 2019, 2020, 2022, 2023 PISM Authors
 *
 * This file is part of PISM.
 *
 * PISM is free software; you can redistribute it and/or modify it under the
 * terms of the GNU General Public License as published by the Free Software
 * Foundation; either version 3 of the License, or (at your option) any later
 * version.
 *
 * PISM is distributed in the hope that it will be useful, but WITHOUT ANY
 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
 * FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
 * details.
 *
 * You should have received a copy of the GNU General Public License
 * along with PISM; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 */
 
#include "pism/regional/EnthalpyModel_Regional.hh"
 
namespace pism {
namespace energy {
 
EnthalpyModel_Regional::EnthalpyModel_Regional(
    std::shared_ptr<const Grid> grid,
    std::shared_ptr<const stressbalance::StressBalance> stress_balance)
    : EnthalpyModel(grid, stress_balance), m_basal_melt_rate_stored(m_grid, "bmr_stored") {
  // Note that the name of this variable (bmr_stored) does not matter: it is
  // *never* read or written. We make a copy of basal_melt_rate_grounded instead.
  m_basal_melt_rate_stored.metadata(0)
      .long_name("time-independent basal melt rate in the no-model-strip")
      .units("m s^-1");
}
EnthalpyModel_Regional::EnthalpyModel_Regional( {…}
 
void EnthalpyModel_Regional::restart_impl(const File &input_file, int record) {
  EnthalpyModel::restart_impl(input_file, record);
 
  m_basal_melt_rate_stored.copy_from(m_basal_melt_rate);
}
void EnthalpyModel_Regional::restart_impl(const File &input_file, int record) {…}
 
void EnthalpyModel_Regional::bootstrap_impl(const File &input_file,
                                            const array::Scalar &ice_thickness,
                                            const array::Scalar &surface_temperature,
                                            const array::Scalar &climatic_mass_balance,
                                            const array::Scalar &basal_heat_flux) {
 
  EnthalpyModel::bootstrap_impl(input_file, ice_thickness, surface_temperature,
                               climatic_mass_balance, basal_heat_flux);
 
  m_basal_melt_rate_stored.copy_from(m_basal_melt_rate);
}
void EnthalpyModel_Regional::bootstrap_impl(const File &input_file, {…}
 
void EnthalpyModel_Regional::initialize_impl(const array::Scalar &basal_melt_rate,
                                             const array::Scalar &ice_thickness,
                                             const array::Scalar &surface_temperature,
                                             const array::Scalar &climatic_mass_balance,
                                             const array::Scalar &basal_heat_flux) {
 
  EnthalpyModel::initialize_impl(basal_melt_rate,
                                 ice_thickness,
                                 surface_temperature,
                                 climatic_mass_balance,
                                 basal_heat_flux);
 
  m_basal_melt_rate_stored.copy_from(m_basal_melt_rate);
}
void EnthalpyModel_Regional::initialize_impl(const array::Scalar &basal_melt_rate, {…}
 
 
void EnthalpyModel_Regional::update_impl(double t, double dt,
                                         const Inputs &inputs) {
 
  unsigned int Mz = m_grid->Mz();
 
  EnthalpyModel::update_impl(t, dt, inputs);
 
  const array::Scalar &no_model_mask = *inputs.no_model_mask;
 
  // The update_impl() call above sets m_work; ghosts are communicated
  // later (in EnergyModel::update()).
  array::AccessScope list{&no_model_mask, &m_work, &m_ice_enthalpy,
      &m_basal_melt_rate, &m_basal_melt_rate_stored};
 
  for (auto p = m_grid->points(); p; p.next()) {
    const int i = p.i(), j = p.j();
 
    if (no_model_mask(i, j) > 0.5) {
      double *new_enthalpy = m_work.get_column(i, j);
      double *old_enthalpy = m_ice_enthalpy.get_column(i, j);
 
      // enthalpy
      for (unsigned int k = 0; k < Mz; ++k) {
        new_enthalpy[k] = old_enthalpy[k];
      }
 
      // basal melt rate
      m_basal_melt_rate(i, j) = m_basal_melt_rate_stored(i, j);
    }
  }
}
void EnthalpyModel_Regional::update_impl(double t, double dt, {…}
 
} // end of namespace energy
} // end of namespace pism