, and plot again:</p><pre class="codeinput">ocean = tmd_interp(<span class="string">'Gr1kmTM_v1.nc'</span>,<span class="string">'mask'</span>,Lat,Lon);

wct(~ocean) = NaN;

figure
pcolor(lon,lat,wct)
shading <span class="string">interp</span>
cb = colorbar;
ylabel(cb,<span class="string">'water column thickness (m)'</span>)
cmocean <span class="string">deep</span> <span class="comment">% optional colormap</span>
</pre><img vspace="5" hspace="5" src="tmd_getting_started_06.png" alt=""> <p>Sometimes you don't