% grd = roms_get_grid(grd_file);

value=0
lon_psi(mask_psi==value) = [];
lat_psi(mask_psi==value) = [];
lon_rho(mask_rho==value) = [];
lat_rho(mask_rho==value) = [];
lon_u(mask_u==value) = [];
lat_u(mask_u==value) = [];
lon_v(mask_v==value) = [];
lat_v(mask_v==value) = [];

plot(lon_psi,lat_psi,'k',lon_psi.',lat_psi.','k');
hold on;
plot(lon_psi(:),lat_psi(:),'r.',...
     lon_rho(:),lat_rho(:),'go',...
       lon_u(:),  lat_u(:),'bx',...
       lon_v(:),  lat_v(:),'cs');
legend('psi','rho','u','v')
hold off
xfac=cos(mean(lat_rho(:))*pi/180);
daspect([ 1 xfac 1] );