SUBHIT
Latest: SUBHIT v2.00 (27NOV18) [Download]		08JAN24

Introduction

For example, this could be useful for merging different molecules from different datasets, starting with this program to extract particular molecules and then using the MRGHIT program to reassemble subsets.

Installing SUBHIT

Download the source code: [subhit.f90]

Then compile with any generic Fortran compiler, eg

gfortran subhit.f90 -o subhit

subhit R-SUBHIT: Running SUBHIT v2.00 Input HITRAN file: HITRAN2012.par Wavenumber range (cm-1) [<CR>=all]: 800 900 [1] List of Molec# to include: <CR> [2] List of Molec# to exclude: <CR> [3] Output HITRAN subset file: example.par I-SUBHIT: Record# 100000 Wavenumber= 829.368699 I-SUBHIT: Record# 200000 Wavenumber= 840.987280 I-SUBHIT: Record# 300000 Wavenumber= 850.284961 I-SUBHIT: Record# 400000 Wavenumber= 859.502340 I-SUBHIT: Record# 500000 Wavenumber= 874.569640 I-SUBHIT: Record# 600000 Wavenumber= 897.545360 I-SUBHIT: Number of records written = 609653 STOP R-SUBHIT: Successful completion

subhit R-SUBHIT: Running SUBHIT v2.00 Input HITRAN file: HITRAN2012.par Wavenumber range (cm-1) [<CR>=all]: <CR> List of Molec# to include: 2.9 2.10 2.11 2.12 Output HITRAN subset file: co2_iso.par I-SUBHIT: Number of records written = 15125 STOP R-SUBHIT: Successful completion

Notes

1. Wavenumber range: just typing <CR> (ie Return/Enter) would reproduce the full spectral range of the original file
2. Included molecules: list of molecules, and only these molecules to included in the output.
   • Specify molecules by their HITRAN numbers separated by spaces (max 200 characters in total), eg 1 2 3 for H2O, CO2 and O3
   • Use numbers after a decimal point to indicate isotopologues, following the HITRAN numbering, eg 1 3 2.1 2.2 would select all H2O and O3 lines, but only isotopologues 1 (626) and 2 (636) for CO2
   • Note that 2.10, 2.11 and 2.12 are used for CO2 838, 837 and 737 rather than the HITRAN designation 0, A and B.
   Typing <CR> means that you don't want to specify a list.
3. Excluded molecules (prompt only appears if user typed <CR> in response to previous). A list of molecules, isotopologues and only these molecules to be excluded from the output, same format as previous. Type <CR> again if all molecules to be included.
4. The program prints a series of progress messages after every 100 000 records of output data are written, and ends with a summary of the number of records written. If the user specifically excluded particular molecules/isotopes there will be an additional summary line giving the number of records excluded.

• RFM home page