%% Heat flux Qi
kN = 2.22;
%-------------- GM ---------------
%(P,J)=(2,1)
kT_Qi_GM_32 = ...
[...
% 13. 1.5e+2 1.6e+1;...%192x96x16x3x2 kymin=0.02
% 11. 5.1e+2 3.5e+2;...%192x96x16x3x2 kymin=0.02
9.0 9.6e+1 3.0e+1;...%192x96x16x3x2 kymin=0.05
7.0 5.0e+1 6.6e+0;...%192x96x16x3x2 kymin=0.05
6.0 3.0e+1 4.8e+0;...%192x96x16x3x2 kymin=0.05
5.0 1.1e+1 9.4e-1;...%192x96x16x3x2 kymin=0.05
4.5 9.2e+0 1.6e+0;...%192x96x16x3x2 kymin=0.05
];
%(P,J)=(4,2)
kT_Qi_GM_53 = ...
[...
% 13. 2.0e+2 1.2e+1;...%96x64x16x3x2 kymin=0.02 (large box)
% 13. 1.1e+2 2.0e+1;...%128x64x16x5x3 kymin=0.02 (large box)
% 13. 1.3e+2 3.5e+1;...%128x64x16x5x3 kymin=0.05
% 11. 1.2e+2 1.6e+1;...%96x64x16x3x2 kymin=0.02 (large box)
% 11. 1.6e+2 1.8e+1;...%128x64x16x5x3 kymin=0.02
% 11. 9.7e+1 2.2e+1;...%128x64x16x5x3 kymin=0.05
9.0 8.3e+1 2.2e+1;...%128x64x16x5x3 kymin=0.05
% 9.0 7.6e+1 2.3e+1;...%256x128x16x3x2 kymin=0.05 (high res)
7.0 4.6e+1 2.3e+0;...%128x64x16x5x3 kymin=0.05
6.0 3.7e+1 6.9e+0;...%128x64x16x5x3 kymin=0.05
5.3 1.9e+1 2.0e+0;...%128x64x16x5x3 kymin=0.05
5.0 1.3e+1 3.3e+0;...%128x64x16x5x3 kymin=0.05
4.5 9.3e+0 1.0e+0;...%128x64x16x5x3 kymin=0.05
];
%(P,J)=(12,2) or higher
kT_Qi_GM_HD = ...
[...
7.0 4.1e+1 6.6e+0;...%192x96x24x13x7 kymin=0.05
5.3 7.5e+0 1.8e+0;...%192x96x24x13x7 kymin=0.05
4.5 1.1e+0 4.0e-1;...%192x96x24x13x7 kymin=0.05
4.5 9.6e-1 1.5e-1;...%128x64x16x13x2 kymin=0.05
4.5 7.9e-1 1.8e-1;...%128x64x16x13x7 kymin=0.05
4.5 1.2e+0 5.4e-1;...%192x64x24x6x4 kymin=0.05 ! Lx is too small... (weird oscillations)
];
%-------------- GENE ---------------
kT_Qi_GENE_SR = ...
[...
13. 2.7e+2 2.2e+1;...%128x64x16x24x12 kymin=0.02 (large box)
% 13. 2.0e+2 6.6e+1;...%128x64x16x24x12 kymin=0.05
11. 1.9e+2 1.7e+1;...%128x64x16x24x12 kymin=0.02 (large box)
% 11. 3.3e+2 1.6e+2;...%128x64x16x24x12 kymin=0.05
9.0 1.1e+2 4.2e+1;...%128x64x16x24x12 kymin=0.05
7.0 3.5e+1 4.6e+0;...%128x64x16x24x12 kymin=0.05
5.3 9.7e+0 6.8e+0;...%128x64x16x24x12 kymin=0.05
4.5 2.3e-1 5.0e-2;...%128x64x16x24x12 kymin=0.05
];
kT_Qi_GENE_HR = ...
[...
5.3 3.8e-1 1.7e-1;...%128x64x16x32x16 Nexc=5 kymin=0.05
];
%% Heat conductivity Xi [Ln/rhoi^2/cs] computed as Xi = Qi/kT/kN
%init
kT_Xi_GM_32 = kT_Qi_GM_32;
kT_Xi_GM_53 = kT_Qi_GM_53;
kT_Xi_GM_HD = kT_Qi_GM_HD;
kT_Xi_GENE = kT_Qi_GENE_SR;
%scale
for i = 2:3
kT_Xi_GM_32 (:,i) = kT_Qi_GM_32 (:,i)./kT_Qi_GM_32 (:,1)./kN;
kT_Xi_GM_53 (:,i) = kT_Qi_GM_53 (:,i)./kT_Qi_GM_53 (:,1)./kN;
kT_Xi_GM_HD (:,i) = kT_Qi_GM_HD (:,i)./kT_Qi_GM_HD (:,1)./kN;
kT_Xi_GENE (:,i) = kT_Qi_GENE_SR (:,i)./kT_Qi_GENE_SR (:,1)./kN;
end
%% Dimits fig 3 data
KT_DIM = [4.0 4.5 5.0 6.0 7.0 9.0 12. 14. 16. 18.];
LLNL_GK_DIM = [5.0 0.0; ...
7.0 2.5;...
9.0 5.0;...
12. 8.0;...
14. 9.0;...
16. 9.5;...
18. 10.];
UCOL_GK_DIM = [5.0 0.5;...
6.0 1.0;...
7.0 1.5];
GFL__97_DIM = [4.0 4.0;...
4.5 5.0;...
5.0 6.5;...
7.0 8.0;...
9.0 11.];
GFL__98_DIM = [5.0 1.5;...
7.0 5.0;...
9.0 7.5];
%% Plot
msz = 8.0; lwt = 2.0;
figure;
if 1
xye = kT_Xi_GM_32;
% errorbar(xye(:,1), xye(:,2),xye(:,3),...
plot(xye(:,1), xye(:,2),...
'>-','DisplayName','192x96x16x3x2',...
'MarkerSize',msz,'LineWidth',lwt); hold on
xye = kT_Xi_GM_53;
% errorbar(xye(:,1), xye(:,2),xye(:,3),...
plot(xye(:,1), xye(:,2),...
'<-','DisplayName','128x64x16x5x3',...
'MarkerSize',msz,'LineWidth',lwt); hold on
xye = kT_Xi_GM_HD;
% errorbar(xye(:,1), xye(:,2),xye(:,3),...
plot(xye(:,1), xye(:,2),...
'^-','DisplayName','128x64x16x13x3',...
'MarkerSize',msz,'LineWidth',lwt); hold on
xye = kT_Xi_GENE;
% errorbar(xye(:,1), xye(:,2),xye(:,3),...
plot(xye(:,1), xye(:,2),...
'+-.k','DisplayName','GENE 128x64x16x24x12',...
'MarkerSize',msz,'LineWidth',lwt); hold on
end
if 1
plot(LLNL_GK_DIM(:,1),LLNL_GK_DIM(:,2),'dk--','DisplayName','Dimits GK, LLNL'); hold on
plot(UCOL_GK_DIM(:,1),UCOL_GK_DIM(:,2),'*k','DisplayName','Dimits PIC, U.COL');
plot(GFL__97_DIM(:,1),GFL__97_DIM(:,2),'+k','DisplayName','Dimits GFL 97');
plot(GFL__98_DIM(:,1),GFL__98_DIM(:,2),'sk','DisplayName','Dimits GFL 98');
end
plot([6.96 6.96],[0 12],'-.r','DisplayName','CBC');
plot([4.0 4.00],[0 12],'-.b','DisplayName','$\kappa_T^{crit}$');
xlabel('$\kappa_T$'); ylabel('$\chi_i$[$L_n/\rho_i^2 v_{thi}$]');
xlim([0 20]);
ylim([0 12]);
title('Dimits et al. 2000, Fig. 3');
legend('show'); grid on;