[cig-commits] [commit] master: Add comment (b728252)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Tue May 20 11:26:39 PDT 2014


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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/19fcaae1dae23206c4e6f1f4e9965965f0128967...6204d1b7337a0741858189f8e4ee9cf2f447ebd6

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

commit b728252f230b4d87fb13cc0440cda232b3c39e4b
Author: anne-glerum <a.c.glerum at uu.nl>
Date:   Tue May 20 00:06:34 2014 +0200

    Add comment


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

b728252f230b4d87fb13cc0440cda232b3c39e4b
 source/postprocess/spherical_velocity_statistics.cc | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/source/postprocess/spherical_velocity_statistics.cc b/source/postprocess/spherical_velocity_statistics.cc
index 592c47d..5797f37 100644
--- a/source/postprocess/spherical_velocity_statistics.cc
+++ b/source/postprocess/spherical_velocity_statistics.cc
@@ -86,12 +86,13 @@ namespace aspect
                                                    fe_values.JxW(q);
               }
           }
-
+      // compute the global sums
       const double global_rad_velocity_square_integral
         = Utilities::MPI::sum (local_rad_velocity_square_integral, this->get_mpi_communicator());
       const double global_tan_velocity_square_integral
         = Utilities::MPI::sum (local_tan_velocity_square_integral, this->get_mpi_communicator());
 
+      // compute the final output by dividing by the volume over which we integrated
       const double rad_vrms = std::sqrt(global_rad_velocity_square_integral) /
                               std::sqrt(this->get_volume());
       const double tan_vrms = std::sqrt(global_tan_velocity_square_integral) /



More information about the CIG-COMMITS mailing list