0001 function fcomp_sync(Ec,tn,f_nosync,tRR_nosync,tnorm_nosync,f_sync,tRR_sync,tnorm_sync,l_opt,graph,p_opt,id_simul)
0002
0003 ntn = length(tn);
0004
0005 if l_opt == 0,
0006 xlim = [0,tn(end)];
0007 mask = 1:ntn;
0008 else
0009 xlim = [tn(2),tn(end)];
0010 mask = 2:ntn;
0011 end
0012
0013 siz = 20+14i;
0014 tit = '';
0015
0016 figure(1),clf
0017
0018 xlab = 't/\tau_c';
0019 ylab = 'n/n_e';
0020 if isfield(graph,'ylim_norm'),
0021 ylim = graph.ylim_norm;
0022 else
0023 ylim = 10.^[-9,1];
0024 graph.ytick_norm = 10.^(-8:2:0);
0025 end
0026
0027 leg = {'ext. RE','int. RE','bulk','total'};
0028
0029 graph1D_jd(tn(mask),tnorm_nosync.re(mask),l_opt > 0,1,'','','',NaN,xlim,ylim,'-','none','r',0.5,siz,gca);
0030 graph1D_jd(tn(mask),tnorm_nosync.ri(mask),l_opt > 0,1,'','','',NaN,xlim,ylim,'-','none','b',0.5,siz,gca);
0031 graph1D_jd(tn(mask),tnorm_nosync.b(mask),l_opt > 0,1,'','','',NaN,xlim,ylim,'-','none',[0,0.5,0],0.5,siz,gca);
0032 graph1D_jd(tn(mask),tnorm_nosync.tot(mask),l_opt > 0,1,xlab,ylab,tit,leg,xlim,ylim,'-','none','k',0.5,siz,gca,0.9,0.7,0.7);
0033
0034 graph1D_jd(tn(mask),tnorm_sync.re(mask),l_opt > 0,1,'','','',NaN,xlim,ylim,'--','none','r',2,siz,gca);
0035 graph1D_jd(tn(mask),tnorm_sync.ri(mask),l_opt > 0,1,'','','',NaN,xlim,ylim,'--','none','b',2,siz,gca);
0036 graph1D_jd(tn(mask),tnorm_sync.b(mask),l_opt > 0,1,'','','',NaN,xlim,ylim,'--','none',[0,0.5,0],2,siz,gca);
0037 graph1D_jd(tn(mask),tnorm_sync.tot(mask),l_opt > 0,1,'','','',NaN,xlim,ylim,'--','none','k',2,siz,gca);
0038
0039 if isfield(graph,'ytick_norm'),
0040 set(gca,'ytick',graph.ytick_norm)
0041 end
0042 set(gca,'YMinorGrid','off')
0043 set(gca,'YMinorTick','on')
0044
0045 figure(2),clf
0046
0047 xlab = 't/\tau_c';
0048 ylab = '\Gamma_R\tau_c';
0049 if isfield(graph,'ylim_gamma'),
0050 ylim = graph.ylim_gamma;
0051 else
0052 ylim = 10.^[-12,-2];
0053 graph.ytick_gamma = 10.^(-12:2:-2);
0054 end
0055
0056 leg = {'ext. RR','int. RR','total RR'};
0057
0058 graph1D_jd(tn(2:ntn),tRR_nosync.re,l_opt > 0,1,'','','',NaN,xlim,ylim,'-','none','r',0.5,siz,gca);
0059 graph1D_jd(tn(2:ntn),tRR_nosync.ri,l_opt > 0,1,'','','',NaN,xlim,ylim,'-','none','b',0.5,siz,gca);
0060 graph1D_jd(tn(2:ntn),tRR_nosync.re + tRR_nosync.ri,l_opt > 0,1,xlab,ylab,tit,leg,xlim,ylim,'-','none','k',0.5,siz,gca,0.9,0.7,0.7);
0061
0062 graph1D_jd(tn(2:ntn),tRR_sync.re,l_opt > 0,1,'','','',NaN,xlim,ylim,'--','none','r',2,siz,gca);
0063 graph1D_jd(tn(2:ntn),tRR_sync.ri,l_opt > 0,1,'','','',NaN,xlim,ylim,'--','none','b',2,siz,gca);
0064 graph1D_jd(tn(2:ntn),tRR_sync.re + tRR_sync.ri,l_opt > 0,1,'','','',NaN,xlim,ylim,'--','none','k',2,siz,gca);
0065
0066 if isfield(graph,'ytick_gamma'),
0067 set(gca,'ytick',graph.ytick_gamma)
0068 end
0069 set(gca,'YMinorGrid','off')
0070 set(gca,'YMinorTick','on')
0071
0072 figure(3),clf
0073
0074 fM = f_nosync{1}(:,end);
0075 f_nosync = f_nosync{end}(:,end);
0076 f_sync = f_sync{end}(:,end);
0077
0078 xlab = 'E_c (MeV)';
0079 ylab = 'f(E_c,\xi = 1)';
0080 xlim = [0,max(Ec)];
0081 if isfield(graph,'ylim_f'),
0082 ylim = graph.ylim_f;
0083 else
0084 ylim = [eps,max(fM)];
0085 graph.ytick_f = 10.^(-18:2:2);
0086 end
0087 tit = '';
0088 siz = 20+14i;
0089
0090 leg = {'f_M','f (no synch.)','f (with synch.)'};
0091 graph1D_jd(Ec,fM,0,1,'','','',NaN,xlim,ylim,'-','none','b',0.5,siz);
0092 graph1D_jd(Ec,f_nosync,0,1,'','','',NaN,xlim,ylim,'-','none','r',0.5,siz);
0093 graph1D_jd(Ec,f_sync,0,1,xlab,ylab,tit,leg,xlim,ylim,'--','none',[0,0.5,0],2,siz,gca,0.9,0.7,0.7);
0094
0095 if isfield(graph,'ytick_f'),
0096 set(gca,'ytick',graph.ytick_f)
0097 end
0098 set(gca,'YMinorGrid','off')
0099 set(gca,'YMinorTick','on')
0100
0101 print_jd(p_opt,['fig_',id_simul,'_nr'],'./figures',1)
0102 print_jd(p_opt,['fig_',id_simul,'_rr'],'./figures',2)
0103 print_jd(p_opt,['fig_',id_simul,'_f1'],'./figures',3)
0104 