Skip to content
Snippets Groups Projects
Commit 52e22443 authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

save scripts

parent 5dd82b38
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
matlab/plot/plot_cosol_mat.m
+ 17
− 11
View file @ 52e22443
......@@ -103,6 +103,7 @@ mfns = {...
'/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_12_P_4_J_2_N_50_kpm_4.0.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landau_P10_J5_dk_5e-2_km_2.0_NFLR_12.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landau_P11_J7_dk_5e-2_km_2.0_NFLR_16.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landau_P12_J6_dk_5e-2_km_1.0_NFLR_16.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landauii_P16_J9_dk_5e-2_km_2.0_NFLR_8.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landau_P16_J9_dk_5e-2_km_2.0_NFLR_8.h5';...
% '/home/ahoffman/gyacomo/iCa/gk.hacked_sugama_P_10_J_5_N_150_kpm_8.0.h5';...
......@@ -115,6 +116,7 @@ CONAME_A = {...
'FC 4 2 NFLR 12'; ...
'LD 10 5 NFLR 12'; ...
'LD 11 7 NFLR 16'; ...
'LD 12 6 NFLR 16'; ...
'LDii 16 9 NFLR 8'; ...
'LD 16 9 NFLR 8'; ...
% 'Hacked SG A';...
......@@ -131,7 +133,7 @@ for j_ = 1:numel(mfns)
matidx = sprintf('%5.5i',idx_);
MAT = h5read(mat_file_name,['/',matidx,'/Caapj/Ciipj']);
gmax(idx_+1) = max((real(eig(MAT))));
wmax(idx_+1) = max(abs(imag(eig(MAT))));
wmax(idx_+1) = max((imag(eig(MAT))));
end
subplot(121)
plot(kp_a,gmax,'DisplayName',CONAME_A{j_}); hold on;
......@@ -140,11 +142,13 @@ for j_ = 1:numel(mfns)
end
subplot(121)
legend('show'); grid on;
ylim([0,100]);
ylim([-10,10]);
xlabel('$k_\perp$'); ylabel('$\gamma_{max}$ from Eig(iCa)')
subplot(122)
legend('show'); grid on;
xlabel('$k_\perp$'); ylabel('$\omega_{max}$ from Eig(iCa)')
ylim([-10,10]);
end
%% Van Kampen plot
......@@ -157,20 +161,22 @@ mfns = {...
'/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_6_P_4_J_2_N_50_kpm_4.0.h5';...
'/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_12_P_4_J_2_N_50_kpm_4.0.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landau_P10_J5_dk_5e-2_km_2.0_NFLR_12.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landau_P12_J6_dk_5e-2_km_1.0_NFLR_16.h5';...
'/home/ahoffman/gyacomo/iCa/gk_landauii_P16_J9_dk_5e-2_km_2.0_NFLR_8.h5';...
'/home/ahoffman/gyacomo/iCa/gk_sugama_tau1e-3_P11_J7_dk_5e-2_km_5.0_NFLR_12.h5';...
'/home/ahoffman/gyacomo/iCa/gk_sugama_tau1e-3_P4_J2_dk_5e-2_km_5.0_NFLR_5.h5';...
'/home/ahoffman/gyacomo/iCa/gk_sugama_tau1e-2_P4_J2_dk_5e-2_km_5.0_NFLR_5.h5';...
% '/home/ahoffman/gyacomo/iCa/gk_sugama_tau1e-3_P11_J7_dk_5e-2_km_5.0_NFLR_12.h5';...
% '/home/ahoffman/gyacomo/iCa/gk_sugama_tau1e-3_P4_J2_dk_5e-2_km_5.0_NFLR_5.h5';...
% '/home/ahoffman/gyacomo/iCa/gk_sugama_tau1e-2_P4_J2_dk_5e-2_km_5.0_NFLR_5.h5';...
};
CONAME_A = {'SG 20 10';...
'PA 20 10';...
'FC 4 2 NFLR 6';...
'FC 4 2 NFLR 12';...
'LD 10 5 NFLR 12';...
'LD 12 6 NFLR 12';...
'LD 16 9 NFLR 16';...
'SG 11 7 NFLR 12, tau 1e-3'; ...
'SG 4 2 NFLR 5, tau 1e-3'; ...
'SG 4 2 NFLR 5, tau 1e-2'; ...
% 'SG 11 7 NFLR 12, tau 1e-3'; ...
% 'SG 4 2 NFLR 5, tau 1e-3'; ...
% 'SG 4 2 NFLR 5, tau 1e-2'; ...
};
TAU_A = [1;...
1;...
......@@ -178,9 +184,9 @@ TAU_A = [1;...
1;...
1;...
1;...
1e-3;...
1e-3;...
1e-2;...
% 1e-3;...
% 1e-3;...
% 1e-2;...
];
grow = {};
puls = {};
......
wk/analysis_heat_flux_paper2.m 0 → 100644
+ 22
− 0
View file @ 52e22443
figure
PJ = '5x3';
% PJ = '9x5';
% PJ = '11x6';
% PJ = '17x9';
% RES = '128x64x24';
RES = '192x96x24';
NU = '0.005';
PJxRES = [PJ,'x',RES];
for CO = {'DGGK','SGGK','LDGK'}
folder = ['/misc/gyacomo23_outputs/paper_2_GYAC23/collision_study/nu',CO{1},'_scan_kT_5.3/',PJxRES,'/nu_',NU];
[t, ~, Qxi] = read_flux_out_XX(folder);
plot(t,Qxi,'.','DisplayName',CO{1}); hold on
end
legend('show');
title([PJ,'x',RES,' nu=',NU]);
ylim([0 25])
\ No newline at end of file
wk/load_metadata_scan.m
+ 5
− 1
View file @ 52e22443
......@@ -10,7 +10,11 @@ addpath(genpath([gyacomodir,'matlab/load'])) % ... add% EXECNAME = 'gyacomo_1.0'
% datafname = 'lin_KBM/12x24_ky_0.05_0.75_P_2_16_DGGK_0.01_be_0.03.mat';
% datafname = 'lin_AUG_scan/12x24_ky_0.05_0.75_P_2_16_DGGK_0.01_be_0.000152.mat';
% datafname = 'lin_Entropy_scan/2x1_ky_0.05_0.75_P_2_8_DGDK_0_be_0.mat';
datafname = 'lin_DTT_AB_rho85_PT_scan/16x24_ky_0.1_1.5_P_2_8_DGGK_0.05_be_0.0034.mat';
% datafname = 'lin_DTT_AB_rho85_PT_scan/16x24_ky_0.1_0.5_P_2_8_DGGK_0.05_be_0.0034.mat';
% datafname = 'lin_DTT_AB_rho85_PT_scan/16x24_ky_0.1_0.5_P_2_8_DGGK_0.05_be_0.0039.mat';
% datafname = 'lin_DTT_AB_rho85_PT_scan/16x24_ky_0.1_1_P_2_8_DGGK_0.05_be_0.0039.mat';
datafname = 'lin_DTT_AB_rho85_PT_scan/16x24_ky_0.1_0.9_P_2_8_kN_1.33_DGGK_0.56901_be_0.0039.mat';
% datafname = 'lin_DTT_AB_rho85_PT_scan/8x24_ky_0.05_1.5_P_2_8_kN_1.33_DGGK_0.1_be_0.0039.mat';
%% Chose if we filter gamma>0.05
FILTERGAMMA = 0;
......
wk/paper_2_scripts_and_results/nonlin_kT_scan_analysis.m
+ 3
− 1
View file @ 52e22443
......@@ -12,6 +12,8 @@ else
rootdir = '/misc/gene_results/kT_scan_nu0';
mrkstyl='d';
end
nustr = '5e-3'; mrkstyl='o';
rootdir = '/misc/gyacomo23_outputs/paper_2_GYAC23/kT_scan_LDGK_nu_5e-3';
Mmax = 6;
msz = 10; lwt = 2.0;
xname = '$\kappa_T (\kappa_N=2.22)$';
......@@ -182,7 +184,7 @@ end
ylabel('$\chi L_N/\rho_s^2 c_s$');
xlabel('$R/L_T$');
title(['$\nu_{DGDK}=$',nustr])
title(['$\nu_{',data.inputs.CONAME,'}=$',nustr])
legend('show');
legend('Location','northwest')
xlim([3.5 10]);
......
wk/paper_2_scripts_and_results/nonlin_nu_scan_analysis.m
+ 19
− 13
View file @ 52e22443
......@@ -2,20 +2,20 @@ kN=2.22;
figure
ERRBAR = 1; LOGSCALE = 0; AU = 0;
resstr={};
kN = 2.22;
msz = 10; lwt = 2.0;
% CO = 'DGGK'; mrkstyl='d';
% CO = 'SGGK'; mrkstyl='s';
CO = 'LDGK'; mrkstyl='o';
CO = 'SGGK'; mrkstyl='s';
% CO = 'LDGK'; mrkstyl='o';
% GRAD = 'kT_7.0'; kT = 6.96;
GRAD = 'kT_5.3'; kT = 5.3;
% GRAD = 'kT_4.5'; kT = 4.5;
xname = ['$\nu_{',CO,'}$ '];
titlename = [CO,', ',GRAD];
scanvarname = 'nu';
rootdir = ['/misc/gyacomo23_outputs/paper_2_GYAC23/collision_study/nu',CO,'_scan_',GRAD];
scanval = {'0.005' '0.01' '0.02' '0.05' '0.1' '0.2' '0.5'};
naming = @(s) num2str(s);
J0 = 00; J1 = 20;
% Get all directories
system(['ls -d ',rootdir,'/*/ > list.txt']);
fid = fopen('list.txt');
......@@ -39,6 +39,7 @@ clrs_ = lines(numel(directories));
for j = 1:numel(directories)
% Get all subdirectories
disp(directories{j});
system(['ls -d ',directories{j},'*/ > list.txt']);
fid = fopen('list.txt');
tline = fgetl(fid); i_ = 1; nus=[]; subdirectories={};
......@@ -56,7 +57,7 @@ for j = 1:numel(directories)
subdirectories = subdirectories(ids); nus = nus(ids);
naming = @(s) sprintf('%1.1f',s); clr_ = clrs_(j,:);
titlename = [CO,', ',GRAD,', ',resstr{j}];
titlename = [CO,', kT=',num2str(kT)];
N = numel(subdirectories);
x = 1:N;
Qx_avg = 1:N;
......@@ -65,9 +66,12 @@ for j = 1:numel(directories)
Chi_std = 1:N;
data = {};
for i = 1:N
subdir = subdirectories{i};
data = compile_results_low_mem(data,subdir,00,20);
subdir = subdirectories{i};
% data = compile_results_low_mem(data,subdir,J0,J1);
try
[t, ~, Q] = read_flux_out_XX(subdir);
data.Ts0D = t;
data.HFLUX_X = Q;
Trange = data.Ts0D(end)*[0.5 0.75];
catch % if data does not exist put 0 everywhere
data.Ts0D = 0;
......@@ -80,16 +84,16 @@ for j = 1:numel(directories)
data.inputs.NU = nus(i);
end
Trange = data.Ts0D(end)*[0.25 0.75];
% Trange = [200 400];
% Trange = [100 400];
%
[~,it0] = min(abs(Trange(1) -data.Ts0D));
[~,it1] = min(abs(Trange(end)-data.Ts0D));
%
Qx_avg(i) = mean(data.HFLUX_X(it0:it1));
Qx_std(i) = std(data.HFLUX_X(it0:it1));
Chi_avg(i) = Qx_avg(i)./data.inputs.K_T/data.inputs.K_N;
Chi_std(i) = Qx_std(i)./data.inputs.K_T/data.inputs.K_N;
x(i) = data.inputs.NU;
Chi_avg(i) = Qx_avg(i)./kT/kN;
Chi_std(i) = Qx_std(i)./kT/kN;
x(i) = nus(i);
subplot(N,2,2*i-1)
hold on;
Qx = data.HFLUX_X;
......@@ -102,19 +106,21 @@ for j = 1:numel(directories)
'Color',clr_); hold on
% plot([T(it0) T(end)],Qx_avg(i)*[1 1],'--k','DisplayName',...
% ['$Q_{avg}=',sprintf('%2.2f',Qx_avg(i)),'\pm',sprintf('%2.2f',Qx_std(i)),'$']);
ylim([0 40]);
end
% plot chi vs nu
subplot(122)
hold on;
if ERRBAR
errorbar(x,Chi_avg,Chi_std,'DisplayName',...
['(',num2str(data.inputs.PMAX),',',num2str(data.inputs.JMAX),')'],...
['(',num2str(Ps(j)-1),',',num2str(Js(j)-1),') ',resstr{j}],...
'color',clr_,'MarkerFaceColor','k','MarkerSize',7,'LineWidth',2); hold on;
else
plot(x,Chi_avg,'DisplayName',...
['(',num2str(data.inputs.PMAX),',',num2str(data.inputs.JMAX),')'],...
['(',num2str(Ps(j)-1),',',num2str(Js(j)-1),') ',resstr{j}],...
'color',clr_,'Marker',mrkstyl);
end
ylim([0 8])
end
% Formatting
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment