World Builder  1.1.0-pre
A geodynamic initial conditions generator
smooth.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2018-2024 by the authors of the World Builder code.
3 
4  This file is part of the World Builder.
5 
6  This program is free software: you can redistribute it and/or modify
7  it under the terms of the GNU Lesser General Public License as published
8  by the Free Software Foundation, either version 2 of the License, or
9  (at your option) any later version.
10 
11  This program 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 Lesser General Public License for more details.
15 
16  You should have received a copy of the GNU Lesser General Public License
17  along with this program. If not, see <https://www.gnu.org/licenses/>.
18 */
19 
20 #ifndef _world_builder_features_subducting_plate_composition_smooth_h
21 #define _world_builder_features_subducting_plate_composition_smooth_h
22 
25 #include <world_builder/world.h>
26 
27 
28 namespace WorldBuilder
29 {
30  namespace Features
31  {
32  using namespace FeatureUtilities;
33  namespace SubductingPlateModels
34  {
35  namespace Composition
36  {
43  class Smooth: public Interface
44  {
45  public:
50 
54  ~Smooth();
55 
59  static
60  void declare_entries(Parameters &prm, const std::string &parent_name = "");
61 
65  void parse_entries(Parameters &prm) override final;
66 
67 
72  double get_composition(const Point<3> &position,
73  const double depth,
74  const unsigned int composition_number,
75  double composition,
76  const double feature_min_depth,
77  const double feature_max_depth,
79  const AdditionalParameters &additional_parameters) const override final;
80 
81  private:
82  // linear temperature submodule parameters
83  double min_distance;
84  double max_distance;
85  double side_distance;
86  std::vector<double> top_fraction;
87  // currently not using the side composition, but maybe usefu if you want another composition towards the end
88  std::vector<double> bottom_fraction;
89  std::vector<unsigned int> compositions;
91 
92  };
93  }
94  }
95  }
96 }
97 
98 #endif