import matplotlib.pyplot as plt
import netCDF4 as nc
import numpy as np
import datetime
from salishsea_tools import viz_tools
import subprocess
%matplotlib inline
from IPython.display import HTML
HTML('''<script>
code_show=true;
function code_toggle() {
if (code_show){
$('div.input').hide();
} else {
$('div.input').show();
}
code_show = !code_show
}
$( document ).ready(code_toggle);
</script>
<form action="javascript:code_toggle()"><input type="submit" value="Click here to toggle on/off the raw code."></form>''')
f = nc.Dataset('/data/vdo/MEOPAR/completed-runs/swiftsure/13apr18tideonly/SalishSea_1h_20180413_20180416_grid_U.nc')
g = nc.Dataset('/data/vdo/MEOPAR/completed-runs/swiftsure/13apr18tideonly/SalishSea_1h_20180413_20180416_grid_V.nc')
mesh = nc.Dataset('/data/vdo/MEOPAR/grid/mesh_mask201702.nc')
fig, ax = plt.subplots(figsize = (8,8))
viz_tools.set_aspect(ax)
ax.set_ylim(200, 500)
ax.set_xlim(0, 250)
ax.pcolormesh(mesh.variables['tmask'][0,0,...]);
grid = nc.Dataset('/data/vdo/MEOPAR/grid/bathymetry_201702.nc')
ufull = f.variables['vozocrtx'][-1,0,200:500:5, 0:250:5]
vfull = g.variables['vomecrty'][-1,0,200:500:5, 0:250:5]
fig, ax = plt.subplots(figsize = (12,12))
ax.pcolormesh(mesh.variables['tmask'][0,0,200:500:5,0:250:5],
cmap = 'Blues', vmin=0, vmax = 3)
Q = ax.quiver(np.ma.masked_array(ufull,
mask = 1 -
mesh.variables['tmask'][0,0,200:500:5, 0:250:5]),
np.ma.masked_array(vfull,
mask = 1 -
mesh.variables['tmask'][0,0,200:500:5,0:250:5]),
scale=17,
width=0.002, headwidth=4
)
Q_key = plt.quiverkey(Q, 0.5, 0.93, 0.5, '0.5 m/s')
ax.set_title(nc.num2date(f.variables['time_counter'][-1],
f.variables['time_counter'].units))
viz_tools.set_aspect(ax);
f = nc.Dataset('/data/vdo/MEOPAR/completed-runs/swiftsure/18may18tideonly/SalishSea_1h_20180518_20180521_grid_U.nc')
g = nc.Dataset('/data/vdo/MEOPAR/completed-runs/swiftsure/18may18tideonly/SalishSea_1h_20180518_20180521_grid_V.nc')
ufull = f.variables['vozocrtx'][-1,0,200:500:5, 0:250:5]
vfull = g.variables['vomecrty'][-1,0,200:500:5, 0:250:5]
fig, ax = plt.subplots(figsize = (12,12))
ax.pcolormesh(mesh.variables['tmask'][0,0,200:500:5,0:250:5],
cmap = 'Blues', vmin=0, vmax = 3)
Q = ax.quiver(np.ma.masked_array(ufull,
mask = 1 -
mesh.variables['tmask'][0,0,200:500:5, 0:250:5]),
np.ma.masked_array(vfull,
mask = 1 -
mesh.variables['tmask'][0,0,200:500:5,0:250:5]),
scale=17,
width=0.002, headwidth=4
)
Q_key = plt.quiverkey(Q, 0.5, 0.93, 0.5, '0.5 m/s')
ax.set_title(nc.num2date(f.variables['time_counter'][-1],
f.variables['time_counter'].units))
viz_tools.set_aspect(ax);
f = nc.Dataset('/data/vdo/MEOPAR/completed-runs/swiftsure/25may18tideonly/SalishSea_1h_20180525_20180528_grid_U.nc')
g = nc.Dataset('/data/vdo/MEOPAR/completed-runs/swiftsure/25may18tideonly/SalishSea_1h_20180525_20180528_grid_V.nc')
ufull = f.variables['vozocrtx'][-1,0,200:500:5, 0:250:5]
vfull = g.variables['vomecrty'][-1,0,200:500:5, 0:250:5]
fig, ax = plt.subplots(figsize = (12,12))
ax.pcolormesh(mesh.variables['tmask'][0,0,200:500:5,0:250:5],
cmap = 'Blues', vmin=0, vmax = 3)
Q = ax.quiver(np.ma.masked_array(ufull,
mask = 1 -
mesh.variables['tmask'][0,0,200:500:5, 0:250:5]),
np.ma.masked_array(vfull,
mask = 1 -
mesh.variables['tmask'][0,0,200:500:5,0:250:5]),
scale=17,
width=0.002, headwidth=4
)
Q_key = plt.quiverkey(Q, 0.5, 0.93, 0.5, '0.5 m/s')
ax.set_title(nc.num2date(f.variables['time_counter'][-1],
f.variables['time_counter'].units))
viz_tools.set_aspect(ax);