diff --git a/.gitignore b/.gitignore
index da9eb16af31b6cb6fd02ce2ea438592d60848e5e..cf060981ea967104df70c1dfd585510bad86c430 100644
--- a/.gitignore
+++ b/.gitignore
@@ -12,6 +12,7 @@
 *.emf
 *.pptx
 *.jpg
+*.gif
 *.jpeg
 *.mat
 *.dat
@@ -34,10 +35,10 @@ iCa/
 *.out
 src/srcinfo.h
 src/srcinfo/srcinfo.h
-*/fort*.90
 local/
 */*.sh
 Gallery/
 .vscode/settings.json
 *figure*
 *results*
+out*
diff --git a/Makefile b/Makefile
index 9112278d95986de1ce17654219df6c1fc245ed03..fd9aa229c40184ff3b9b294aee34764211229b87 100644
--- a/Makefile
+++ b/Makefile
@@ -1,9 +1,9 @@
 include local/dirs.inc
 include local/make.inc
 
-EXEC = $(BINDIR)/helaz3
-EFST = $(BINDIR)/helaz3_fst
-EDBG = $(BINDIR)/helaz3_dbg
+EXEC = $(BINDIR)/gyacomo
+EFST = $(BINDIR)/gyacomo_fst
+EDBG = $(BINDIR)/gyacomo_dbg
 
 F90 = mpiifort
 # #F90 = ftn #for piz-daint cluster
diff --git a/local/dirs.inc b/local/dirs.inc
index 3c5dca40e6a23ce33f5bf82929e62e097501c62a..2695b83f22c3574509d9d1ef3c4a4b800f7f6616 100644
--- a/local/dirs.inc
+++ b/local/dirs.inc
@@ -1,5 +1,5 @@
 #Local code, binaries, pputils library
-PREFIX   = $(HOME)/HeLaZ
+PREFIX   = $(HOME)/gyacomo
 SRCDIR   = $(PREFIX)/src
 BINDIR   = $(PREFIX)/bin
 OBJDIR   = $(PREFIX)/obj
diff --git a/wk/analysis_HeLaZ.m b/wk/analysis_gyacomo.m
similarity index 84%
rename from wk/analysis_HeLaZ.m
rename to wk/analysis_gyacomo.m
index 370c61dd64bb1df607e3ec3109728a6ced6a39d7..70c4f92d238e890828fbdbcec728897f51597abd 100644
--- a/wk/analysis_HeLaZ.m
+++ b/wk/analysis_gyacomo.m
@@ -1,11 +1,16 @@
-addpath(genpath([helazdir,'matlab'])) % ... add
-addpath(genpath([helazdir,'matlab/plot'])) % ... add
-addpath(genpath([helazdir,'matlab/compute'])) % ... add
-addpath(genpath([helazdir,'matlab/load'])) % ... add
+%% UNCOMMENT FOR TUTORIAL
+% gyacomodir = pwd; gyacomodir = gyacomodir(1:end-2); % get code directory
+% resdir = '.'; %Name of the directory where the results are located
+% JOBNUMMIN = 00; JOBNUMMAX = 10;
+%%
+addpath(genpath([gyacomodir,'matlab'])) % ... add
+addpath(genpath([gyacomodir,'matlab/plot'])) % ... add
+addpath(genpath([gyacomodir,'matlab/compute'])) % ... add
+addpath(genpath([gyacomodir,'matlab/load'])) % ... add
 
 %% Load the results
-LOCALDIR  = [helazdir,'results/',outfile,'/'];
-MISCDIR   = ['/misc/HeLaZ_outputs/results/',outfile,'/'];
+LOCALDIR  = [gyacomodir,resdir,'/'];
+MISCDIR   = ['/misc/gyacomo_outputs/',resdir,'/']; %For long term storage
 system(['mkdir -p ',MISCDIR]);
 system(['mkdir -p ',LOCALDIR]);
 CMD = ['rsync ', LOCALDIR,'outputs* ',MISCDIR]; disp(CMD);
@@ -28,8 +33,8 @@ options.TAVG_1   = 1000;%0.9*data.Ts3D(end); % Averaging times duration
 options.NCUT     = 4;              % Number of cuts for averaging and error estimation
 options.NMVA     = 1;              % Moving average for time traces
 % options.ST_FIELD = '\Gamma_x';   % chose your field to plot in spacetime diag (e.g \phi,v_x,G_x)
-options.ST_FIELD = '\phi';          % chose your field to plot in spacetime diag (e.g \phi,v_x,G_x)
-options.INTERP   = 1;
+% options.ST_FIELD = '\phi';          % chose your field to plot in spacetime diag (e.g \phi,v_x,G_x)
+% options.INTERP   = 1;
 fig = plot_radial_transport_and_spacetime(data,options);
 save_figure(data,fig,'.png')
 end
@@ -45,8 +50,8 @@ if 0
 % Options
 options.INTERP    = 1;
 options.POLARPLOT = 0;
-options.NAME      = '\phi';
-% options.NAME      = 'N_i^{00}';
+% options.NAME      = '\phi';
+options.NAME      = 'N_i^{00}';
 % options.NAME      = 'v_y';
 % options.NAME      = 'n_i^{NZ}';
 % options.NAME      = '\Gamma_x';
@@ -57,7 +62,7 @@ options.PLAN      = 'xy';
 options.COMP      = 'avg';
 % options.TIME      = data.Ts5D(end-30:end);
 % options.TIME      =  data.Ts3D;
-options.TIME      = [1:0.2:500];
+options.TIME      = [0:500];
 data.EPS          = 0.1;
 data.a = data.EPS * 2000;
 create_film(data,options,'.gif')