[cig-commits] [commit] master: disable pressure normalization for now (7da25c9)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/76d4275352ef2cae5de9a073acd1c03a92c2670c...4f3d06fd1f3754419813db37ec9ef7f0f6f3cb15

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

commit 7da25c96c42a2d5986428cac92af31f07eaadaeb
Author: Timo Heister <timo.heister at gmail.com>
Date:   Mon May 12 18:03:05 2014 -0400

    disable pressure normalization for now


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

7da25c96c42a2d5986428cac92af31f07eaadaeb
 source/simulator/helper_functions.cc | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/source/simulator/helper_functions.cc b/source/simulator/helper_functions.cc
index 51abcc0..3a2859a 100644
--- a/source/simulator/helper_functions.cc
+++ b/source/simulator/helper_functions.cc
@@ -535,6 +535,10 @@ namespace aspect
     if (parameters.pressure_normalization == "no")
       return;
 
+    // TODO: I need to think about how to implement this. Choose "no" for now.
+    Assert(introspection.block_indices.velocities != introspection.block_indices.pressure,
+        ExcNotImplemented());
+
     double my_pressure = 0.0;
     double my_area = 0.0;
     if (parameters.pressure_normalization == "surface")