example_3IncAngles.m

% Script to create the plots TB H-pol for three incidence angles e.g. 20, 30, 40°
% Region min to max coordinates [47.5 7.5 50 10]: Case study Neckar Catchment.
% This script has been used to produce Figure 2 for JORS paper.
% More information: https://github.com/pablosaa/Process_SMOS_L1C
%
% Copiright 2017, Pablo Saavedra G.
%

% Geographical limits:
lat_lon_box = [47.5 7.5 50 10];  % Limits used in Figure 2 JORS paper
latlim = lat_lon_box([1,3]);
lonlim = lat_lon_box([2,4]);
% selecting DBL file to process and plot:
file_example = ['SM_OPER_MIR_SCLF1C_20150702T042618_20150702T051937_620_001_1.DBL'];
if exist(file_example,'file'),
  TSF = Process_SMOSxL1C(file_example,lat_lon_box);  % no MAT-file output!
else
  % Open file browser GUI:
  TSF = Process_SMOSxL1C;
end

idx_pol =1; % 1 for H, 2 for V
inc_ang = [35 40 45]; % Incidence Angles indexes (index=1 is 0 deg.)
N_ang = length(inc_ang);
figure;
set(gcf,'PaperPositionMode','auto','Position',[117 500 1326 395]);
for idx=1:N_ang,
  ax(idx) = subplot(1,N_ang,idx);
  scatter(TSF.GridPoint_Longitude,TSF.GridPoint_Latitude,250,...
	  TSF.TB_Fixed_IncAngle(:,inc_ang(idx)+1,idx_pol),'o','filled');
  tl=title(sprintf('TB_H [K] @ \\theta_{inc} = %2d^\\circ',...
			   inc_ang(idx)),...
	   'FontName','utopia','FontSize',15);
  tmp=get(ax(idx),'XTickLabel');
  set(ax(idx),'XTickLabel',[tmp repmat(' °E',size(tmp,1),1)]);
  tmp=get(ax(idx),'YTickLabel');
  set(ax(idx),'YTickLabel',[tmp repmat(' °N',size(tmp,1),1)]);
  if idx==2, xlabel('LONGITUDE','FontSize',15); end
  if idx==1, ylabel('LATITUDE','FontSize',15); end
end
set(ax,'FontName','utopia','FontSize',15,'LineWidth',1,'Color','none',...
    'Box','on','GridLineStyle','--','XLim',lonlim,'YLim',latlim,...
    'TickDir','out','XMinorTick','on','YMinorTick','on','Layer','top',...
    'LineWidth',1.5,'TickLength',[0.02 0.05],'CLim',[205 285]);

colormap(jet);
hb = colorbar('FontSize',15,'FontName','utopia');
set(hb,'Position',[0.9178 0.1164 0.0105 0.7570]);
return;

% end of script