In [1]:

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

In [2]:

f=nc.Dataset('/data/eolson/results/MEOPAR/SS36runs/CedarRuns/t15r15DiatAlpha0/SalishSea_1h_20150220_20150410_ptrc_T_20150312-20150321.nc')

In [3]:

mesh=nc.Dataset('/ocean/eolson/MEOPAR/NEMO-forcing/grid/mesh_mask201702_noLPE.nc')

In [4]:

plt.pcolormesh(mesh.variables['tmask'][0,0,:,:])

Out[4]:

<matplotlib.collections.QuadMesh at 0x7f3f51ff0c50>

In [5]:

plt.pcolormesh(np.mean(f.variables['diatoms'][:,:,200:400,100],0))

Out[5]:

<matplotlib.collections.QuadMesh at 0x7f3f5162a588>

In [8]:

plt.pcolormesh(np.mean(f.variables['diatoms'][:,0,:,:],0))
plt.colorbar()

Out[8]:

<matplotlib.colorbar.Colorbar at 0x7f3f425303c8>

In [9]:

plt.pcolormesh(np.mean(f.variables['flagellates'][:,0,:,:],0))
plt.colorbar()

Out[9]:

<matplotlib.colorbar.Colorbar at 0x7f3f419f4198>

In [ ]: