In [1]:
import numpy as np
import netCDF4 as nc
import matplotlib.pyplot as plt
from salishsea_tools import viz_tools, places,geo_tools
%matplotlib inline
In [3]:
mesh=nc.Dataset('/data/eolson/results/MEOPAR/NEMO-forcing-new/grid/mesh_mask201702.nc')
In [4]:
bathy=nc.Dataset('/data/eolson/results/MEOPAR/NEMO-forcing-new/grid/bathymetry_201702.nc')
In [6]:
mesh.variables.keys()
Out[6]:
dict_keys(['nav_lon', 'nav_lat', 'time_counter', 'tmask', 'umask', 'vmask', 'fmask', 'tmaskutil', 'umaskutil', 'vmaskutil', 'fmaskutil', 'glamt', 'glamu', 'glamv', 'glamf', 'gphit', 'gphiu', 'gphiv', 'gphif', 'e1t', 'e1u', 'e1v', 'e1f', 'e2t', 'e2u', 'e2v', 'e2f', 'ff', 'mbathy', 'misf', 'isfdraft', 'e3t_0', 'e3u_0', 'e3v_0', 'e3w_0', 'gdept_0', 'gdepu', 'gdepv', 'gdepw_0', 'gdept_1d', 'gdepw_1d', 'e3t_1d', 'e3w_1d'])
In [7]:
bathy.variables.keys()
Out[7]:
dict_keys(['nav_lon', 'nav_lat', 'Bathymetry'])
In [16]:
latQU39=50.0307
lonQU39=-125.0992
In [17]:
fig,ax=plt.subplots(1,1,figsize=(8,5))
ax.pcolormesh(bathy['nav_lon'],bathy['nav_lat'],bathy['Bathymetry'])
ax.set_xlim(-125.5,-124)
ax.set_ylim(49.5,50.3)
ax.plot(lonQU39,latQU39,'r*')
viz_tools.set_aspect(ax,coords='map')
<ipython-input-17-473eb2e09dcc>:2: MatplotlibDeprecationWarning: shading='flat' when X and Y have the same dimensions as C is deprecated since 3.3. Either specify the corners of the quadrilaterals with X and Y, or pass shading='auto', 'nearest' or 'gouraud', or set rcParams['pcolor.shading']. This will become an error two minor releases later. ax.pcolormesh(bathy['nav_lon'],bathy['nav_lat'],bathy['Bathymetry'])
Out[17]:
1.1363636363636362
In [14]:
places.PLACES
Out[14]:
{'Campbell River': {'lon lat': (-125.24, 50.04),
'stn number': 8074,
'mean sea lvl': 2.916,
'hist max sea lvl': 5.35,
'wind grid ji': (190, 102),
'NEMO grid ji': (747, 125),
'ww3 grid ji': (453, 109)},
'Cherry Point': {'lon lat': (-122.766667, 48.866667),
'stn number': 9449424,
'mean sea lvl': 3.543,
'hist max sea lvl': 5.846,
'wind grid ji': (122, 166),
'NEMO grid ji': (343, 342),
'ww3 grid ji': (193, 462)},
'Friday Harbor': {'lon lat': (-123.016667, 48.55),
'stn number': 9449880,
'mean sea lvl': 2.561,
'hist max sea lvl': 4.572,
'wind grid ji': (108, 155),
'NEMO grid ji': (300, 267),
'ww3 grid ji': (124, 427)},
'Halfmoon Bay': {'lon lat': (-123.912, 49.511),
'stn number': 7830,
'NEMO grid ji': (549, 254),
'wind grid ji': (158, 136),
'ww3 grid ji': (331, 297),
'mean sea lvl': 3.14,
'hist max sea lvl': 5.61},
'Nanaimo': {'lon lat': (-123.93, 49.16),
'stn number': 7917,
'mean sea lvl': 3.08,
'hist max sea lvl': 5.47,
'wind grid ji': (142, 133),
'NEMO grid ji': (484, 208),
'ww3 grid ji': (261, 298)},
'Neah Bay': {'lon lat': (-124.6, 48.4),
'stn number': 9443090,
'mean sea lvl': 1.925,
'hist max sea lvl': 4.359,
'wind grid ji': (111, 105),
'NEMO grid ji': (384, 15),
'ww3 grid ji': (89, 200)},
'New Westminster': {'lon lat': (-122.90535, 49.203683),
'stn number': 7654,
'mean sea lvl': 1.3,
'hist max sea lvl': 4.66,
'NEMO grid ji': (423, 363),
'wind grid ji': (138, 164)},
'Patricia Bay': {'lon lat': (-123.4515, 48.6536),
'stn number': 7277,
'mean sea lvl': 2.256,
'hist max sea lvl': 4.38,
'NEMO grid ji': (351, 214),
'wind grid ji': (115, 143),
'ww3 grid ji': (145, 363)},
'Point Atkinson': {'lon lat': (-123.25, 49.33),
'stn number': 7795,
'mean sea lvl': 3.09,
'hist max sea lvl': 5.61,
'wind grid ji': (146, 155),
'NEMO grid ji': (468, 329),
'ww3 grid ji': (296, 393)},
'Port Renfrew': {'lon lat': (-124.421, 48.555),
'stn number': 8525,
'mean sea lvl': 1.937,
'hist max sea lvl': 4.359,
'NEMO grid ji': (401, 61),
'wind grid ji': (117, 112),
'ww3 grid ji': (123, 226)},
'Sandy Cove': {'lon lat': (-123.23, 49.34),
'stn number': 7786,
'NEMO grid ji': (468, 333),
'wind grid ji': (146, 155),
'ww3 grid ji': (294, 396),
'mean sea lvl': 3.1,
'hist max sea lvl': 5.61},
'Squamish': {'lon lat': (-123.155, 49.694),
'stn number': 7811,
'NEMO grid ji': (532, 389),
'wind grid ji': (162, 160),
'ww3 grid ji': (370, 404),
'mean sea lvl': 3.14,
'hist max sea lvl': 5.61},
'Victoria': {'lon lat': (-123.3707, 48.424666),
'stn number': 7120,
'mean sea lvl': 1.881,
'hist max sea lvl': 3.76,
'wind grid ji': (104, 144),
'NEMO grid ji': (302, 196),
'ww3 grid ji': (90, 374)},
'Woodwards Landing': {'lon lat': (-123.0754, 49.1251),
'stn number': 7610,
'hist max sea lvl': 4.66,
'mean sea lvl': 1.84,
'NEMO grid ji': (414, 329),
'wind grid ji': (135, 138)},
'Boundary Bay': {'lon lat': (-122.925, 49.0),
'stn number': None,
'hist max sea lvl': 5.61,
'mean sea lvl': 3.09,
'NEMO grid ji': (380, 335),
'wind grid ji': (129, 162),
'ww3 grid ji': (222, 439)},
'Calamity Point': {'lon lat': (-123.1276, 49.31262),
'stn number': 7724,
'mean sea lvl': 3.001,
'NEMO grid ji': None,
'wind grid ji': (456, 344),
'ww3 grid ji': None},
'Vancouver Harbour': {'lon lat': (-123.1069, 49.28937),
'stn number': 7735,
'mean sea lvl': 3.001,
'NEMO grid ji': None,
'wind grid ji': (143, 159),
'ww3 grid ji': None},
'Port Moody': {'lon lat': (-122.8658, 49.28814),
'stn number': 7755,
'mean sea lvl': 3.143,
'NEMO grid ji': None,
'wind grid ji': (142, 166),
'ww3 grid ji': None},
'Indian Arm Head': {'lon lat': (-122.8864, 49.4615),
'stn number': 7774,
'mean sea lvl': 3.052,
'NEMO grid ji': None,
'wind grid ji': (150, 167),
'ww3 grid ji': None},
'2nd Narrows Rail Bridge': {'lon lat': (-123.0247222, 49.2938889),
'stn number': 3160171,
'mean sea lvl': None,
'NEMO grid ji': None,
'wind grid ji': (143, 161),
'ww3 grid ji': None},
'Tsawwassen': {'lon lat': (-123.132722, 49.006165),
'stn number': None,
'mean sea lvl': None,
'NEMO grid ji': (396, 305),
'wind grid ji': (130, 155),
'ww3 grid ji': None,
'in berth radius': 0.0015},
'Duke Pt.': {'lon lat': (-123.89095676900132, 49.16340592936349),
'stn number': None,
'mean sea lvl': None,
'NEMO grid ji': (481, 213),
'wind grid ji': (142, 134),
'ww3 grid ji': None,
'in berth radius': 0.002},
'Horseshoe Bay': {'lon lat': (-123.2728, 49.3742),
'stn number': None,
'mean sea lvl': None,
'NEMO grid ji': (478, 331),
'wind grid ji': (148, 154),
'ww3 grid ji': None},
'Departure Bay': {'lon lat': (-123.8909, 49.1632),
'stn number': None,
'mean sea lvl': None,
'NEMO grid ji': (481, 213),
'wind grid ji': (142, 134),
'ww3 grid ji': None},
'Swartz Bay': {'lon lat': (-123.4102, 48.6882),
'stn number': None,
'mean sea lvl': None,
'NEMO grid ji': (354, 225),
'wind grid ji': (117, 144),
'ww3 grid ji': None},
'Vancouver': {'lon lat': (-123.1207, 49.2827)},
'Hope': {'lon lat': (-121.4419, 49.3858)},
'British Columbia': {'lon lat': (-123.6, 49.9)},
'Washington State': {'lon lat': (-123.8, 47.8)},
'Pacific Ocean': {'lon lat': (-125.6, 48.1)},
'Juan de Fuca Strait': {'lon lat': (-124.7, 48.47)},
'Puget Sound': {'lon lat': (-122.67, 48)},
'Strait of Georgia': {'lon lat': (-123.8, 49.3)},
'Central SJDF': {'lon lat': (-123.9534, 48.281677),
'NEMO grid ji': (315, 95),
'GEM2.5 grid ji': (101, 124)},
'Baynes Sound': {'lon lat': (-124.86022, 49.60356),
'NEMO grid ji': (635, 126)},
'S3': {'lon lat': (-123.558, 49.125),
'NEMO grid ji': (450, 258),
'GEM2.5 grid ji': (138, 144)},
'Cluster_1': {'NEMO grid ji': (241, 212),
'lon lat': (48.215, -123.099),
'Vector Stn': '64'},
'Cluster_2': {'NEMO grid ji': (294, 127),
'lon lat': (48.261, -123.717),
'Vector Stn': '69'},
'Cluster_3': {'NEMO grid ji': (376, 291),
'lon lat': (48.899, -123.138),
'Vector Stn': '45'},
'Cluster_4': {'NEMO grid ji': (282, 305),
'lon lat': (48.555, -122.75),
'Vector Stn': '53'},
'Cluster_5': {'NEMO grid ji': (344, 271),
'lon lat': (48.735, -123.135),
'Vector Stn': '57'},
'Cluster_6': {'NEMO grid ji': (320, 68),
'lon lat': (48.249, -124.11),
'Vector Stn': '73'},
'Cluster_7': {'NEMO grid ji': (504, 246),
'lon lat': (49.317, -123.801),
'Vector Stn': '27'},
'Cluster_8': {'NEMO grid ji': (646, 168),
'lon lat': (49.726, -124.679),
'Vector Stn': '12'},
'Cluster_9': {'NEMO grid ji': (423, 300), 'lon lat': (49.101, -123.249)},
'Central node': {'lon lat': (-123.425825, 49.040066666),
'depth': 294,
'NEMO grid ji': (424, 266),
'wind grid ji': (133, 147),
'ONC stationCode': 'SCVIP'},
'Delta BBL node': {'lon lat': (-123.339633, 49.074766),
'depth': 143,
'NEMO grid ji': (424, 283),
'wind grid ji': (134, 150),
'ONC stationCode': 'LSBBL'},
'Delta DDL node': {'lon lat': (-123.32972, 49.08495),
'depth': 107,
'NEMO grid ji': (426, 286),
'wind grid ji': (135, 150),
'ONC stationCode': 'USDDL'},
'East node': {'lon lat': (-123.316836666, 49.04316),
'depth': 164,
'NEMO grid ji': (417, 283),
'wind grid ji': (133, 150),
'ONC stationCode': 'SEVIP'},
'Ballenas Islands': {'lon lat': (-124.16, 49.35),
'NEMO grid ji': (536, 197),
'wind grid ji': (152, 127)},
'Discovery Island': {'lon lat': (-123.226, 48.425),
'NEMO grid ji': (291, 219),
'wind grid ji': (104, 148)},
'Entrance Island': {'lon lat': (-123.811, 49.209),
'NEMO grid ji': (484, 231),
'wind grid ji': (143, 137)},
'Race Rocks': {'lon lat': (-123.531, 48.298),
'NEMO grid ji': (288, 159),
'wind grid ji': (99, 137)},
'Sand Heads': {'lon lat': (-123.3, 49.1),
'stn number': 7594,
'NEMO grid ji': (426, 293),
'mean sea lvl': 2.875,
'hist max sea lvl': 5.3950000000000005,
'GEM2.5 grid ji': (135, 151),
'wind grid ji': (135, 151),
'ww3 grid ji': (246, 385)},
'Saturna Island': {'lon lat': (-123.045, 48.784),
'NEMO grid ji': (347, 290),
'wind grid ji': (119, 156)},
'Sisters Islet': {'lon lat': (-124.43, 49.49),
'NEMO grid ji': (582, 175),
'GEM2.5 grid ji': (160, 120),
'wind grid ji': (160, 120)},
'Esquimalt': {'lon lat': (-123.439, 48.432),
'NEMO grid ji': (307, 189),
'wind grid ji': (105, 141)},
'Pam Rocks': {'lon lat': (-123.299, 49.488),
'NEMO grid ji': (502, 341),
'wind grid ji': (153, 154)},
'Halibut Bank': {'lon lat': (-123.72, 49.34),
'NEMO grid ji': (503, 261),
'GEM2.5 grid ji': (149, 141),
'wind grid ji': (149, 141),
'EC buoy number': 46146},
'Sentry Shoal': {'lon lat': (-125.0, 49.92),
'NEMO grid ji': (707, 145),
'GEM2.5 grid ji': (183, 107),
'wind grid ji': (183, 107),
'EC buoy number': 46131},
'Egmont': {'lon lat': (-123.93, 49.75), 'NEMO grid ji': (598, 282)},
'Comox Airport': {'lon lat': (-124.9, 49.717),
'NEMO grid ji': (660, 134),
'wind grid ji': (173, 108)},
'Squamish Airport': {'lon lat': (-123.161, 49.783),
'wind grid ji': (166, 161)},
'YVR': {'lon lat': (-123.184, 49.195),
'GEM2.5 grid ji': (139, 155),
'wind grid ji': (139, 155)},
'Hoodsport': {'lon lat': (-123.1126, 47.4218), 'NEMO grid ji': (89, 114)},
'Twanoh': {'lon lat': (-123.0083, 47.375), 'NEMO grid ji': (72, 123)},
'DabobBay': {'lon lat': (-122.8029, 47.8031), 'NEMO grid ji': (141, 205)},
'PointWells': {'lon lat': (-122.3972, 47.7612), 'NEMO grid ji': (104, 259)},
'CarrInlet': {'lon lat': (-122.73, 47.28), 'NEMO grid ji': (34, 152)},
'Hansville': {'lon lat': (-122.627, 47.90733), 'NEMO grid ji': (148, 243)},
'Dockton': {'lon lat': (-122.45722, 47.37611), 'NEMO grid ji': (34, 204)},
'PointWilliams': {'lon lat': (-122.40612, 47.53716),
'NEMO grid ji': (62, 231)},
'Sandheads': {'lon lat': (-123.3, 49.1),
'stn number': 7594,
'NEMO grid ji': (426, 293),
'mean sea lvl': 2.875,
'hist max sea lvl': 5.3950000000000005,
'GEM2.5 grid ji': (135, 151),
'wind grid ji': (135, 151),
'ww3 grid ji': (246, 385)}}
In [21]:
# NEMO grid
j,i=geo_tools.find_closest_model_point(lonQU39,latQU39,mesh.variables['nav_lon'][:,:],mesh.variables['nav_lat'][:,:])
print(j,i)
736 144
In [22]:
# ops grid
ops=nc.Dataset('/results/forcing/atmospheric/GEM2.5/operational/ops_y2015m08d16.nc')
ops.variables.keys()
Out[22]:
dict_keys(['atmpres', 'nav_lat', 'nav_lon', 'precip', 'qair', 'solar', 'tair', 'therm_rad', 'time_counter', 'u_wind', 'v_wind', 'x', 'y'])
In [30]:
jo,io=geo_tools.find_closest_model_point(lonQU39,latQU39,ops.variables['nav_lon'][:,:]-360,ops.variables['nav_lat'][:,:],grid='GEM2.5',)
print(jo,io)
189 106
In [ ]: