The vel2fits.py

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

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:

python vel2fits.py

Example

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

See RTgallery for velocity integrated emission maps produced with RT-synth