The vel2fits.py
===============

Convert PPV output to FITS format. Edit the first lines:

.. code-block:: python

   data = np.loadtxt('RT_vel_SILCC64_+z_CO10.dat')  # Input file
   ipix = 64                                        # Resolution
   csize = 62.5                                     # Cloud size (pc)
   outputname = 'SILCC64_vel_CO10.fits'            # Output name
   channels = 201                                   # Number of channels

Run with:

.. code-block:: bash

   python vel2fits.py

Example
~~~~~~~

The figure below shows an example FITS file generated with ``vel2fits.py`` and viewed with ``ds9``.

.. figure:: ../images/ds9example.png
   :width: 98%

See :ref:`RTgallery` for velocity integrated emission maps produced with **RT-synth**