scripts update

wk/analysis_3D.m

@@ -10,9 +10,9 @@ system(['mkdir -p ',MISCDIR]);
 CMD = ['rsync ', LOCALDIR,'outputs* ',MISCDIR]; disp(CMD);
 system(CMD);
 % Load outputs from jobnummin up to jobnummax
-JOBNUMMIN = 07; JOBNUMMAX = 20;
+JOBNUMMIN = 10; JOBNUMMAX = 11;
 data = compile_results(MISCDIR,JOBNUMMIN,JOBNUMMAX); %Compile the results from first output found to JOBNUMMAX if existing
-
+data.localdir = LOCALDIR;
 
 %% PLOTS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 default_plots_options
@@ -40,7 +40,7 @@ end
 if 0
 %% MOVIES %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % Options
-options.INTERP    = 1;
+options.INTERP    = 0;
 options.POLARPLOT = 0;
 options.NAME      = '\phi';
 % options.NAME      = 'N_i^{00}';
@@ -48,12 +48,12 @@ options.NAME      = '\phi';
 % options.NAME      = 'n_i^{NZ}';
 % options.NAME      = '\Gamma_x';
 % options.NAME      = 'n_i';
-options.PLAN      = 'xz';
+options.PLAN      = 'yz';
 % options.NAME      = 'f_e';
 % options.PLAN      = 'sx';
 options.COMP      = 'avg';
 % options.TIME      = dat.Ts5D;
-options.TIME      = 1100:0.1:1170;
+options.TIME      = 1250:0.1:1370;
 data.EPS          = 0.1;
 data.a = data.EPS * 2000;
 create_film(data,options,'.gif')
@@ -65,17 +65,17 @@ if 0
 options.INTERP    = 0;
 options.POLARPLOT = 0;
 options.AXISEQUAL = 1;
-options.NAME      = '\phi';
+% options.NAME      = '\phi';
 % options.NAME      = 'n_i';
-% options.NAME      = 'N_i^{00}';
+options.NAME      = 'N_i^{00}';
 % options.NAME      = 'T_i';
 % options.NAME      = '\Gamma_x';
 % options.NAME      = 'k^2n_e';
-options.PLAN      = 'xz';
+options.PLAN      = 'yz';
 % options.NAME      = 'f_i';
 % options.PLAN      = 'sx';
-options.COMP      = 9;
-options.TIME      = [1100 1150];
+options.COMP      = 'avg';
+options.TIME      = [1300 1350];
 data.a = data.EPS * 2e3;
 fig = photomaton(data,options);
 save_figure(data,fig)

wk/analysis_header.m

@@ -5,8 +5,8 @@ helazdir = '/home/ahoffman/HeLaZ/';
 outfile ='';
 outfile ='';
 outfile ='';
-outfile ='shearless_cyclone/128x128x16_12_3_L_120_CBC_1.0';
-% outfile ='shearless_cyclone/128x128x16xdmax_L_120_CBC_1.0';
+outfile ='shearless_cyclone/128x128x16_CBC_SG_study';
+% outfile ='shearless_cyclone/128x128x16xdmax_6_L_120_CBC_1.0';
 % outfile ='shearless_cyclone/linear_CBC_100_D_16';
 % outfile ='quick_run/CLOS_1_64x64_5x3_L_120_kN_2.0_kT_0.5_nu_1e-01_SGGK';
 % outfile ='pedestal/64x64x16x2x1_L_300_LnT_20_nu_0.1';

wk/debug_script.m

-system('cd ../results/debug/test_diag; mpirun -np 8 ./helaz3_dbg 2 2 2; cd $HOME/HeLaZ/wk');
+% system('cd ../results/dev/test_diag; ./helaz3_dbg; cd $HOME/HeLaZ/wk');
+% system('cd ../results/dev/test_diag; mpirun -np 2 ./helaz3_dbg 1 2 1; cd $HOME/HeLaZ/wk');
+system('cd ../results/dev/test_diag; mpirun -np 8 ./helaz3_dbg 2 2 2; cd $HOME/HeLaZ/wk');
 
 
-filename = '../results/debug/test_diag/outputs_00.h5';
+filename = '../results/dev/test_diag/outputs_00.h5';
 
+% test phi
 f_        = h5read(filename,'/data/var3d/phi/000004/');
 f_gatherv = h5read(filename,'/data/var3d/phi_gatherv/000004/');
 f_        = f_.real        + 1i*f_.imaginary;
@@ -10,9 +13,18 @@ f_gatherv = f_gatherv.real + 1i*f_gatherv.imaginary;
 err = sum(sum(sum(abs(f_-f_gatherv))))/sum(sum(sum(abs(f_))));
 disp(['error_phi = ',sprintf('%2.2e',err)]);    
 
+% test Ni00
 f_        = h5read(filename,'/data/var3d/Ni00/000004/');
 f_gatherv = h5read(filename,'/data/var3d/Ni00_gatherv/000004/');
 f_        = f_.real        + 1i*f_.imaginary;
 f_gatherv = f_gatherv.real + 1i*f_gatherv.imaginary;
 err = sum(sum(sum(abs(f_-f_gatherv))))/sum(sum(sum(abs(f_))));
 disp(['error_Ni00 = ',sprintf('%2.2e',err)]);    
+
+% test Nipj
+f_        = h5read(filename,'/data/var5d/moments_i/000004/');
+f_gatherv = h5read(filename,'/data/var5d/moments_i_gatherv/000004/');
+f_        = f_.real        + 1i*f_.imaginary;
+f_gatherv = f_gatherv.real + 1i*f_gatherv.imaginary;
+err = sum(sum(sum(sum(sum(abs(f_-f_gatherv))))))/sum(sum(sum(sum(sum(abs(f_))))));
+disp(['error_Nipj = ',sprintf('%2.2e',err)]);    

wk/quick_run.m

@@ -113,7 +113,7 @@ end
 
 if 0
 %% Ballooning plot
-options.time_2_plot = data.Ts3D(end);
+options.time_2_plot = [0.9 1]*data.Ts3D(end);
 options.kymodes     = [0.5];
 options.normalized  = 1;
 options.field       = 'phi';