diff --git a/src/moments_eq_rhs_mod.F90 b/src/moments_eq_rhs_mod.F90
index 66d4287da98f9453c634bd7cc57da513a9d6b27f..aad64e460ed8064f396e086d0a1a596fb197087e 100644
--- a/src/moments_eq_rhs_mod.F90
+++ b/src/moments_eq_rhs_mod.F90
@@ -12,7 +12,7 @@ SUBROUTINE compute_moments_eq_rhs
   USE prec_const
   USE collision
   USE time_integration
-  USE geometry, ONLY: gradz_coeff, dBdz, Ckxky, Gamma_NL, Gamma_phipar
+  USE geometry, ONLY: gradz_coeff, dBdz, Ckxky, Gamma_NL!, Gamma_phipar
   USE calculus, ONLY : interp_z, grad_z, grad_z2
   IMPLICIT NONE
 
diff --git a/src/processing_mod.F90 b/src/processing_mod.F90
index 6440fbef5e4050e7b84062781723871e29a106fe..aceca8249b08f257faef0a77c49a96afe2ba2f70 100644
--- a/src/processing_mod.F90
+++ b/src/processing_mod.F90
@@ -406,7 +406,7 @@ SUBROUTINE compute_nadiab_moments_z_gradients_and_interp
   USE fields,           ONLY : moments_i, moments_e, phi, psi
   USE array,            ONLY : kernel_e, kernel_i, nadiab_moments_e, nadiab_moments_i, &
                                ddz_nepj, ddzND_nepj, interp_nepj,&
-                               ddz_nipj, ddzND_nipj, interp_nipj, ddz_phi
+                               ddz_nipj, ddzND_nipj, interp_nipj!, ddz_phi
   USE time_integration, ONLY : updatetlevel
   USE model,            ONLY : qe_taue, qi_taui,q_o_sqrt_tau_sigma_e, q_o_sqrt_tau_sigma_i, &
                                KIN_E, CLOS, beta
@@ -517,12 +517,12 @@ SUBROUTINE compute_nadiab_moments_z_gradients_and_interp
     ENDDO
   ENDDO
 
-  ! Phi parallel gradient (experimental, should be negligible)
-  DO ikx = ikxs,ikxe
-    DO iky = ikys,ikye
-      CALL grad_z(0,phi(iky,ikx,izgs:izge), ddz_phi(iky,ikx,izs:ize))
-    ENDDO
-  ENDDO
+  ! Phi parallel gradient (not implemented fully, should be negligible)
+  ! DO ikx = ikxs,ikxe
+  !   DO iky = ikys,ikye
+  !     CALL grad_z(0,phi(iky,ikx,izgs:izge), ddz_phi(iky,ikx,izs:ize))
+  !   ENDDO
+  ! ENDDO
 
   ! Execution time end
   CALL cpu_time(t1_process)