save scripts

matlab/compile_results.m

@@ -9,6 +9,7 @@ function [DATA] = compile_results(DIRECTORY,JOBNUMMIN,JOBNUMMAX)
     DATA.MUy_EVOL  = []; % evolution of parameter mu between jobs
     DATA.MUz_EVOL  = []; % evolution of parameter mu between jobs
     DATA.K_N_EVOL = []; %
+    DATA.K_T_EVOL = []; %
     DATA.L_EVOL   = []; % 
     DATA.DT_EVOL  = []; %
     % FIELDS
@@ -238,6 +239,7 @@ function [DATA] = compile_results(DIRECTORY,JOBNUMMIN,JOBNUMMAX)
             DATA.MUy_EVOL  = [DATA.MUy_EVOL  DATA.MUy    DATA.MUy];
             DATA.MUz_EVOL  = [DATA.MUz_EVOL  DATA.MUz    DATA.MUz];
             DATA.K_N_EVOL  = [DATA.K_N_EVOL DATA.K_N   DATA.K_N];
+            DATA.K_T_EVOL  = [DATA.K_T_EVOL DATA.K_T   DATA.K_T];
             DATA.L_EVOL    = [DATA.L_EVOL    DATA.L      DATA.L];
             DATA.DT_EVOL   = [DATA.DT_EVOL   DATA.DT_SIM DATA.DT_SIM];
             

matlab/plot/statistical_transport_averaging.m

@@ -44,7 +44,7 @@ if options.NPLOTS > 0
     legend(['$Q_x^\infty=$',sprintf('%2.2e',transp_seg_avg(end))])
 end   
 res.time_seg = time_seg;
-res.Qx_t     = 
+res.Qx_t     = transp_seg_avg;
 res.Gx_avg = mean(gamma);
 res.Gx_std = std (gamma);
 disp(['G_x=',sprintf('%2.2e',res.Gx_avg),'+-',sprintf('%2.2e',res.Gx_std)]);

wk/analysis_gyacomo.m

@@ -37,7 +37,7 @@ if 1
 % options.TAVG_0   = data.TJOB_SE(2*jid_+1);%0.4*data.Ts3D(end);
 % options.TAVG_1   = data.TJOB_SE(2*(jid_+1));%0.9*data.Ts3D(end); % Averaging times duration
 options.TAVG_0   = 100;
-options.TAVG_1   = 1000;
+options.TAVG_1   = 500;
 options.NCUT     = 5;              % 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)
@@ -52,8 +52,8 @@ if 0
 %% statistical transport averaging
 Gavg =[]; Gstd = [];
 Qavg =[]; Qstd = [];
-figure; hold on;
-plot(data.Ts0D,data.Qx);
+% figure; hold on;
+% plot(data.Ts0D,data.HFLUX_X);
 for i_ = 1:2:numel(data.TJOB_SE) 
 % i_ = 3; 
 disp([num2str(data.TJOB_SE(i_)),' ',num2str(data.TJOB_SE(i_+1))])
@@ -66,6 +66,11 @@ Qavg = [Qavg res.Qx_avg]; Qstd = [Qstd res.Qx_std];
 end
 % disp(Gavg); disp(Gstd);
 disp(Qavg); disp(Qstd);
+figure
+% errorbar(data.NU_EVOL(2:2:end),Qavg,Qstd,'--s');
+% xlabel('$\nu$'); ylabel('$Q_x^\infty$');
+errorbar(data.K_T_EVOL(2:2:end),Qavg,Qstd,'--s','DisplayName',data.paramshort);
+xlabel('$\kappa_T$'); ylabel('$Q_x^\infty$');
 end
 if 0
 %% MOVIES %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -93,7 +98,7 @@ options.RESOLUTION = 256;
 create_film(data,options,'.gif')
 end
 
-if 1
+if 0
 %% fields snapshots
 % Options
 options.INTERP    = 0;

wk/header_3D_results.m

@@ -46,7 +46,7 @@ PARTITION  = '/misc/gyacomo_outputs/';
 % resdir = 'paper_2_nonlinear/kT_6.96_nu_0.01/5x3x192x96x24';
 
 %% kT=5.3 results
-resdir = 'paper_2_nonlinear/kT_5.3/5x3x128x64x64';
+% resdir = 'paper_2_nonlinear/kT_5.3/5x3x128x64x64';
 % resdir = 'paper_2_nonlinear/kT_5.3/5x3x128x64x24';
 % resdir = 'paper_2_nonlinear/kT_5.3/7x4x128x64x24';
 % resdir = 'paper_2_nonlinear/kT_5.3/7x4x128x64x24_MUxy_0';
@@ -59,8 +59,22 @@ resdir = 'paper_2_nonlinear/kT_5.3/5x3x128x64x64';
 % resdir = 'paper_2_nonlinear/kT_5.3/11x6x128x64x24';
 % resdir = 'paper_2_nonlinear/kT_5.3/11x6x128x64x64';
 
-%% Scan nu
-resdir = 'paper_2_nonlinear/nu_scan_kT_5.3/FCGK_5x3x128x64x24_dp';
+%% kT=6.96 results
+% resdir = 'paper_2_nonlinear/kT_6.96_dp/3x2x128x64x24';
+% resdir = 'paper_2_nonlinear/kT_6.96_dp/3x2x192x96x32';
+% resdir = 'paper_2_nonlinear/kT_6.96_dp/5x3x192x96x32';
+resdir = 'paper_2_nonlinear/kT_6.96_dp/7x4x192x96x32';
+
+%% Scan kT
+% resdir = 'paper_2_nonlinear/kT_scan_nu_1e-3/5x3x128x64x24_dp';
+% resdir = 'paper_2_nonlinear/kT_scan_nu_1e-3/5x3x192x96x32_dp';
+% resdir = 'paper_2_nonlinear/kT_scan_nu_1e-3/7x4x128x64x24_dp';
+% resdir = 'paper_2_nonlinear/kT_scan_nu_1e-3/7x4x192x96x32_dp';
+
+%% Scan nu, kT = 5.3
+% resdir = 'paper_2_nonlinear/nu_scan_kT_5.3/FCGK_5x3x128x64x24_dp';
+% resdir = 'paper_2_nonlinear/nu_scan_kT_5.3/DGGK_7x4x128x64x24_dp';
+% resdir = 'paper_2_nonlinear/nu_scan_kT_5.3/SGGK_7x4x128x64x24_dp';
 
 %% Miller
 % resdir = 'paper_2_nonlinear/kT_4.15_miller/7x4x128x64x24';