comp_luke_REregion_jd

0001 clear
0002 p_opt = -1;
0003 c_opt = true;
0004 xi_opt = true;
0005 pmax = 20;
0006 
0007 lukedata_1 = load('DKE_RESULTS_Runaway_alpha3_nosync.mat');
0008 lukedata_2 = load('DKE_RESULTS_Runaway_alpha3_sync_0p6.mat');
0009 
0010 xlim = [-3,9];
0011 xtick = -3:3:9;
0012 ylim = [0,6];
0013 ytick = [0:2:6];
0014 %
0015 
0016 if c_opt,
0017     colors = {'r','b','r','b','r'};
0018     styles = {'-','--','-','--','-.'};
0019     sizes = [2,2,0.5,0.5,0.5];
0020     colstr = '_col';
0021 else
0022     colors = {'k','k','k','k','k'};
0023     styles = {'-','--','-','--','-.'};
0024     sizes = [2,2,0.5,0.5,0.5];
0025     colstr = '';
0026 end
0027 
0028 p_1 = lukedata_1.fgrid.pn*sqrt(lukedata_1.betath2);
0029 p_2 = lukedata_2.fgrid.pn*sqrt(lukedata_2.betath2);
0030 
0031 if ~isfield(lukedata_1,'pmhucrit') || isempty(lukedata_1.pmhucrit) || all(isnan(lukedata_1.pmhucrit)),
0032     pmhucrit_1 = NaN(1,length(p_1));
0033     for ip = 1:length(p_1),
0034         pmhucrit_1(ip) = min([1,lukedata_1.fgrid.mhu(lukedata_1.tacc{end}(ip,:,:))]);
0035     end
0036 else
0037     pmhucrit_1 = lukedata_1.pmhucrit;
0038 end
0039 
0040 if ~isfield(lukedata_2,'pmhucrit') || isempty(lukedata_2.pmhucrit) || all(isnan(lukedata_2.pmhucrit)),
0041     pmhucrit_2 = NaN(1,length(p_2));
0042     for ip = 1:length(p_2),
0043         pmhucrit_2(ip) = min([1,lukedata_2.fgrid.mhu(lukedata_2.tacc{end}(ip,:,:))]);
0044     end
0045 else
0046     pmhucrit_2 = lukedata_2.pmhucrit;
0047 end
0048 
0049 ppar_1 = p_1.*pmhucrit_1;
0050 pperp_1 = p_1.*sqrt(1 - pmhucrit_1.^2);
0051 
0052 ppar_2 = p_1.*pmhucrit_2;
0053 pperp_2 = p_1.*sqrt(1 - pmhucrit_2.^2);
0054 
0055 pc = 1/sqrt(lukedata_1.alpha - 1);
0056 pperp0 = sqrt((lukedata_2.alpha - 1)/lukedata_2.taurinv);
0057 
0058 ppar_xi = p_1/lukedata_1.alpha;
0059 pperp_xi = p_1*sqrt(1 - 1/lukedata_1.alpha^2);
0060 
0061 % PLOT
0062 
0063 figure(1),clf
0064 
0065 xlab = 'p_{||}';
0066 ylab = 'p_{\perp}';
0067 tit = '';
0068 leg = {'without ALD','with ALD'};
0069 siz = 20+14i;
0070 
0071 graph1D_jd(ppar_1,pperp_1,0,0,'','','',NaN,xlim,ylim,styles{1},'none',colors{1},sizes(1),siz,gca);
0072 graph1D_jd(ppar_2,pperp_2,0,0,xlab,ylab,tit,leg,xlim,ylim,styles{2},'none',colors{2},sizes(2),siz,gca);
0073 graph1D_jd([pc,pc],ylim,0,0,'','','',NaN,xlim,ylim,styles{3},'none',colors{3},sizes(3),siz,gca);
0074 graph1D_jd(xlim,[pperp0,pperp0],0,0,'','','',NaN,xlim,ylim,styles{4},'none',colors{4},sizes(4),siz,gca);
0075 if xi_opt,
0076     graph1D_jd(ppar_xi,pperp_xi,0,0,'','','',NaN,xlim,ylim,styles{5},'none',colors{5},sizes(5),siz,gca);
0077     text(0,ylim(2)/2,'\xi=E_{||}^{-1}','fontsize',20)
0078     xistr = '_xiopt';
0079 else
0080     xistr = '';
0081 end
0082 axis('equal');axis([xlim,ylim]);
0083 text(pc,ylim(2)/2,'p_{||}=p_c','fontsize',20)
0084 text(xlim(1)/2,pperp0,'p_{\perp}=p_{\perp0}','fontsize',20)
0085 set(gca,'xtick',xtick)
0086 set(gca,'ytick',ytick)
0087 
0088 print_jd(p_opt,['fig_luke_Reregion',xistr,colstr],'figures/',1)