PISM, A Parallel Ice Sheet Model  stable v2.0.4 committed by Constantine Khrulev on 2022-05-25 12:02:27 -0800
PatersonBuddCold.hh
Go to the documentation of this file.
1 /* Copyright (C) 2015, 2016, 2018, 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 
20 #ifndef _PATERSONBUDDCOLD_H_
21 #define _PATERSONBUDDCOLD_H_
22 
23 #include "PatersonBudd.hh"
24 
25 namespace pism {
26 namespace rheology {
27 
28 //! Cold case of Paterson-Budd
30 public:
31  PatersonBuddCold(const std::string &prefix, const Config &config, EnthalpyConverter::Ptr EC);
32  virtual ~PatersonBuddCold() = default;
33 
34  //! Return the temperature T corresponding to a given value A=A(T).
35  double tempFromSoftness(double A) const;
36 
37 protected:
38  // takes care of hardness...
39  double softness_from_temp(double T_pa) const;
40 
41  // ignores pressure and uses non-pressure-adjusted temperature
42  double flow_from_temp(double stress, double temp,
43  double , double) const;
44 };
45 
46 bool FlowLawIsPatersonBuddCold(const FlowLaw &flow_law,
48 
49 } // end of namespace rheology
50 } // end of namespace pism
51 
52 
53 #endif /* _PATERSONBUDDCOLD_H_ */
A class for storing and accessing PISM configuration flags and parameters.
std::shared_ptr< EnthalpyConverter > Ptr
EnthalpyConverter::Ptr EC() const
Definition: FlowLaw.cc:71
virtual ~PatersonBuddCold()=default
double softness_from_temp(double T_pa) const
double tempFromSoftness(double A) const
Return the temperature T corresponding to a given value A=A(T).
PatersonBuddCold(const std::string &prefix, const Config &config, EnthalpyConverter::Ptr EC)
double flow_from_temp(double stress, double temp, double, double) const
The flow law (temperature-dependent version).
Cold case of Paterson-Budd.
Derived class of FlowLaw for Paterson-Budd (1982)-Glen ice.
Definition: PatersonBudd.hh:29
bool FlowLawIsPatersonBuddCold(const FlowLaw &flow_law, const Config &config, EnthalpyConverter::Ptr EC)