commented lines for a possible parallel drive implementation

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
 

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)