import marvin

my_cube = marvin.tools.Cube('8485-1901')
spectrum = my_cube[15, 15].flux
ax = spectrum.plot(show_std=True)
ax.set_xlim(6000, 8000)