clear all
lonm=load('/homedata/dkonsta/modele/grilles_lmdz/longitude_2x2');
latm=load('/homedata/dkonsta/modele/grilles_lmdz/latitude_2x2');
latm=-latm;
lonmax=length(lonm);
latmax=length(latm);


chemin='/bdd/CFMIP/GOCCP/Dimitra/CRef/'
list=dir('/bdd/CFMIP/GOCCP/Dimitra/CRef/CRef*nc')
for ifile=1:length(list)
   ifile
   ncid=netcdf([chemin list(ifile).name])
   cref_mod2(ifile,:,:)=ncid{'CRef_mod250'}(:,:);
   cref_mod1(ifile,:,:)=ncid{'CRef_mod1km'}(:,:);
   cref_par(ifile,:,:)=ncid{'CRef_par'}(:,:);   
end 


cref_mod2(cref_mod2<-100)=NaN;
cref_mod2=squeeze(nanmean(cref_mod2,1));


%for ilat=1:latmax
%  if (latm(ilat)>30 | latm(ilat)<-30) 
%      cref_mod2(:,ilat)=NaN;
%  end
%end

a=[0.75 0 1; 0.5 0 1; 0 0 1; 0 0.5 1; 0 1 1; 0 1 0.5; 0 1 0; 0.5 1 0; 1 1 0; 1 0.5 0; 1 0 0; 1 0 0.5];

lonm(1)=-180;
lonm(180)=180;

figure
m_proj('hammer-aitoff','clongitude',-150)
m_pcolor(lonm,latm,cref_mod2')
hold on
m_pcolor(lonm-360,latm,cref_mod2')
shading flat
colormap(a)
m_coast('color','k')
m_grid('xaxis','middle')
caxis([0 0.6])
title('CLOUD REFLECTANCE MODIS 250M')
colorbar
