tIsothermalGlen.cc - 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
---
tIsothermalGlen.cc (1739B)
---
1 /* Copyright (C) 2015, 2016, 2017 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 #include <cmath>
21
22 #include "IsothermalGlen.hh"
23 #include "pism/util/ConfigInterface.hh"
24
25 namespace pism {
26 namespace rheology {
27
28 IsothermalGlen::IsothermalGlen(const std::string &prefix,
29 const Config &config, EnthalpyConverter::Ptr ec)
30 : PatersonBudd(prefix, config, ec) {
31 m_name = "isothermal Glen";
32
33 m_softness_A = config.get_number("flow_law.isothermal_Glen.ice_softness");
34 m_hardness_B = pow(m_softness_A, m_hardness_power);
35 }
36
37 double IsothermalGlen::flow_impl(double stress, double, double, double) const {
38 return m_softness_A * pow(stress, m_n-1);
39 }
40
41 double IsothermalGlen::softness_impl(double, double) const {
42 return m_softness_A;
43 }
44
45 double IsothermalGlen::hardness_impl(double, double) const {
46 return m_hardness_B;
47 }
48
49 double IsothermalGlen::flow_from_temp(double stress, double, double, double) const {
50 return m_softness_A * pow(stress,m_n-1);
51 }
52
53 } // end of namespace rheology
54 } // end of namespace pism