diff --git a/local/dirs.inc b/local/dirs.inc
index 5d0bbe8dcb029e29f502ec997eaa3cedb75e42b9..c7b2c0ea07373e4f2c0e3d6be20165cf8db8e11d 100644
--- a/local/dirs.inc
+++ b/local/dirs.inc
@@ -1,17 +1,13 @@
 #Local code, binaries, pputils library
-PREFIX   = $(HOME)/HeLaZ#write the path to the root of your HeLaZ distrib here
+PREFIX   = $(HOME)/HeLaZ
 SRCDIR   = $(PREFIX)/src
 BINDIR   = $(PREFIX)/bin
 OBJDIR   = $(PREFIX)/obj
 LIBDIR   = $(PREFIX)/lib
 MODDIR   = $(PREFIX)/mod
 
-#FMDIR      = $(PREFIX)/FM
-#FFTW3DIR   = /home/ahoffman/lib/fftw-3.3.8
-#FUTILS_DIR = /home/ahoffman/local/spcpc/intel
-
 FMDIR      = $(HOME)/lib/FM
 FFTW3DIR   = $(HOME)/lib/fftw-3.3.8
-FUTILS_DIR = $(HOME)/lib/futils
+FUTILS_DIR = $(HOME)/lib/futils/src
 
 # Naming ideas : HeLaZ, MoNoLiT (Moment Non Linear Torroidal)
diff --git a/matlab/write_sbash.m b/matlab/write_sbash.m
new file mode 100644
index 0000000000000000000000000000000000000000..573e6302d0a8ecf4713077742277378a96ce8289
--- /dev/null
+++ b/matlab/write_sbash.m
@@ -0,0 +1,49 @@
+% Write the input script "fort.90" with desired parameters
+INPUT = 'setup_and_run.sh';
+fid = fopen(INPUT,'wt');
+
+fprintf(fid,[...
+'#!/bin/bash\n',...
+'mkdir -p $CINECA_SCRATCH/HeLaZ/wk\n',...
+'mkdir -p $CINECA_SCRATCH/HeLaZ/bin\n\n',...
+...
+'cd $CINECA_SCRATCH/HeLaZ/wk/\n',...
+'cp $HOME/HeLaZ/wk/fort.90 .\n',...
+'cp $HOME/HeLaZ/wk/batch_script.sh .\n',...
+...%'cp $HOME/HeLaZ/bin/helaz ../bin\n\n',...
+'cp -r $HOME/HeLaZ/iCa ..\n\n',...
+...
+'mkdir -p ', BASIC.RESDIR,'\n',...
+...
+'sbatch batch_script.sh\n',...
+'echo $CINECA_SCRATCH/HeLaZ/results/',BASIC.SIMID,'/',BASIC.PARAMS,'/out.txt']);
+
+fclose(fid);
+system(['cp setup_and_run.sh ',BASIC.RESDIR,'/.']);
+
+% Write the sbatch script
+INPUT = 'batch_script.sh';
+fid = fopen(INPUT,'wt');
+
+fprintf(fid,[...
+'#!/bin/bash\n',...
+'#SBATCH --time=', CLUSTER.TIME,'\n',...
+'#SBATCH --nodes=', CLUSTER.NODES,'\n',...
+'#SBATCH --cpus-per-task=', CLUSTER.CPUPT,'\n',...
+'#SBATCH --ntasks-per-node=', CLUSTER.NTPN,'\n',...
+'#SBATCH --mem=', CLUSTER.MEM,'\n',...
+'#SBATCH --error=', BASIC.RESDIR,'err.txt\n',...
+'#SBATCH --output=', BASIC.RESDIR,'out.txt\n',...
+'#SBATCH --account=FUA34_GBSedge\n',...
+'#SBATCH --partition=skl_fua_',CLUSTER.PART,'\n\n',...
+...% '#SBATCH --job-name=',PARAMS,'\n\n',...
+'module load intel\n',...
+'module load intelmpi\n',...
+'module load autoload hdf5/1.10.4--intelmpi--2018--binary\n',...
+'module load fftw\n',...
+'srun --cpu-bind=cores ./../bin/helaz']);
+
+fclose(fid);
+system(['cp batch_script.sh ',BASIC.RESDIR,'/.']);
+
+system('scp {fort.90,setup_and_run.sh,batch_script.sh} ahoffman@login.marconi.cineca.it:/marconi/home/userexternal/ahoffman/HeLaZ/wk');