ASPECT
sphere.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2014 - 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 #ifndef _aspect_geometry_model_sphere_h
21 #define _aspect_geometry_model_sphere_h
22 
26 
27 namespace aspect
28 {
29  namespace GeometryModel
30  {
31  using namespace dealii;
32 
33  template <int dim>
34  class Sphere : public Interface<dim>, public SimulatorAccess<dim>
35  {
36  public:
40  Sphere();
41 
45  void create_coarse_mesh (parallel::distributed::Triangulation<dim> &coarse_grid) const override;
46 
54  double length_scale () const override;
55 
69  double depth(const Point<dim> &position) const override;
70 
75  double height_above_reference_surface(const Point<dim> &position) const override;
76 
80  Point<dim> representative_point(const double depth) const override;
81 
85  double maximal_depth() const override;
86 
94  std::set<types::boundary_id>
95  get_used_boundary_indicators () const override;
96 
102  std::map<std::string,types::boundary_id>
103  get_symbolic_boundary_names_map () const override;
104 
110  bool
111  has_curved_elements() const override;
112 
117  aspect::Utilities::Coordinates::CoordinateSystem natural_coordinate_system() const override;
118 
124  bool
125  point_is_in_domain(const Point<dim> &point) const override;
126 
132  std::array<double,dim> cartesian_to_natural_coordinates(const Point<dim> &position) const override;
133 
139  Point<dim> natural_to_cartesian_coordinates(const std::array<double,dim> &position) const override;
140 
141 
142  static
143  void
145 
146  void
147  parse_parameters (ParameterHandler &prm) override;
148 
152  double
153  radius () const;
154 
155  private:
159  double R;
160 
165  };
166  }
167 }
168 
169 
170 #endif
float depth
void declare_parameters(ParameterHandler &prm)
Point< spacedim > point(const gp_Pnt &p, const double tolerance=1e-10)
const SphericalManifold< dim > spherical_manifold
Definition: sphere.h:164