[cig-commits] [commit] master: change depth_average zones from 100 to 10 (d75f6a6)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/925399ce4de0e953aa20de00940677b11013e5bc...327ebd4b084f4d740908b7fe9ca583f6420d2874

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

commit d75f6a617872de113eba838c4ba853f865f1c91f
Author: Timo Heister <timo.heister at gmail.com>
Date:   Sun Jan 18 17:43:07 2015 -0500

    change depth_average zones from 100 to 10
    
    In each cell we place 10 quadrature points, which means 100 zones
    require probably 50 cells in depth direction, which typical computations
    do not have.


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

d75f6a617872de113eba838c4ba853f865f1c91f
 doc/modules/changes.h               | 6 ++++++
 source/postprocess/depth_average.cc | 4 ++--
 2 files changed, 8 insertions(+), 2 deletions(-)

diff --git a/doc/modules/changes.h b/doc/modules/changes.h
index d1c691f..2d8f84a 100644
--- a/doc/modules/changes.h
+++ b/doc/modules/changes.h
@@ -6,6 +6,12 @@
  *
  *
  * <ol>
+ * <li> The default for the number of zones in "depth average" got reduced
+ * from 100 to 10. A large number of zones causes erroneous output on coarse
+ * meshes.
+ * <br>
+ * (Timo Heister, 2015/01/19)
+ *
  * <li> New: ASPECT now uses meshes for the spherical shell where each node is
  * placed on concentric shells. Previously, upon mesh refinement, new vertices
  * were placed at an averaged location of the vertices of the mother cell,
diff --git a/source/postprocess/depth_average.cc b/source/postprocess/depth_average.cc
index 072e46f..4edd05f 100644
--- a/source/postprocess/depth_average.cc
+++ b/source/postprocess/depth_average.cc
@@ -197,11 +197,11 @@ namespace aspect
                              "Units: years if the "
                              "'Use years in output instead of seconds' parameter is set; "
                              "seconds otherwise.");
-          prm.declare_entry ("Number of zones", "100",
+          prm.declare_entry ("Number of zones", "10",
                              Patterns::Integer (1),
                              "The number of zones in depth direction within which we "
                              "are to compute averages. By default, we subdivide the entire "
-                             "domain into 100 depth zones and compute temperature and other "
+                             "domain into 10 depth zones and compute temperature and other "
                              "averages within each of these zones. However, if you have a "
                              "very coarse mesh, it may not make much sense to subdivide "
                              "the domain into so many zones and you may wish to choose "