From 9bf8acffde347e271764eb460a0895b5be0b9a5e Mon Sep 17 00:00:00 2001
From: asjeb <151923475+asjeb@users.noreply.github.com>
Date: Wed, 12 Jun 2024 10:23:50 +0200
Subject: [PATCH] Update smash/fcore/operator/md_gr_operator.f90

simple precision

Co-authored-by: Francois Colleoni <110899888+inoelloc@users.noreply.github.com>
---
 smash/fcore/operator/md_gr_operator.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/smash/fcore/operator/md_gr_operator.f90 b/smash/fcore/operator/md_gr_operator.f90
index 99fee787..721d4ee9 100644
--- 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