From 5c727b907e8ad0bd8f436ebf79a72d2c1c55ceea Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Mon, 11 Jan 2021 10:33:48 +0100
Subject: [PATCH] Debugging 1D linear
---
src/moments_eq_rhs.F90 | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/src/moments_eq_rhs.F90 b/src/moments_eq_rhs.F90
index 8e25ed8f..bdbfc06e 100644
--- a/src/moments_eq_rhs.F90
+++ b/src/moments_eq_rhs.F90
@@ -127,6 +127,8 @@ SUBROUTINE moments_eq_rhs
kzloope : DO ikz = ikzs,ikze
kr = krarray(ikr) ! Poloidal wavevector
kz = kzarray(ikz) ! Toroidal wavevector
+ IF (Nkz .EQ. 1) kz = krarray(ikr) ! If 1D simulation we put kr as kz
+
kperp2 = kr**2 + kz**2 ! perpendicular wavevector
be_2 = kperp2 * sigmae2_taue_o2 ! Kernel argument
@@ -338,6 +340,8 @@ SUBROUTINE moments_eq_rhs
kzloopi : DO ikz = ikzs,ikze
kr = krarray(ikr) ! Poloidal wavevector
kz = kzarray(ikz) ! Toroidal wavevector
+ IF (Nkz .EQ. 1) kz = krarray(ikr) ! If 1D simulation we put kr as kz
+
kperp2 = kr**2 + kz**2 ! perpendicular wavevector
bi_2 = kperp2 * sigmai2_taui_o2 ! Kernel argument
--
GitLab