Merge pull request #2620 from mitchcs/casda-docs-cutout-channel-example

Casda docs cutout channel example

Author: bsipocz

docs/casda/casda.rst

@@ -158,7 +158,7 @@ below:
     >>> url_list = casda.cutout(subset[:1], coordinates=centre, radius=14*u.arcmin)
     >>> filelist = casda.download_files(url_list, savedir='/tmp')
 
-An example script to download a 3D cutout from the WALLABY Pre-Pilot Eridanus cube at a specified position and velocity
+An example script to download a 3D cutout from the WALLABY Pre-Pilot Eridanus cube at a specified position and velocity range
 is shown below:
 
 .. doctest-skip::
@@ -176,6 +176,23 @@ is shown below:
     >>> url_list = casda.cutout(eridanus_cube, coordinates=centre, radius=9*u.arcmin, band=freq)
     >>> filelist = casda.download_files(url_list, savedir='/tmp')
 
+An example script to download a 3D cutout of a spectral channel range from the WALLABY Pre-Pilot Eridanus cube 
+is shown below:
+
+.. doctest-skip::
+
+    >>> from astropy import coordinates, units as u, wcs
+    >>> from astroquery.casda import Casda
+    >>> centre = coordinates.SkyCoord.from_name('NGC 1371')
+    >>> casda = Casda()
+    >>> casda.login(username='[email protected]')
+    >>> result = Casda.query_region(centre, radius=30*u.arcmin)
+    >>> public_data = Casda.filter_out_unreleased(result)
+    >>> eridanus_cube = public_data[public_data['filename'] == 'Eridanus_full_image_V3.fits']
+    >>> channel = (5, 10)
+    >>> url_list = casda.cutout(eridanus_cube, channel=channel)
+    >>> filelist = casda.download_files(url_list, savedir='/tmp')
+
 
 Reference/API
 =============