[cig-commits] r1331 - trunk/aspect/source/postprocess

[heister at dealii.org]

Author: heister
Date: 2012-10-25 09:32:57 -0600 (Thu, 25 Oct 2012)
New Revision: 1331

Modified:
   trunk/aspect/source/postprocess/table_velocity_statistics.cc
   trunk/aspect/source/postprocess/velocity_statistics.cc
Log:
fix output of gravity in velocity statistics

Modified: trunk/aspect/source/postprocess/table_velocity_statistics.cc
===================================================================
--- trunk/aspect/source/postprocess/table_velocity_statistics.cc	2012-10-25 15:20:16 UTC (rev 1330)
+++ trunk/aspect/source/postprocess/table_velocity_statistics.cc	2012-10-25 15:32:57 UTC (rev 1331)
@@ -168,7 +168,7 @@
           // we do not compute the compositions but give the functions below the value 0.0 instead
           std::vector<double> composition_values(this->n_compositional_fields(),0.0);
 
-          Point<dim> representative_point = Point<dim>::unit_vector(dim-1);
+          Point<dim> representative_point = this->get_geometry_model().representative_point(0);
           const double gravity = this->get_gravity_model().gravity_vector(representative_point).norm();
           const double Ra = material_model.reference_density()*
                             gravity*

Modified: trunk/aspect/source/postprocess/velocity_statistics.cc
===================================================================
--- trunk/aspect/source/postprocess/velocity_statistics.cc	2012-10-25 15:20:16 UTC (rev 1330)
+++ trunk/aspect/source/postprocess/velocity_statistics.cc	2012-10-25 15:32:57 UTC (rev 1331)
@@ -151,7 +151,7 @@
               // we do not compute the compositions but give the functions below the value 0.0 instead
               std::vector<double> composition_values(this->n_compositional_fields(),0.0);
 
-              Point<dim> representative_point = Point<dim>::unit_vector(dim-1);
+              Point<dim> representative_point = this->get_geometry_model().representative_point(0);
               const double gravity = this->get_gravity_model().gravity_vector(representative_point).norm();
               const double Ra = material_model.reference_density()*
                                 gravity*