In [1]:
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy
import xarray
import cmocean as cm

Load high resolution COSIMA ocean sea surface temperature and convert to degrees celcius

In [2]:
ds = xarray.open_dataset('https://dapds00.nci.org.au/thredds/dodsC/cj50/access-om2/cf-compliant/access-om2-01/v20191202/jra55v13_iaf/ocean/surface-temp/surface-temp_access-om2-01_198501_198512.nc')
surface_temp = ds.surface_temp -273.15

Plot data using a Nearside Perspective projection using cartopy

In [3]:
plt.figure(figsize=(24,12))

ax = plt.axes(projection=ccrs.NearsidePerspective(central_longitude=145.0, 
                                                  central_latitude=-20.0, 
                                                  satellite_height=8500000))
ax.background_img(name='bluemarble', resolution='extrahigh')

surface_temp.isel(time=8).plot(ax=ax, 
                  transform=ccrs.PlateCarree(), 
                  x='xt_ocean', 
                  y='yt_ocean',
                  vmin=-2, 
                  vmax=30,
                  cmap=cm.cm.thermal);
In [ ]: