source: ETALON/CLIO/Analysis/GFW/Corrections.m @ 715

close all; clear all; clc
load CLIOSpec.mat

phi=0:0.5:6;
for c=2:(size(bin,2)-2)
   bin2(:,c-1)=bin_single(:,c);
%    bin2(:,c-1)=bin(:,c);

end
 SEY=[flip(bin2(:,2:end)')',bin2]';

% figure
%  [hc hc]=contourf(thet,[ -flip(phi(2:end)),phi],SEY)
%   set(hc,'LineStyle','none');
phiNew=-6:0.1:6;
 for c=1:size(SEY,2)
     SEYi(:,c)=pchip([ -flip(phi(2:end)),phi],SEY(:,c),phiNew);
 end
 
%  figure
%  [hc hc]=contourf(thet,phiNew,SEYi)
%   set(hc,'LineStyle','none');
%   
 Thet=30:0.1:140;

  for c=1:size(SEYi,1)
     SEYij(c,:)=pchip(thet,SEYi(c,:), Thet);
  end
   figure
 [hc hc]=contourf(Thet,phiNew,SEYij)
  set(hc,'LineStyle','none');
  TH=48:7:125
  
  Aed_mm=25;%Aperture entrance diameter,  mm 
lgr_mm=180;%Grating length, mm
el_mm=6;

NormKV=[1.3458    1.2208    1.1325    1.0711    1.0305    1.0072    1.0002    1.0072 1.0305    1.0711    1.1325    1.2208];%180 length grating, 25 input mirror, 48:7:125
lgr_e_mm = abs(Aed_mm./sin(Thet/180*pi));           % variable effective grating length, determined basically by the cone entrance diameter of XXXmm*/
lgr_e_mm(lgr_e_mm>lgr_mm)=lgr_mm;
CR=0.1.*lgr_e_mm(:);%J;


SEYij=SEYij.*repmat(CR,[1 size(SEYij,1)])';

 for  c=1:12
  
load(['map' num2str(TH(c)) '.mat']) 

MB=Map(:,41:end-40)'./max(max(Map));

CSEY=SEYij(:,[find(Thet>=TH(c)-10&Thet<=TH(c)+10)]);
  figure
  contourf(MB.*CSEY)
  Integral(c)=sum(sum(MB.*CSEY))*(0.1/180*pi)^2.*NormKV(c);%*
% SEYij(:,find(Thet>=38&Thet<=58))
 end
 figure
 hold on
 plot(TH,Integral)
Spectr=Integral
 
%   save SPRP15.mat Integral
%    save SPRP10.mat  Spectr