[cig-commits] commit 2406 by bangerth to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 2406

Leave 2 TODOs where we currently use the wrong formula.

U   trunk/aspect/source/postprocess/dynamic_topography.cc
U   trunk/aspect/source/postprocess/visualization/dynamic_topography.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2406&peg=2406

Diff:
Modified: trunk/aspect/source/postprocess/dynamic_topography.cc
===================================================================
--- trunk/aspect/source/postprocess/dynamic_topography.cc	2014-04-04 03:28:25 UTC (rev 2405)
+++ trunk/aspect/source/postprocess/dynamic_topography.cc	2014-04-04 03:40:38 UTC (rev 2406)
@@ -99,6 +99,7 @@
                   const double viscosity = out.viscosities[q];
                   const double density   = out.densities[q];
 
+//TODO: We need to subtract 2/3*div(u) from the stress here in the compressible case
                   const SymmetricTensor<2,dim> stress = 2 * viscosity * in.strain_rate[q];
 
                   const Tensor<1,dim> gravity = this->get_gravity_model().gravity_vector(location);

Modified: trunk/aspect/source/postprocess/visualization/dynamic_topography.cc
===================================================================
--- trunk/aspect/source/postprocess/visualization/dynamic_topography.cc	2014-04-04 03:28:25 UTC (rev 2405)
+++ trunk/aspect/source/postprocess/visualization/dynamic_topography.cc	2014-04-04 03:40:38 UTC (rev 2406)
@@ -91,6 +91,7 @@
             const double viscosity = out.viscosities[q];
             const double density   = out.densities[q];
 
+//TODO: We need to subtract 2/3*div(u) from the stress here in the compressible case
             const SymmetricTensor<2,dim> stress = 2 * viscosity * in.strain_rate[q];
 
             const Tensor<1,dim> gravity = this->get_gravity_model().gravity_vector(location);