comp_hxr_jd

0001 function comp_hxr_jd(hxr_path,filename,p_opt,n1,n2)
0002 %
0003 if nargin < 5,
0004     n1 = 3;
0005     n2 = 3;
0006 end
0007 if nargin < 3,
0008     p_opt = -1;
0009 end
0010 %
0011 %load([hxr_path,'HXR_RESULTS_',id,'.mat'],'Zbremdiag','hxrexp');
0012 %
0013 if strfind(filename,'HXR_RESULTS'),
0014     id = filename(13:strfind(filename,'.mat')-1);
0015     load([hxr_path,filename],'Zbremdiag','hxrexp');
0016     %
0017     kfit = Zbremdiag.kfit;
0018     brem_bd_diag = Zbremdiag.brem_bd_diag;
0019     %
0020 elseif strfind(filename,'PROC_LUKE'),
0021     id = filename(19:strfind(filename,'.mat')-1);
0022     load([hxr_path,filename],'data_proc');
0023     %
0024     kfit = data_proc.hxr.kfit;
0025     brem_bd_diag = data_proc.hxr.brem_bd_diag;
0026     hxrexp = data_proc.hxr.hxrexp.';
0027     %
0028 end
0029 %
0030 nc = size(brem_bd_diag,1);
0031 leg = {'sim.','exp.'};
0032 %
0033 figure(1), clf
0034 %
0035 siz = 14+10*1i;
0036 %
0037 for ik = 1:length(kfit),
0038     %
0039     %figure(ik), clf
0040     subplot(n1,n2,ik)
0041     %
0042     graph1D_jd(1:nc,brem_bd_diag(:,ik),0,0,'Chord number','dN_{pulse}/dt [s^{-1}]',['Pulse count rate, E=',num2str(kfit(ik)),' keV'],NaN,'','','-','none','b',2,siz,gca,0.9);
0043     graph1D_jd(1:nc,hxrexp(ik,1:nc),0,0,'','','',leg,'','','-','none','r',2,siz);
0044     %
0045     %keyboard
0046     %
0047 end
0048 %
0049 print_jd(p_opt,['Fig_hxrcomp_',id],[hxr_path,'figures'],1);
0050 %
0051