In [1]:

    

%matplotlib notebook


    

In [2]:

    

import pylab as pl


    

In [3]:

    

import pyspeckit

twotwo = pyspeckit.Cube('gcnh3_11.fits')
oneone = pyspeckit.Cube('gcnh3_22.fits')
oneone.xarr.velocity_convention = 'radio'
twotwo.xarr.velocity_convention = 'radio'
stack = pyspeckit.CubeStack([oneone,twotwo])

stack.plot_special = pyspeckit.wrappers.fitnh3.plotter_override
stack.plot_special_kwargs = {'fignum':3, 'vrange':[-30,135]}


    

    




WARNING: AstropyDeprecationWarning: astropy.extern.six will be removed in 4.0, use the six module directly if it is still needed [astropy.extern.six]
Invalid 'BLANK' keyword in header.  The 'BLANK' keyword is only applicable to integer data, and will be ignored in this HDU.
Invalid 'BLANK' keyword in header.  The 'BLANK' keyword is only applicable to integer data, and will be ignored in this HDU.






    




INFO: Creating Cube Stack [pyspeckit.cubes.SpectralCube]
INFO: Concatenating data [pyspeckit.cubes.SpectralCube]

In [4]:

    

stack.mapplot(estimator=np.nanmax)
pl.figure(3)
stack.plot_spectrum(5,5)

pl.show()