Changeset 491 in ETALON for CLIO

Timestamp:

Apr 18, 2016, 2:47:58 PM (8 years ago)

Author:

hodnevuc

Message:

 

File:

• CLIO/CTR+gfw/GFWnew/GFWnew.m (modified) (7 diffs)

CLIO/CTR+gfw/GFWnew/GFWnew.m

@@ -54 +54 @@
 phiDiv=3;%steps in phi
 saveflag=0;%set to 1 if you want to save model parameters
+basefig=0;%base number of figures
+icolor='r';%color of lines on plot
 if (~isempty(var))
 strread(var,'%s','delimiter',';')
@@ -134 +136 @@
 SpecInter=pchip( FreqHalf(1:Nmax),SpectrumHalf(1:Nmax),FreqInter);%interpolation
 if PlotFunc(1)==1
-figure(1)
-plot(FreqInter,SpecInter,'-k','linew',2)
+figure(basefig+1)
+plot(FreqInter,SpecInter,icolor,'linew',2)
 hold on
 grid on
@@ -298 +300 @@
 
 if PlotFunc(1)==1
-figure(1)
+figure(basefig+1)
 hold on
-plot(c./(lam*1e-3)*1e-12,FFact,'or','linew',2)
+plot(c./(lam*1e-3)*1e-12,FFact,'ok','linew',2)
 legend('Form factor','Used points')
 end
@@ -382 +384 @@
 %% Plots
 if PlotFunc(2)==1
- figure(2)
- plot(thet,bin_single(:,phiDiv+2),'k','linew',2)
+ figure(basefig+2)
+ plot(thet,bin_single(:,phiDiv+2),icolor,'linew',2)
  set(gca,'fontsize',16)
 ylabel('dI_1/d\Omega, [J sr^{-1}cm^{-1}]')
@@ -391 +393 @@
 end
 if PlotFunc(3)==1
-figure(3)
- plot(thet,sum_tr,'r','linew',2)
+figure(basefig+3)
+ plot(thet,sum_tr,icolor,'linew',2)
 % thet_rad=thet/180*pi;
 % polar(thet_rad(:),sum_tr(:));
@@ -406 +408 @@
 
 if PlotFunc(4)==1
-figure(4)
-plot(time,Profile,'r','linew',2)
+figure(basefig+4)
+plot(time,Profile,icolor,'linew',2)
  set(gca,'fontsize',16)
 ylabel('Amplitude, [a.u.]')
@@ -416 +418 @@
 
 if PlotFunc(5)==1
-figure(5)
-semilogy(nu_Hz*1e-12,SP_spec_1,'r','linew',2)
+figure(basefig+5)
+semilogy(nu_Hz*1e-12,SP_spec_1,icolor,'linew',2)
  hold on
 %  load CTR.mat