ASPECT
SAVANI_perturbation.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2011 - 2019 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_SAVANI_perturbation_h
23 #define _aspect_initial_temperature_SAVANI_perturbation_h
24 
26 #include <aspect/utilities.h>
27 
28 namespace aspect
29 {
30  namespace InitialTemperature
31  {
32  using namespace dealii;
33 
34  namespace internal
35  {
36  namespace SAVANI
37  {
38  class SphericalHarmonicsLookup;
39  class SplineDepthsLookup;
40  }
41  }
42 
52  template <int dim>
53  class SAVANIPerturbation : public Interface<dim>, public ::aspect::SimulatorAccess<dim>
54  {
55  public:
60 
61 
66  void
67  initialize () override;
68 
72  //TODO
73  virtual
74  double get_Vs (const Point<dim> &position) const;
75 
79  double initial_temperature (const Point<dim> &position) const override;
80 
84  static
85  void
86  declare_parameters (ParameterHandler &prm);
87 
91  void
92  parse_parameters (ParameterHandler &prm) override;
93 
94 
95  private:
96 
101  {
103  constant
104  };
105 
110 
114  std::string data_directory;
116 
122 
137 
144 
150 
155  unsigned int max_order;
156 
163 
168  std::unique_ptr<internal::SAVANI::SphericalHarmonicsLookup> spherical_harmonics_lookup;
169 
174  std::unique_ptr<internal::SAVANI::SplineDepthsLookup> spline_depths_lookup;
175 
180 
184  unsigned int vs_to_density_index;
185 
190 
191  };
192 
193  }
194 }
195 
196 #endif
aspect::Utilities::AsciiDataProfile< dim > profile
std::unique_ptr< internal::SAVANI::SphericalHarmonicsLookup > spherical_harmonics_lookup
void declare_parameters(ParameterHandler &prm)
std::unique_ptr< internal::SAVANI::SplineDepthsLookup > spline_depths_lookup
Definition: compat.h:88