MORSE Input Files

Last Updated: 09FEB05

When MORSE is run, certain input files are always required: the driver table (morse.drv) and files associated with mandatory sections of the driver table. Additional files may be required by optional sections of the driver table. Click on File column in the table below for details of a particular file.

MORSE Input Files
Mandatory Sections
Section File Description
morse.drv Driver Table
*MWL *.lst Microwindow List
*MDB MW*.DAT Microwindow Database
*L1C L1C* Level 1C Data
*CLI *.atm Climatological Data
*ILS *.ils (Apodised) Instrument Line Shape
*FOV *.fov Field of View
Optional Sections
Section File Description
*ACV *.acv A Priori Covariance
*HIT *.bin HITRAN Binary File
*L2P *.atm Previously Retrieved Profiles
*LUT CS*DAT Look-Up Tables
*NTE *.nte Vibrational Temperatures (non-LTE)
*XSC .xsc Molecular Cross-Sections

Notes

  1. Most input files start with an arbitrary number of optional comment records, identified by the first character being `!'
  2. Blank records are ignored
  3. Numerical data records are read free format unless otherwise specifed
The rest of this document describes each input file in detail

morse.drv (mandatory)

Description
Driver Table
Controls operation of MORSE

Specified
Fixed name (see Note 1)

Format
See driver.html

Notes
  1. The filename morse.drv is hard-wired in subroutine msedrv.for.

Examples
morse_pt.drv for pT retrieval
morse_h2o.drv for H2O retrieval

*.lst (mandatory)

Description
List(s) of microwindows to be used in retrieval.

Specified
*MWL section of the driver table

Format
Records up to C*120
! [comment]
Any record with first character ! is ignored
SEQ# MW_Label WnoL WnoU AltL AltU [Abs(1) Abs(2) ...etc]
SEQ# I*4 Sequence Number (Dummy) - see Note 1 eg 1
MW_Label C*8 Microwindow Label - see Note 2 eg PT__0001
WnoL R*8 MW Lower Wavenumber Boundary [cm-1] eg 686.400
WnoU R*8 MW Upper Wavenumber Boundary [cm-1] eg 689.400
AltL R*4 MW Lower Tangent Altitude [km] eg 30
AltU R*4 MW Upper Tangent Altitude [km] eg 68
Abs(i) I*4 (Optional) HITRAN Index of absorber#i - see Note 3 eg 2 (=CO2)

Notes
  1. The first 4 characters of a non-comment record normally contain the number of the microwindow in order of selection, however these are not used by MORSE. Microwindows are processed in the order in which they appear in the .lst file.
  2. The Microwindow Label is assumed to occupy positions 5:8 in the record. All subsequent data is read free-format.
  3. If any absorber indices are specified in the .lst file these will override the list of absorbers in the Microwindow Database File
  4. Unlike MORSE, the L1C program requires a single list file in order to extract microwindows from L1B data. However this can be created from different microwindow lists (eg different species) simply by concatenating the individual files, eg
    cat pt_mws.lst h2o_mws.lst > all_mws.lst

Examples
pt_000.lst Microwindows for pT retrieval
h2o_000.lst Microwindows for H2O retrieval

MW*DAT (mandatory-ish)

Description
Microwindow Database file
File containing, amongst other things, the spectral masks, irregular grids and absorber lists associated with each microwindow.

Specified
*MDB section of the driver table

Format
These are ASCII files whose format is too complicated to detail here, but is described in ESA documents
The Revised Microwindow Database Format
Tech Note PO-TN-IMK-GS-0002 (April, 1999)
ASCII Input Data Interface Control Document (DD53)
PO-IF-DOG-GS-0002 (April, 1999)

Notes
  1. The *.imk microwindow files generated by mwmake program are in fact single microwindow versions of the MW*DAT files, minus the irregular grids which are derived by a separate process - see PO-MA-OXF-GS-0019.pdf

Examples
Directory ./MW contains Microwindow Database Files for various species and spectral resolutions.
See ./MW/README.txt for further details.

L1C* (mandatory)

Description
MIPAS L1C Data
Apodised spectra for selected microwindows

Specified
*L1C section of the driver table

Format
Described at http://www.atm.ox.ac.uk/group/mipas/L1C

Notes
  1. L1C data can be generated from ESA L1B spectra by the IDL program l1c.pro, or its predecessor lv2pp.pro.

Example
L1C_02081_SCAN_13 (Scan#13, Orbit 2081)

*.atm (mandatory)

Description
Climatological profiles
Used as both a priori and initial guess.

Specified
*CLI section of the driver table

Format
RFM .atm file format
! [comments]
Any record starting with ! is ignored
NLev
NLev I*4 Number of profile levels eg 121
*HGT
*HGT C*4 Indicates Height Profile to follow fixed
Hgt(1) Hgt(2) ... Hgt(NLev)
Hgt(i) R*4 Altitude [km] of profile level i eg, 0,1 ... 120
*Species
*Species Up to C*8 Indicates Species Profile to follow eg *H2O
Prf(1) Prf(2) ... Prf(NLev)
Prf(i) R*4 Profile value at level i eg 936.8 ...
*Species etc ...
Repeat for next Species, etc ..
*END
*END C*4 End of file marker fixed

Notes
  1. Comment records may be inserted anywhere except between the *Species record and the subsequent records containing profile values
  2. The *Species record is read as far as the first space (C*8 max, so maximum length of the species name is C*7), anything else on the same record is ignored. Also case insensitive (everything converted to lower case internally).
  3. Apart from *HGT, special *Species names are
    • *TEM Temperature [K]
    • *PRE Pressure [mb]
    • *AEROSOL Extinction [/km]
    Anything else is compared to the list of absorbing species and, if matched, assumed to be the concentration in [ppmv]. Anything not matched is ignored.

Examples
ESA `Initial Guess' Climatology, as created by U.Leicester are in directory ./IG - see ./IG/README.txt
Also http://www.atm.ox.ac.uk/RFM/rfm_downloads.html#atm for other `atmospheres'.
aerosol.atm `zero' aerosol profile.

*.ils (mandatory)

Description
(Apodised) Instrument Line Shape file(s)

Specified
*ILS section of the driver table

Format
RFM .ils file format.
! [comments]
Initial records starting with ! are ignored
Npt Pt1 PtD [Wno1 Wno2]
Npt I*4 Number of tabulated points eg 701
Pt1 R*8 Relative wavenumber [cm-1] of 1st point eg -0.175
PtD R*8 Tabulation wavenumber increment [cm-1] eg 0.0005
Wno1 R*8 (Optional) Lower wavenumber limit [cm-1] for application eg 685.0
Wno2 R*8 (Optional) Upper wavenumber limit [cm-1] for application eg 970.0
ILS(1) ILS(2) ... ILS(Npt)
ILS(i) R*8 ILS function tabulated at point i eg 8.7317463e-07

Notes
  1. A single ILS file can be used for the complete spectrum (in which case omit Wno1, Wno2 in header), or individual files for different spectral ranges
  2. Where both a general and a spectrally limited ILS file can be used for a particular microwindow, the spectrally limited file will be used in preference. Any other ambiguities (eg 2 general files, 2 spectrally limited files applicable to the same microwindow) will result in a fatal error message.
  3. ILS files can be simply concatenated
  4. ILS functions are renormalised so may be scaled by any arbitrary factor.
  5. ILS functions are interpolated to the internal wavenumber grid so changing Pt1 is a simple way to introduce a frequency shift and change PtD to change the AILS width.

Examples
ofm.ils Idealised AILS for 0.025cm-1 sampling, Norton Beer Strong Apodisation, applicable over whole range.
mipas_a.ails Pre-launch computation of the MIPAS A-band AILS.

*.fov (mandatory)

Description
Field of View file

Specified
*FOV section of the driver table

Format
RFM .fov file format.
! [comments]
Initial records starting with ! are ignored
Npt
Npt I*4 Number of tabulated points eg 5
Alt(1) Alt(2) ... Alt(Npt)
Alt(i) R*4 Relative altitude [km] of point i eg -1.5
FOV(1) FOV(2) ... FOV(Npt)
FOV(i) R*4 Tabulated FOV Response at point i eg 0.2798

Notes
  1. The end two tabulated points FOV(1), FOV(Npt) must be zero (in effect, it the positions Alt(1), Alt(Npt) rather than the values of these two points define the edges of the FOV)
  2. The FOV function is renormalised internally, so any arbitrary scaling can be applied to the tabulated values.

Examples
rfm_1km5.fov Field of View represented as 5 points with 1.5km spacing

*.acv (optional)

Description
A Priori Covariance Matrix

Specified
*ACV section of the driver table

Format
Same as morse.cov output file.
! [comments]
Initial records with first character `!' are ignored
N M
N I*4 First dimension of matrix eg, 36
M I*4 Second dimension of matrix (See Note 2) eg, 36
SA(1,1) SA(2,1) ... SA(N,M)
SA(i,j) R*8 Element i,j of the Covariance Matrix eg, 0.151308903009384

Notes
  1. Although covariance matrices are always square, the file format is intended to be used for any matrix, so both dimensions M, N are specified (it is checked that these are both the same when the file is read).

Examples
morse.cov

*.bin (optional)

Description
HITRAN (-type) Binary File

Specified
*HIT section of the driver table

Format
HITRAN ASCII data converted into a direct-access binary file using the RFM program hitbin.for
See www.hitran.com

Notes
  1. The same binary files can be used on most unix and linux systems, although for Windows and Sun machines it may be necessary to regenerate the file using hitbin since the byte-ordering appears to be different.
  2. There is a MIPAS-specific database in HITRAN ASCII format - see http://www.ifac.cnr.it/retrieval/database.html. Once you have obtained the necessary permission we can give you the binary version for use with MORSE.

Example
File hitran_dummy.bin in directory ./bin (just H2O lines from 791-808cm-1)

*.atm (optional)

Description
Previously retrieved profiles from the same L1B scan

Specified
*L2P section of the driver table

Format
RFM .atm file format.

Notes
  1. Same format as input files for *CLI section except that a *SWEEP profile is also accepted.

Examples
morse_pt.atm
morse_h2o.atm

CS*DAT (optional)

Description
SVD-compressed Look-Up Tables of absorption coefficients

Specified
Directories specified in *LUT section of the driver table
Filenames are constructed in the form:
CS_[MW_Label]_[ii].DAT
[MW_Label] A8 Microwindow Label eg H2O_0001
[ii] I2.2 HITRAN Index of the absorber eg 02 (=CO2).

Format
See PT-TN-OXF-GS_0011.txt

Notes
  1. The SVD-compressed LUTs are those used in the ESA retrieval, generated as part of the microwindow processing suite mwproc (see PO-MA-OXF-GS-0019.pdf)
  2. In theory, instead of the SVD-compressed LUTs, the LUTs generated directly by the RFM using the TAB option can also be used (provided they have the same filenames), however this hasn't been tested.

Examples
In directory ./CS (see ./CS/README.txt)

*.nte (optional)

Description
Vibrational Temperatures (non-LTE calculations)

Specified
*NTE section of the driver table

Format
RFM .nte file format.

Examples

.xsc (optional)

Description
Heavy Molecule Cross-Section data

Specified
Directories specified in *XSC section of the driver table
Filenames are constructed in the form:
[gas].xsc
[gas] up to C*7 name of the molecule (See Note 1) eg f12

Format
RFM .xsc file format, except restricted to one molecule per file. Basically the same format as the heavy molecule cross section data distributed with HITRAN except for the addition of standardised header records.

Notes
  1. Molecule names are expected to be lower case and using the F.. codes for CFCs eg f11 for CFCl3

Examples
See http://www.atm.ox.ac.uk/RFM/rfm_downloads.html#xsc
aerosol.xsc flat aerosol spectrum.
Back to MORSE Index