def plot_model(model):
    cube = cubes[model]
    lon = cube.mesh.nodes[:, 0]
    lat = cube.mesh.nodes[:, 1]
    ind = -1 # just take the last time index for now
    zcube = cube[ind]
    triang = tris[model]

    fig, ax = plt.subplots(figsize=(7, 7),
                           subplot_kw=dict(projection=ccrs.PlateCarree()))
    ax.set_extent([lon.min(), lon.max(), lat.min(), lat.max()])
    ax.coastlines()
    levs = np.arange(-1, 5, 0.2)
    cs = ax.tricontourf(triang, zcube.data, levels=levs)
    fig.colorbar(cs)
    ax.tricontour(triang, zcube.data, colors='k',levels=levs)
    tvar = cube.coord('time')
    tstr = tvar.units.num2date(tvar.points[ind])
    gl = ax.gridlines(draw_labels=True)
    gl.xlabels_top = gl.ylabels_right = False
    title = ax.set_title('%s: Elevation (m): %s' % (zcube.attributes['title'], tstr))
    return fig, ax