GEIHIT
Latest: GEIHIT v4.00 (09AUG18)		15NOV24
This needs updating for HITRAN 2020!

GEISA 2015 data is available from

• cds-espri.ipsl.upmc.fr/etherTypo/?id=1732

Since most of the essential parameters for radiative transfer modelling are the same in GEISA and HITRAN the conversion doesn't require anything too complicated [Details...], but the following adjustments are made by GEIHIT

• GEISA contains molecules and isotopologues which are not defined in HITRAN 2016 - these are discarded by default, but can also be included (see GEISA v HITRAN Molecules)
• The following GEISA lines are excluded
  • Initial 10 lines of H2CO which have wavenumber 0.0
  • Many lines (particularly NH3 and SO2, but also PH3 and H2S) without any defined air-broadened half-width
  • H2O lines without a temperature coefficient for the half-width
  A summary of excluded lines is printed to the terminal at the end of the run.
• The GEISA values sometimes have an extra significant figure compared to HITRAN, in which case the value is truncated to fit the HITRAN format.
• GEISA vibrational levels are described by 25 characters compared to 15 for HITRAN. In most cases it is possible to rewrite these in the HITRAN form (subroutine GQFORM and, specifically for NH3, NH3FORM) but the GEISA notation for C4H2 and C2HD (3rd isotopologue of C2H2) is anomalous and so is left unchanged.
• GEISA accuracies are given in physical units, while HITRAN uses integer codes for accuracy ranges - the conversion is performed by function IDXACC.
• There is no GEISA equivalent for the HITRAN Reference fields (identifying the original source of each transition), the line mixing flags, or the upper and lower statistical weights. These are set to blank or 0.0 in the output.

Download the source code: [geihit.f90]

Compile with any generic FORTRAN compiler, eg

gfortran geihit.f90 -o geihit

(Note: GEIHIT makes use of the intrinsic MOVE_ALLOC subroutine introduced with FORTRAN 2003, but any generic FORTRAN compiler should cope with this.)

Running GEIHIT

geihit R-GEIHIT: Running program GEIHIT v4.00 [1] GEISA 2015 input file: line_GEISA2015_asc_gs08_v1.0_all [2] HITRAN output file: geisa.par Include non-HITRAN mol/iso? (Y/N) [=N]: Y I-HITMOL: Assign HITRAN Mol,Iso ID=62 1 to GEISA Mol,Iso=HNC 142 [3] I-HITMOL: Assign HITRAN Mol,Iso ID=23 4 to GEISA Mol,Iso=HCN 224 I-GEIHIT: Line# 100000 6.996648 cm-1 I-GEIHIT: Line# 200000 10.781763 cm-1 I-GEIHIT: Line# 300000 14.669723 cm-1 ... I-HITMOL: Assign HITRAN Mol,Iso ID=51 1 to GEISA Mol,Iso=C3H4 341 ... I-HITMOL: Assign HITRAN Mol,Iso ID=63 1 to GEISA Mol,Iso=C6H6 266 ... I-HITMOL: Assign HITRAN Mol,Iso ID=61 1 to GEISA Mol,Iso=C3H8 221 ... I-HITMOL: Assign HITRAN Mol,Iso ID= 4 8 to GEISA Mol,Iso=N2O 556 ... I-GEIHIT: Line#3400000 1792.480720 cm-1 I-HITMOL: Assign HITRAN Mol,Iso ID=60 1 to GEISA Mol,Iso=GeH4 411 I-HITMOL: Assign HITRAN Mol,Iso ID=19 6 to GEISA Mol,Iso=OCS 634 I-GEIHIT: Line#3500000 1982.129400 cm-1 I-GEIHIT: Line#3600000 2111.540140 cm-1 I-HITMOL: Assign HITRAN Mol,Iso ID= 4 6 to GEISA Mol,Iso=N2O 458 I-HITMOL: Assign HITRAN Mol,Iso ID= 4 7 to GEISA Mol,Iso=N2O 548 I-GEIHIT: Line#3700000 2270.603842 cm-1 ... I-GEIHIT: Line#5000000 11437.821300 cm-1 W-GEIHIT: Summary of excluded lines ... H2CO 10 lines from 0.000 - 0.000 cm-1: invalid wavenumber NH3 12281 lines from 0.327 - 1515.168 cm-1: no air half-width data H2O 504 lines from 247.760 - 6676.146 cm-1: no T.Coeff data PH3 264 lines from 866.468 - 3490.690 cm-1: no air half-width data H2S 153 lines from 1100.120 - 4195.623 cm-1: no air half-width data SO2 5861 lines from 1301.624 - 1397.007 cm-1: no air half-width data STOP R-GEIHIT: Successful completion [4]

Notes

1. The first line tells you which version (4.00) of the program is being executed
2. Then follows the dialogue which asks for
   1. GEISA filename.
   2. HITRAN-type output filename
   3. Whether you want to include additional molecules and isotopologues in the output which are not standard in HITRAN 2016 (type 'Y' or 'y' if 'yes'; anything else, including <CR> in this case, interpreted as 'no').
3. The program then prints a series of progress messages which will be one of the following

   I-GEIHIT: Line# ...

   a progress message every 100 000 input lines, with the associated wavenumber

   I-HITMOL: Assign HITRAN Mol,Iso ID= ... (only if non-HITRAN mol/iso selected for inclusion)

   Notification of assignment of non-standard HITRAN molecular,isotopic ID on first occurrence

   W-GEIHIT: Summary of excluded lines ...

   Then follows a summary for each molecule of number of lines excluded, wavenumber range spanned by excluded lines, and reason for exclusion.

4. The program should end with the statement: STOP R-GEIHIT: Successful completion

• HITRAN2016 defines 49 molecules with line-transition parameters [List], of which three (#30 SF6, #35 ClONO2 and #42 CF4) are generally not included in line lists since the tabulated cross-section data are more complete/accurate.
• GEISA 2015 defines 52 molecules [List] of which 3 (CH3D, C2HD and HDO) are deuterium isotopologues, leaving 49 different molecules.
• 44 molecules are common to both HITRAN2016 and GEISA 2015
• The 5 GEISA molecules not included in HITRAN are: HNC (NB distinct from HCN), GeH4, C3H4 (Propyne), C6H6 (Benzene) and C3H8 (Propane)
• The 5 HITRAN 2016 molecules not included in GEISA 2015 are: O (monatomic oxygen), HOBr, H2, CS, COCl2 (Phosgene)
• GEISA has some additional isotopologues which are not present in HITRAN2016 for N2O, OCS and HCN.

If you type 'Y' to the prompt for non-HITRAN species the new assignments are listed in the table on the right (as well as in the I-HITMOL messages generated through the run). The new molecule# assignments are somewhat arbitary: 51 was chosen for C3H4 since it is the number assigned to this number in the new Gamache TIPS data (which extends the HITRAN molecular numbering up to 53), and the other molecules, which are not in the TIPS data, just assigned #60 onwards. Isotopologue assignments are more meaningful and represent decreasing order of abundance (#1 = most common isotopologue). For the new molecules, GEISA data only contains the main isotopologue. Edit the subroutine HITMOL if you want to change these assignments.

List GEISA2015 Molecules/Isotopes not defined in HITRAN2016 GEISA Assigned Molec Iso Molec Iso N2O 458 4 6 N2O 548 4 7 N2O 556 4 8 OCS 634 19 6 HCN 224 23 4 C3H4 341 51 1 GeH4 411 60 1 C3H8 221 61 1 HNC 142 62 1 C6H6 266 63 1

Mapping from GEISA to HITRAN

• RFM home page