In [2]:

    

import iris
import iris.plot as iplt
import matplotlib.pyplot as plt


    

    




/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/importlib/_bootstrap.py:219: RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
  return f(*args, **kwds)
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/importlib/_bootstrap.py:219: RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
  return f(*args, **kwds)
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/importlib/_bootstrap.py:219: RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
  return f(*args, **kwds)

In [3]:

    

%matplotlib inline


    

In [4]:

    

my_basin_file = '/g/data/r87/dbi599/DRSv2/CMIP5/CSIRO-Mk3-6-0/historical/fx/ocean/r0i0p0/basin/latest/basin_fx_CSIRO-Mk3-6-0_historical_r0i0p0.nc'
cmip_basin_file = '/g/data/ua6/DRSv2/CMIP5/CSIRO-Mk3-6-0/historical/fx/ocean/r0i0p0/basin/latest/basin_fx_CSIRO-Mk3-6-0_historical_r0i0p0.nc'


    

In [9]:

    

my_basin_cube = iris.load_cube(my_basin_file)
#my_basin_cube


    

In [6]:

    

cmip_basin_cube = iris.load_cube(cmip_basin_file)
#cmip_basin_cube


    

    




/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/site-packages/iris/fileformats/cf.py:798: UserWarning: Missing CF-netCDF measure variable 'areacello', referenced by netCDF variable 'basin'
  warnings.warn(message % (variable_name, nc_var_name))

In [7]:

    

iplt.pcolormesh(cmip_basin_cube)
plt.show()


    

    




/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/site-packages/iris/fileformats/netcdf.py:395: UserWarning: WARNING: missing_value not used since it
cannot be safely cast to variable data type
  var = variable[keys]
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/site-packages/iris/fileformats/netcdf.py:395: UserWarning: WARNING: _FillValue not used since it
cannot be safely cast to variable data type
  var = variable[keys]






    

In [10]:

    

iplt.pcolormesh(my_basin_cube[0, ::])
plt.gca().coastlines()
plt.show()


    

In [ ]: