From ef1f55d585fa510b08e9fe78c0d14a0062d1a6bb Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 24 Jun 2022 09:16:33 +0200
Subject: [PATCH] electron fluxes update
---
src/processing_mod.F90 | 18 +++++++++---------
1 file changed, 9 insertions(+), 9 deletions(-)
diff --git a/src/processing_mod.F90 b/src/processing_mod.F90
index 100496ec..f1b8eaad 100644
--- a/src/processing_mod.F90
+++ b/src/processing_mod.F90
@@ -166,11 +166,11 @@ END SUBROUTINE compute_radial_electron_transport
! 1D diagnostic to compute the average radial particle transport <T_i v_ExB_x>_xyz
SUBROUTINE compute_radial_ion_heatflux
- USE fields, ONLY : moments_i, moments_e, phi
- USE array, ONLY : kernel_e, kernel_i, HF_phi_correction_operator
+ USE fields, ONLY : moments_i, phi
+ USE array, ONLY : kernel_i, HF_phi_correction_operator
USE geometry, ONLY : Jacobian, iInt_Jacobian
USE time_integration, ONLY : updatetlevel
- USE model, ONLY : qe_taue, qi_taui, KIN_E, tau_i, tau_e
+ USE model, ONLY : qi_taui, KIN_E, tau_i
USE calculus, ONLY : simpson_rule_z
IMPLICIT NONE
COMPLEX(dp) :: hflux_local, integral
@@ -229,11 +229,11 @@ END SUBROUTINE compute_radial_ion_heatflux
! 1D diagnostic to compute the average radial particle transport <T_e v_ExB_x>_xyz
SUBROUTINE compute_radial_electron_heatflux
- USE fields, ONLY : moments_e, moments_e, phi
- USE array, ONLY : kernel_e, kernel_e, HF_phi_correction_operator
+ USE fields, ONLY : moments_e, phi
+ USE array, ONLY : kernel_e, HF_phi_correction_operator
USE geometry, ONLY : Jacobian, iInt_Jacobian
USE time_integration, ONLY : updatetlevel
- USE model, ONLY : qe_taue, qi_taui, KIN_E, tau_e, tau_e
+ USE model, ONLY : qe_taue, KIN_E, tau_e
USE calculus, ONLY : simpson_rule_z
IMPLICIT NONE
COMPLEX(dp) :: hflux_local, integral
@@ -258,9 +258,9 @@ SUBROUTINE compute_radial_electron_heatflux
-(n_dp+1._dp)*moments_e(ip0_e,in+1,iky,ikx,izs:ize,updatetlevel)&
-n_dp*moments_e(ip0_e,in-1,iky,ikx,izs:ize,updatetlevel))
ENDDO
- ! Add polarisation contribution
- integrant(izs:ize) = integrant(izs:ize) + Jacobian(izs:ize,0)*tau_e*imagu*ky_&
- *CONJG(phi(iky,ikx,izs:ize))*phi(iky,ikx,izs:ize) * HF_phi_correction_operator(iky,ikx,izs:ize)
+ ! ! Add polarisation contribution
+ ! integrant(izs:ize) = integrant(izs:ize) + Jacobian(izs:ize,0)*tau_e*imagu*ky_&
+ ! *CONJG(phi(iky,ikx,izs:ize))*phi(iky,ikx,izs:ize) * HF_phi_correction_operator(iky,ikx,izs:ize)
! Integrate over z
call simpson_rule_z(integrant,integral)
hflux_local = hflux_local + integral*iInt_Jacobian
--
GitLab