w vectors at our indices and plot, overlay with the vector representation of the new velocities

figure;

pcolorjw(g.lon, g.lat, Uobj.magnitude(itime, klev, :, :));
colorbar;
arrows(g.lon(1:end,1:end), g.lat(1:end,1:end),...
    Uobj.vectors(itime, klev, 1:end,1:end), 0.08, 'black');
title(datestr(g.time));
dasp(44);

%% Subset to jj,ii range: New Plot / New Vectors

figure;
% [ U, g ] = cgrid_uv2rho(nc, uname, vname, hname, aname, itime, klev, 1:58, 1:70);
pcolorjw(g.lon(1:58, 1:70), g.lat(1