~, ~] = size(hycom.(currvar).data);
0136             <span class="keyword">if</span> iscolumn(hycom.Depth.data)
0137                 hdepth = permute(repmat(-hycom.Depth.data, [1, hx, hy]), [2, 3, 1]);
0138             <span class="keyword">else</span>
0139                 hdepth = -hycom.Depth.data;
0140             <span class="keyword">end</span>
0141 
0142             <span class="comment">% Make a land mask of the HYCOM domain (based on the surface</span>
0143             <span class="comment">% layer 