PISM, A Parallel Ice Sheet Model  stable v2.0.4 committed by Constantine Khrulev on 2022-05-25 12:02:27 -0800
MohrCoulombYieldStress.hh
Go to the documentation of this file.
1 // Copyright (C) 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019, 2021 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 "YieldStress.hh"
23 
24 #include "pism/util/iceModelVec.hh"
25 
26 namespace pism {
27 
28 class IceModelVec2T;
29 
30 class IceModelVec2CellType;
31 
32 //! @brief PISM's default basal yield stress model which applies the
33 //! Mohr-Coulomb model of deformable, pressurized till.
35 public:
37  virtual ~MohrCoulombYieldStress() = default;
38 
39  void set_till_friction_angle(const IceModelVec2S &input);
40 protected:
41  void restart_impl(const File &input_file, int record);
42  void bootstrap_impl(const File &input_file, const YieldStressInputs &inputs);
43  void init_impl(const YieldStressInputs &inputs);
44 
45  void define_model_state_impl(const File &output) const;
46  void write_model_state_impl(const File &output) const;
47 
49 
50  MaxTimestep max_timestep_impl(double t) const;
51  void update_impl(const YieldStressInputs &inputs, double t, double dt);
52 
53  void finish_initialization(const YieldStressInputs &inputs);
54 
56 
57  std::shared_ptr<IceModelVec2T> m_delta;
58 private:
59  void till_friction_angle(const IceModelVec2S &bed_topography,
60  IceModelVec2S &result);
61 
62  void till_friction_angle(const IceModelVec2S &basal_yield_stress,
63  const IceModelVec2S &till_water_thickness,
64  const IceModelVec2S &ice_thickness,
65  const IceModelVec2CellType &cell_type,
66  IceModelVec2S &result);
67 };
68 
69 } // end of namespace pism
70 
71 #endif /* _PISMMOHRCOULOMBYIELDSTRESS_H_ */
High-level PISM I/O class.
Definition: File.hh:51
std::shared_ptr< const IceGrid > ConstPtr
Definition: IceGrid.hh:233
"Cell type" mask. Adds convenience methods to IceModelVec2Int.
Combines the max. time step with the flag indicating if a restriction is active. Makes is possible to...
Definition: MaxTimestep.hh:31
std::shared_ptr< IceModelVec2T > m_delta
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 write_model_state_impl(const File &output) const
The default (empty implementation).
void set_till_friction_angle(const IceModelVec2S &input)
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
void till_friction_angle(const IceModelVec2S &bed_topography, IceModelVec2S &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.
MohrCoulombYieldStress(IceGrid::ConstPtr g)
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)
Definition: YieldStress.hh:44
static const double g
Definition: exactTestP.cc:39
std::map< std::string, Diagnostic::Ptr > DiagnosticList
Definition: Diagnostic.hh:117