{ "cells": [ { "cell_type": "code", "execution_count": 103, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import netCDF4 as nc\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from salishsea_tools import geo_tools, tidetools, viz_tools, loadDataFRP\n", "import matplotlib.cm as cm\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "base=nc.Dataset('/ocean/vdo/MEOPAR/completed-runs/threemonthbase/testD/SalishSea_1h_20170406_20170415_grid_T_20170410-20170410.nc')\n", "testa=nc.Dataset('/ocean/vdo/MEOPAR/completed-runs/threemonthsa/test44/SalishSea_1h_20170406_20170415_grid_T_20170410-20170410.nc')\n", "testb=nc.Dataset('/ocean/vdo/MEOPAR/completed-runs/threemonthsb/testd/SalishSea_1h_20170406_20170415_grid_T_20170410-20170410.nc')" ] }, { "cell_type": "code", "execution_count": 62, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Bathymetry = nc.Dataset('/data/vdo/MEOPAR/NEMO-forcing/grid/bathymetry_201702.nc')\n", "bathy, X, Y = tidetools.get_bathy_data(Bathymetry)" ] }, { "cell_type": "code", "execution_count": 65, "metadata": { "collapsed": true }, "outputs": [], "source": [ "mesh = nc.Dataset('/data/vdo/MEOPAR/NEMO-forcing/grid/mesh_mask201702.nc')\n", "tmask = mesh.variables['tmask'][:]" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "(40, 898, 398) 415 337 1.0\n", "(40, 898, 398) 415 337 2.0\n", "(40, 898, 398) 442 259 3.0\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/vdo/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:3883: FutureWarning: comparison to `None` will result in an elementwise object comparison in the future.\n", " check = self.filled(0).__eq__(other)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "(40, 898, 398) 438 268 4.0\n", "(40, 898, 398) 434 278 5.0\n", "(40, 898, 398) 432 281 6.0\n", "(40, 898, 398) 430 285 7.0\n", "(40, 898, 398) 428 289 8.0\n", "(40, 898, 398) 427 291 9.0\n", "(40, 898, 398) 412 291 10.0\n", "(40, 898, 398) 443 258 11.0\n", "(40, 898, 398) 438 268 12.0\n", "(40, 898, 398) 434 278 13.0\n", "(40, 898, 398) 432 281 14.1\n", "(40, 898, 398) 432 281 14.2\n", "(40, 898, 398) 432 287 15.0\n", "(40, 898, 398) 432 292 16.0\n", "(40, 898, 398) 427 291 17.0\n", "(40, 898, 398) 415 337 18.0\n" ] } ], "source": [ "stationdata, casts = loadDataFRP.loadDataFRP_SSGrid()" ] }, { "cell_type": "code", "execution_count": 27, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "
\n", " | Station | \n", "Date_UTC | \n", "Time_UTC_hhmmss | \n", "
---|---|---|---|
0 | \n", "1.0 | \n", "20170410 | \n", "17:54:17 | \n", "
1 | \n", "2.0 | \n", "20170410 | \n", "18:05:11 | \n", "
2 | \n", "3.0 | \n", "20170410 | \n", "19:44:22 | \n", "
3 | \n", "4.0 | \n", "20170410 | \n", "20:25:40 | \n", "
4 | \n", "5.0 | \n", "20170410 | \n", "21:05:12 | \n", "
5 | \n", "6.0 | \n", "20170410 | \n", "21:40:15 | \n", "
6 | \n", "7.0 | \n", "20170410 | \n", "21:58:48 | \n", "
7 | \n", "8.0 | \n", "20170410 | \n", "22:30:56 | \n", "
8 | \n", "9.0 | \n", "20170410 | \n", "22:45:20 | \n", "
9 | \n", "10.0 | \n", "20170531 | \n", "17:19:23 | \n", "
10 | \n", "11.0 | \n", "20170531 | \n", "18:13:05 | \n", "
11 | \n", "12.0 | \n", "20170531 | \n", "18:51:36 | \n", "
12 | \n", "13.0 | \n", "20170531 | \n", "19:24:38 | \n", "
13 | \n", "14.1 | \n", "20170531 | \n", "19:50:40 | \n", "
14 | \n", "14.2 | \n", "20170531 | \n", "19:53:25 | \n", "
15 | \n", "15.0 | \n", "20170531 | \n", "20:12:26 | \n", "
16 | \n", "16.0 | \n", "20170531 | \n", "20:41:47 | \n", "
17 | \n", "17.0 | \n", "20170531 | \n", "21:01:03 | \n", "
18 | \n", "18.0 | \n", "20170531 | \n", "22:05:25 | \n", "