diff --git a/Makefile_wsl_debian b/Makefile_wsl_debian
new file mode 100644
index 0000000000000000000000000000000000000000..36c1754b3998fd63b9ad059ee6c52d9eb03cb386
--- /dev/null
+++ b/Makefile_wsl_debian
@@ -0,0 +1,290 @@
+include local/dirs.inc
+include local/make.inc
+#Different namings depending on the make input
+EXEC = $(BINDIR)/gyacomo				#all
+EFST = $(BINDIR)/gyacomo_fast		#fast
+EDBG = $(BINDIR)/gyacomo_debug	#debug
+EALP = $(BINDIR)/gyacomo_alpha  #alpha
+
+# F90 = mpiifort
+F90 = mpif90
+# #F90 = ftn #for piz-daint cluster
+# # Add Multiple-Precision Library
+# EXTLIBS += -L$(FMDIR)/lib
+# EXTINC += -I$(FMDIR)/mod
+# # Add local fftw dir
+# EXTLIBS += -L$(FFTWDIR)/lib
+# EXTINC += -I$(FFTWDIR)/include
+# # Add lapack
+# EXTLIBS += -L$(LAPACKDIR)/lib
+# EXTINC += -I$(LAPACKDIR)/mod
+# Standard version with optimized compilation
+all: dirs src/srcinfo.h
+# all: F90FLAGS = -O3 -xHOST
+all: F90FLAGS = -std=legacy -ffree-line-length-0 -O3
+all: $(EXEC)
+# Fast compilation
+fast: dirs src/srcinfo.h
+fast: F90FLAGS = -fast
+fast: $(EFST)
+# Debug version with all flags
+debug: dirs src/srcinfo.h
+debug: F90FLAGS = -g -std=legacy -ffree-line-length-0 #-traceback -CB -ftrapuv -warn all -debug all
+debug: $(EDBG)
+# Alpha version, optimized as all but creates another binary
+alpha: dirs src/srcinfo.h
+alpha: F90FLAGS = -O3 -xHOST
+alpha: $(EALP)
+
+install: dirs src/srcinfo.h $(EXEC) mvmod
+
+run: all
+	(cd wk; $(EXEC);)
+
+dirs:
+	mkdir -p $(BINDIR)
+	mkdir -p $(OBJDIR)
+	mkdir -p $(MODDIR)
+
+src/srcinfo.h:
+	( cd src/srcinfo; $(MAKE);)
+
+clean: cleanobj cleanmod
+	@rm -f src/srcinfo.h
+	@rm -f src/srcinfo/srcinfo.h
+
+cleanobj:
+	@rm -f $(OBJDIR)/*o
+
+cleanmod:
+	@rm -f $(MODDIR)/*mod
+	@rm -f *.mod
+
+cleanbin:
+	@rm -f $(EXEC)
+
+mvmod:
+	mv *.mod mod/.
+
+$(OBJDIR)/diagnose.o : src/srcinfo.h
+
+FOBJ=$(OBJDIR)/advance_field_mod.o $(OBJDIR)/array_mod.o $(OBJDIR)/auxval.o \
+$(OBJDIR)/basic_mod.o $(OBJDIR)/coeff_mod.o $(OBJDIR)/closure_mod.o \
+$(OBJDIR)/collision_mod.o $(OBJDIR)/nonlinear_mod.o $(OBJDIR)/control.o \
+$(OBJDIR)/diagnose.o $(OBJDIR)/diagnostics_par_mod.o $(OBJDIR)/endrun.o \
+$(OBJDIR)/fields_mod.o $(OBJDIR)/fourier_mod.o $(OBJDIR)/geometry_mod.o \
+$(OBJDIR)/ghosts_mod.o $(OBJDIR)/grid_mod.o $(OBJDIR)/inital.o \
+$(OBJDIR)/initial_par_mod.o $(OBJDIR)/lag_interp_mod.o $(OBJDIR)/main.o \
+$(OBJDIR)/memory.o $(OBJDIR)/miller_mod.o $(OBJDIR)/model_mod.o \
+$(OBJDIR)/moments_eq_rhs_mod.o $(OBJDIR)/numerics_mod.o $(OBJDIR)/parallel_mod.o \
+$(OBJDIR)/ppexit.o $(OBJDIR)/ppinit.o $(OBJDIR)/prec_const_mod.o \
+$(OBJDIR)/processing_mod.o $(OBJDIR)/readinputs.o $(OBJDIR)/restarts_mod.o \
+$(OBJDIR)/solve_EM_fields.o $(OBJDIR)/stepon.o $(OBJDIR)/tesend.o \
+$(OBJDIR)/time_integration_mod.o $(OBJDIR)/utility_mod.o
+
+ $(EXEC): $(FOBJ)
+	$(F90) $(LDFLAGS) $(OBJDIR)/*.o $(EXTMOD) $(EXTINC) $(EXTLIBS) -o $@
+
+ $(EFST): $(FOBJ)
+	$(F90) $(LDFLAGS) $(OBJDIR)/*.o $(EXTMOD) $(EXTINC) $(EXTLIBS) -o $@
+
+ $(EDBG): $(FOBJ)
+	$(F90) $(LDFLAGS) $(OBJDIR)/*.o $(EXTMOD) $(EXTINC) $(EXTLIBS) -o $@
+
+ $(EALP): $(FOBJ)
+	$(F90) $(LDFLAGS) $(OBJDIR)/*.o $(EXTMOD) $(EXTINC) $(EXTLIBS) -o $@
+
+
+ $(OBJDIR)/advance_field_mod.o : src/advance_field_mod.F90 \
+   $(OBJDIR)/grid_mod.o $(OBJDIR)/array_mod.o $(OBJDIR)/initial_par_mod.o \
+	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/time_integration_mod.o $(OBJDIR)/basic_mod.o \
+	 $(OBJDIR)/fields_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/advance_field_mod.F90 -o $@
+
+ $(OBJDIR)/array_mod.o : src/array_mod.F90 \
+	 $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/array_mod.F90 -o $@
+
+ $(OBJDIR)/auxval.o : src/auxval.F90 \
+ 	 $(OBJDIR)/fourier_mod.o $(OBJDIR)/memory.o $(OBJDIR)/model_mod.o \
+	 $(OBJDIR)/geometry_mod.o  $(OBJDIR)/grid_mod.o $(OBJDIR)/numerics_mod.o \
+	 $(OBJDIR)/parallel_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/auxval.F90 -o $@
+
+ $(OBJDIR)/basic_mod.o : src/basic_mod.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/basic_mod.F90 -o $@
+
+ $(OBJDIR)/calculus_mod.o : src/calculus_mod.F90  \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o \
+	 $(OBJDIR)/parallel_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/calculus_mod.F90 -o $@
+
+ $(OBJDIR)/coeff_mod.o : src/coeff_mod.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/model_mod.o \
+	 $(OBJDIR)/basic_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/coeff_mod.F90 -o $@
+
+ $(OBJDIR)/closure_mod.o : src/closure_mod.F90 \
+ 	 $(OBJDIR)/model_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/grid_mod.o \
+	 $(OBJDIR)/array_mod.o $(OBJDIR)/fields_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/closure_mod.F90 -o $@
+
+ $(OBJDIR)/collision_mod.o : src/collision_mod.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/fields_mod.o \
+   $(OBJDIR)/grid_mod.o $(OBJDIR)/model_mod.o $(OBJDIR)/prec_const_mod.o \
+	 $(OBJDIR)/time_integration_mod.o $(OBJDIR)/utility_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/collision_mod.F90 -o $@
+
+ $(OBJDIR)/control.o : src/control.F90 \
+	 $(OBJDIR)/auxval.o $(OBJDIR)/geometry_mod.o $(OBJDIR)/prec_const_mod.o \
+   $(OBJDIR)/basic_mod.o $(OBJDIR)/ppexit.o $(OBJDIR)/ppinit.o \
+	 $(OBJDIR)/readinputs.o $(OBJDIR)/tesend.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/control.F90 -o $@
+
+ $(OBJDIR)/diagnose.o : src/diagnose.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/processing_mod.o $(OBJDIR)/array_mod.o \
+   $(OBJDIR)/basic_mod.o $(OBJDIR)/diagnostics_par_mod.o $(OBJDIR)/fields_mod.o \
+	 $(OBJDIR)/grid_mod.o $(OBJDIR)/initial_par_mod.o $(OBJDIR)/model_mod.o \
+	 $(OBJDIR)/time_integration_mod.o\
+	 $(OBJDIR)/parallel_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/diagnose.F90 -o $@
+
+ $(OBJDIR)/diagnostics_par_mod.o : src/diagnostics_par_mod.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/basic_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/diagnostics_par_mod.F90 -o $@
+
+ $(OBJDIR)/endrun.o : src/endrun.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/basic_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/endrun.F90 -o $@
+
+ $(OBJDIR)/fields_mod.o : src/fields_mod.F90 \
+   $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/fields_mod.F90 -o $@
+
+ $(OBJDIR)/fourier_mod.o : src/fourier_mod.F90 \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/fourier_mod.F90 -o $@
+
+ $(OBJDIR)/geometry_mod.o : src/geometry_mod.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/calculus_mod.o $(OBJDIR)/miller_mod.o \
+   $(OBJDIR)/grid_mod.o $(OBJDIR)/model_mod.o $(OBJDIR)/prec_const_mod.o \
+	 $(OBJDIR)/utility_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/geometry_mod.F90 -o $@
+
+ $(OBJDIR)/ghosts_mod.o : src/ghosts_mod.F90 \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/fields_mod.o $(OBJDIR)/grid_mod.o\
+   $(OBJDIR)/geometry_mod.o $(OBJDIR)/ppinit.o  $(OBJDIR)/time_integration_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/ghosts_mod.F90 -o $@
+
+ $(OBJDIR)/grid_mod.o : src/grid_mod.F90 \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/model_mod.o $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/grid_mod.F90 -o $@
+
+ $(OBJDIR)/inital.o : src/inital.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/fields_mod.o \
+	 $(OBJDIR)/initial_par_mod.o $(OBJDIR)/model_mod.o $(OBJDIR)/numerics_mod.o \
+	 $(OBJDIR)/solve_EM_fields.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/ghosts_mod.o \
+	 $(OBJDIR)/grid_mod.o \
+	 $(OBJDIR)/restarts_mod.o $(OBJDIR)/time_integration_mod.o $(OBJDIR)/utility_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/inital.F90 -o $@
+
+ $(OBJDIR)/initial_par_mod.o : src/initial_par_mod.F90 \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/initial_par_mod.F90 -o $@
+
+ $(OBJDIR)/lag_interp_mod.o : src/lag_interp_mod.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/lag_interp_mod.F90 -o $@
+
+ $(OBJDIR)/main.o : src/main.F90 $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/main.F90 -o $@
+
+ $(OBJDIR)/memory.o : src/memory.F90 $ \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/collision_mod.o\
+   $(OBJDIR)/fields_mod.o $(OBJDIR)/model_mod.o $(OBJDIR)/time_integration_mod.o \
+	 $(OBJDIR)/grid_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/memory.F90 -o $@
+
+ $(OBJDIR)/miller_mod.o : src/miller_mod.F90 \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/model_mod.o $(OBJDIR)/lag_interp_mod.o \
+   $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/miller_mod.F90 -o $@
+
+ $(OBJDIR)/model_mod.o : src/model_mod.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/model_mod.F90 -o $@
+
+ $(OBJDIR)/moments_eq_rhs_mod.o : src/moments_eq_rhs_mod.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/calculus_mod.o $(OBJDIR)/fields_mod.o \
+	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o $(OBJDIR)/model_mod.o \
+	 $(OBJDIR)/time_integration_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/moments_eq_rhs_mod.F90 -o $@
+
+ $(OBJDIR)/nonlinear_mod.o : src/nonlinear_mod.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/fourier_mod.o \
+	 $(OBJDIR)/fields_mod.o $(OBJDIR)/grid_mod.o $(OBJDIR)/model_mod.o\
+	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/time_integration_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/nonlinear_mod.F90 -o $@
+
+ $(OBJDIR)/numerics_mod.o : src/numerics_mod.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o  $(OBJDIR)/basic_mod.o $(OBJDIR)/coeff_mod.o \
+	 $(OBJDIR)/utility_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/numerics_mod.F90 -o $@
+
+ $(OBJDIR)/parallel_mod.o : src/parallel_mod.F90 \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/parallel_mod.F90 -o $@
+
+ $(OBJDIR)/ppexit.o : src/ppexit.F90 \
+ 	 $(OBJDIR)/prec_const_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/coeff_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/ppexit.F90 -o $@
+
+ $(OBJDIR)/ppinit.o : src/ppinit.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o\
+   $(OBJDIR)/fields_mod.o $(OBJDIR)/array_mod.o $(OBJDIR)/time_integration_mod.o \
+	 $(OBJDIR)/basic_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/ppinit.F90 -o $@
+
+ $(OBJDIR)/prec_const_mod.o : src/prec_const_mod.F90
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/prec_const_mod.F90 -o $@
+
+ $(OBJDIR)/processing_mod.o : src/processing_mod.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o \
+   $(OBJDIR)/fields_mod.o $(OBJDIR)/basic_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/processing_mod.F90 -o $@
+
+ $(OBJDIR)/readinputs.o : src/readinputs.F90  \
+   $(OBJDIR)/diagnostics_par_mod.o $(OBJDIR)/initial_par_mod.o $(OBJDIR)/model_mod.o  \
+  $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o $(OBJDIR)/time_integration_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/readinputs.F90 -o $@
+
+ $(OBJDIR)/restarts_mod.o : src/restarts_mod.F90   \
+	 $(OBJDIR)/diagnostics_par_mod.o $(OBJDIR)/grid_mod.o $(OBJDIR)/time_integration_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/restarts_mod.F90 -o $@
+
+ $(OBJDIR)/solve_EM_fields.o : src/solve_EM_fields.F90 \
+ 	 $(OBJDIR)/array_mod.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/grid_mod.o \
+	 $(OBJDIR)/ghosts_mod.o $(OBJDIR)/fields_mod.o $(OBJDIR)/array_mod.o \
+	 $(OBJDIR)/time_integration_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/parallel_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/solve_EM_fields.F90 -o $@
+
+ $(OBJDIR)/stepon.o : src/stepon.F90 \
+ 	 $(OBJDIR)/initial_par_mod.o $(OBJDIR)/prec_const_mod.o $(OBJDIR)/advance_field_mod.o \
+   $(OBJDIR)/basic_mod.o $(OBJDIR)/nonlinear_mod.o $(OBJDIR)/grid_mod.o \
+	 $(OBJDIR)/array_mod.o $(OBJDIR)/numerics_mod.o $(OBJDIR)/fields_mod.o \
+	 $(OBJDIR)/ghosts_mod.o $(OBJDIR)/moments_eq_rhs_mod.o $(OBJDIR)/solve_EM_fields.o\
+	 $(OBJDIR)/utility_mod.o $(OBJDIR)/model_mod.o $(OBJDIR)/time_integration_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/stepon.F90 -o $@
+
+ $(OBJDIR)/tesend.o : src/tesend.F90 \
+ 	 $(OBJDIR)/basic_mod.o $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/tesend.F90 -o $@
+
+ $(OBJDIR)/time_integration_mod.o : src/time_integration_mod.F90 \
+   $(OBJDIR)/basic_mod.o  $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/time_integration_mod.F90 -o $@
+
+ $(OBJDIR)/utility_mod.o : src/utility_mod.F90  \
+   $(OBJDIR)/grid_mod.o $(OBJDIR)/basic_mod.o $(OBJDIR)/prec_const_mod.o
+	$(F90) -c $(F90FLAGS) $(FPPFLAGS) $(EXTMOD) $(EXTINC) src/utility_mod.F90 -o $@
diff --git a/local/dirs_wsl_debian.inc b/local/dirs_wsl_debian.inc
new file mode 100644
index 0000000000000000000000000000000000000000..1391f01cfad5de94f76f47f1158d7845f05a65bd
--- /dev/null
+++ b/local/dirs_wsl_debian.inc
@@ -0,0 +1,16 @@
+#Local code, binaries, pputils library
+PREFIX   = $(HOME)/gyacomo
+SRCDIR   = $(PREFIX)/src
+BINDIR   = $(PREFIX)/bin
+OBJDIR   = $(PREFIX)/obj
+LIBDIR   = $(PREFIX)/lib
+MODDIR   = $(PREFIX)/mod
+
+HDF5_LIB   = /usr/local/hdf5-1.8.22/lib
+FMDIR      = $(HOME)/lib/FM
+FFTW3DIR   = $(HOME)/lib/fftw-3.3.10
+#FFTW3DIR  = $(FFTW_HOME)#for Marconi
+FUTILS_DIR = $(HOME)/lib/futils/src
+#LAPACKDIR = $(HOME)/lib/lapack
+
+# Naming ideas : HeLaZ, MoNoLiT (Moment Non Linear Torroidal)
diff --git a/testcases/cyclone_example/fort.90 b/testcases/cyclone_example/fort.90
index e595b41a0a5c39dee1064a074fc77283af24b09c..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391 100644
Binary files a/testcases/cyclone_example/fort.90 and b/testcases/cyclone_example/fort.90 differ
diff --git a/testcases/cyclone_example/fort_00.90 b/testcases/cyclone_example/fort_00.90
new file mode 100644
index 0000000000000000000000000000000000000000..fd926cba93a53a4bd2ee432dc0fa030c9dfc9c6c
--- /dev/null
+++ b/testcases/cyclone_example/fort_00.90
@@ -0,0 +1,81 @@
+&BASIC
+  nrun   = 100000000
+  dt     = 0.01
+  tmax   = 50
+  maxruntime = 356400
+/
+&GRID
+  pmaxe  = 4
+  jmaxe  = 2
+  pmaxi  = 4
+  jmaxi  = 2
+  Nx     = 128
+  Lx     = 120
+  Ny     = 64
+  Ly     = 160
+  Nz     = 16
+  Nexc   = 0
+  SG     = .true.
+/
+&GEOMETRY
+  geom   = 's-alpha'
+  q0     = 1.4
+  shear  = 0.8
+  eps    = 0.18
+  parallel_BC = 'dirichlet'
+/
+&OUTPUT_PAR
+  nsave_0d = 50
+  nsave_1d = -1
+  nsave_2d = -1
+  nsave_3d = 100
+  nsave_5d = 500
+  write_doubleprecision = .false.
+  write_gamma = .true.
+  write_hf    = .true.
+  write_phi   = .true.
+  write_Na00  = .false.
+  write_Napj  = .true.
+  write_Sapj  = .false.
+  write_dens  = .true.
+  write_temp  = .true.
+  job2load    = -1
+/
+&MODEL_PAR
+  ! Collisionality
+  CLOS    = 0
+  NL_CLOS = 0
+  LINEARITY = 'nonlinear'
+  KIN_E   = .false.
+  mu_x    = 0.1
+  mu_y    = 0.1
+  N_HD    = 2
+  mu_z    = 2.0
+  mu_p    = 0
+  mu_j    = 0
+  nu      = 0.05
+  tau_e   = 1
+  tau_i   = 1
+  sigma_e = 0.023338
+  sigma_i = 1
+  q_e     = -1
+  q_i     = 1
+  K_Ni    = 2.22
+  K_Ti    = 6.92
+/
+&COLLISION_PAR
+  collision_model = 'DG'
+  gyrokin_CO      = .false.
+  interspecies    = .true.
+  mat_file        = 'null'
+/
+&INITIAL_CON
+  INIT_OPT      = 'blob'
+  ACT_ON_MODES  = 'donothing'
+  init_background  = 0
+  init_noiselvl = 1e-3
+  iseed         = 42
+/
+&TIME_INTEGRATION_PAR
+  numerical_scheme = 'RK4'
+/
diff --git a/testcases/zpinch_example/fort_00.90 b/testcases/zpinch_example/fort_00.90
index cd43ea92c69b2958c5c2163cfc85f351b68edffd..67312e5472dc8d8c4ae88ed959ac2e9744dbd936 100644
--- a/testcases/zpinch_example/fort_00.90
+++ b/testcases/zpinch_example/fort_00.90
@@ -1,7 +1,7 @@
 &BASIC
   nrun   = 100000000
   dt     = 0.01
-  tmax   = 200
+  tmax   = 50
   maxruntime = 356400
 /
 &GRID
@@ -33,7 +33,7 @@
   write_gamma = .t.
   write_hf    = .t.
   write_phi   = .t.
-  write_Na00  = .t.
+  write_Na00  = .f.
   write_Napj  = .t.
   write_Sapj  = .f.
   write_dens  = .t.