From d35c4f2064f70ebbaa9e39edad229a16999987f7 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Thu, 2 Jun 2022 15:27:31 +0200
Subject: [PATCH] scripts update
---
wk/analysis_3D.m | 20 ++++++++++----------
wk/analysis_header.m | 4 ++--
wk/debug_script.m | 16 ++++++++++++++--
wk/quick_run.m | 2 +-
4 files changed, 27 insertions(+), 15 deletions(-)
diff --git a/wk/analysis_3D.m b/wk/analysis_3D.m
index feea3ac3..5b429ee4 100644
--- a/wk/analysis_3D.m
+++ b/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)
diff --git a/wk/analysis_header.m b/wk/analysis_header.m
index 55ad7a20..e347fe84 100644
--- a/wk/analysis_header.m
+++ b/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';
diff --git a/wk/debug_script.m b/wk/debug_script.m
index 3954d020..1d53f6cb 100644
--- a/wk/debug_script.m
+++ b/wk/debug_script.m
@@ -1,8 +1,11 @@
-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)]);
diff --git a/wk/quick_run.m b/wk/quick_run.m
index 05284974..6435cd4b 100644
--- a/wk/quick_run.m
+++ b/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';
--
GitLab