grd=roms_get_grid('../frc2008/swf_grid6.nc');

roms_plot_mesh(grd,1)
roms_plot_bathy(grd,'zebra(2,64)',[1 2 3 4 5 6:3:50])

%shell pt indices
i=129;
j=163;

hold on

plot3(grd.lon_rho(i),grd.lat_rho(j),100,'k*')


