From f6c793cbb60eaae844dd4b8075790fe905ec159a Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Mon, 14 Nov 2022 10:23:05 +0100
Subject: [PATCH] doubt on the NL geometric prefactor
---
src/geometry_mod.F90 | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/src/geometry_mod.F90 b/src/geometry_mod.F90
index 84bac735..83a972fd 100644
--- a/src/geometry_mod.F90
+++ b/src/geometry_mod.F90
@@ -139,9 +139,6 @@ CONTAINS
ENDDO
! coefficient in the front of parallel derivative
gradz_coeff(iz,eo) = 1._dp / jacobian(iz,eo) / hatB(iz,eo)
- ! Factor in front of the nonlinear term
- hatB_NL(iz,eo) = Jacobian(iz,eo)&
- *(gxx(iz,eo)*gyy(iz,eo) - gxy(iz,eo)**2)/hatB(iz,eo)
ENDDO
ENDDO
@@ -213,6 +210,8 @@ CONTAINS
ENDDO
! coefficient in the front of parallel derivative
gradz_coeff(iz,eo) = 1._dp / Jacobian(iz,eo) / hatB(iz,eo)
+ ! Nonlinear term prefactor
+ hatB_NL(iz,eo) = 1._dp !Jacobian(iz,eo)*(gxx(iz,eo)*gyy(iz,eo) - gxy(iz,eo)**2)/hatB(iz,eo)
ENDDO zloop
ENDDO parity
--
GitLab