
   clear all; close all;

   mini=17;
   delta=0.5;
   maxi=28;

   dur_plot='2003-2007';

   file='/loceanfs/pulsation/cholod/data/reynolds/sst_reynolds_weekly-mean_1981-10-29_2010-12-26.pc.nc';

   name='REYNOLDS';
   name_title=regexprep(name,'_','\\_');

%%% Lecture du fichier REYNOLDS
   disp(['lecture sst dans ' file ' ...']);
   nc=netcdf(file,'r'); 
   ntime=length(nc('time_counter'));
   lat=nc{'lat'}(:);
   lon=nc{'lon'}(:)-360;
   tos=nc{'sst'}(:,:,:);
   close(nc);

%  scale_factor=0.01
   tos=tos .* 0.01;

%  950  :  2 Jan 2000
%  1106 : 29 Dec 2002
%  1367 : 30 Dec 2007
   tos2=tos(1106:1367,:,:);

%  moyenne sur la première dimension et conserve le nombre de dimension
   mtos=squeeze(mean(tos2(:,:,:),1));

   opengl neverselect % pour pb cartes graphiques

%  lon(10) = -100.5000
%  lon(11) = -99.5000
%  lon(40) = -70.5000
%  lon(41) = -69.5000
%  lat(15) = 10.5000
%  lat(16) = 9.5000
%  lat(55) = -29.5000
%  lat(56) = -30.5000

%%%  Plot 1
   figure;

   pcolor(lon(10:41),lat(15:56),mtos(15:56,10:41));
   shading flat;

%  contourf(lon(10:41),lat(15:56),mtos(15:56,10:41),[mini:delta:maxi]);
%  hold on
%  contour(lon(10:41),lat(15:56),mtos(15:56,10:41),[mini:delta:maxi]);

   set(gca,'FontS',14); set(gca,'PlotBoxAspectRatio',[0.7 1 1]); box on;
   hcc=colorbar; set(hcc,'FontS',14);
   caxis ([mini maxi])
   title(['sst mean (C) - ' name_title ' (' dur_plot ')'],'FontSize',18);
   figname1=['map_sst_mean_' name '_' dur_plot];
   print ('-dpng',figname1)



