diff --git a/wk/HD_study.m b/wk/HD_study.m
index baa0b2b00d1b4fe4f2ff215a15befe8f248bf481..57920a911780f81152fc5281a8033c62d1d55773 100644
--- a/wk/HD_study.m
+++ b/wk/HD_study.m
@@ -1,11 +1,13 @@
% bursts simulations
b_=[...
+ 5e-1, 5e+0;...
1e-1, 1e-1;...
1e-1, 3e-2;...
1e-1, 2e-2;...
1e-1, 1e-2;...
- 1e-2, 3e-3;...
5e-2, 3e-3;...
+ 5e-2, 2e-3;...
+ 1e-2, 3e-3;...
1e-3, 2e-2;...
];
% converged plateau simulations
@@ -26,29 +28,23 @@ dp_=[...
1e-3, 5e-3;...
1e-3, 2.5e-3,...
];
-% resolution issue plateau?
+% not sure
rp_=[...
5e-2, 1e-3;...
+ 1e-2, 1e-3;...
1e-2, 1e-4;...
];
figure; set(gcf, 'Position', [100, 100, 900, 400])
grid on; xlim([5e-4,5e0]); ylim([5e-5,5e+0]);
set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
xlabel('$\nu_{DGGK}$'); ylabel('$\mu_{HD}$'); hold on;
-
-scatter( b_(:,1), b_(:,2),'o',...
- 'MarkerFaceColor',[0.6350 0.0780 0.1840],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
- 'DisplayName','Bursts');
-scatter(cp_(:,1),cp_(:,2),'s',...
- 'MarkerFaceColor',[0.4660 0.6740 0.1880],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
- 'DisplayName','Converged Plateau');
-scatter(dp_(:,1),dp_(:,2),'v',...
- 'MarkerFaceColor',[0.4660 0.6740 0.1880],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
- 'DisplayName','Moving Plateau');
-scatter(rp_(:,1),rp_(:,2),'h',...
- 'MarkerFaceColor',[0.9290 0.6940 0.1250],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
- 'DisplayName','Moving Plateau');
-% plot([1 1],[5e-5 5e-1],'--','Color',[0.4660 0.6740 0.1880]);
+
+% Trajectory of simulations
+% HD_study/100x50_L_50_P_2_J_1_eta_0.6_nu_1e+00_DGGK_CLOS_0_mu_1e-02/
+mu_ = [0 1 2 5 5];
+nu_ = [1 1 1 1 0.5];
+plot(nu_,mu_,'x-','DisplayName','100x50, L=050, P,J=2,1,');
+
% Trajectory of simulations
% HD_study/100x50_L_50_P_2_J_1_eta_0.6_nu_1e-01_DGGK_CLOS_0_mu_1e-02/
@@ -69,4 +65,18 @@ plot(nu_,mu_,'x-','DisplayName','150x75, L=100, P,J=2,1,');
% HD_study/150x75_L_100_P_2_J_1_eta_0.6_nu_5e-02_DGGK_CLOS_0_mu_3e-02/
mu_ = [3e-3 1e-3];
nu_ = [5e-2 5e-2];
-plot(nu_,mu_,'x-','DisplayName','150x75, L=100, P,J=2,1,');
\ No newline at end of file
+plot(nu_,mu_,'x-','DisplayName','150x75, L=100, P,J=2,1,');
+
+scatter( b_(:,1), b_(:,2),'o',...
+ 'MarkerFaceColor',[0.6350 0.0780 0.1840],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
+ 'DisplayName','Bursts');
+scatter(cp_(:,1),cp_(:,2),'s',...
+ 'MarkerFaceColor',[0.4660 0.6740 0.1880],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
+ 'DisplayName','Converged Plateau');
+scatter(dp_(:,1),dp_(:,2),'v',...
+ 'MarkerFaceColor',[0.4660 0.6740 0.1880],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
+ 'DisplayName','Moving Plateau');
+scatter(rp_(:,1),rp_(:,2),'h',...
+ 'MarkerFaceColor',[0.9290 0.6940 0.1250],'MarkerEdgeColor',[0 0 0],'SizeData',50,...
+ 'DisplayName','not sure');
+plot([1 1],[5e-5 5e-1],'--','Color',[0.4660 0.6740 0.1880]);
\ No newline at end of file
diff --git a/wk/analysis_2D.m b/wk/analysis_2D.m
index 6626eaa054314f7c18944e7cadd28be64d9c6305..aacdf5c5ac53793f1704e10fea44497e52f8e5b3 100644
--- a/wk/analysis_2D.m
+++ b/wk/analysis_2D.m
@@ -2,13 +2,17 @@ addpath(genpath('../matlab')) % ... add
for i_ = 1
% for ETA_ =[0.6:0.1:0.9]
%% Load results
-if 0% Local results
-% outfile ='HD_study/100x50_L_50_P_2qe_J_1_eta_0.6_nu_1e-01_DGGK_CLOS_0_mu_1e-02';
+if 1% Local results
+outfile ='';
+outfile ='';
+outfile ='';
+outfile ='';
+outfile ='';
+outfile ='';
+outfile ='HD_study/300x150_L_100_P_2_J_1_eta_0.6_nu_1e-01_DGGK_CLOS_0_mu_1e-03';
+% outfile ='HD_study/200x100_L_200_P_2_J_1_eta_0.6_nu_1e+00_DGGK_CLOS_0_mu_0e+00';
% outfile ='HD_study/150x75_L_100_P_2_J_1_eta_0.6_nu_5e-02_DGGK_CLOS_0_mu_3e-03';
-% outfile ='HD_study/200x100_L_100_P_2_J_1_eta_0.6_nu_1e-01_DGGK_CLOS_0_mu_3e-02';
-% outfile ='HD_study/300x150_L_100_P_2_J_1_eta_0.6_nu_1e-01_DGGK_CLOS_0_mu_1e-03';
% outfile ='HD_study/150x75_L_100_P_2_J_1_eta_0.6_nu_1e-02_DGGK_CLOS_0_mu_3e-03';
-outfile ='HD_study/200x100_L_200_P_2_J_1_eta_0.6_nu_1e+00_DGGK_CLOS_0_mu_0e+00';
BASIC.RESDIR = ['../results/',outfile,'/'];
BASIC.MISCDIR = ['/misc/HeLaZ_outputs/results/',outfile,'/'];
CMD = ['cp ', BASIC.RESDIR,'outputs* ',BASIC.MISCDIR]; disp(CMD);
@@ -18,7 +22,6 @@ outfile ='';
outfile ='';
outfile ='';
outfile ='';
-outfile ='';
outfile ='/marconi_scratch/userexternal/ahoffman/HeLaZ/results/kobayashi/300x150_L_100_P_10_J_5_eta_0.71429_nu_1e-02_PAGK_CLOS_0_mu_0e+00/out.txt';
% outfile = outcl{i_};
% load_marconi(outfile);
@@ -404,7 +407,7 @@ end
if 1
%% |phi_k|^2 spectra (Kobayashi 2015 fig 3)
-tstart = 3000; tend = 4000;
+tstart = 3000; tend = 5000;
[~,itstart] = min(abs(Ts2D-tstart));
[~,itend] = min(abs(Ts2D-tend));
trange = itstart:itend;
@@ -445,7 +448,7 @@ clear Z_rth phi_kp ni_kp Ti_kp
end
%%
-t0 =2250;
+t0 =000;
[~, it02D] = min(abs(Ts2D-t0));
[~, it05D] = min(abs(Ts5D-t0));
skip_ = 5;
diff --git a/wk/continue_multiple_runs_marconi.m b/wk/continue_multiple_runs_marconi.m
index 50852ff3a09e5fa084ae9a5b9bfa0be71e3d453d..0e701e36ee5c8c0123b432a56cf4c076e6e44a29 100644
--- a/wk/continue_multiple_runs_marconi.m
+++ b/wk/continue_multiple_runs_marconi.m
@@ -62,7 +62,7 @@ function [] = continue_run(outfilename)
line_= A{43}; mu_old = str2num(line_(13:end));
A{43} = [' mu = ',num2str(0*mu_old)];
% change N
- line_= A{13}; N_old = str2num(line_(10:end));
+ line_= A{13}; N_old = str2num(line_(8:end));
A{13} = [' Nr = ',num2str(N_old)];
A{15} = [' Nz = ',num2str(N_old)];
% change L