ain: 

ocean = tmd_interp('Gr1kmTM_v1.nc','mask',Lat,Lon); 

wct(~ocean) = NaN; 

figure
pcolor(lon,lat,wct)
shading interp 
cb = colorbar; 
ylabel(cb,'water column thickness (m)')
cmocean deep % optional colormap 

%%
% Sometimes you don't want to interpolate, and you just want the raw data. 
% Use <tmd_data_documentation.html |tmd_data|> to get things like height
% coefficients: 

% Load complex height coefficients:
[h,x,y,cons] = tmd_data('Gr1kmTM_v1.nc','h'); 

% Examine outputs: 
whos h x y cons

%%  
