tring">'Gr1kmTM_v1.nc'</span>,<span class="string">'hAm'</span>,<span class="string">'constituents'</span>,<span class="string">'m2'</span>);

<span class="comment">% Load the global ocean mask:</span>
ocean = tmd_data(<span class="string">'Gr1kmTM_v1.nc'</span>,<span class="string">'mask'</span>);

figure
h=imagesc(x,y,M2);
h.AlphaData = ocean; <span class="comment">% makes land transparent</span>
axis <span class="string">xy</span> <span class="string">image</span> <span class="comment">% sets axes equal 