Update smash/fcore/operator/md_gr_operator.f90

simple precision Co-authored-by: Francois Colleoni <[email protected]>
DassHydro · Jun 12, 2024 · 9bf8acf · 9bf8acf
1 parent 96d63fb
commit 9bf8acf
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/smash/fcore/operator/md_gr_operator.f90 b/smash/fcore/operator/md_gr_operator.f90
@@ -94,7 +94,7 @@ subroutine hortonian_production(pn, en, cp, beta, alpha1, hp, pr, perc, dt)
 
         inv_cp = 1._sp / cp
         pr = 0._sp
-        gam = 1 - exp(-pn * alpha1)
+        gam = 1._sp - exp(-pn * alpha1)
         lambda = sqrt(1 - gam)
         inv_lambda = 1._sp / lambda