Loading [MathJax]/extensions/tex2jax.js
PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
BlatterMod.hh
Go to the documentation of this file.
1/* Copyright (C) 2020, 2021, 2022, 2023 PISM Authors
2 *
3 * This file is part of PISM.
4 *
5 * PISM is free software; you can redistribute it and/or modify it under the
6 * terms of the GNU General Public License as published by the Free Software
7 * Foundation; either version 3 of the License, or (at your option) any later
8 * version.
9 *
10 * PISM is distributed in the hope that it will be useful, but WITHOUT ANY
11 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
12 * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
13 * details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with PISM; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
18 */
19
20#ifndef PISM_BLATTER_MOD_H
21#define PISM_BLATTER_MOD_H
22
23#include <memory> // std::shared_ptr
24
25#include "pism/stressbalance/blatter/Blatter.hh"
26#include "pism/stressbalance/SSB_Modifier.hh"
27
28namespace pism {
29namespace stressbalance {
30
31/*!
32 * The "modifier" post-processing the velocity field computed using `Blatter`.
33 */
34class BlatterMod : public SSB_Modifier {
35public:
36 BlatterMod(std::shared_ptr<Blatter> solver);
37 virtual ~BlatterMod() = default;
38
39 void init();
40
41 void update(const array::Vector &sliding_velocity,
42 const Inputs &inputs,
43 bool full_update);
44private:
45 std::shared_ptr<Blatter> m_solver;
46
47 void transfer(const array::Scalar &ice_thickness);
48
49 void compute_max_diffusivity(const array::Vector &velocity,
50 const array::Scalar &ice_thickness,
51 const array::Scalar1 &surface);
52};
53
54} // end of namespace stressbalance
55} // end of namespace pism
56
57#endif /* PISM_BLATTER_MOD_H */
std::shared_ptr< Blatter > m_solver
Definition BlatterMod.hh:45
void update(const array::Vector &sliding_velocity, const Inputs &inputs, bool full_update)
Definition BlatterMod.cc:53
void transfer(const array::Scalar &ice_thickness)
Definition BlatterMod.cc:75
void compute_max_diffusivity(const array::Vector &velocity, const array::Scalar &ice_thickness, const array::Scalar1 &surface)
Shallow stress balance modifier (such as the non-sliding SIA).