diff --git a/src/processing_mod.F90 b/src/processing_mod.F90
index 7a4b2783b0bef6e32636279a27e3344605ef1237..2c34792828e2c1564c921dee25d05a124c74b2a9 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)