From 33cad120c25be111fe9b3a0c79a716a8d7fc99bb Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 11 Nov 2022 14:01:45 +0100
Subject: [PATCH] typo for electron kernels. now working fine
---
src/numerics_mod.F90 | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/src/numerics_mod.F90 b/src/numerics_mod.F90
index 09fed922..1c769cae 100644
--- a/src/numerics_mod.F90
+++ b/src/numerics_mod.F90
@@ -69,8 +69,7 @@ DO iz = izgs,izge
j_int = jarray_e(ij)
j_dp = REAL(j_int,dp)
y_ = sigmae2_taue_o2 * kparray(iky,ikx,iz,eo)**2
- IF(j_int .LE. 0) THEN
- factj = 1._dp
+ IF(j_int .LT. 0) THEN
kernel_e(ij,iky,ikx,iz,eo) = 0._dp
ELSE
factj = GAMMA(j_dp+1._dp)
--
GitLab