From f7125b3d4df45799e48a8244476cc9d4cbc8e0d0 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Wed, 19 Oct 2022 18:42:52 +0200
Subject: [PATCH] comment the correction operator for heat flux
---
src/processing_mod.F90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/processing_mod.F90 b/src/processing_mod.F90
index 7a4b2783..2c347928 100644
--- a/src/processing_mod.F90
+++ b/src/processing_mod.F90
@@ -370,8 +370,8 @@ SUBROUTINE compute_radial_electron_heatflux
ENDDO
ENDIF
! Add polarisation contribution
- integrant(izs:ize) = integrant(izs:ize) + tau_e*imagu*ky_&
- *CONJG(phi(iky,ikx,izs:ize))*phi(iky,ikx,izs:ize) * HF_phi_correction_operator(iky,ikx,izs:ize)
+ ! integrant(izs:ize) = integrant(izs:ize) + 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
integrant(izgs:izge) = Jacobian(izgs:izge,0)*integrant(izgs:izge)
call simpson_rule_z(integrant,integral)
--
GitLab