cngf-pf

commit 06114430ffc3496532a3c7466bdf65c7073a7baa
parent 86bd3afebaef909928eafbecd1d52c892ae40116
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon, 30 Aug 2021 09:09:58 +0200

use magnitude of shear strain rate for computing inertia number

Diffstat:
M
simulation.c
 | 
2
+-

1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/simulation.c b/simulation.c
@@ -421,7 +421,7 @@ lithostatic_pressure_distribution(struct simulation *sim)
 inline static double
 inertia_number(double gamma_dot_p, double d, double sigma_n_eff, double rho_s)
 {
-	return gamma_dot_p * d / sqrt(sigma_n_eff / rho_s);
+	return fabs(gamma_dot_p) * d / sqrt(sigma_n_eff / rho_s);
 }
 
 void