[cig-commits] [commit] master: adjust a default value to account for the use of a different temperature in the original source of the melting model (623ff74)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/7a1bf75b781206bfe69ef7443f651e398e0c1522...9b2e4324e5cc6300267d41b2bfa35b70f9a44cd4

>---------------------------------------------------------------

commit 623ff746c27ea5164a19cbbbef9b528ab1604cd6
Author: Juliane Dannberg <dannberg at gfz-potsdam.de>
Date:   Fri Sep 19 18:52:25 2014 +0200

    adjust a default value to account for the use of a different temperature in the original source of the melting model


>---------------------------------------------------------------

623ff746c27ea5164a19cbbbef9b528ab1604cd6
 source/material_model/latent_heat_melt.cc         | 6 +++++-
 source/postprocess/visualization/melt_fraction.cc | 6 +++++-
 2 files changed, 10 insertions(+), 2 deletions(-)

diff --git a/source/material_model/latent_heat_melt.cc b/source/material_model/latent_heat_melt.cc
index b48cb67..14a933c 100644
--- a/source/material_model/latent_heat_melt.cc
+++ b/source/material_model/latent_heat_melt.cc
@@ -607,11 +607,15 @@ namespace aspect
                              "function that approximates the solidus "
                              "of pyroxenite. "
                              "Units: $°C$.");
-          prm.declare_entry ("D2", "1.23e-7",
+          prm.declare_entry ("D2", "1.329e-7",
                              Patterns::Double (),
                              "Prefactor of the linear pressure term "
                              "in the quadratic function that approximates "
                              "the solidus of pyroxenite. "
+                             "Note that this factor is different from the "
+                             "value given in Sobolev, 2011, because they use "
+                             "the potential temperature whereas we use the "
+                             "absolute temperature. "
                              "Units: $°C/Pa$.");
           prm.declare_entry ("D3", "-5.1e-18",
                              Patterns::Double (),
diff --git a/source/postprocess/visualization/melt_fraction.cc b/source/postprocess/visualization/melt_fraction.cc
index d055bba..60d1b07 100644
--- a/source/postprocess/visualization/melt_fraction.cc
+++ b/source/postprocess/visualization/melt_fraction.cc
@@ -221,11 +221,15 @@ namespace aspect
                                  "function that approximates the solidus "
                                  "of pyroxenite. "
                                  "Units: $°C$.");
-              prm.declare_entry ("D2", "1.23e-7",
+              prm.declare_entry ("D2", "1.329e-7",
                                  Patterns::Double (),
                                  "Prefactor of the linear pressure term "
                                  "in the quadratic function that approximates "
                                  "the solidus of pyroxenite. "
+                                 "Note that this factor is different from the "
+                                 "value given in Sobolev, 2011, because they use "
+                                 "the potential temperature whereas we use the "
+                                 "absolute temperature. "
                                  "Units: $°C/Pa$.");
               prm.declare_entry ("D3", "-5.1e-18",
                                  Patterns::Double (),