ASPECT
initial_lithostatic_pressure.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_boundary_traction_initial_lithospheric_pressure_h
23 #define _aspect_boundary_traction_initial_lithospheric_pressure_h
24 
27 
28 
29 namespace aspect
30 {
31  namespace BoundaryTraction
32  {
33  using namespace dealii;
34 
41  template <int dim>
42  class InitialLithostaticPressure : public Interface<dim>, public SimulatorAccess<dim>
43  {
44  public:
45 
52  void initialize () override;
53 
54 
60  Tensor<1,dim>
61  boundary_traction (const types::boundary_id boundary_indicator,
62  const Point<dim> &position,
63  const Tensor<1,dim> &normal_vector) const override;
64 
65 
69  static
70  void
71  declare_parameters (ParameterHandler &prm);
72 
73 
77  void
78  parse_parameters (ParameterHandler &prm) override;
79 
80  private:
81 
85  unsigned int n_points;
86 
90  double delta_z;
91 
92  /*
93  * The user-specified point where to calculate the pressure profile.
94  * The vertical coordinate/radius is ignored.
95  */
97 
101  std::vector<double> pressure;
102 
107  double interpolate_pressure (const Point<dim> &p) const;
108 
112  types::boundary_id bottom_boundary_id;
113 
121  };
122  }
123 }
124 
125 
126 #endif
void declare_parameters(ParameterHandler &prm)
Definition: compat.h:42