Forecast template

Salish Sea NEMO Model Daily Forecast Figures for Sea Surface Height¶

Set-up¶

In [1]:

from __future__ import division

import datetime
from glob import glob
import os

from IPython.core.display import HTML
import netCDF4 as nc

from salishsea_tools.nowcast import figures

%matplotlib inline

In [2]:

def results_dataset(period, grid, results_dir):
    """Return the results dataset for period (e.g. 1h or 1d)
    and grid (e.g. grid_T, grid_U) from results_dir.
    """
    filename_pattern = 'SalishSea_{period}_*_{grid}.nc'
    filepaths = glob(os.path.join(results_dir, filename_pattern.format(period=period, grid=grid)))
    return nc.Dataset(filepaths[0])

In [3]:

run_date = datetime.datetime(2014,12,7)

# Results dataset location
results_home = '/ocean/sallen/allen/research/MEOPAR/SalishSea/forecast/'
results_dir = os.path.join(results_home, run_date.strftime('%d%b%y').lower())
# model winds
model_path='/ocean/sallen/allen/research/MEOPAR/Operational/fcst/'

Load the results datasets:

In [4]:

grid_T_hr = results_dataset('1h', 'grid_T', results_dir)

In [5]:

bathy = nc.Dataset('/data/nsoontie/MEOPAR/NEMO-forcing/grid/bathy_meter_SalishSea2.nc')

Display the figures:

In [6]:

HTML('<h2>{:%d%b%y} Forecast Figures</h2>'.format(run_date))

Out[6]:

07Dec14 Forecast Figures

Point Atkinson water levels¶

In [7]:

fig = figures.PA_tidal_predictions(grid_T_hr)

Tidal predictions calculated with t_tide: http://www.eos.ubc.ca/~rich/#T_Tide

In [8]:

fig = figures.compare_tidalpredictions_maxSSH(grid_T_hr, bathy, model_path, name='Victoria')

In [9]:

fig = figures.compare_tidalpredictions_maxSSH(grid_T_hr, bathy, model_path)

In [10]:

fig = figures.compare_tidalpredictions_maxSSH(grid_T_hr, bathy, model_path, name='Campbell River')

Water levels at NOAA stations¶

Observed water levels and tidal predictions are from NOAA: http://tidesandcurrents.noaa.gov/stations.html?type=Water+Levels

In [11]:

fig = figures.compare_water_levels(grid_T_hr, bathy,PST=1)

In [12]:

fig = figures.plot_thresholds_all(grid_T_hr, bathy, model_path)

Weather¶

Observations from Environment Canada data. http://climate.weather.gc.ca/

Modelled winds are from the High Resolution Deterministic Prediction System of Environment Canada. https://weather.gc.ca/grib/grib2_HRDPS_HR_e.html

In [13]:

fig = figures.Sandheads_winds(grid_T_hr,bathy,model_path)

In [14]:

reload(figures)
fig = figures.average_winds_at_station(grid_T_hr, bathy, model_path, station = 'all')

In [15]:

reload(figures)
fig = figures.winds_at_max_ssh(grid_T_hr, bathy, model_path, station = 'all')

In [15]: