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_05.png" alt=""> <p>Above, the water column thickness is zero wherever there's land. Perhaps you'd prefer to set the land grid cells to NaN rather than zero. Here's how: Just interpolate the ocean mask, set the non-ocean pixels to NaN