From 072b720f20385d396fde39efde881b6fbe06f15d Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 16 Feb 2024 16:07:54 +0100
Subject: [PATCH] corrrect adiabatic ion model
---
src/solve_EM_fields.F90 | 5 +----
1 file changed, 1 insertion(+), 4 deletions(-)
diff --git a/src/solve_EM_fields.F90 b/src/solve_EM_fields.F90
index d28a4ac..dfbc9d2 100644
--- a/src/solve_EM_fields.F90
+++ b/src/solve_EM_fields.F90
@@ -69,10 +69,7 @@ CONTAINS
!!!!!!!!!!!!!!! adiabatic ion model
! Candy et al. 2007, rho_i = -q_i/tau_i phi
IF (ADIAB_I) THEN
- DO iz = 1,local_nz
- izi = iz+ngz/2
- rho(iz) = rho(iz) - 0*q_i/tau_i * phi(iky,ikx,izi)
- ENDDO
+ ! No ion contribution
ENDIF
!!!!!!!!!!!!!!! Inverting the poisson equation
--
GitLab