%matplotlib inline
import matplotlib.pyplot as plt
import netCDF4 as nc
import numpy as np
from salishsea_tools import bathy_tools

bathy_orig = nc.Dataset('../../nemo-forcing/grid/bathy_meter_SalishSea2.nc', 'r')
lats = bathy_orig.variables['nav_lat']
lons = bathy_orig.variables['nav_lon']

fig = bathy_tools.plot_colourmesh(
    bathy_orig, 'VENUS, o-model-T-grid, *-ADCP', 
    axis_limits=(-123.6, -122.9, 48.9, 49.2))
plt.plot(lons[416,282],lats[416,282],'ow')
plt.plot(lons[424,267],lats[416,282],'ow')
plt.plot(-123-(19.055/60.),49+(2.531/60.),'*r')
plt.plot(-123-(25.5287/60.),49+(2.4083/60.),'*r')