% Read MODIS TERRA or AQUA data distribuidos no formato NETCDF

file='A20180012018365.L3m_YR_SST_sst_9km.nc';

% O comando ncdisp mostra na tela toda a informacao nos metadados.
% Verifique essas informacoes para descobrir os nomes e tamanhos da
% variáveis lon, lat e SST.
% Apos descobrir essas informacoes, comente o comando para nao ficar mostrando
% todas essas informacoes na tela.

%   ncdisp(file)

info=ncinfo(file);
% Informações sobre os dados podem ser acessados como Atributos, seus nomes e valores.
% Por exemplo, o valor 62 do info.Attributes.Value: info.Attributes(62).Value mostra
% o valor máximo do dado de temperatura na superfície.
% O info.Attributes.Name mostra o que o valor significa.

maxsst=info.Attributes(62).Value; % valor máximo da SST

% Ler os vetores de Lat e Lon

lon=ncread(file,'lon');
lat=ncread(file,'lat');

% Ler a matrix de dados
data=ncread(file,'sst'); 
sst=double(data');

% Remover valores maiores que o maxsst
L=sst>maxsst;
sst(L)=nan;

imagesc(lon,lat,sst);axis('xy');c=colorbar;
cc=get(c,'title');
set(cc,'string','^oC');
xlabel('longitude');
ylabel('latitude');

print -djpeg modis_sst