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
OptTillphiYieldStress.hh
Go to the documentation of this file.
1// Copyright (C) 2011--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 _PISMOPTTILLPHIYIELDSTRESS_H_
20#define _PISMOPTTILLPHIYIELDSTRESS_H_
21
22#include "pism/basalstrength/MohrCoulombYieldStress.hh"
23
24namespace pism {
25
26//! Iterative optimization of the till friction angle.
28public:
29 OptTillphiYieldStress(std::shared_ptr<const Grid> g);
30 virtual ~OptTillphiYieldStress() = default;
31
32private:
34
35 void update_tillphi(const array::Scalar &ice_surface_elevation,
36 const array::Scalar &bed_topography,
37 const array::CellType &mask);
38
39 void init_t_last(const File &input_file);
40 void init_usurf_target(const File &input_file);
41
42 void init_impl(const YieldStressInputs &inputs);
43 void bootstrap_impl(const File &input_file, const YieldStressInputs &inputs);
44 void restart_impl(const File &input_file, int record);
45 void update_impl(const YieldStressInputs &inputs, double t, double dt);
46
47 MaxTimestep max_timestep_impl(double t) const;
48
49 void define_model_state_impl(const File &output) const;
50 void write_model_state_impl(const File &output) const;
51
55
57
58 // convergence threshold:
59 double m_dhdt_min;
60
61 // lower and upper bounds of a tillphi adjustment:
62 double m_dphi_min;
63 double m_dphi_max;
64
65 // constants defining the lower bound of tillphi:
66 double m_phi0_min;
67 double m_phi0_max;
68 double m_topg_min;
69 double m_topg_max;
70
71 // the upper bound of tillphi:
72 double m_phi_max;
73
74 //! time of the last till friction angle update
75 double m_t_last;
76 //! Update interval in seconds
78 //! Temporal resolution to use when checking whether it's time to update
79 double m_t_eps;
80 //! Name of the variable used to store the last update time.
81 std::string m_time_name;
82};
83
84} // end of namespace pism
85
86#endif /* _PISMOPTTILLPHIYIELDSTRESS_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...
PISM's default basal yield stress model which applies the Mohr-Coulomb model of deformable,...
void define_model_state_impl(const File &output) const
MaxTimestep max_timestep_impl(double t) const
void bootstrap_impl(const File &input_file, const YieldStressInputs &inputs)
void init_t_last(const File &input_file)
void init_impl(const YieldStressInputs &inputs)
double m_update_interval
Update interval in seconds.
DiagnosticList diagnostics_impl() const
double m_t_last
time of the last till friction angle update
void update_tillphi(const array::Scalar &ice_surface_elevation, const array::Scalar &bed_topography, const array::CellType &mask)
virtual ~OptTillphiYieldStress()=default
void restart_impl(const File &input_file, int record)
std::string m_time_name
Name of the variable used to store the last update time.
double m_t_eps
Temporal resolution to use when checking whether it's time to update.
void write_model_state_impl(const File &output) const
The default (empty implementation).
void update_impl(const YieldStressInputs &inputs, double t, double dt)
void init_usurf_target(const File &input_file)
Iterative optimization of the till friction angle.
"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