fix the units for the darcy velocity postprocessor

geodynamics · Nov 4, 2024 · ac71f79 · ac71f79
1 parent 4aa1942
commit ac71f79
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/source/postprocess/visualization/darcy_velocity.cc b/source/postprocess/visualization/darcy_velocity.cc
@@ -87,9 +87,9 @@ namespace aspect
             const double fluid_density = fluid_out->fluid_densities[q];
             const double permeability = fluid_out->permeabilities[q];
             const Tensor<1,dim> solid_velocity = in.velocity[q];
-            const Tensor<1,dim> darcy_velocity = solid_velocity -
+            const Tensor<1,dim> darcy_velocity = (solid_velocity -
                                                  permeability / fluid_viscosity / porosity * gravity *
-                                                 (solid_density - fluid_density) * velocity_scaling_factor;
+                                                 (solid_density - fluid_density)) * velocity_scaling_factor;
 
             for (unsigned int k=0; k<dim; ++k)
               computed_quantities[q](k) = darcy_velocity[k];