tMohrCoulombYieldStress.hh - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
HTML git clone git://src.adamsgaard.dk/pism
DIR Log
DIR Files
DIR Refs
DIR LICENSE
---
tMohrCoulombYieldStress.hh (2492B)
---
1 // Copyright (C) 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019 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 #include "pism/util/iceModelVec2T.hh"
26
27 namespace pism {
28
29 class IceModelVec2CellType;
30
31 //! @brief PISM's default basal yield stress model which applies the
32 //! Mohr-Coulomb model of deformable, pressurized till.
33 class MohrCoulombYieldStress : public YieldStress {
34 public:
35 MohrCoulombYieldStress(IceGrid::ConstPtr g);
36 virtual ~MohrCoulombYieldStress();
37
38 void set_till_friction_angle(const IceModelVec2S &input);
39 private:
40 void restart_impl(const File &input_file, int record);
41 void bootstrap_impl(const File &input_file, const YieldStressInputs &inputs);
42 void init_impl(const YieldStressInputs &inputs);
43
44 void define_model_state_impl(const File &output) const;
45 void write_model_state_impl(const File &output) const;
46
47 DiagnosticList diagnostics_impl() const;
48
49 MaxTimestep max_timestep_impl(double t) const;
50 void update_impl(const YieldStressInputs &inputs, double t, double dt);
51
52 void finish_initialization(const YieldStressInputs &inputs);
53 private:
54 void till_friction_angle(const IceModelVec2S &bed_topography,
55 IceModelVec2S &result);
56
57 void till_friction_angle(const IceModelVec2S &basal_yield_stress,
58 const IceModelVec2S &till_water_thickness,
59 const IceModelVec2S &ice_thickness,
60 const IceModelVec2CellType &cell_type,
61 IceModelVec2S &result);
62
63 IceModelVec2S m_till_phi;
64
65 IceModelVec2T::Ptr m_delta;
66 };
67
68 } // end of namespace pism
69
70 #endif /* _PISMMOHRCOULOMBYIELDSTRESS_H_ */