% show header
ncdisp('http://icdc.cen.uni-hamburg.de/thredds/dodsC/ecco2-temp-all')
% get potential temperature variable at time step 170, depth level 1
theta=ncread('http://icdc.cen.uni-hamburg.de/thredds/dodsC/ecco2-temp-all', 'THETA',[1 1 1 170], [1440 720 1 1]);
% get lat and lon 
lat=ncread('http://icdc.cen.uni-hamburg.de/thredds/dodsC/ecco2-temp-all', 'LATITUDE_T',[1], [720]);
lon=ncread('http://icdc.cen.uni-hamburg.de/thredds/dodsC/ecco2-temp-all', 'LONGITUDE_T',[1], [1440]);
% set missing value to NaN
theta(theta<-1000)=NaN
% plot map
imagesc(lon,lat,theta')
axis xy