tHooke.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
       ---
       tHooke.cc (1637B)
       ---
            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 "Hooke.hh"
           23 #include "pism/util/ConfigInterface.hh"
           24 
           25 namespace pism {
           26 namespace rheology {
           27 
           28 // Hooke
           29 
           30 Hooke::Hooke(const std::string &prefix,
           31              const Config &config, EnthalpyConverter::Ptr ec)
           32   : PatersonBudd(prefix, config, ec) {
           33   m_name = "Hooke";
           34 
           35   m_Q_Hooke  = config.get_number("flow_law.Hooke.Q");
           36   m_A_Hooke  = config.get_number("flow_law.Hooke.A");
           37   m_C_Hooke  = config.get_number("flow_law.Hooke.C");
           38   m_K_Hooke  = config.get_number("flow_law.Hooke.k");
           39   m_Tr_Hooke = config.get_number("flow_law.Hooke.Tr");
           40 }
           41 
           42 Hooke::~Hooke() {
           43   // empty
           44 }
           45 
           46 double Hooke::softness_from_temp(double T_pa) const {
           47   return m_A_Hooke * exp(-m_Q_Hooke/(m_ideal_gas_constant * T_pa)
           48                          + 3.0 * m_C_Hooke * pow(m_Tr_Hooke - T_pa, -m_K_Hooke));
           49 }
           50 
           51 } // end of namespace rheology
           52 } // end of namespace pism