ASPECT
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules Pages
compute_entropy_profile.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2016 - 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_adiabatic_conditions_compute_entropy_profile_h
23 #define _aspect_adiabatic_conditions_compute_entropy_profile_h
24 
25 
28 
29 namespace aspect
30 {
31  namespace AdiabaticConditions
32  {
43  template <int dim>
44  class ComputeEntropyProfile : public Interface<dim>
45  {
46  public:
51 
59  void initialize () override;
60 
71  bool is_initialized() const override;
72 
76  double temperature (const Point<dim> &p) const override;
77 
81  double pressure (const Point<dim> &p) const override;
82 
86  double density (const Point<dim> &p) const override;
87 
92  double density_derivative (const Point<dim> &p) const override;
93 
94  static
95  void
96  declare_parameters (ParameterHandler &prm);
97 
98  void
99  parse_parameters (ParameterHandler &prm) override;
100 
101  private:
109 
113  unsigned int n_points;
114 
119 
125  std::vector<double> temperatures;
126  std::vector<double> pressures;
127  std::vector<double> densities;
128 
133  double delta_z;
134 
141  std::shared_ptr<const aspect::InitialComposition::Manager<dim>> initial_composition_manager;
142 
147  double get_property (const Point<dim> &p,
148  const std::vector<double> &property) const;
149  };
150  }
151 }
152 
153 
154 #endif
double pressure(const Point< dim > &p) const override
double get_property(const Point< dim > &p, const std::vector< double > &property) const
double density(const Point< dim > &p) const override
static void declare_parameters(ParameterHandler &prm)
std::shared_ptr< const aspect::InitialComposition::Manager< dim > > initial_composition_manager
double temperature(const Point< dim > &p) const override
void parse_parameters(ParameterHandler &prm) override
double density_derivative(const Point< dim > &p) const override