ASPECT
ascii_reference_profile.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2016 - 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 
22 #ifndef _aspect_model_ascii_reference_profile_h
23 #define _aspect_model_ascii_reference_profile_h
24 
27 #include <aspect/utilities.h>
28 namespace aspect
29 {
30  namespace MaterialModel
31  {
32  using namespace dealii;
33 
45  template <int dim>
47  {
48  public:
53 
59  void
60  initialize () override;
61 
62  void evaluate(const MaterialModel::MaterialModelInputs<dim> &in,
63  MaterialModel::MaterialModelOutputs<dim> &out) const override;
64 
73  bool is_compressible () const override;
85  static
86  void
87  declare_parameters (ParameterHandler &prm);
88 
92  void
93  parse_parameters (ParameterHandler &prm) override;
94 
102  void
103  create_additional_named_outputs (MaterialModel::MaterialModelOutputs<dim> &out) const override;
104 
105 
106  private:
110  bool tala;
111 
115  double viscosity;
116 
122 
129  std::vector<double> viscosity_prefactors;
130 
136  std::vector<double> transition_depths;
137 
139 
144 
149  unsigned int density_index;
151  unsigned int specific_heat_index;
152  unsigned int compressibility_index;
153 
158  unsigned int seismic_vp_index;
159  unsigned int seismic_vs_index;
160  unsigned int seismic_dvp_dT_index;
161  unsigned int seismic_dvs_dT_index;
162  };
163 
164  }
165 }
166 
167 #endif
void declare_parameters(ParameterHandler &prm)
aspect::Utilities::AsciiDataProfile< dim > profile
Definition: compat.h:42