Where shall we evaluate the storm surge performance of the model? The station needs to have water level data for the period of 2002-2010 (this is when we have wind data).
Point Atkinson (49.34,-123.25)
Victoria Harbour (48.42,-123.37)
Patricia Bay (48.65,-123.45)
Campbell River (50.04,-125.25)
All other sites are not in our domain or do not have water level data for these points
from salishsea_tools import tidetools
bathy, X, Y = tidetools.get_SS_bathy_data()
x1, y1 = tidetools.find_closest_model_point(-123.25,49.34,X,Y,bathy)
print('Point Atkinson: i='+str(x1)+', j='+str(y1))
x1, y1 = tidetools.find_closest_model_point(-123.37,48.42,X,Y,bathy)
print('Victoria Harbour: i='+str(x1)+', j='+str(y1))
x1, y1 = tidetools.find_closest_model_point(-123.45,48.65,X,Y,bathy)
print('Patricia Bay: i='+str(x1)+', j='+str(y1))
x1, y1 = tidetools.find_closest_model_point(-125.25,50.04,X,Y,bathy)
print('Campbell River: i='+str(x1)+', j='+str(y1))
Point Atkinson: i=468, j=328 Victoria Harbour: i=298, j=195 Patricia Bay: i=351, j=213 Campbell River: i=747, j=123