% read woa13 data downloaded from 
% http://www.nodc.noaa.gov/
% Read the readme file.
%
% these data are distributed as netcdf files, in a 180x360
% maps for 102 levels (5500m) at extended standard depths for T. 
%
gap=[.01,.05]; marg_h=[.01,.01]; marg_w=[.03,.01];

% Construct the standard depths, latitude and longitude vectors
dep=ncread('woa18_decav_t00_01.nc','depth');
lat=...;
lon=...;

% Read the temperature obtained by objective analysis
tem=...;

% Use permute to change the orientation of the map
tem=permute(tem,[2 1 3]);

%  Cut the study region: Mediterraneo

ila=find(lat>=30&lat<=45);
ilo=find(lon>=-5&lon<=40);
rla=lat(ila);
rlo=lon(ilo);
rtem=tem(ila,ilo,:);

% plot the SST map
subtightplot(1,1,1,gap,marg_h,marg_w);
imagesc(rlo,rla,squeeze(rtem(:,:,1)));axis('xy','equal','tight');t=colorbar;
title(['Temperatura da superfície do mar'])
xlabel('longitude')
ylabel('latitude')
tt=get(t,'title');
set(tt,'string','(^oC)')

%print -djpeg sst00