hfs_nh3_blanking

Back to the main HfS page

hfs_nh3_blanking is a procedure to flag the fit values in the output out and FITS files from hfs_nh3_cube, according to any line parameter value, error, or |value|/error.

The procedure reads the values and errors of the line fit parameters from the out files, and you can choose which line parameter to use for blanking the fit. Several passes can be done, using different line fit parameters or blanking criteria.
The output files are the blanked files <parfile>_comp#_blank.out and <parfile>_<fit_par>_comp#_blank.fits, without the fit values flagged out by hfs_blanking.
Example of a hfs_nh3_blanking run:

$ hfs_nh3_blanking i22134.par ________________________________________________________ HfS_NH3_blanking. Blanking fit values from hfs_nh3_cube Robert Estalella, 2017/05 _________________________________________________________ Par file: i22134.par dim= [512, 512, 16, 9] Storage allocated successfully for arrays Reading out file i22134_comp1.out Reading out file i22134_comp2.out
Number of pixels Total: 49 Fitted: 49 Comp 1: 46 Comp 2: 49 HfS NH3 Blanking Menu_________________________________ 0. Quit without saving 1. Save blanked files and exit 2. Delta_V 3. V_lsr1 4. A*tau1_m 5. tau1_m 6. V_lsr2 7. A*tau2_m Choose option or parameter to use for blanking (0-7): 2
Delta_V Value Error |Val|/Err Comp Min Max Min Max Min Max 1 8.68E-02 2.34E+00 1.89E-01 4.95E+00 1.02E-01 1.09E+00 2 9.73E-02 1.53E+00 1.44E-01 8.76E+00 1.15E-01 3.25E+00
Blanking Criterion Menu_______________________________ 0. Back 1. Blanking for Error > cutoff 2. Blanking for Error < cutoff 3. Blanking for |Value|/Error < cutoff 4. Blanking for Value > cutoff 5. Blanking for Value < cutoff Choose option (0-5): 1 Blanking (1). Enter cutoff: 1.0
Blanking pass 1 Blanked Left Pixels: 4 45 Comp 1: 16 30 Comp 2: 10 39
Delta_V Value Error |Val|/Err Comp Min Max Min Max Min Max 1 8.68E-02 5.23E-01 1.89E-01 9.52E-01 1.26E-01 5.49E-01 2 9.73E-02 1.36E+00 1.44E-01 9.98E-01 1.59E-01 3.25E+00
Blanking Criterion Menu_______________________________ 0. Back 1. Blanking for Error > cutoff 2. Blanking for Error < cutoff 3. Blanking for |Value|/Error < cutoff 4. Blanking for Value > cutoff 5. Blanking for Value < cutoff Choose option (0-5): 5 Blanking (5). Enter cutoff: 0.1
Blanking pass 2 Blanked Left Pixels: 2 43 Comp 1: 5 25 Comp 2: 2 37
Delta_V Value Error |Val|/Err Comp Min Max Min Max Min Max 1 1.00E-01 5.23E-01 2.39E-01 9.52E-01 1.26E-01 5.49E-01 2 1.08E-01 1.36E+00 1.44E-01 9.98E-01 1.59E-01 3.25E+00
Blanking Criterion Menu_______________________________ 0. Back 1. Blanking for Error > cutoff 2. Blanking for Error < cutoff 3. Blanking for |Value|/Error < cutoff 4. Blanking for Value > cutoff 5. Blanking for Value < cutoff Choose option (0-5): 0 HfS NH3 Blanking Menu_________________________________ 0. Quit without saving 1. Save blanked files and exit 2. Delta_V 3. V_lsr1 4. A*tau1_m 5. tau1_m 6. V_lsr2 7. A*tau2_m Choose option or parameter to use for blanking (0-7): 1
Writing blanked out files i22134_comp1_blank.out i22134_comp2_blank.out
Writing blanked fits files maps/i22134_dvline_comp1_blank.fits maps/i22134_vlsr1_comp1_blank.fits maps/i22134_atau1m_comp1_blank.fits maps/i22134_tau1m_comp1_blank.fits maps/i22134_vlsr2_comp1_blank.fits maps/i22134_atau2m_comp1_blank.fits maps/i22134_tau2m_comp1_blank.fits maps/i22134_Tex_comp1_blank.fits maps/i22134_Trot_comp1_blank.fits maps/i22134_Tk_comp1_blank.fits maps/i22134_f0N11_comp1_blank.fits maps/i22134_f1N11_comp1_blank.fits maps/i22134_f0N22_comp1_blank.fits maps/i22134_f1N22_comp1_blank.fits maps/i22134_f0Nnh3_comp1_blank.fits maps/i22134_f1Nnh3_comp1_blank.fits maps/i22134_dvline_comp2_blank.fits maps/i22134_vlsr1_comp2_blank.fits maps/i22134_atau1m_comp2_blank.fits maps/i22134_tau1m_comp2_blank.fits maps/i22134_vlsr2_comp2_blank.fits maps/i22134_atau2m_comp2_blank.fits maps/i22134_tau2m_comp2_blank.fits maps/i22134_Tex_comp2_blank.fits maps/i22134_Trot_comp2_blank.fits maps/i22134_Tk_comp2_blank.fits maps/i22134_f0N11_comp2_blank.fits maps/i22134_f1N11_comp2_blank.fits maps/i22134_f0N22_comp2_blank.fits maps/i22134_f1N22_comp2_blank.fits maps/i22134_f0Nnh3_comp2_blank.fits maps/i22134_f1Nnh3_comp2_blank.fits

Input

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

Output

<parfile>_comp#_blank.out: an ASCII file for each velocity component, without the fit values flagged out by hfs_blanking.
FITS files for each component in folder maps, with non-flagged pixels, and NaN for pixels flagged out by hfs_blanking. Each FITS file has two planes, the first plane with the values of the parameter, and the second plane with the uncertainties:
maps/<parfile>_dvline_comp#_blank.fits
maps/<parfile>_vlsr1_comp#_blank.fits
maps/<parfile>_atau1m_comp#_blank.fits
maps/<parfile>_tau1m_comp#_blank.fits
maps/<parfile>_vlsr2_comp#_blank.fits
maps/<parfile>_atau2m_comp#_blank.fits
maps/<parfile>_tau2m_comp#_blank.fits
maps/<parfile>_Tex_comp#_blank.fits
maps/<parfile>_Trot_comp#_blank.fits
maps/<parfile>_Tk_comp#_blank.fits
maps/<parfile>_f0N11_comp#_blank.fits
maps/<parfile>_f1N11_comp#_blank.fits
maps/<parfile>_f0N22_comp#_blank.fits
maps/<parfile>_f1N22_comp#_blank.fits
maps/<parfile>_f0Nnh3_comp#_blank.fits
maps/<parfile>_f1Nnh3_comp#_blank.fits