In [1]:

import numpy as np
import matplotlib.pyplot as plt
import netCDF4 as nc

%matplotlib inline

In [2]:

fbathy=nc.Dataset('/data/eolson/results/MEOPAR/NEMO-forcing-new/grid/bathymetry_201702.nc')

In [3]:

fphy=nc.Dataset('/results/SalishSea/nowcast-green.201812/01jun16/SalishSea_1d_20160601_20160601_ptrc_T.nc')

In [4]:

fbathy.variables.keys()

Out[4]:

dict_keys(['nav_lon', 'nav_lat', 'Bathymetry'])

In [5]:

fphy.variables.keys()

Out[5]:

dict_keys(['nav_lat', 'nav_lon', 'bounds_lon', 'bounds_lat', 'area', 'deptht', 'deptht_bounds', 'nitrate', 'time_centered', 'time_centered_bounds', 'time_counter', 'time_counter_bounds', 'ammonium', 'silicon', 'diatoms', 'flagellates', 'ciliates', 'microzooplankton', 'dissolved_organic_nitrogen', 'particulate_organic_nitrogen', 'biogenic_silicon', 'Fraser_tracer', 'mesozooplankton'])

In [6]:

fbathy.variables['Bathymetry']

Out[6]:

<class 'netCDF4._netCDF4.Variable'>
float64 Bathymetry(y, x)
    _FillValue: 0.0
    least_significant_digit: 0.1
    units: metres
    long_name: sea_floor_depth
unlimited dimensions: 
current shape = (898, 398)
filling on

In [7]:

fig,ax=plt.subplots(1,1,figsize=(5,10))
ax.pcolormesh(fphy.variables['mesozooplankton'][0,0,:,:])
ax.contour(fbathy.variables['Bathymetry'][:,:],[100,],colors='w',linewidths=1)

Out[7]:

<matplotlib.contour.QuadContourSet at 0x7f0280092520>

In [8]:

fig,ax=plt.subplots(1,1,figsize=(5,10))
ax.contourf(fphy.variables['mesozooplankton'][0,0,:,:],20)
ax.contour(fbathy.variables['Bathymetry'][:,:],[100,],colors='w',linewidths=1)

Out[8]:

<matplotlib.contour.QuadContourSet at 0x7f027ff88550>

In [9]:

fbathy.close()
fphy.close()

In [ ]: