ASPECT
adiabatic.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2012 - 2024 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_initial_temperature_adiabatic_h
23 #define _aspect_initial_temperature_adiabatic_h
24 
27 #include <aspect/utilities.h>
28 
29 #include <deal.II/base/parsed_function.h>
30 
31 namespace aspect
32 {
33  namespace InitialTemperature
34  {
35  using namespace dealii;
36 
45  namespace BoundaryLayerAgeModel
46  {
47  enum Kind
48  {
51  function
52  };
53 
60  parse (const std::string &parameter_name,
61  const ParameterHandler &prm);
62  }
63 
72  template <int dim>
73  class Adiabatic : public Interface<dim>, public Utilities::AsciiDataBoundary<dim>
74  {
75  public:
79  Adiabatic ();
80 
81  void initialize () override;
82 
83  // avoid -Woverloaded-virtual:
85 
89  double initial_temperature (const Point<dim> &position) const override;
90 
94  static
95  void
96  declare_parameters (ParameterHandler &prm);
97 
101  void
102  parse_parameters (ParameterHandler &prm) override;
103 
104  private:
109  types::boundary_id surface_boundary_id;
110 
117 
118  /* Age of the lower thermal boundary layer. */
120 
125  double radius;
130  double amplitude;
131  /*
132  * Position of the initial temperature perturbation (in the
133  * center or at the boundary of the model domain).
134  */
136 
137  /*
138  * Deviation from adiabaticity in a subadiabatic mantle
139  * temperature profile. 0 for an adiabatic temperature
140  * profile.
141  */
143 
148 
149  /*
150  * Whether to use the half space cooling model, or the plate cooling
151  * model
152  */
153  std::string cooling_model;
154 
155  /*
156  * Depth to the base of the lithosphere for plate cooling model, in m
157  */
159 
165  std::unique_ptr<Functions::ParsedFunction<1>> function;
166 
170  Functions::ParsedFunction<dim> age_function;
171 
177 
181  double top_boundary_layer_age(const Point<dim> &position) const;
182  };
183  }
184 }
185 
186 
187 #endif
BoundaryLayerAgeModel::Kind top_boundary_layer_age_model
Definition: adiabatic.h:147
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:59
Functions::ParsedFunction< dim > age_function
Definition: adiabatic.h:170
Utilities::Coordinates::CoordinateSystem coordinate_system
Definition: adiabatic.h:176
BoundaryLayerAgeModel::Kind parse(const std::string &parameter_name, const ParameterHandler &prm)
Definition: compat.h:42
types::boundary_id surface_boundary_id
Definition: adiabatic.h:109