up

wk/HD_study.m

 % 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

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; 

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