diff --git a/src/inital.F90 b/src/inital.F90
index 612980513ea9c61e1b8d737b089e6c0947afab4f..011e3bef94aae8112edf6de08a2a6dff8279b365 100644
--- a/src/inital.F90
+++ b/src/inital.F90
@@ -10,7 +10,6 @@ SUBROUTINE inital
   USE closure,          ONLY: apply_closure_model
   USE ghosts,           ONLY: update_ghosts_moments, update_ghosts_EM
   USE restarts,         ONLY: load_moments, job2load
-  USE numerics,         ONLY: play_with_modes, save_EM_ZF_modes
   USE processing,       ONLY: compute_fluid_moments
   USE model,            ONLY: KIN_E, LINEARITY
   USE nonlinear,        ONLY: compute_Sapj, nonlinear_init
@@ -73,11 +72,6 @@ SUBROUTINE inital
   CALL update_ghosts_EM
   !! End of phi and moments initialization
 
-  ! Save (kx,0) and (0,ky) modes for num exp
-  CALL save_EM_ZF_modes
-  ! Freeze/Wipe some selected modes (entropy,zonal,turbulent)
-  CALL play_with_modes
-
   ! Load the COSOlver collision operator coefficients
   IF(cosolver_coll) &
   CALL load_COSOlver_mat