replaced hard coded damper index

geodynamics · Aug 11, 2024 · c1f4134 · c1f4134
1 parent ada9989
commit c1f4134
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/source/material_model/rheology/composite_visco_plastic.cc b/source/material_model/rheology/composite_visco_plastic.cc
@@ -657,7 +657,7 @@ namespace aspect
         // arranged in parallel with the viscoplastic elements.
         // The total stress is equal to the sum of the viscoplastic stress and
         // damper stress.
-        const double damper_stress = 2. * damper_viscosity * (edot_ii - partial_strain_rates[4]);
+        const double damper_stress = 2. * damper_viscosity * (edot_ii - partial_strain_rates[isostrain_damper_strain_rate_index]);
         const double total_stress = viscoplastic_stress + damper_stress;
 
         // 6) Return the effective creep viscosity using the total stress