clear all; close all;

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

   dur_plot='2003-2007';

   file=['/loceanfs/pulsation/cholod/data/TMI_AMSRE/TMIAMSRE_' dur_plot '_Peru.nc'];
   namerun='TMI_AMSRE';
   namerun_title='TMI\_AMSRE';

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

%%%  Plot 1
   figure;

   pcolor(lon(:,:),lat(:,:),mtos(:,:));
   shading flat;
%  contourf(lon(:,:),lat(:,:),mtos(:,:),[mini:delta:maxi]);
%  hold on
%  contour (lon(:,:),lat(:,:),mtos(:,:),[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) - ' namerun_title ' (' dur_plot ')'],'FontSize',18);
   figname1=['map_sst_mean_' namerun '_' dur_plot];
   print ('-dpng',figname1)


%quit