ASPECT
spherical_shell.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2011 - 2023 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_spherical_shell_h
23 #define _aspect_initial_temperature_spherical_shell_h
24 
27 
28 
29 namespace aspect
30 {
31  namespace InitialTemperature
32  {
33  using namespace dealii;
34 
41  template <int dim>
42  class SphericalHexagonalPerturbation : public Interface<dim>, public SimulatorAccess<dim>
43  {
44  public:
48  double initial_temperature (const Point<dim> &position) const override;
49 
53  static
54  void
55  declare_parameters (ParameterHandler &prm);
56 
60  void
61  parse_parameters (ParameterHandler &prm) override;
62 
63  private:
70  unsigned int angular_mode;
71 
80 
84  double R1;
85 
86  };
87 
88 
95  template <int dim>
96  class SphericalGaussianPerturbation : public Interface<dim>, public SimulatorAccess<dim>
97  {
98  public:
99 
104 
108  double initial_temperature (const Point<dim> &position) const override;
109 
113  static
114  void
115  declare_parameters (ParameterHandler &prm);
116 
120  void
121  parse_parameters (ParameterHandler &prm) override;
122 
123  private:
124  double angle;
125  double depth;
126  double amplitude;
127  double sigma;
128  double sign;
129  unsigned int npoint;
130 
131  std::vector<double> radial_position;
132  std::vector<double> geotherm;
133 
137  double R0;
138 
142  double R1;
143  };
144  }
145 }
146 
147 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:59
Definition: compat.h:42