update scripts

matlab/compile_results.m

@@ -298,7 +298,7 @@ function [DATA] = compile_results(DIRECTORY,JOBNUMMIN,JOBNUMMAX)
         end
         if W_NA00
         DATA.Ni00   = zeros(Nky,Nkx,Nz,numel(Ts3D_));   DATA.Ni00(:,:,1:Nz,:)   = Ni00_;
-        DATA.Nipjz  = zeros(Nky,Nkx,Nz,numel(Ts3D_));  DATA.Nipjz(:,:,1:Nz,:)   = Nipjz_;
+        DATA.Nipjz  = zeros(Pi_new,Ji_new,Nz,numel(Ts3D_));  DATA.Nipjz(:,:,1:Nz,:)   = Nipjz_;
         end
         if W_DENS
         DATA.DENS_I = zeros(Nky,Nkx,Nz,numel(Ts3D_)); DATA.DENS_I(:,:,1:Nz,:)   = DENS_I_;

wk/analysis_gyacomo.m

 %% UNCOMMENT FOR TUTORIAL
 gyacomodir = pwd; gyacomodir = gyacomodir(1:end-2); % get code directory
-resdir = '../testcases/zpinch_example'; %Name of the directory where the results are located
-PARTITION = '';
-JOBNUMMIN = 00; JOBNUMMAX = 10;
+% resdir = '../testcases/zpinch_example'; %Name of the directory where the results are located
+% PARTITION = '';
+% JOBNUMMIN = 00; JOBNUMMAX = 10;
 %%
 addpath(genpath([gyacomodir,'matlab'])) % ... add
 addpath(genpath([gyacomodir,'matlab/plot'])) % ... add

wk/header_3D_results.m

@@ -59,9 +59,8 @@ resdir = 'paper_2_nonlinear/kT_5.3/5x3x128x64x64';
 % resdir = 'paper_2_nonlinear/kT_5.3/11x6x128x64x24';
 % resdir = 'paper_2_nonlinear/kT_5.3/11x6x128x64x64';
 
-%% Old stuff
-% resdir = 'CBC/kT_4.5_128x64x16x13x7/';
-% resdir = 'CBC/kT_5.3_192x96x24x13x7/';
+%% Scan nu
+resdir = 'paper_2_nonlinear/nu_scan_kT_5.3/FCGK_5x3x128x64x24_dp';
 
 %% Miller
 % resdir = 'paper_2_nonlinear/kT_4.15_miller/7x4x128x64x24';