diff --git a/DOUAR.f90 b/DOUAR.f90
index 7b335de6051142c913d103db776c1c548e546bbf..f28dcbf8c076d335dc1a449c018323db948ddb89 100644
--- a/DOUAR.f90
+++ b/DOUAR.f90
@@ -729,7 +729,7 @@ do while (istep.le.params%nstep)
             call show_time (total,step,inc,1,'write global output$')
             call write_global_output(params,istep,iter,current_time,osolve,ov,vo,surface,cl,'debug')
             call mpi_barrier (mpi_comm_world,ierr)
-            call slices (params,osolve,ov,vo,sections,istep,iter,iter_nl)
+!            call slices (params,osolve,ov,vo,sections,istep,iter,iter_nl)
          end if
 
          !--------------------------------------------------------------------------------------
@@ -768,7 +768,7 @@ do while (istep.le.params%nstep)
       !--------------------------------------------------------------------------------------
 
       call show_time (total,step,inc,1,'slicing the cube$')
-      call slices(params,osolve,ov,vo,sections,istep,iter,iter_nl)
+!      call slices(params,osolve,ov,vo,sections,istep,iter,iter_nl)
 
       call heap (threadinfo,'osolve%eviscosity','main',size(osolve%eviscosity),'dp',-1) ; deallocate (osolve%eviscosity)
       call heap (threadinfo,'osolve%q','main',size(osolve%q),'dp',-1)                   ; deallocate (osolve%q)