Loading [MathJax]/extensions/tex2jax.js
PISM, A Parallel Ice Sheet Model 2.2.1-cd005eec8 committed by Constantine Khrulev on 2025-03-07
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
MohrCoulombYieldStress.hh
Go to the documentation of this file.
1// Copyright (C) 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019, 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#ifndef _PISMMOHRCOULOMBYIELDSTRESS_H_
20#define _PISMMOHRCOULOMBYIELDSTRESS_H_
21
22#include "pism/basalstrength/YieldStress.hh"
23
24namespace pism {
25
26//! @brief PISM's default basal yield stress model which applies the
27//! Mohr-Coulomb model of deformable, pressurized till.
29public:
30 MohrCoulombYieldStress(std::shared_ptr<const Grid> g);
31 virtual ~MohrCoulombYieldStress() = default;
32
33 void set_till_friction_angle(const array::Scalar &input);
34protected:
35 void restart_impl(const File &input_file, int record);
36 void bootstrap_impl(const File &input_file, const YieldStressInputs &inputs);
37 void init_impl(const YieldStressInputs &inputs);
38
39 void define_model_state_impl(const File &output) const;
40 void write_model_state_impl(const File &output) const;
41
43
44 MaxTimestep max_timestep_impl(double t) const;
45 void update_impl(const YieldStressInputs &inputs, double t, double dt);
46
47 void finish_initialization(const YieldStressInputs &inputs);
48
50
51 std::shared_ptr<array::Forcing> m_delta;
52private:
53 void till_friction_angle(const array::Scalar &bed_topography,
54 array::Scalar &result);
55
56 void till_friction_angle(const array::Scalar &basal_yield_stress,
57 const array::Scalar &till_water_thickness,
58 const array::Scalar &ice_thickness,
59 const array::CellType &cell_type,
60 array::Scalar &result);
61};
62
63} // end of namespace pism
64
65#endif /* _PISMMOHRCOULOMBYIELDSTRESS_H_ */
High-level PISM I/O class.
Definition File.hh:55
Combines the max. time step with the flag indicating if a restriction is active. Makes is possible to...
void init_impl(const YieldStressInputs &inputs)
void restart_impl(const File &input_file, int record)
void update_impl(const YieldStressInputs &inputs, double t, double dt)
DiagnosticList diagnostics_impl() const
void set_till_friction_angle(const array::Scalar &input)
void write_model_state_impl(const File &output) const
The default (empty implementation).
void define_model_state_impl(const File &output) const
void finish_initialization(const YieldStressInputs &inputs)
MaxTimestep max_timestep_impl(double t) const
virtual ~MohrCoulombYieldStress()=default
std::shared_ptr< array::Forcing > m_delta
void till_friction_angle(const array::Scalar &bed_topography, array::Scalar &result)
Computes the till friction angle phi as a piecewise linear function of bed elevation,...
void bootstrap_impl(const File &input_file, const YieldStressInputs &inputs)
Initialize the pseudo-plastic till mechanical model.
PISM's default basal yield stress model which applies the Mohr-Coulomb model of deformable,...
The PISM basal yield stress model interface (virtual base class)
"Cell type" mask. Adds convenience methods to array::Scalar.
Definition CellType.hh:30
static const double g
Definition exactTestP.cc:36
std::map< std::string, Diagnostic::Ptr > DiagnosticList