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
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
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);