clear all;close all

d=dir('dt_global_twosat_msla_h_201*.nc');

nt=length(d);

t=zeros(1,nt);
eta=zeros(nt,1440);

for i=1:length(d)
display(['Processing file ' num2str(i) ' of ' num2str(nt) '.'])

file=d(i).name;
    
% Coloque aqui os comandos para ler tempo, x (lon), y (lat),
% e z (anomalia da altura da superfície).
  
time=
x=
y=
z=
    
% O que faz essa parte?
  
iy=find(y==9.875);
t(i)=time+datenum(1950,1,1)-1;
eta(i,:)=z(iy,:);

end

figure
imagesc(x,t,eta);