RFM Driver Table Sections

*ATM Section: Atmospheric Profiles

25MAR19

Type
Primary Section#5

Description
Atmospheric profile data.
Usually one or more .atm format files describing the atmosphere to be modelled.

Format
Multiple fields, order is significant.
FieldTypeDescription
FILATM C200 Name of an .atm file
PARAM=VALUE C=R Explicit assignment of a constant value to a profile

QualifierTypeDescriptionRange
(PSI) R Horizontal angle [deg] -90:+90
Type: I=Integer; R=Real; D=Double Precision; Cn=character string, length n.

Notes
  1. The first profile in the first file should be either *HGT (height, km) or *PRE (pressure, hPa), and this defines the internal vertical grid to be used by the RFM, the first (lowest) grid point (not necessarily 0km), defining the surface and the last (highest) grid point defining the top of the atmosphere. For this reason you may find it convenient to have this first file contain a single *HGT or *PRE profile just to define the grid.
  2. Having established whether a height or pressure grid is used, all subsequent files must also contain a profile of the same type (although not necessarily as the first profile in the file) to allow the RFM to interpolate to the internal grid.
  3. For the vertical interpolation, Temperature, height/ln(p) and ln(VMR) are assumed to vary linearly with the vertical coordinate ln(p)/height. The LIN Flag forces VMR rather than ln(VMR) to vary linearly with vertical coordinate.
  4. All data in subsequent files replaces any previously loaded data, e.g. if *TEM (temperature) profiles appear in the first and second files, the value in the second file will be used. Hence the order matters.
  5. If subsequent .atm files include profiles specified over smaller vertical ranges than the internal grid, the end values of the short profiles are duplicated, i.e., the short profiles are extended over the full vertical range assuming the end values remain constant. The exception is pressure, which is linearly extrapolated on ln(p) from the pair of end values.
  6. Apart from the initial user-defined grid, extra profile levels may be inserted by use of the JAC, LAY, LEV or OBS Flags.
  7. The PRF Flag can be used to write out the actual profile that ends up being used in the RFM.
  8. The header record of each accepted profile, plus *END markers, are written to the rfm.log file, together with notification if this replaces an existing profile.
  9. After all .atm files are read, there is a check that profiles of Temperature, Pressure and VMR for absorber specified in the *GAS section have been found (& extinction profile for aerosol), otherwise a fatal error message results.
  10. A pressure profile is always required but, with a pressure grid, altitude is only required for limb-viewing calculations, if the HYD flag is used, or to identify levels for perturbations (JAC flag) or output levels (LEV or FLX flags).
  11. The PARAM=VALUE format, where PARAM represents any parameter that might be found in an .atm file (without the initial '*' character, e.g. TEM=280) can be used as a shorthand for setting a constant profile value. This is also useful for setting homogeneous path conditions (HOM Flag).
  12. Isotopes: if isotopes are specified in any input files (eg by having *CH4(1) as one of the profiles) then the RFM treats these as separate gases - one gas for each specified isotopic profile plus one gas for the remaining isotopes (even if profiles have been supplied for all isotopes). Any subsequent profiles for just *CH4 will not overwrite the isotopic profile.

    2D Atmospheres

  13. The default is to assume 1-dimensional atmospheres, i.e. no variation in atmospheric parameters along profile surfaces. The GRA flag enables 2D atmospheres, and profile files then require an extra coordinate, PSI to define the horizontal location and (usually) the OBS flag to specify the observer's horizontal location (as PSIOBS — see *OBS section).
  14. Angles PSI are defined as subtended at the local centre of curvature, which is approximately the same as the spacing of degrees of latitude, i.e., PSI=1 corresponds to a distance of around 111 km on the surface of the earth. Direction is defined so that the observer is at high values looking towards low values, i.e. so that profiles with lower values of PSI are further away. Tolerance for distinguishing different PSI angles (fixed by PSITOL parameter in prfgra.for) is 1E-3 deg, approximately 100m along the surface. [simple geometry also tells you that any variation in elevation angle of the observation will also move the tangent point location, PSI, by the same angle]
  15. Absolute PSI values are defined in one of two ways, depending on whether the OBS flag enabled
    • OBS flag enabled (usual case): PSI values are defined relative to the observer location PSIOBS in the *OBS section.
    • OBS flag disabled: PSI=0 is defined as the location of all the tangent points either refracted or geometric according to the specification of the *TAN section. In this case, multiple tangent paths implies multiple observer locations.
  16. When the GRA Flag is enabled, a two dimensional atmosphere is established so that profiles are written (or overwritten) at each location PSI independently, sharing only the vertical grid established by the first .atm file in the *ATM section. However, it is not necessary to supply profiles for all absorbers at all locations since, after all profiles have been read, there is an interpolation step to fill in any missing profiles from those that have been supplied. Suggested use is that the atmosphere is fully defined for PSI=0, and only the horizontally varying profiles are supplied at other locations.
  17. If no profile is explicitly supplied at PSI=0, then a reference profile at PSI=0 is created by interpolating the supplied profiles.

Examples
Zero-dimensional atmosphere (i.e., homogeneous path)
*ATM
  TEM=296
  pre=1013.25
  H2O=100                  ! 100ppmv of water vapour        
One-dimensional atmosphere
*ATM                
     hgtgrd.atm     ! Define altitude grid for calculation (z only)
     std.atm        ! US Std Atmosphere (z,p,T,H2O,CO2,O3,N2O,CO,CH4,O2)
     minor.atm      ! Minor species (HNO3, etc)
     co2_2005.atm   ! Replace US Std. CO2 vmrs with 2005 concentrations
Two-dimensional atmosphere, defined at 3 locations
*ATM
     std.atm(0)     ! US Std Atm, all species, at PSI=0 
     t_near.atm(2)  ! "near" side temperature profile at PSI=2 deg
     p_near.atm(2)  ! Pressure profile, also at PSI=2 deg
     t_far.atm(-1)  ! "far" side temperature profile, at PSI=-1 deg
     minor.atm(2)   ! Minor species, initially placed at PSI=2
                    ! All std and minor species apart from p,T will be 
                    ! duplicated at PSI=-1, 0 and +2. 
                    ! The std p profile at PSI=0 will be duplicated at PSI=-1

Bugs
Bug#9 (Fixed v5.02)
Bug#109 (Fixed v4.34)
Bug#82 (Fixed v4.31)