added stress output

geodynamics · Nov 10, 2023 · 0cfc1ab · 0cfc1ab
1 parent 3391ba8
commit 0cfc1ab
Show file tree

Hide file tree

Showing 5 changed files with 45,013 additions and 45,015 deletions.
diff --git a/source/postprocess/visualization/shear_stress.cc b/source/postprocess/visualization/shear_stress.cc
@@ -68,11 +68,10 @@ namespace aspect
         for (unsigned int q=0; q<n_quadrature_points; ++q)
           {
             // Compressive stress is negative by the sign convention
-            // used by the engineering community, and internally by ASPECT.
+            // used by the engineering community, and as input and used
+            // internally by ASPECT.
             // Here, we change the sign of the stress to match the
-            // sign convention used by the geoscience community,
-            // which is also the sign convention for stress expected
-            // in ASPECT parameter files.
+            // sign convention used by the geoscience community.
             SymmetricTensor<2,dim> shear_stress;
 
             // If elasticity is enabled, the deviatoric stress is stored

diff --git a/source/postprocess/visualization/stress.cc b/source/postprocess/visualization/stress.cc
@@ -68,11 +68,10 @@ namespace aspect
         for (unsigned int q=0; q<n_quadrature_points; ++q)
           {
             // Compressive stress is negative by the sign convention
-            // used by the engineering community, and internally by ASPECT.
+            // used by the engineering community, and as input and used
+            // internally by ASPECT.
             // Here, we change the sign of the stress to match the
-            // sign convention used by the geoscience community,
-            // which is also the sign convention for stress expected
-            // in ASPECT parameter files.
+            // sign convention used by the geoscience community.
             SymmetricTensor<2,dim> stress = in.pressure[q] * unit_symmetric_tensor<dim>();
 
             // If elasticity is enabled, the deviatoric stress is stored

diff --git a/tests/viscoelastic_stress_build-up.prm b/tests/viscoelastic_stress_build-up.prm
@@ -77,7 +77,7 @@ end
 subsection Postprocess
   subsection Visualization
     set Time between graphical output = 1e3
-    set List of output variables = shear stress, stress second invariant
+    set List of output variables = stress, shear stress, stress second invariant
     set Output format = gnuplot
   end
 end