ASPECT
depth_dependent.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2014 - 2022 by the authors of the ASPECT code.
3 
4  This file is part of ASPECT.
5 
6  ASPECT is free software; you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation; either version 2, or (at your option)
9  any later version.
10 
11  ASPECT is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with ASPECT; see the file LICENSE. If not see
18  <http://www.gnu.org/licenses/>.
19 */
20 
21 #ifndef _aspect_material_model_depth_dependent_h
22 #define _aspect_material_model_depth_dependent_h
23 
24 #include <deal.II/base/function_lib.h>
27 #include <deal.II/base/parsed_function.h>
29 
30 #include <boost/lexical_cast.hpp>
31 
32 namespace aspect
33 {
34  namespace MaterialModel
35  {
36  using namespace dealii;
37 
47  template <int dim>
49  {
50  public:
54  void initialize() override;
55 
60  void update() override;
61 
66  void
67  evaluate (const typename Interface<dim>::MaterialModelInputs &in,
68  typename Interface<dim>::MaterialModelOutputs &out) const override;
72  static void
73  declare_parameters (ParameterHandler &prm);
74 
78  void
79  parse_parameters (ParameterHandler &prm) override;
80 
85  bool is_compressible () const override;
86 
87  private:
88 
93  {
97  None
98  };
99 
100  // The viscosity that will be used as a reference for scaling the depth-dependent
101  // prefactor.
103 
108 
113  std::unique_ptr<Functions::InterpolatedTensorProductGridData<1>> viscosity_file_function;
114 
118  double
119  viscosity_from_list(const double &depth) const;
120 
124  double
125  viscosity_from_file(const double &depth) const;
126 
131  double
132  calculate_depth_dependent_prefactor(const double &depth) const;
133 
137  std::vector<double> depth_values;
138  std::vector<double> viscosity_values;
139 
144  Functions::ParsedFunction<1> viscosity_function;
145 
149  std::unique_ptr<MaterialModel::Interface<dim>> base_model;
150 
155  std::unique_ptr<Rheology::AsciiDepthProfile<dim>> depth_dependent_rheology;
156  };
157  }
158 }
159 
160 #endif
std::unique_ptr< Rheology::AsciiDepthProfile< dim > > depth_dependent_rheology
std::unique_ptr< Functions::InterpolatedTensorProductGridData< 1 > > viscosity_file_function
Functions::ParsedFunction< 1 > viscosity_function
void declare_parameters(ParameterHandler &prm)
std::unique_ptr< MaterialModel::Interface< dim > > base_model
Definition: compat.h:88