save scripts

wk/analysis_gene.m

@@ -22,12 +22,13 @@ addpath(genpath([helazdir,'matlab/load'])) % ... add
 % folder = '/misc/gene_results/ZP_kn_2.5_large_box/';
 % folder = '/misc/gene_results/CBC/128x64x16x24x12/';
 % folder = '/misc/gene_results/CBC/196x96x20x32x16_02/';
-folder = '/misc/gene_results/CBC/128x64x16x6x4/';
+% folder = '/misc/gene_results/CBC/128x64x16x6x4/';
 % folder = '/misc/gene_results/CBC/KT_5.3_128x64x16x24x12_01/';
 % folder = '/misc/gene_results/CBC/KT_4.5_128x64x16x24x12_01/';
 % folder = '/misc/gene_results/CBC/KT_9_128x64x16x24x12/';
 % folder = '/misc/gene_results/CBC/KT_13_large_box_128x64x16x24x12/';
 % folder = '/misc/gene_results/CBC/Lapillone_Fig6/';
+folder = '/misc/gene_results/Z-pinch/HP_kN_1.6_adapt_mu/';
 gene_data = load_gene_data(folder);
 gene_data = invert_kxky_to_kykx_gene_results(gene_data);
 if 1

wk/analysis_gyacomo.m

@@ -28,10 +28,11 @@ FMT = '.fig';
 if 1
 %% Space time diagramm (fig 11 Ivanov 2020)
 % data.scale = 1;%/(data.Nx*data.Ny)^2;
-options.TAVG_0   = 350;%0.4*data.Ts3D(end);
-options.TAVG_1   = 1000;%0.9*data.Ts3D(end); % Averaging times duration
+i_ = 19; disp([num2str(data.NU_EVOL(i_)),' ',num2str(data.NU_EVOL(i_+1))])
+options.TAVG_0   = data.TJOB_SE(i_);%0.4*data.Ts3D(end);
+options.TAVG_1   = data.TJOB_SE(i_+1);%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.NMVA     = 100;              % 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;
@@ -51,7 +52,8 @@ if 0
 options.INTERP    = 1;
 options.POLARPLOT = 0;
 % options.NAME      = '\phi';
-options.NAME      = 'N_i^{00}';
+options.NAME      = '\omega_z';
+% options.NAME      = 'N_i^{00}';
 % options.NAME      = 'v_y';
 % options.NAME      = 'n_i^{NZ}';
 % options.NAME      = '\Gamma_x';
@@ -59,33 +61,35 @@ options.NAME      = 'N_i^{00}';
 options.PLAN      = 'xy';
 % options.NAME      = 'f_i';
 % options.PLAN      = 'sx';
-options.COMP      = 'avg';
+% options.COMP      = 'avg';
 % options.TIME      = data.Ts5D(end-30:end);
 % options.TIME      =  data.Ts3D;
-options.TIME      = [0:500];
+options.TIME      = [100:1:300];
 data.EPS          = 0.1;
 data.a = data.EPS * 2000;
-create_film(data,options,'.gif')
+options.RESOLUTION = 256;
+create_film(data,options,'.avi')
 end
 
 if 1
 %% 2D snapshots
 % Options
-options.INTERP    = 0;
+options.INTERP    = 1;
 options.POLARPLOT = 0;
-options.AXISEQUAL = 0;
-% options.NAME      = '\phi';
+options.AXISEQUAL = 1;
+options.NAME      = '\phi';
 % options.NAME      = '\psi';
-% options.NAME      = 'n_e';
-options.NAME      = 'N_i^{00}';
+% options.NAME      = 'n_i';
+% options.NAME      = 'N_i^{00}';
 % options.NAME      = 'T_i';
 % options.NAME      = '\Gamma_x';
 % options.NAME      = 'k^2n_e';
-options.PLAN      = 'kxky';
+options.PLAN      = 'xy';
 % options.NAME      'f_i';
 % options.PLAN      = 'sx';
 options.COMP      = 'avg';
-options.TIME      = [1000 1100 1200];
+options.TIME      = [200];
+
 data.a = data.EPS * 2e3;
 fig = photomaton(data,options);
 % save_figure(data,fig)
@@ -95,10 +99,10 @@ if 0
 %% 3D plot on the geometry
 options.INTERP    = 0;
 options.NAME      = '\phi';
-options.PLANES    = [1:1:16];
+options.PLANES    = [1];
 options.TIME      = [30];
 options.PLT_MTOPO = 1;
-options.PLT_FTUBE = 1;
+options.PLT_FTUBE = 0;
 data.EPS = 0.4;
 data.rho_o_R      = 3e-3; % Sound larmor radius over Machine size ratio
 fig = show_geometry(data,options);

wk/header_3D_results.m

@@ -26,16 +26,7 @@ helazdir = '/home/ahoffman/HeLaZ/';
 % outfile ='shearless_cyclone/64x32x16x9x5_CBC_120';
 
 % outfile = 'shearless_cyclone/64x32x16x5x3_CBC_CO/64x32x16x5x3_CBC_LRGK';
-%% ZPINCH
-% outfile ='Zpinch_rerun/Kn_2.5_200x48x5x3';
-% outfile ='Zpinch_rerun/Kn_2.5_256x128x5x3';
-% outfile ='Zpinch_rerun/Kn_2.5_312x196x5x3_Lx_400_Ly_200';
-% outfile ='Zpinch_rerun/Kn_2.5_256x64x5x3';
-% outfile ='Zpinch_rerun/Kn_2.0_200x48x9x5_large_box';
-% outfile ='Zpinch_rerun/Kn_2.0_256x64x9x5_Lx_240_Ly_120';
-% outfile ='Zpinch_rerun/Kn_1.6_256x128x7x4';
-% outfile ='Zpinch_rerun/Kn_1.6_200x48x11x6';
-% outfile ='Zpinch_rerun/Kn_1.6_256x128x21x11';
+
 
 %% CBC
 % outfile = 'CBC/64x32x16x5x3';