ASPECT
melt_fraction.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2013 - 2021 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_postprocess_visualization_melt_fraction_h
23 #define _aspect_postprocess_visualization_melt_fraction_h
24 
27 
28 #include <deal.II/numerics/data_postprocessor.h>
29 
30 
31 namespace aspect
32 {
33  namespace Postprocess
34  {
35  namespace VisualizationPostprocessors
36  {
45  template <int dim>
47  : public DataPostprocessorScalar<dim>,
48  public SimulatorAccess<dim>,
49  public Interface<dim>
50  {
51  public:
52  MeltFraction ();
53 
54  void
55  evaluate_vector_field(const DataPostprocessorInputs::Vector<dim> &input_data,
56  std::vector<Vector<double>> &computed_quantities) const override;
57 
61  static
62  void
63  declare_parameters (ParameterHandler &prm);
64 
68  void
69  parse_parameters (ParameterHandler &prm) override;
70 
71  private:
76  // for the solidus temperature
77  double A1; // °C
78  double A2; // °C/Pa
79  double A3; // °C/(Pa^2)
80 
81  // for the lherzolite liquidus temperature
82  double B1; // °C
83  double B2; // °C/Pa
84  double B3; // °C/(Pa^2)
85 
86  // for the liquidus temperature
87  double C1; // °C
88  double C2; // °C/Pa
89  double C3; // °C/(Pa^2)
90 
91  // for the reaction coefficient of pyroxene
92  double r1; // cpx/melt
93  double r2; // cpx/melt/GPa
94  double M_cpx; // mass fraction of pyroxenite
95 
96  // melt fraction exponent
97  double beta;
98 
103  // for the melting temperature
104  double D1; // °C
105  double D2; // °C/Pa
106  double D3; // °C/(Pa^2)
107 
108  // for the melt-fraction dependence of productivity
109  double E1;
110  double E2;
111  };
112  }
113  }
114 }
115 
116 #endif
void evaluate_vector_field(const DataPostprocessorInputs::Vector< dim > &input_data, std::vector< Vector< double >> &computed_quantities) const override
Definition: compat.h:59