%%  the hope that  it  will  be  useful,  but  without  any  %
%%  warranty.  You  can  redistribute  it and/or modify  it  %
%%  under  the  terms  of  the  GNU  General Public License  %
%%  as  published  by  the   Free  Software  Foundation  at  %
%%       http://www.gnu.org/copyleft/gpl.html                %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Octave script plot SST
%
% 24/10/2015: andrea.doglioli@univ-amu.fr
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function plot_SST(fname,sst_interval)

load(fname);

sst_min=sst_interval(1);
sst_max=sst_interval(2);

sst(sst<sst_min)=sst_min;
sst(sst>sst_max)=sst_max;

%[C,hc]=m_contourf(lon,lat,sst);
%set(hc,'EdgeColor','none') 

m_pcolor(lon,lat,sst)
shading flat

caxis([sst_min sst_max])

%h.LevelList=[16:0.1:20]
%h.LevelListMode='manual'

%cMap = jet; cMap = [cMap; 1 1 1] ; colormap(cMap); 
colormap(jet);
colorbar

title(fname,'interpreter','none')


return
