Add cbf heat flux

doc/make_cite_html.py

@@ -77,7 +77,8 @@
                "dannberg:etal:2022" : "entropy",
                "dannberg:grain:size" : "grainsize",
                "dannberg:etal:2021" : "boukaremelt",
-               "neuharth:etal:2022" : "fastscape"
+               "neuharth:etal:2022" : "fastscape",
+               "dannberg:gassmoeller:etal:2024" : "cbfheatflux"
 }
 
 want = []
@@ -114,7 +115,8 @@
         "dannberg:etal:2022" : "Dannberg, J., Gassmöller, R., Li, R., Lithgow-Bertelloni, C. and Stixrude, L. 2022. An entropy method for geodynamic modelling of phase transitions: capturing sharp and broad transitions in a multiphase assemblage. Geophysical Journal International, 231(3), pp.1833-1849. doi:10.1093/gji/ggac293. https://doi.org/10.1093/gji/ggac293.",
         "dannberg:grain:size" : "Dannberg, J., Eilon, Z., Faul, U., Gassmöller, R., Moulik, P. and Myhill, R. 2017. The importance of grain size to mantle dynamics and seismological observations. Geochemistry, Geophysics, Geosystems, 18(8), pp.3034-3061. doi:10.1002/2017GC006944. https://doi.org/10.1002/2017GC006944.",
         "dannberg:etal:2021" : "Dannberg, J., Myhill, R., Gassmöller, R. and Cottaar, S., 2021. The morphology, evolution and seismic visibility of partial melt at the core–mantle boundary: implications for ULVZs. Geophysical Journal International, 227(2), pp.1028-1059.",
-        "neuharth:etal:2022" : "Neuharth, D., Brune, S., Wrona, T., Glerum, A., Braun, J. and Yuan, X. 2022. Evolution of rift systems and their fault networks in response to surface processes. Tectonics, 41(3), e2021TC007166. doi:10.1029/2021TC007166. https://doi.org/10.1029/2021TC007166."
+        "neuharth:etal:2022" : "Neuharth, D., Brune, S., Wrona, T., Glerum, A., Braun, J. and Yuan, X. 2022. Evolution of rift systems and their fault networks in response to surface processes. Tectonics, 41(3), e2021TC007166. doi:10.1029/2021TC007166. https://doi.org/10.1029/2021TC007166.",
+        "dannberg:gassmoeller:etal:2024" : "Dannberg, J., Gassmöller, R., Thallner, D., LaCombe, F. and Sprain, C. 2024. Changes in core–mantle boundary heat flux patterns throughout the supercontinent cycle. Geophysical Journal International, 237(3), pp.1251-1274. doi:10.1093/gji/ggae075. https://doi.org/10.1093/gji/ggae075."
 }
 
 

doc/sphinx/references.bib

@@ -11986,6 +11986,21 @@ @article{Liu2019
   journal = {Geophysical Journal International}
 }
 
+@article{dannberg:gassmoeller:etal:2024,
+author = {Dannberg, Juliane and Gassmöller, Rene and Thallner, Daniele and LaCombe, Frederick and Sprain, Courtney},
+title = "{Changes in core-mantle boundary heat flux patterns throughout the supercontinent cycle}",
+journal = {Geophysical Journal International},
+volume = {237},
+number = {3},
+pages = {1251-1274},
+year = {2024},
+month = {02},
+issn = {0956-540X},
+doi = {10.1093/gji/ggae075},
+url = {https://doi.org/10.1093/gji/ggae075},
+eprint = {https://academic.oup.com/gji/article-pdf/237/3/1251/57214885/ggae075.pdf}
+}
+
 @Article{Fraters2019c,
 AUTHOR = {Fraters, M. and Thieulot, C. and van den Berg, A. and Spakman, W.},
 TITLE = {The Geodynamic World Builder: a solution for complex initial conditions in numerical modeling},

include/aspect/postprocess/heat_flux_densities.h

@@ -41,6 +41,12 @@ namespace aspect
     class HeatFluxDensities : public Interface<dim>, public ::aspect::SimulatorAccess<dim>
     {
       public:
+        /**
+         * Initialize the postprocessor.
+         */
+        void
+        initialize() override;
+
         /**
          * Evaluate.
          */

include/aspect/postprocess/heat_flux_map.h

@@ -41,6 +41,13 @@ namespace aspect
        * The consistent Galerkin FEM for computing derived boundary quantities in thermal and or fluids
        * problems. International Journal for Numerical Methods in Fluids, 7(4), 371-394.
        *
+       * The implementation of the method is benchmarked in
+       *
+       * Juliane Dannberg, Rene Gassmöller, Daniele Thallner, Frederick LaCombe, Courtney Sprain (2024).
+       * Changes in core–mantle boundary heat flux patterns throughout the supercontinent cycle,
+       * Geophysical Journal International, Volume 237, Issue 3, June 2024, Pages 1251–1274,
+       * https://doi.org/10.1093/gji/ggae075.
+       *
        * In summary, the method solves the temperature equation again on the boundary faces, with known
        * temperatures and solving for the boundary fluxes that satisfy the equation. Since the
        * equation is only formed on the faces and it can be solved using only diagonal matrices,
@@ -83,6 +90,12 @@ namespace aspect
     class HeatFluxMap : public Interface<dim>, public ::aspect::SimulatorAccess<dim>
     {
       public:
+        /**
+         * Initialize the postprocessor.
+         */
+        void
+        initialize() override;
+
         /**
          * Evaluate the solution for the heat flux.
          */

include/aspect/postprocess/heat_flux_statistics.h

@@ -40,6 +40,12 @@ namespace aspect
     class HeatFluxStatistics : public Interface<dim>, public ::aspect::SimulatorAccess<dim>
     {
       public:
+        /**
+         * Initialize the postprocessor.
+         */
+        void
+        initialize() override;
+
         /**
          * Evaluate the solution for some heat_flux statistics.
          */

include/aspect/postprocess/visualization/heat_flux_map.h

@@ -44,8 +44,17 @@ namespace aspect
           public Interface<dim>
       {
         public:
+          /**
+           * Constructor.
+           */
           HeatFluxMap();
 
+          /**
+           * Initialize the postprocessor.
+           */
+          void
+          initialize() override;
+
           /**
            * Fill the temporary storage variables with the
            * heat flux for the current time step.

source/postprocess/dynamic_topography.cc

@@ -505,6 +505,8 @@ namespace aspect
     void
     DynamicTopography<dim>::parse_parameters (ParameterHandler &prm)
     {
+      CitationInfo::add("geoid");
+
       prm.enter_subsection("Postprocess");
       {
         prm.enter_subsection("Dynamic topography");

source/postprocess/heat_flux_densities.cc

@@ -33,6 +33,15 @@ namespace aspect
 {
   namespace Postprocess
   {
+    template <int dim>
+    void
+    HeatFluxDensities<dim>::initialize ()
+    {
+      CitationInfo::add("cbfheatflux");
+    }
+
+
+
     template <int dim>
     std::pair<std::string,std::string>
     HeatFluxDensities<dim>::execute (TableHandler &statistics)

source/postprocess/heat_flux_map.cc

@@ -452,6 +452,13 @@ namespace aspect
       }
     }
 
+    template <int dim>
+    void
+    HeatFluxMap<dim>::initialize ()
+    {
+      CitationInfo::add("cbfheatflux");
+    }
+
 
 
     template <int dim>

source/postprocess/heat_flux_statistics.cc

@@ -33,6 +33,15 @@ namespace aspect
 {
   namespace Postprocess
   {
+    template <int dim>
+    void
+    HeatFluxStatistics<dim>::initialize ()
+    {
+      CitationInfo::add("cbfheatflux");
+    }
+
+
+
     template <int dim>
     std::pair<std::string,std::string>
     HeatFluxStatistics<dim>::execute (TableHandler &statistics)

source/postprocess/visualization/heat_flux_map.cc

@@ -44,6 +44,15 @@ namespace aspect
 
 
 
+      template <int dim>
+      void
+      HeatFluxMap<dim>::initialize ()
+      {
+        CitationInfo::add("cbfheatflux");
+      }
+
+
+
       template <int dim>
       void
       HeatFluxMap<dim>::update ()