From e80c03b0efcfdcf34c0c3f1476c572b07ddd345d Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Tue, 21 Jun 2022 17:08:08 +0200
Subject: [PATCH] typo for kinetic electrons nonlinear term
---
src/nonlinear_mod.F90 | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/nonlinear_mod.F90 b/src/nonlinear_mod.F90
index 0fd98081..90fe637b 100644
--- a/src/nonlinear_mod.F90
+++ b/src/nonlinear_mod.F90
@@ -86,13 +86,13 @@ SUBROUTINE compute_nonlinear
bracket_sum_r = 0._dp ! initialize sum over real nonlinear term
nloope: DO in = 1,nmax+1 ! Loop over laguerre for the sum
! First convolution terms
- F_cmpx(ikxs:ikxe,ikxs:ikxe) = phi(ikxs:ikxe,ikxs:ikxe,iz) * kernel_e(in, ikxs:ikxe,ikxs:ikxe, iz, eo)
+ F_cmpx(ikys:ikye,ikxs:ikxe) = phi(ikys:ikye,ikxs:ikxe,iz) * kernel_e(in, ikys:ikye,ikxs:ikxe, iz, eo)
! Second convolution terms
- G_cmpx(ikxs:ikxe,ikxs:ikxe) = 0._dp ! initialization of the sum
+ G_cmpx(ikys:ikye,ikxs:ikxe) = 0._dp ! initialization of the sum
smax = MIN( (in-1)+(ij-1), Jmaxe );
DO is = 1, smax+1 ! sum truncation on number of moments
- G_cmpx(ikxs:ikxe,ikxs:ikxe) = G_cmpx(ikxs:ikxe,ikxs:ikxe) + &
- dnjs(in,ij,is) * moments_e(ip,is,ikxs:ikxe,ikxs:ikxe,iz,updatetlevel)
+ G_cmpx(ikys:ikye,ikxs:ikxe) = G_cmpx(ikys:ikye,ikxs:ikxe) + &
+ dnjs(in,ij,is) * moments_e(ip,is,ikys:ikye,ikxs:ikxe,iz,updatetlevel)
ENDDO
!/!\ this function add its result to bracket_sum_r (hard to read sorry) /!\
CALL poisson_bracket_and_sum(F_cmpx,G_cmpx)
--
GitLab