from marvin.tools import Maps
maps = Maps('8485-1901')

# must use Maps.getSpaxel() to get cube and modelcube
spax = maps.getSpaxel(x=17, y=17, xyorig='lower', cube=True, modelcube=True)

# mask out pixels lacking model fit
no_fit = ~spax.full_fit.masked.mask

# extra arguments to plot are passed to the matplotlib routine
ax = spax.flux.plot(label='observed')
ax.plot(spax.full_fit.wavelength[no_fit], spax.full_fit.value[no_fit], label='model')
ax.legend()