In this sense it is similar to the RFM, although in the RFM most optional sections require the relevant 3-character flag in the *FLG section to 'enable' the option, while MORSE automatically uses any optional section in the driver table, without any need for an associated flag.
Structure of MORSE Driver Table | ||
---|---|---|
Mandatory sections (in this order) | ||
1. | Comment record written to output file headers | |
2. | Option flags | |
3. | Input spectra | |
4. | | Vertical grid as Altitude levels, or |
| Pressure levels, or | |
| Measurement Tangent Point levels | |
5. | | Target quantities to be retrieved |
6. | | Microwindow list |
7. | | Atmospheric profiles |
8. | | Instrument Line Shape files, or |
| Spectral Response Function files | |
9. | | Field of View Shape (omitted if nadir-viewing) |
Optional sections (in any order) | ||
| Accuracy settings | |
| A Priori Covariance | |
| A Priori Standard Deviation | |
| Cloud Detection Criteria | |
| Convergence Criteria | |
| Internal spectral fine grid | |
| Irregular Spectral Grid | |
| HITRAN Database | |
| Look Up Table Filenames | |
| Vibrational Temperatures for non-LTE forward model | |
| Specify Output Files | |
| Limit range of pixels to be retrieved | |
| Surface parameters | |
| Specify line shapes for line-by-line calculations | |
| SVD-compressed Look Up Tables | |
| Directories for Molecular Cross Sections | |
Termination Record (last record that will be read) | ||
where ABC is some code (case-insensitive) defining the contents of the section. These 4 characters have to occupy the first 4 positions in the record.*ABC
The driver table is terminated with
(remember to add a 'carriage return' <CR> after the *END)*END
The first 8 sections (9 for limb-viewing) are mandatory and the sequence is fixed. Thereafter sections are optional and in any order, terminated with the *END record.
Only the first 4 characters of any record starting with * are read, so comments may be added to these records without any exclamation mark if you need to remind yourself what each does.
Type Description Storage
I Integer 4 bytes
R Real 4 bytes
D Double Precision 8 bytes
Cn Character string, length n n bytes
The rest of this document describes the contents
of each section.
which will appear in the second record in output files (eg for orbit 12345) as
*HDR Section (#1)
Field Type Description
C79 Text for output file
headers
*FLG Section (#2)
Field Type Value Description
C3
ACC
More stringent default accuracy/convergence criteria
AHY
Suppress hydrostatic equilibrium
APR
Output A Priori information
in .rtv files
CLD
Use cloud detection
CTM
Use molecular (eg H2O) continuum data
DFS
Include Degrees of Freedom for Signal
in .rtv files
FGD
Use full, as opposed to irregular, grid
FIX
Assume fixed relative altitudes (eg detector array)
FMW
use full microwindow (ie set all masks = TRUE)
FVZ
Suppress FOV distortion due to refraction
GEO
Use geometric (ie non-refractive) ray paths
GRA
Use 2D atmospheric model
INF
Include Information Content [bits]
in .rtv files
LIN
Assume VMR varies linearly with
altitude rather than logarithmically
MIX
Use CO2 line mixing
MWO
Write outputs after each microwindow
in .rtv files
QAD
Use simple quadratic fit to line wings
rather than inverse quadratic
RES
Include vertical resolution factor
in .rtv files
RSM
Allow run to resume from last pixel
*L1C Section (#3)
Field Type Description Requirement
C200
Name of .l1c file
Mandatory
C200
Name of .noi file
Mandatory for Nadir
*ALT, *PRE *TAN Section (#4)
Label Grid
*ALT Altitude grid
*PRE Pressure grid
*TAN Tangent Height grid
[2]
Field Type
Description
C200
Name of
.atm File
*RTV Section (#5)
Section structure
Retrieval Grid, single record
Retrieval Parameters
[and retrieval grid subsets], multiple records
Field
Type
Value Description Units Restrictions
Cn
*
Retrieve on all (NLev) grid levels
Only item in record
0
Column retrievals
Only item in record
i
Retrieve on level i
1 : NLev
i:j
Retrieve on all levels i:j inclusive
1 : NLev
i:j:k
Retrieve on levels i:j with step k
1 : NLev
C14
TEM
Temperature [2]
K
PRE
Pressure [2]
hPa
Not on *PRE levels
CTM
Continuum [3]
km-1
AEROSOL
Aerosol [4]
km-1
SFCTEM
Surface Temperature
K
OFF
Radiometric Offset [3]
nW/(cm2.sr.cm-1)
[gas]
Molecule VMR[5]
ppmv
from
List of gases
*MIC Section (#6)
Field
Type
Description
Cn
Name of Microwindow File
eg nh3.mic
Section Structure - Nadir Viewing
LAB
WNOMIN
WNOMAX
NABS
[NCHN]
Rec#1: Microwindow Header
ABS(1)
ABS(2) Rec#2: Microwindow absorbers
[CHN(1)
CHN(2) Rec#3+: Microwindow channels
Repeat for further microwindows
ISEQ LAB WNOL WNOU ALTL ALTU [ABS(1), ABS(2) ...]
[repeat for further microwindows]
Field
Type
Description Range
ISEQ
C4
Microwindow sequence number (not used)
eg ' 1 '
LAB
C8
Microwindow Label
eg PT__0001
WNOL
D
MW Lower Wavenumber Boundary [cm-1]
eg 686.400
WNOU
D
MW Upper Wavenumber Boundary [cm-1]
eg 689.400
ALTL
R
MW Lower Tangent Altitude [km]
eg 30
ALTU
R
MW Upper Tangent Altitude [km]
eg 68
ABS(i)
I
RFM/HITRAN
Index of absorber#i (Optional)
eg 2 (=CO2)
MDBFIL=MDBFIL
Field
Type
Description Range
MDBFIL=MDBFIL
C*
(Optional) name of microwindow database file (MW*DAT)
eg MW_PT__040.DAT
*MIC
NH3_0001 930 940 4 5 ! use 930-940cm-1 range, 4 absorbers, 5 channels
nh3 co2 h2o o3
930.75, 931.0, 931.25, 931.5 931.75
*MIC
pt_016.lst ! pT microwindows
MDBFIL=MW_PT__510.DAT ! pT FR mode database
*ATM Section (#7)
Field Type Description
ATMFIL C200
Name of
.atm,
.orb or
.glo file
Qualifier Type Description
(-QAL1;QAL2;...) C200
List of species to exclude from file
(QAL1;QAL2;...) C200
List of only species to include from file
*ILS Section (#8)
Field Type Description
ILSFIL C200
Name(s) of .ils file(s)
*ILS
iasi.ils ! IASI ILS file
*SRF Section (#8)
*FOV Section (#9) (Limb Only)
Field Type Description
FOVFIL C200
Name of .fov file
*FOV
mipas_5pt.fov ! 5 pt representation of MIPAS FOV
*ACC Section (Optional)
Field Type Description Default
NGJACC=
NGJACC
I
Number of gases for which Jacobians are computed
No. rtvd species
PGCACC=
PGCACC
R
Fractional change in pressure requiring new C-G path calc
0.01
TGCACC=
TGCACC
R
Change in temperature [K] requiring new C-G path calc
0.1
JCPACC=
JCPACC
R
Fraction of Jacobian perturbation requiring additional path
0.01
JCCACC=
JCCACC
R
Fraction of Jacobian perturbation requiring new C-G path calc
0.5
*ACC
NGJACC=4 TGCACC=0.2
*ASD Section (Optional)
Field
Type
Description Units Range Default
CORREL=
CORREL R
Correlation length [km] 0:100 50km
TEMASD=
TEMASD R
Temperature SD [K] 0.001:1000 10K
SFCTEMASD=
SFCTEMASD R
Surface Temperature SD [K] 0.001:1000 10K
PREASD=
PREASD R
Pressure SD [%] 0.001:1000 50%
[gas]ASD=
VMRASD R
VMR SD [%] 0.001:1000 100%
EXTASD=
EXTASD R
Extinction SD [km-1] 0.001:1000 0.005
POIASD=
POIASD R
Pointing SD [km] 0.001:1000 0.15
*ASD
CORREL=0 ! Uncorrelated A Priori
CH4=10 ! Change CH4 to 10% uncertainty
*CLD Section (Optional)
Field
Type Description Default
PIXCLD=PIXCLD
R
Maximum Pixel Cloud percentage (0:100) 0.0 [%]
Field
Type Description Default
CIXCLD=CIXCLD
R
Maximum Cloud Index (-1 or positive) 1.8
RADCLD=RADCLD
R
Minimum Cloud Radiance (-1 or positive) 125.0
[nW/(cm2.sr.cm-1)]
TOPCLD=TOPCLD
R
Maximum Cloud Altitude (-1 or positive) 30.0 [km]
*CLD
RADCLD=-1 ! Remove Cloud Radiance threshold test
CIXCLD=2.2 ! Set more stringent Cloud Index threshold
*CNV Section (Optional)
Record Type
Description Default
MAXITR=
MAXITR
I
Maximum number of iterations (GE 1) 10
CHILIM=
CHILIM
R
Minimum value of ChiSq for convergence 1.0
DELCHI=
DELCHI
R
Minimum Change in ChiSq for convergence 0.1
GAMMAX=
GAMMAX
R
Maximum value for Levenberg-Marquardt parameter 1000.0
*CNV
CHILIM=2.0
MAXITR=5
*GRD Section (Optional)
Field Type Description
GRDFIL C200
Name of irregular grid file
*HIT Section (Optional)
Field Type Description
HITFIL C200
Name of binary spectroscopic
database file
*HIT
../morse_files/hitran_mipas.bin
*LUT Section (Optional)
Field Type Description
LUTFIL C200
Name of LUT file
*LUT
[need to add some usable examples]
*NTE Section (Optional)
Field Type Description
NTEFIL C200
Name of
.nte file
*NTE
co2_day_200km.nte ! mid-lat daytime vib.temps
*OUT Section (Optional)
Field Type Description Default
OUTDIR=OUTDIR C193
Output directory
[local directory]
ATMFIL=ATMFIL C193
Retrieved Atmosphere
no output
DIAFIL=DIAFIL C193
Diagnostics
no output
RESFIL=RESFIL C193
Residual Spectra
no output
RSMFIL=RSMFIL C193
Resume file
morse.rsm
RTVFIL=RTVFIL C193
Retrieval Vector
morse.rtv
SWPFIL=SWPFIL C193
Sweep Diagnostics
no output
*OUT
RTVFIL=morse_#ORB.rtv
ATMFIL=morse_#ORB.atm
outdir=../Morse_Outputs/
*PIX Section (Optional)
NPIX = number of pixels in L1C file
Field Type Description Range
IPIX1 I
1st pixel to be retrieved 1:NPIX
[IPIX2] I
(Optional) last pixel to be retrieved IPIX1:NPIX
*PIX
1 10 ! Process just the first 10 pixels
*SHP Section (Optional)
Field Type Description
SHPNAM
C3 Lineshape to be used
GAS
C7 Molecule to which lineshape is applied,
or '*'
Field Description
VOI
Voigt Lineshape
LOR
Lorentz Lineshape
DOP
Doppler Lineshape
VVW or
VAN
Van Vleck Weisskopf Lineshape
CHI
Voigt plus Chi-factor
*SHP
Chi-factor CO2 ! Use CO2 sub-Lorentzian wings
LORENTZ H2O CH4 ! Use Lorentz shape for H2O and CH4
LORENTZ N2O ! Also N2O
VOI * ! Use Voigt shape for all other absorbers [default anyway]
DOPP CO ! Doppler shape for CO
*SVD Section (Optional)
Field Type Description
SVDFIL C200
Name of SVD-compressed LUT file
CS*DAT
*LUT
./examples/CS_PT__0511_01.DAT ! pT MW#501 H2O LUT
./examples/CS_PT__0511_02.DAT ! pT MW#501 CO2 LUT
CS_*DAT ! all other PT MWs & absorber LUTs
*XSC Section (Optional)
Field Type Description
XSCFIL C200
Filenames for .xsc files
*XSC
./examples/aerosol.xsc ! just aerosol.xsc in current directory
../xscfiles/*.xsc ! directory for all other .xsc files