hfs_nh3_view

Back to the main HfS page

hfs_nh3_view is an interactive graphic procedure for displaying spectra from a FITS data cube and the corresponding synthetic spectra. The data and synthetic spectra at any pixel can be extracted in ASCII files.

It is a menu-driven program, and uses the PGplot Graphics Subroutine Library to display the spectral data, fit, and residual.
Example of a hfs_nh3_view run:

$ hfs_nh3_view test.par______________________________________________ HfS_NH3_view. View data and synthetic spectra Robert Estalella, 2015/12

______________________________________________                                                          naxis= [128, 128, 124, 2]                                                                              Storage allocated successfully for arrays

Fits file: test1.fits (1) Cols (2) Rows (3) Chan CTYPEi= RA---SIN DEC--SIN VELO-LSR NAXISi= 128 128 124 CRPIXi= 64.000 65.000 64.000 CRVALi= 333.784 58.819 -18300.000 CDELTi= -1.917E-04 1.917E-04 -3.089E+02 CROTAi= 0.000E+00 0.000E+00 0.000E+00 EQUINOX= 2000.0 Sorting velocities in ascending order
Fits file: test2.fits (1) Cols (2) Rows (3) Chan CTYPEi= RA---SIN DEC--SIN VELO-LSR NAXISi= 128 128 120 CRPIXi= 64.000 65.000 61.000 CRVALi= 333.784 58.819 -18300.000 CDELTi= -1.917E-04 1.917E-04 -3.085E+02 CROTAi= 0.000E+00 0.000E+00 0.000E+00 EQUINOX= 2000.0 Sorting velocities in ascending order Reading out file test_comp1.out Reading out file test_comp2.out
Delt_x, Delt_y: ( 3, 3) ( -2.07, +2.07) TRC: (0033, 0043) (+0021.39, -0015.18) BLC: (0030, 0040) (+0023.46, -0017.25)
In NH3(1,1) mom0 map (UPPER PANEL): Click LEFT button (A) to select POSITION Click MIDDLE button (D) to switch COMPONENT Click RIGHT button (X) to EXIT Press P key to PRINT spectrum and fit

A plot of the integrated intensity for the NH₃ (1,1) and (2,2) lines is shown, and you can select with the mouse the position in the (1,1) map for which the the spectra (data and synthetic) are shown.
You can select to show the integrated intensity for all channels, or for the channel ranges of each velocity component.
You can print the spectra (data and synthetic) for the pixel selected.

Input

<parfile>.par, the same parameter file used as input to hfs_nh3_cube. From the information in this file, hfs_nh3_view reads the corresponding FITS data cubes and the <parfile>_comp#.out file for each velocity component (created by hfs_nh3_cube), with the parameters of the fitted spectra.

Output

<source_11>_<xoffset>_<yoffset>.spec, <source_22>_<xoffset>_<yoffset>.spec: data spectra for the pixel selected, with the hanning filtering and boxcar smoothing given in the <parfile>.par file.
<source_11>_<xoffset>_<yoffset>.synt, <source_22>_<xoffset>_<yoffset>.synt: synthetic spectra fitted for the pixel selected.

Report abuse