diff --git a/src/array_mod.F90 b/src/array_mod.F90
index acb0f640f16de12e7a0ff73540ca7f5956e22785..062c9debd3ea2fda9e8528d7b42556850dd4e52d 100644
--- a/src/array_mod.F90
+++ b/src/array_mod.F90
@@ -27,7 +27,7 @@ MODULE array
   ! non-zonal modes (kx=0)
   COMPLEX(dp), DIMENSION(:,:,:,:), ALLOCATABLE :: moments_e_NZ
   COMPLEX(dp), DIMENSION(:,:,:,:), ALLOCATABLE :: moments_i_NZ
-  COMPLEX(dp), DIMENSION(:,:),     ALLOCATABLE :: phi_EM
+  COMPLEX(dp), DIMENSION(:,:),     ALLOCATABLE :: phi_NZ
 
   ! Non linear term array (ip,ij,iky,ikx,iz)
   COMPLEX(dp), DIMENSION(:,:,:,:,:), ALLOCATABLE :: Sepj ! electron
diff --git a/src/memory.F90 b/src/memory.F90
index 86d063e5d8b381a46dbefa50a0d447cc0a952d19..06bf123b970d6c74cf0e556209203e8fd514fca9 100644
--- a/src/memory.F90
+++ b/src/memory.F90
@@ -16,7 +16,7 @@ SUBROUTINE memory
   CALL allocate_array(           phi, ikys,ikye, ikxs,ikxe, izgs,izge)
   CALL allocate_array(           psi, ikys,ikye, ikxs,ikxe, izgs,izge)
   CALL allocate_array(        phi_ZF, ikxs,ikxe, izs,ize)
-  CALL allocate_array(        phi_EM, ikys,ikye, izs,ize)
+  CALL allocate_array(        phi_NZ, ikys,ikye, izs,ize)
   CALL allocate_array(inv_poisson_op, ikys,ikye, ikxs,ikxe, izs,ize)
   CALL allocate_array( inv_ampere_op, ikys,ikye, ikxs,ikxe, izs,ize)
   CALL allocate_array(   inv_pol_ion, ikys,ikye, ikxs,ikxe, izs,ize)
diff --git a/src/numerical_experiments_mod.F90 b/src/numerical_experiments_mod.F90
index ce2b3ba32e78064905b0235e8186d914074dc116..c316a123293ff0215f067bc2becc3549f6a1bd5f 100644
--- a/src/numerical_experiments_mod.F90
+++ b/src/numerical_experiments_mod.F90
@@ -17,7 +17,7 @@ CONTAINS
 !******************************************************************************!
 SUBROUTINE save_EM_ZF_modes
   USE fields
-  USE array, ONLY : moments_e_ZF, moments_i_ZF, phi_ZF, moments_e_NZ,moments_i_NZ,phi_EM
+  USE array, ONLY : moments_e_ZF, moments_i_ZF, phi_ZF, moments_e_NZ,moments_i_NZ,phi_NZ
   USE grid
   USE time_integration, ONLY: updatetlevel
   USE model, ONLY: KIN_E
@@ -38,18 +38,18 @@ SUBROUTINE save_EM_ZF_modes
     IF(KIN_E) &
     moments_e_NZ(ips_e:ipe_e,ijs_e:ije_e,ikys:ikye,izs:ize) = moments_e(ips_e:ipe_e,ijs_e:ije_e,ikys:ikye,ikx_0,izs:ize,updatetlevel)
     moments_i_NZ(ips_i:ipe_i,ijs_i:ije_i,ikys:ikye,izs:ize) = moments_i(ips_i:ipe_i,ijs_i:ije_i,ikys:ikye,ikx_0,izs:ize,updatetlevel)
-    phi_EM(ikys:ikye,izs:ize) = phi(ikys:ikye,ikx_0,izs:ize)
+    phi_NZ(ikys:ikye,izs:ize) = phi(ikys:ikye,ikx_0,izs:ize)
   ELSE
     IF(KIN_E) &
     moments_e_NZ(ips_e:ipe_e,ijs_e:ije_e,ikys:ikye,izs:ize) = 0._dp
     moments_i_NZ(ips_e:ipe_e,ijs_i:ije_i,ikys:ikye,izs:ize) = 0._dp
-    phi_EM(ikys:ikye,izs:ize) = 0._dp
+    phi_NZ(ikys:ikye,izs:ize) = 0._dp
   ENDIF
 END SUBROUTINE
 
 SUBROUTINE play_with_modes
   USE fields
-  USE array, ONLY : moments_e_ZF, moments_i_ZF, phi_ZF, moments_e_NZ,moments_i_NZ,phi_EM
+  USE array, ONLY : moments_e_ZF, moments_i_ZF, phi_ZF, moments_e_NZ,moments_i_NZ,phi_NZ
   USE grid
   USE time_integration, ONLY: updatetlevel
   USE initial_par, ONLY: ACT_ON_MODES
@@ -100,7 +100,7 @@ SUBROUTINE play_with_modes
       IF(KIN_E) &
       moments_e(ips_e:ipe_e,ijs_e:ije_e,ikys:ikye,ikx_0,izs:ize,updatetlevel) = moments_e_NZ(ips_e:ipe_e,ijs_e:ije_e,ikys:ikye,izs:ize)
       moments_i(ips_i:ipe_i,ijs_i:ije_i,ikys:ikye,ikx_0,izs:ize,updatetlevel) = moments_i_NZ(ips_i:ipe_i,ijs_i:ije_i,ikys:ikye,izs:ize)
-      phi(ikys:ikye,ikx_0,izs:ize) = phi_EM(ikys:ikye,izs:ize)
+      phi(ikys:ikye,ikx_0,izs:ize) = phi_NZ(ikys:ikye,izs:ize)
     ENDIF
   CASE('amplify_zonal')
     IF(KIN_E) &