Changeset 541 in ETALON for CLIO/CTR+gfw

Timestamp:

Apr 22, 2016, 9:59:42 PM (8 years ago)

Author:

delerue

Message:

CLIO paper updated

Location:

CLIO/CTR+gfw

Files:

• CTR.m (modified) (1 diff)
• GFWnew/GFWnew.m (copied) (copied from CLIO/CTR+gfw/GFWnew/GFWnew.m) (9 diffs)

CLIO/CTR+gfw/CTR.m

@@ -126 +126 @@
  title(['CTR spectrum'])
  grid on
+ xlim([0 5])
 %  xlim([min(c./nu*1e+3) 2])
 end

CLIO/CTR+gfw/GFWnew/GFWnew.m

@@ -63 +63 @@
     load([P_file]);
   else
-      disp('Set correct file name!')
+      error('Set correct file name!')
   end
 else
@@ -70 +70 @@
 sigma_t= pulse_len_ps/(sqrt(2*log(2))*(1+epsl))*1e-3;  %The sigma of the leading edge    ns
 
-nu=(1:1e+10:3e+13).*1e-9;%Ghz
+nu=(1:1e+10:3e+16).*1e-9;%Ghz
 
 sigma_L= sigma_t*epsl;
@@ -220 +220 @@
 FreqHalf=Frequency(1:round(length(Frequency)/2));%half of frequency range
 SpectrumHalf=SpectrumHalf./max(SpectrumHalf);%normalization of Spectrum
-Nmax=find(FreqHalf>=2.5,1);%2.5 THz is more than enoght
+Nmax=find(FreqHalf>=60,1);%2.5 THz is more than enoght
 % Interpolation is applyed to decrease error in calculation (freq step in FFT is to wide for this calculation ) 
 FreqInter=FreqHalf(1):((FreqHalf(2)-FreqHalf(1))*1e-3):FreqHalf(Nmax);%interpolation diapason
@@ -226 +226 @@
 if PlotFunc(1)==1
 figure(1)
+clf
 plot(FreqInter,SpecInter,'-k','linew',2)
 hold on
@@ -252 +253 @@
  if isempty(s_coh_long)
      s_coh_long=0;%ACHTUNG!!!
-     disp('fail to find Form Factro value')
+     error('fail to find Form Factor value')
  end
 else
@@ -367 +368 @@
 if PlotFunc(2)==1
  figure(2)
+ clf
  plot(thet,bin_single(:,phiDiv+2),'k','linew',2)
  set(gca,'fontsize',16)
@@ -376 +378 @@
 if PlotFunc(3)==1
 figure(3)
+clf
  plot(thet,sum_tr,'r','linew',2)
 % thet_rad=thet/180*pi;
@@ -391 +394 @@
 if PlotFunc(4)==1
 figure(4)
+clf
 plot(time,Profile,'r','linew',2)
  set(gca,'fontsize',16)
@@ -401 +405 @@
 if PlotFunc(5)==1
 figure(5)
+clf
 semilogy(nu_Hz*1e-12,SP_spec_1,'r','linew',2)
  hold on