GEIHIT | ![]() | ![]() |
Latest: GEIHIT v4.00 (09AUG18) | 22FEB25 | |
This needs updating for HITRAN 2020! |
GEISA 2015 data is available from
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
Download the source code: [geihit.f90]
Compile with any generic FORTRAN compiler, eg
In this example the gfortran compiler is used [More...] and the executable is given the name geihit.gfortran geihit.f90 -o geihit
(Note:
A typical run might be
(
The program takes a few minutes to complete.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
In most cases the GEISA
format contains the equivalent
HITRAN molecular and isotopologue ID, with values of (-9, 0) for (molecules,
isotopologues) which are not recognised by HITRAN.
However GEISA 2015
doesn't know about the new additions for HITRAN 2016 so the
Note the use of '0', 'A' and 'B' in the HITRAN .par format for isotopologue#10, 11 and 12. |
|
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. |
|
Parameter | Format | GEISA | Missing | Notes |
---|---|---|---|---|
Molecule# | I2 | Field K | Set in HITMOL | See GEISA v HITRAN Molecules |
Isotopologue# | A1 | Field L | Set in HITMOL | See GEISA v HITRAN Molecules |
Wavenumber | F12.6 | Field A | n/a | Exclude line if value = 0.0 |
Intensity | 1PE10.3 | Field B | n/a | Truncate by 1 digit |
Ein. A Coeff | 1PE10.3 | Field M | Set to 0.0 | |
Air Half-Width | F5.3/F5.4 | Field C | Exclude | Truncate by 1 digit if ≥ 1.0 |
Self Half-Width | F5.4/F5.3/F5.2 | Field N | Set to Air Half-Width | Truncate by 1,2 digits if ≥ 1,10 |
Lower Energy | F10.4 | Field D | Set to -1.0 | |
T-Coeff | F4.2 | Field F | Exclude | |
Pressure Shift | F8.6/F8.5 | Field O | Set to 0.0 | Truncate by 1 digit if ≥ 10 |
Upper Vib | A15 | Field E1 | Converted in GQFORM,NH3FORM | |
Lower Vib | A15 | Field E2 | Converted in GQFORM,NH3FORM | |
Upper Rot | A15 | Field E3 | Copied | |
Lower Rot | A15 | Field E4 | Copied | |
Error Codes | 6I1 | Fields A',B',C',N',F',O' | Set to 0 | Converted in IDXACC |
Ref Codes | 12X | undefined | n/a | set blank |
Line Mixing | 1X | undefined | n/a | set blank |
Upper Wgt | F7.1 | undefined | n/a | set 0.0 |
Lower Wgt | F7.1 | undefined | n/a | set 0.0 |