clear all;
close all;


ncload('CERES_EBAF_TOA_Terra_Edition1A_200003-200510.nc');

meansw=squeeze(nanmean(swup,1));   % total cloud vector
meanlw=squeeze(nanmean(lwup,1));   % total cloud vector
meantot=squeeze(nanmean(net,1));   % total cloud vector



% Planisphere (miller graph)
figure
subplot(3,1,1)
m_proj('miller','clongitude',0);
m_pcolor(lon,lat,meansw)


 y(1:length(lon))=0;
m_line(lon,y','color','black')
 y(1:length(lon))=30;
m_line(lon,y','color','black')
 y(1:length(lon))=60;
m_line(lon,y','color','black')
 y(1:length(lon))=-60;
m_line(lon,y','color','black')
 y(1:length(lon))=-30;
m_line(lon,y','color','black')

m_coast('color','black') ; % pour les cotes continents
%caxis([0 1])
m_grid('xtick', [], 'ytick', []);
shading flat;     % pour ne pas mettre de grille
colorbar;
title('CERES-EBAF SW for 200003-200510')


% Planisphere (miller graph)
subplot(3,1,2)
m_proj('miller','clongitude',0);
m_pcolor(lon,lat,meanlw)

 y(1:length(lon))=0;
m_line(lon,y','color','black')
 y(1:length(lon))=30;
m_line(lon,y','color','black')
 y(1:length(lon))=60;
m_line(lon,y','color','black')
 y(1:length(lon))=-60;
m_line(lon,y','color','black')
 y(1:length(lon))=-30;
m_line(lon,y','color','black')

m_coast('color','black') ; % pour les cotes continents
%caxis([0 1])
m_grid('xtick', [], 'ytick', []);
shading flat;     % pour ne pas mettre de grille
colorbar;
title('CERES-EBAF LW for 200003-200510')



% Planisphere (miller graph)
subplot(3,1,3)
m_proj('miller','clongitude',0);
m_pcolor(lon,lat,meantot)

 y(1:length(lon))=0;
m_line(lon,y','color','black')
 y(1:length(lon))=30;
m_line(lon,y','color','black')
 y(1:length(lon))=60;
m_line(lon,y','color','black')
 y(1:length(lon))=-60;
m_line(lon,y','color','black')
 y(1:length(lon))=-30;
m_line(lon,y','color','black')

m_coast('color','black') ; % pour les cotes continents
%caxis([0 1])
m_grid('xtick', [], 'ytick', []);
shading flat;     % pour ne pas mettre de grille
colorbar;
title('CERES-EBAF TOT for 200003-200510')



