{
"cells": [
{
"cell_type": "code",
"execution_count": 100,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"import numpy as np\n",
"import netCDF4 as nc\n",
"import matplotlib.pyplot as plt\n",
"import scipy.io as sio\n",
"import pandas as pd\n",
"import datetime\n",
"import xarray as xr\n",
"from salishsea_tools import tidetools, geo_tools, viz_tools\n",
"from matplotlib.colors import LinearSegmentedColormap\n",
"import os\n",
"import h5py\n",
"import glob\n",
"%matplotlib inline\n",
"\n",
"from IPython.core.display import display, HTML\n",
"display(HTML(\"\"))"
]
},
{
"cell_type": "code",
"execution_count": 41,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
""
],
"text/plain": [
""
]
},
"execution_count": 41,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"from IPython.display import HTML\n",
"\n",
"HTML('''\n",
"\n",
"''')\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# In one SalishSeaCast grid"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [],
"source": [
"f0 = nc.Dataset('/ocean/vdo/MIDOSS/Lagrangian_AKNS_crude-0_InOneSSGrid-0.nc')\n",
"f1 = nc.Dataset('/ocean/vdo/MIDOSS/Lagrangian_AKNS_crude-1_InOneSSGrid-1.nc')\n",
"f2 = nc.Dataset('/ocean/vdo/MIDOSS/Lagrangian_AKNS_crude-2_InOneSSGrid-2.nc')\n",
"f3 = nc.Dataset('/ocean/vdo/MIDOSS/Lagrangian_AKNS_crude-3_InOneSSGrid-3.nc')\n",
"f4 = nc.Dataset('/ocean/vdo/MIDOSS/Lagrangian_AKNS_crude-4_InOneSSGrid-4.nc')"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA3YAAAEzCAYAAAB9p5F7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAY6UlEQVR4nO3df4ykd30f8PeHu4t9yCFHwYjDdmNaB1Cg+EwO12l+CBxLGFpwJaAhUpOAYrlBcZMSEVorUhIjVQpNWyoMoTICaihRbbkQOY5NZUTdgBTOMfTuagtcXRsIBlMTwCYnzEk23/6xj2F97O3O7u3s7md4vaQRM8+P2c+M3178nmefZ2qMEQAAAPp60nYPAAAAwOlR7AAAAJpT7AAAAJpT7AAAAJpT7AAAAJpT7AAAAJpbs9hV1ZlVdVdVHamqe6vq2mn5f6qqv6yqw9PtwLS8quodVXWsqo5W1Yvm/SJgOZmlG5mlI7mlG5ll0e2eYZsTSS4dYxyvqj1JPllVt0/rfmuMcfNJ2788yY9Nt7+f5N3T/8JWkVm6kVk6klu6kVkW2ppH7MaS49PDPdNttW81vyLJB6b9PpVkX1XtP/1RYTYySzcyS0dySzcyy6Kb6Ry7qtpVVYeTPJjkjjHGoWnVv54OTb+9qs6Ylp2T5IvLdr9/WgZbRmbpRmbpSG7pRmZZZLP8KWbGGI8lOVBV+5J8pKpekOSaJF9J8kNJrk/yL5O8NUmt9BQnL6iqq5JctfRg90/UmU/dyPyQJBmPfPWvxxhnf/exzLLDbUVmE7ll85yc2cTvWnY2maWblTK7HjMVu+/+sDEeqqo7k1w+xvi30+ITVfX+JG+eHt+f5Lxlu52b5MsrPNf1WfqXJ0968jPGGc/9J+scHb7n24ff9YWVlsssO9VWZHZ6PrllU5wqs4nftexMMks3q2V2FrNcFfPs6VONVNXeJJcl+dzjf2NcVZXkHye5Z9rlliS/NF1J6JIkD48xHjidIWE9ZJZuZJaO5JZuZJZFN8sRu/1JbqiqXVkqgjeNMW6tqo9X1dlZOkx9OMmvTtvfluQVSY4l+VaSN2z+2LAqmaUbmaUjuaUbmWWhrVnsxhhHk1y0wvJLT7H9SPJrpz8abIzM0o3M0pHc0o3MsuhmuiomAAAAO5diBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0JxiBwAA0Nzu7R4AaOpp55563dfu37o5AIAt842/eOcp19335b9Zdd9Lrrhms8dhGUfsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmvMF5cCGfOOj/2pD+z31xVdv8iQAwE7w3Gf98HaP8APNETsAAIDmFDsAAIDmFDsAAIDm1ix2VXVmVd1VVUeq6t6quvak9ddV1fFlj19fVV+tqsPT7cp5DA6nIrN0I7N0JLd0I7MsulkunnIiyaVjjONVtSfJJ6vq9jHGp6rqYJJ9K+xz4xjDFRLYLjJLNzJLR3JLNzLLQlvziN1Y8vinF3um26iqXUn+IMlb5jgfrJvM0o3M0pHc0o3MsuhmOseuqnZV1eEkDya5Y4xxKMnVSW4ZYzywwi6vrqqjVXVzVZ23ifPCTGSWbmSWjuSWbmSWRTZTsRtjPDbGOJDk3CQXV9XPJnltkutW2PxPkpw/xnhhko8luWGl56yqq6rq7qq6ezz6yMamh1OQWbqZR2YTuWW+/K6lG5llka3rC8rHGA9V1Z1JXprkgiTHqipJnlxVx8YYF4wxvrZsl/ckedspnuv6JNcnyZOe/Iyxgdl/YHzjL955ynW+7Hl1Mks3m5nZ6fnklrnzu3Y+/P///MjsxsnezjXLVTHPrqp90/29SS5L8ukxxjPHGOePMc5P8q0xxgXTNvuX7f6qJJ/d/LHh1GSWbmSWjuSWbmSWRTfLEbv9SW6YTix9UpKbxhi3rrL9r1fVq5I8muTrSV5/2lPC+sgs3cgsHckt3cgsC23NYjfGOJrkojW2OWvZ/WuSXHP6o8HGyCzdyCwdyS3dyCyLbqaLpwAAALBzKXYAAADNKXYAAADNKXYAAADNKXYAAADNresLytkevgiSnUgu6caXPQOwyByxAwAAaE6xAwAAaE6xAwAAaE6xAwAAaE6xAwAAaE6xAwAAaE6xAwAAaE6xAwAAaM4XlAPwA+FbJx7b7hFgXZ764qu3ewSgEUfsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmlPsAAAAmvMF5QD8QDjnp39ju0cAgLlxxA4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKA5xQ4AAKC5NYtdVZ1ZVXdV1ZGqureqrj1p/XVVdXzZ4zOq6saqOlZVh6rq/M0fG05NZulIbulGZulGZll0sxyxO5Hk0jHGhUkOJLm8qi5Jkqo6mGTfSdv/SpJvjDEuSPL2JG/bxHlhFjJLR3JLNzJLNzLLQluz2I0lj396sWe6jaraleQPkrzlpF2uSHLDdP/mJD9XVbVJ88KaZJaO5JZuZJZuZJZFN9M5dlW1q6oOJ3kwyR1jjENJrk5yyxjjgZM2PyfJF5NkjPFokoeTPG3zRoa1ySwdyS3dyCzdyCyLbPcsG40xHktyoKr2JflIVf1sktcmeckKm6/0Scb4vo2qrkpyVZJkz1kzjguzkVk6klu6kVm6kVkW2bquijnGeCjJnUlemuSCJMeq6vNJnlxVx6bN7k9yXpJU1e4kP5Lk6ys81/VjjINjjIO1e++GXwCsRmbpSG7pRmbpRmZZRLNcFfPs6VONVNXeJJcl+fQY45ljjPPHGOcn+dZ0YmmS3JLkl6f7r0ny8THG9326AfMis3Qkt3Qjs3Qjsyy6Wf4Uc3+SG6YTS5+U5KYxxq2rbP/eJB+cPu34epLXnf6YsC4yS0dySzcySzcyy0Jbs9iNMY4muWiNbc5adv/bWfpbZdgWMktHcks3Mks3MsuiW9c5dgAAAOw8ih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzih0AAEBzaxa7qjqzqu6qqiNVdW9VXTstf++07GhV3VxVZ03LX19VX62qw9Ptynm/CFhOZulGZulIbulGZll0u2fY5kSSS8cYx6tqT5JPVtXtSd40xvhmklTVv09ydZLfn/a5cYxx9VwmhrXJLN3ILB3JLd3ILAttzWI3xhhJjk8P90y3sexfgEqyN8mY15CwHjJLNzJLR3JLNzLLopvpHLuq2lVVh5M8mOSOMcahafn7k3wlyfOSXLdsl1cvO5x93mYPDWuRWbqRWTqSW7qRWRbZTMVujPHYGONAknOTXFxVL5iWvyHJs5J8NsnPT5v/SZLzxxgvTPKxJDes9JxVdVVV3V1Vd49HHznNlwFPJLN0M4/MJnLLfPldSzcyyyJb11UxxxgPJbkzyeXLlj2W5MYkr54ef22McWJa/Z4kP3GK57p+jHFwjHGwdu/dwOiwNpmlm83M7LSt3DJ3ftfSjcyyiGa5KubZVbVvur83yWVJ7quqC6ZlleSVST43Pd6/bPdXZemTD9gyMks3MktHcks3Msuim+WqmPuT3FBVu7JUBG9K8qdJPlFVT0lSSY4keeO0/a9X1auSPJrk60lev9lDwxpklm5klo7klm5kloU2y1Uxjya5aIVVP3WK7a9Jcs1pzgUbJrN0I7N0JLd0I7MsunWdYwcAAMDOo9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0p9gBAAA0t2axq6ozq+quqjpSVfdW1bXT8vdOy45W1c1Vdda0/IyqurGqjlXVoao6f74vAZ5IZulIbulGZulGZll0sxyxO5Hk0jHGhUkOJLm8qi5J8qYxxoVjjBcm+askV0/b/0qSb4wxLkjy9iRvm8PcsBqZpSO5pRuZpRuZZaGtWezGkuPTwz3TbYwxvpkkVVVJ9iYZ0zZXJLlhun9zkp+btoEtIbN0JLd0I7N0I7MsupnOsauqXVV1OMmDSe4YYxyalr8/yVeSPC/JddPm5yT5YpKMMR5N8nCSp23y3LAqmaUjuaUbmaUbmWWRzVTsxhiPjTEOJDk3ycVV9YJp+RuSPCvJZ5P8/LT5Sp9kjJMXVNVVVXV3Vd09Hn1kQ8PDqcgsHckt3cgs3cgsi2xdV8UcYzyU5M4kly9b9liSG5O8elp0f5LzkqSqdif5kSRfX+G5rh9jHBxjHKzdezc0PKxFZulIbulGZulGZllEs1wV8+yq2jfd35vksiT3VdUF07JK8sokn5t2uSXJL0/3X5Pk42OM7/t0A+ZFZulIbulGZulGZll0u2fYZn+SG6pqV5aK4E1J/jTJJ6rqKVk6TH0kyRun7d+b5INVdSxLn2q8btOnhtXJLB3JLd3ILN3ILAttzWI3xjia5KIVVv3UKbb/dpLXnuZcsGEyS0dySzcySzcyy6Jb1zl2AAAA7DyKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHOKHQAAQHNrFruqOrOq7qqqI1V1b1VdOy3/UFXdV1X3VNX7qmrPtPwlVfVwVR2ebr8z7xcBy8ks3cgs3cgsHckti273DNucSHLpGOP4FPRPVtXtST6U5J9O2/xRkiuTvHt6/Ikxxj/a9GlhNjJLNzJLNzJLR3LLQluz2I0xRpLj08M9022MMW57fJuquivJuXOZENZJZulGZulGZulIbll0M51jV1W7qupwkgeT3DHGOLRs3Z4kv5jko8t2+cnpMPftVfX8TZ0YZiCzdCOzdCOzdCS3LLKZit0Y47ExxoEsfYJxcVW9YNnqP0zyZ2OMT0yPP5PkR8cYFya5Lskfr/ScVXVVVd1dVXePRx/Z+CuAFcgs3cwjs4ncMj8yS0f++4BFtq6rYo4xHkpyZ5LLk6SqfjfJ2Ul+c9k23xxjHJ/u35ZkT1U9fYXnun6McXCMcbB27934K4BVyCzdbGZmp/Vyy1zJLB357wMW0SxXxTy7qvZN9/cmuSzJ56rqyiQvS/ILY4zvLNv+mVVV0/2Lp5/xtXkMDyuRWbqRWbqRWTqSWxbdLFfF3J/khqralaVA3zTGuLWqHk3yhSR/PmX+w2OMtyZ5TZI3TusfSfK66WRV2CoySzcySzcyS0dyy0Kb5aqYR5NctMLyFfcdY7wzyTtPfzTYGJmlG5mlG5mlI7ll0a3rHDsAAAB2HsUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACgOcUOAACguTWLXVWdWVV3VdWRqrq3qq6dln+oqu6rqnuq6n1VtWdaXlX1jqo6VlVHq+pF834RsJzM0o3M0pHc0o3MsuhmOWJ3IsmlY4wLkxxIcnlVXZLkQ0mel+TvJdmb5Mpp+5cn+bHpdlWSd2/20LAGmaUbmaUjuaUbmWWhrVnsxpLj08M9022MMW6b1o0kdyU5d9rmiiQfmFZ9Ksm+qto/j+FhJTJLNzJLR3JLNzLLopvpHLuq2lVVh5M8mOSOMcahZev2JPnFJB+dFp2T5IvLdr9/WgZbRmbpRmbpSG7pRmZZZLtn2WiM8ViSA1W1L8lHquoFY4x7ptV/mOTPxhifmB7XSk9x8oKquipLh7WT5MS3D7/rnpO32UZPT/LX2z3EMuZZ23OXP5DZHWGnzbTT5pl7ZhO5XSfzrO65Jy/wu3bbmWd1Mrvz/pmYZ3Xfl9n1mKnYPW6M8VBV3Znk8iT3VNXvJjk7yT9bttn9Sc5b9vjcJF9e4bmuT3J9klTV3WOMg+sbfX7Ms7qdNk+yNNNKy2V2++y0mXbiPCst38zMTs8ntzMyz+pOldnE79rtYp7Vyax51rIT5zmd/We5KubZ06caqaq9SS5L8rmqujLJy5L8whjjO8t2uSXJL01XErokycNjjAdOZ0hYD5mlG5mlI7mlG5ll0c1yxG5/khuqaleWiuBNY4xbq+rRJF9I8udVlSQfHmO8NcltSV6R5FiSbyV5w1wmh1OTWbqRWTqSW7qRWRbamsVujHE0yUUrLF9x3+mKQr+2zjmuX+f282ae1e20eZJlM8nsjrHTZtqx82xRZp/wM3cI86xuR8/jd+2OYJ7Vyax51rJQ89RSZgEAAOhqpq87AAAAYOfakmJXVedV1X+vqs9W1b1V9RsnrX9zVY2qevr0uKrqHVV1rKqOVtWLtnmel1TVw1V1eLr9zlbMU1W/V1VfWvZzX7Fsn2um9+e+qnrZds5TVedX1SPLlv/HrZhnWvfPp/fg3qr6N8uWn9b7I7Mbm0dmV59nWiezK88js6vMI7MyK7OzzTOtk9mV55HZVeZZiMyOMeZ+y9LJqi+a7v9wkv+d5Menx+cl+W9ZOmn16dOyVyS5PUvfH3JJkkPbPM9Lkty61e9Pkt9L8uYVtv/xJEeSnJHk2Un+T5Jd2zjP+Unu2Yb356VJPpbkjGndMzbr/ZFZmZVZmZVZmZVZmZVZme2U2S05YjfGeGCM8Znp/t8k+WySc6bVb0/yljzxCx+vSPKBseRTSfZV1f5tnGeu1phnJVck+S9jjBNjjL/M0tWaLt7GeeZqlXnemOT3xxgnpnUPTruc9vsjs6c1z0pkVmZl9vTmmSuZldk5zDNXMiuzc5hnrrYis1t+jl1VnZ+lKxIdqqpXJfnSGOPISZudk+SLyx7fnzn9g5hxniT5yao6UlW3V9Xz5zHLyfNMi66eDt+/r6qeOi3blvdnlXmS5NlV9T+r6n9U1c/MY5YV5nlOkp+pqkPTz33xtNmmvj8yO/s80yKZPfU8MiuzG50nkVmZldlZ5pFZmd3oPEnzzG5psauqs5L81yT/IsmjSX47yUp/31srLNv0TxzWMc9nkvzoGOPCJNcl+ePNnuXkecYY30zy7iR/N8mBJA8k+XePb7rC7nN9f9aY54Ekf3uMcVGS30zyR1X1lC2YZ3eSp2bpzxt+K8lNVVXZxPdHZmefR2ZnmkdmZXaj88iszMrsbPPIrMxudJ72md2yYldVe7L0Ij40xvhwlt7QZyc5UlWfT3Juks9U1TOz1EjPW7b7uUm+vF3zjDG+OcY4niRjjNuS7KnpRNQ5zpMxxv8bYzw2xvhOkvfke4dft+P9OeU80yHir033P52lvwF+zrznydL78OHpTxzuSvKdJE/PJr0/MrvueWR2jXkiszK7wXlkVmZldrZ5IrMyu8F5FiKzY04nCC6/ZalxfiDJf1hlm8/neyd3/sM88WTTu7Z5nmcm3/3Ov4uT/NXjj+c5T5L9y+6/KUt/Z5skz88TT6b8v9nck03XO8/Zj//8JH8nyZeS/K0tmOdXk7x1uv+cLB2urs14f2RWZmVWZmVWZmVWZmVWZjtldtPCtcYL+eksHTo8muTwdHvFKsGrJO/KUlP+X0kObvM8Vye5d3pzP5XkH2zFPEk+OL3+o0luOSmIvz29P/clefl2zpPk1cven88keeUWzfNDSf5zknumn3vpZr0/MiuzMiuzMiuzMiuzMiuznTL7eGsHAACgqS2/KiYAAACbS7EDAABoTrEDAABoTrEDAABoTrEDAABoTrEDAABoTrEDAABoTrEDAABo7v8DtYcLUZ8aU5sAAAAASUVORK5CYII=\n",
"text/plain": [
"