We present here the original observations used in the manuscript 
"SOFIA FEEDBACK Survey: The Pillars of Creation in [C II] and Molecular Lines"
(Karim et al., Astronomical Journal, 2023; hereinafter Karim+2023)

We also present a snapshot copy of the scoby (Spectra from Catalogs of OB
Stars) software repository, some model data necessary for it, and some
examples of how to run it.

If you have any questions or comments, please send email to either
Ramsey Karim: rlkarim@terpmail.umd.edu
Marc Pound: mpound@umd.edu

All images/datacubes are in FITS format (Flexible Image Transport System;
https://fits.gsfc.nasa.gov/fits_standard.html).  

A. Combined Array for Research in Millimeter-wave Astronomy (CARMA)
-------------------------------------------------------------------
1. Spectral line datacubes
Spectral resolution ~ 0.08 km/s 

    M16_hcn.fits  - HCN(J=1-0), equatorial coordinates
    M16_hcop.fits - HCO+(J=1-0), equatorial coordinates
    M16_cs.fits   - CS(J=2-1), equatorial coordinates
    M16_n2hp.fits - N2H+(J=1-0, F1=2-1, F=2-1), equatorial coordinates.

2. Integrated intensity maps

The following are integrated intensity (moment 0) maps created from the
datacubes above.  The HCN and N2H+ maps are integrated over all frequency
components, i.e., including satellite lines.  

    M16_hcn_mom0.fits - HCN(J=1-0), equatorial coordinates
    M16_hcop_mom0.fits- HCO+(J=1-0), equatorial coordinates
    M16_cs_mom0.fits - CS(J=2-1), equatorial coordinates
    M16_n2hp_mom0.fits- N2H+(J=1-0), equatorial coordinates


3. Continuum map
This map is made from the continuum windows of the CARMA data centered
around 92 GHz.

    M16_3mm_continuum.fits - 92 GHz continuum, equatorial coordinates


B. Berkeley-Illinois-Maryland Array (BIMA) 
------------------------------------------
These were originally published by Pound 1998
(https://ui.adsabs.harvard.edu/abs/1998ApJ...493L.113P/abstract)
and Pound 2007 
(https://ui.adsabs.harvard.edu/abs/2007Ap%26SS.307..187P/abstract)

1. Spectral line datacubes
   Spectral resolution ~ 0.25 km/s

    M16_co1-0.fits   - 12CO(J=1-0), equatorial coordinates
    M16_13co1-0.fits - 13CO(J=1-0), equatorial coordinates
    M16_c18O1-0.fits - C18O(J=1-0), equatorial coordinates


2. Integrated intensity maps

These are integrated intensity (moment 0) maps created from the
BIMA datacubes above.  

    M16_co1-0_mom0.fits   - 12CO(J=1-0), equatorial coordinates
    M16_12co1-0_mom0.fits - 13CO(J=1-0), equatorial coordinates
    M16_c18o1-0_mom0.fits - C18O(J=1-0), equatorial coordinates


C. Atacama Pathfinder Experiment Telescope (APEX)
-------------------------------------------------
1. Spectral line data
   Spectral resolution - 0.25 km/s

These were observed by Nicola Schneider (University of Cologne) using the
SEPIA660 instrument on the APEX telescope.

    M16_co6-5.fits - 12CO(J=6-5), equatorial coordinates

The APEX LAsMA CO 3-2 datacubes referenced in the paper are not included here but will be released with a future publication.


D. Stratospheric Observatory For Infrared Astronomy (SOFIA)
------------------------------------------------------------
1. Spectral line data 
   Spectral resolution - 15.4 km/s

This is datacube [O I] 63 micron emission, observed by Nicola Schneider
(University of Cologne) with the GREAT instrument.

    M16_OI_63.fits - [O I] 63 micron, equatorial coordinates

The [C II] datacube is publicly archived with the rest of the FEEDBACK
survey data by the NASA/IPAC Infrared Science Archive (IRSA)  at
https://irsa.ipac.caltech.edu/Missions/sofia.html

E. SOFTWARE
------------
1. scoby
This is a zip file containing the August 2023 version of the scoby
repository. scoby is developed and will remain on GitHub: 
https://github.com/ramseykarim/scoby

scoby is a Python software library intended to be imported as a module into
Python code.  Some configuration is needed to run the code; detailed
instructions are found within the code itself and within the scoby-nb
examples, included here and described below.

    scoby-main.zip - scoby software


2. scoby examples
This is a zip file containing the August 2023 version of the scoby-nb
repository. scoby-nb is developed and will remain on GitHub: 
https://github.com/ramseykarim/scoby-nb/

scoby-nb is a collection of Jupyter notebooks that demonstrate how to use the
scoby software. They do not generate the FUV field maps used in Karim+2023,
but can be used as a guide in recreating these maps.

    scoby-nb-main.zip - scoby example notebooks

3. OB star model spectra
This is a zip file containing the OB star model spectrum grid from the PoWR
(Pottsdam Wolf-Rayet Models) Grid of stellar models (Gräfener, Koesterke,
& Hamann 2002, A&A; Hainich et al. 2019, A&A).  These were retrieved from
the PoWR website in 2020. Other models and more information can be found on
the website: https://www.astro.physik.uni-potsdam.de/~wrh/PoWR/powrgrid1.php

    OB.zip - OB star model spectrum grid