- Timestamp:
- Apr 18, 2016, 2:51:08 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
CLIO/CTR+gfw/CTR.m
r465 r492 24 24 pulse_len_ps=2.5;%Nominal bunch duration (FWHM), ps 25 25 gam=100;%relativistic gamma 26 26 basefig=0;%base number of figures 27 icolor='r';%color of lines on plot 27 28 if (~isempty(var)) 28 29 disp(var); … … 91 92 E_total_1=sum(CTR_spec_N(imin:end))*(omega(2)-omega(1)); 92 93 if PlotFunc(1)==1 93 figure( 10)94 loglog(3e+8./nu*1e+3,FFT_spec_om, 'k','linew',2)94 figure(basefig+1) 95 loglog(3e+8./nu*1e+3,FFT_spec_om,icolor,'linew',2) 95 96 hold on 96 97 % loglog(3e+8./(nu2(1:round(length(nu))))*1e+3,FFT_spec_om,'r','linew',2) … … 106 107 107 108 if PlotFunc(2)==1 108 figure( 20)109 loglog(c./nu*1e+3,CTR_spec_1*1e+18, 'r','linew',2)109 figure(basefig+2) 110 loglog(c./nu*1e+3,CTR_spec_1*1e+18,icolor,'linew',2) 110 111 hold on 111 112 grid on … … 117 118 118 119 if PlotFunc(3)==1 119 figure( 30)120 semilogy(c./nu*1e+3,CTR_spec_N*1e+18, 'r','linew',2)120 figure(basefig+3) 121 semilogy(c./nu*1e+3,CTR_spec_N*1e+18,icolor,'linew',2) 121 122 set(gca,'fontsize',16) 122 123 ylabel('dI/d\omega, [\muJ/THz]') … … 134 135 135 136 if PlotFunc(4)==1 136 figure( 40)137 plot(time*1e+3,Profile./max(Profile), 'r','linew',2)137 figure(basefig+4) 138 plot(time*1e+3,Profile./max(Profile),icolor,'linew',2) 138 139 set(gca,'fontsize',16) 139 140 title(['Profile, fwhm=' num2str(pulse_len_ps) ' ps']) … … 146 147 147 148 if PlotFunc(5)==1 148 figure( 50)149 figure(basefig+5) 149 150 [h h]=contourf(c./nu*1e+3,Theta*1e+3,d2Udisk.*1e+18,40); 150 151 set(h,'LineStyle','none');
Note: See TracChangeset
for help on using the changeset viewer.