RFM Driver Table Sections

*GAS Section: List of Absorbers

01APR19

Type
Primary Section#4

Description
List of absorbing species required for calculations
Otional qualifiers to select particular bands, isotopes and continuum handling.

Format
Multiple fields, arbitrary order.
FieldTypeDescriptionRange
FILGAS C200 Name of a .gas file
GAS C7 Chemical Formula. See List of Absorbers
IDX I HITRAN/RFM index of molecule See List of Absorbers
* C1 All significant molecules See Notes.
Type: I=Integer; R=Real; D=Double Precision; Cn=character string, length n.

Various 'qualifiers', contained in a pair of brackets (...), may be appended immediately after (no spaces) the GAS, IDX or * fields, and their order may be significant.
QualifierTypeDescriptionRange
(ISO) I HITRAN Isotope index ±1:(N.isotopes for molecule)
(LGQ:UGQ) 2I Lower, Upper HITRAN Global Quantum indices 1:999 or '*'
(CTMGAS) C5 Continuum treatment for gas 'CTM' or 'NOCTM'
(LIM) I Threshold for '*' absorbers 0:99
(H2OWNG) I Treatment of H2O 25cm-1 offset 'SUB' or 'NOSUB'

Notes
  1. The type of field is identified by first checking to see if it is an existing file (FILGAS), then if it is a '*' character, then if it can be interpreted as a numerical index (IDX) and finall as a molecule name (GAS).

    Wildcard

  2. The wildcard character '*' can be used on its own or in conjunction with any other gases (in which case it will add any additional absorbers), but may only appear once in the *GAS section, and not in any .gas file. The main use of this is if you do not know in advance which absorbers might be significant in a particular spectral range. The RFM has an internal database (in file optdat_dat.f90) of approximate optical depths for a zenith path through a standard atmosphere of all absorbing species in the range 0–20 000 cm-1.
  3. The gas wildcard character has optional 'optical strength' qualifier (LIM) which can select a cut off for absorbers, eg '*(10)' is a reasonable value to exclude all undetectable absorbers. A value '*(0)' is equivalent to just '*'. The 'optical strength' parameter is notionally 15 + log10(max optical depth), so a value 15 corresponds to optical depth 1 (hence '*(10)' excludes all absorbers with optical depths < 1E-5 in a vertical path).
  4. If the gas wildcard is used, an entry appears in the rfm.log file listing all absorbers found together with their 'optical strengths', i.e. you could just use the RFM to produce an ordered list of significant absorbers within any spectral range (performed by the subroutine gasall_sub.f90).

    Continuum

  5. The (CTM) and (NOCTM) qualifiers are only applicable to molecules with continuum models: CO2, H2O, O2, N2 and the the CTM Flag (*FLG section) is also required to have any effect.
    • (CTM): Only the continuum component will be modelled for this molecule - all HITRAN line transitions will be omitted (so no other qualifiers may be used).
    • (NOCTM): Exclude the continuum component for this molecule. This may be placed anywhere in the qualifier string.
  6. H2OWNG The definition of the standard H2O continuum is that it includes the cumulative absorption of all line wings beyond 25 cm-1 from the line centre. So to combine H2O line and continuum terms the standard procedure would be to calculate the optical depth contribution of each H2O line at ±25 cm-1 from line centre, and subtract this value across the whole line shape between ±25 cm-1. This is all handled automatically within the RFM, depending on whether or not the H2O continuum is included in the calculation.

    An issue arises when trying to model H2O absorption using look-up tables (TAB, LUT flags). Since the H2O continuum term has a strong dependence on p,T and H2O concentration it requires a high density of points to be modelled in a LUT. On the other hand, the continuum part is relatively simple to compute explicitly. So it is usually desirable to have H2O LUTs which contain only the truncated H2O line contributions which can then simply be added to the computed continuum term.

    The RFM anticipates this behaviour so, if calculating LUTs (TAB flag, without the CTM flag) it will subtract the 25 cm-1 offset from all H2O lines, and if using LUTs (LUT flag, with the CTM flag) it will simply add these contributions to any H2O continuum term. The H2OWNG qualifier allows the user to alter this default behaviour, eg

    • H2O(SUB): calculate H2O lines with the 25 cm-1 subtraction but without the continuum. This might be useful for a line-by-line calculation to verify the accuracy of spectra generated with a standard H2O .tab file
    • H2O(NOSUB): override-continuum subtraction with TAB flag. This could be used to force the RFM to generate H2O .tab fileS which include the continuum term.
    Note that the RFM itself has no way of distinguishing whether an H2O LUT has been calculated with or without the continuum term.

    Isotopologues (see RFM Handling of Isotopologues)

  7. The Isotope (ISO) and band (LGQ:UGQ) qualifiers may only be used for 'line' molecules (i.e. not cross-section molecules), those with HITRAN/RFM indices 1:99. The effect is to apply a selection criterion to the HITRAN database so that only transitions from the specified isotopes/transitions are read in.
  8. The isotope number (ISO) refers to the HITRAN nomenclature, ordered such that 1=most abundant. Negative values can also be used, interpreted as excluding particular isotopes.
  9. As special cases, the RFM accepts commonly used abbreviations 'HDO' and 'CH3D' as equivalent to 'H2O(4)' and 'CH4(3)'.

    Vibrational Bands (See RFM Vibrational Level Assignments)

  10. The vibrational bands (LGQ:UGQ) are identified by the RFM indices for each vibrational state, whose definition depends on the form of the molecule but 1=ground state in all cases. Either the lower or upper indicies numbers may be replaced by a wildcard character '*', indicating that all transitions with a particular upper/lower level are to be included. However, these are checked for logical consistency with other selections, eg you can't have (1:*) and (1:4) in the qualifier list for the same molecule.
  11. More complicated selection rules can be applied by combining the (ISO) and (LGQ:UGQ) qualifiers. The general form of the qualifiers is [list_of_isotopes][list_of_bands], and implicitly the list starts with all isotopes and finishes with all bands; eg ...(1)(2)(3:4)(1:6)(3)(2:6)... selects bands (3:4) and (1:6) for isotopes (1) and (2), and band (2:6) for isotope (3). If the first qualifier(s) (or only qualifiers) are a list of bands, then these are selected for all isotopes. If the last qualifier(s) (or only qualifiers) are a list of isotopes, then all bands are selected for these isotopes.

Example
*GAS                       
   H2O CO2 N2O CH4 HNO3 O3 ! MIPAS target species
   O2(CTM) N2(CTM) ! Oxygen and Nitrogen - continua only
   minor.gas       ! a file of minor species
   *(10)                   ! Plus any other significant absorbers 
   CO(1)(1:3) NO(1)(1:*)   ! Only primary isotopes, and specific bands

Bugs
Bug#15 (Fixed v5.02)
Bug#114 (Fixed v4.35)
Bug#111 (Fixed v4.34)
Bug#90 (Fixed v4.32)
Bug#89 (Fixed v4.32)
Bug#87 (Fixed v4.32)