"comment">%Whole domain:</span>
<span class="comment">% [ U, g ] = cgrid_uv2rho(nc, uname, vname, hname, aname, itime, klev);</span>
Uobj = hvar.getvectors(uvar, vvar, avar);
</pre><h2>Grab the data at the indices of interest<a name="4"></a></h2><pre class="codeinput">itime = 3; <span class="comment">% 3rd time step</span>
klev = -1; <span class="comment">% last (top) layer</span>
<span class="comment">% Grab coordinate information at our indices</span>
g = Uobj.grid(itime, klev, :, :);
</pre><h2>Initialize