Mandatory sections | ||
---|---|---|
*HDR | 1. Comment record written to output file headers | |
*FLG | 2. Option flags | |
*L1C | 3. Input spectra | |
*RTV | 4. Target quantities to be retrieved | |
*[LEV] | 5. Retrieval levels (*ALT, *PRE, *SWP or *TAN) | |
*MWL | 6. Microwindow list files | |
*MDB | 7. Microwindow Database files | |
*CLI | 8. Climatological profiles | |
*ILS | 9. Instrument Line Shape files | |
*FOV | 10. Field of View Shape (omitted if nadir-viewing) | |
Optional sections | ||
*ACC | Accuracy settings | |
*ACV | A Priori Covariance | |
*ASD | A Priori Standard Deviation | |
*CLD | Cloud Detection Criteria | |
*CNV | Convergence Criteria | |
*FMZ | Forward Model Altitudes | |
*HIT | HITRAN Database | |
*L2P | Level 2 Input Profiles | |
*LUT | Directories for Look Up Tables | |
*NTE | Vibrational Temperatures for non-LTE forward model | |
*OUT | Rename Output Files | |
*SHP | Specify line shapes for line-by-line calculations | |
*XSC | Directories for Molecular Cross Sections | |
Termination Record | ||
*END |
*ABCwhere 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.
The driver table is terminated with
*END(remember to add a <CR> after the *END)
The first 10 sections (*HDR ... *FOV) are mandatory and the sequence is fixed (or 9 sections for nadir-viewing geometry since FOV is not specified for nadir-viewing). 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
The rest of this document describes the contents
of each section.
*HDR Section (#1)
Field Type Description
HDRREC
C*79 Text for output file
headers
*HDR
O3 Rtvl generated with morse.drv 14NOV04 ! [and this part won't be written]
*FLG Section (#2)
Field Type Value Description
FLG C*3
ACC
More stringent default accuracy/convergence criteria
New in v3
AHY Suppress hydrostatic equilibrium
Only with PRE retrieval
APR Output A Priori information
(morse.apr)
CLD Use cloud detection
COV Output covariance matrix
(morse.cov)
CTM Use molecular (eg H2O) continuum data
DIA Output iteration diagnostics
(morse.dia)
FGD Use full, as opposed to irregular, grid
FIX Assume fixed relative altitudes
(eg detector array)
Only with PRE retrieval New in v3
FMW use full microwindow (ie set all masks = TRUE)
New in v3
FVZ Suppress FOV distortion due to refraction
GEO Use geometric (ie non-refractive) ray paths
LIN Assume VMR varies linearly with
altitude rather than logarithmically
MWO Write
outputs for each microwindow
QAD Use simple quadratic fit to line wings
rather than inverse quadratic
RES Output residual spectra
(morse.res)
*FLG
CLD APR RES ! Use cloud-detection, generate A Priori and Residuals files
*L1C Section (#3)
Field Type Description Requirement
1. L1CFIL C*80
Name of L1C file
Mandatory
2. SWPFIL C*80
Name of .swp file
Optional
*L1C
L1C_02081_SCAN#13 ! Scan#13, orbit 2081
../pt/morse.swp ! Modified Sweep altitudes from pT retrieval output
*RTV Section (#4)
Field Type Value Description Units Restrictions
RTV C*10
TEM Temperature [K]
PRE Pressure [mb]
Not on *PRE levels
CTM Continuum [/km]
AEROSOL Aerosol [/km]
SFT Surface Temperature
[K]
OFF Radiometric Offset
[nW/(cm2.sr.cm-1)]
[gas] Molecule VMR
[ppmv] [gas] from List of gases
[gas]([iso]) Isotopomeric VMR
[ppmv] [iso] from
List of isotopes
[gas]([iso])([vib])
Vibrational Temperature
[K] [vib] is HITRAN Global Quantum Number (1=ground state)
[gas](lin) Linear Retrieval
[ppmv] v13AUG07 onwards
*RTV
TEM PRE CTM O3 H2O(1) H2O(2) ! Joint pT, O3, ctm and main two H2O isotopes
*ALT, *PRE, *SWP or *TAN Section (#5)
*LEV Mandatory
Section Header defines grid type
GRD(i) Mandatory
Repeat for i=1:ngrd retrieval grid levels
SUBSET=[rtv(j)] Optional
Start Repeat for j=1:nrtv retrieved parameters
IDXSUB(k,j) Optional
Repeat for k subset levels of parameter j
End Repeat for j
REFLEV=[rlev] Optional Insert at any point
REFSWP=[rswp] Optional Insert at any point
Field Type Value Description
*LEV C*4
*ALT Altitudes [km]
*PRE Pressures [mb]
*SWP Tangent points,
specified by Sweep indices
*TAN Tangent points,
specified by Nominal altitudes [km]
GRD R*4
Retrieval Grid levels (*ALT, *PRE, *TAN)
I Retrieval Grid Sweep indices (*SWP)
C*1 *
Shorthand to indicate all L1B tangent points/sweeps
v10may07
SUBSET= C*8
[rtv] Define subset for retrieved parameter [rtv]
IDXSUB I
Index of GRD for retrieval subset level
REFLEV= I
[rlev] Index of retrieval level used as reference
REFSWP= I
[rswp] Index of sweep used as reference
*TAN
12 15 18 21 24 27 30 33 36 39 42 47 52 60 68 ! Retrieve 12-68km tan.pts
76 ! plus extra level at 76km
SUBSET=SF6
12 15 18 30 42 ! Retrieve SF6 on fewer levels
REFLEV=3 ! Use nominal 18km sweep as reference level
SUBSET=CTM
12 15 18 21 ! Only retrieve continuum to 21km
*MWL Section (#6)
Field Type Description
FILMWL C*80
Names of .lst files
*MWL
pt_000.lst ! pT microwindows
h2o_000.lst ! plus H2O microwindows for a joint retrieval with pT
*MDB Section (#7)
Field Type Description
MDBFIL C*80
Names of MW*DAT files
*MDB
MW_PT__103.DAT MW_H2O_103.DAT ! pT & H2O microwindows
* ! construct any other MWs from .lst file
*CLI Section (#8)
Field Type Description
CLIFIL C*80
Names of .atm files
*CLI
day.atm ! Mid-latitude dayime profiles
aerosol.atm ! zero extinction aerosol profile
h2o.atm ! H2O retrieved from previous scan (so different alt.grid)
*ILS Section (#9)
Field Type Description
ILSFIL C*80
Names of .ils files
*ILS
mipas_a.ils ! MIPAS A band AILS
mipas_ab.ils ! MIPAS AB band AILS
nbstrong.ils ! Idealised Norton-Beer Strong ILS for other bands
*FOV Section (#10)
Field Type Description
ILSFIL C*80
Name of .fov file
*FOV
rfm_1km5.fov ! 5 pt representation of MIPAS FOV @1.5km spacing
*ACC Section (Optional)
Field Type Description Default
NGJACC=
I
Number of gases for which Jacobians are computed
No. rtvd species
PGCACC= R*4
Fractional change in pressure requiring new C-G path calc
0.01
TGCACC= R*4
Change in temperature [K] requiring new C-G path calc
0.1
JCPACC= R*4
Fraction of Jacobian perturbation requiring additional path
0.01
JCCACC= R*4
Fraction of Jacobian perturbation requiring new C-G path calc
0.5
*ACC
NGJACC=4 TGCACC=0.2
*ACV Section (Optional)
(optional, default=0)
Field Type Description Default
ACVFIL
C*80
A Priori Covariance Matrix
RELAX=
R*4
Relaxation parameter
0
SA' = SA*EXP(-R) + SC*(1-EXP(-R))
where SA' is the A Priori covariance actually used to
start of the retrieval,
SA is the covariance read from the file specified in this section,
and SC is the climatological covariance constructed in the same way as
the default a priori (including modifications in the
*ASD section). A value RELAX=0.0 (which is
the default) means that the specified A Priori covariance is used
exactly as supplied, a value RELAX=0.1 implies relaxing to
climatological covariance with a 1/e time constant equivalent to
10 scans, 0.01 is a hundred scans, etc.
*ACV
morse.cov ! Use previous output covariance
relax=0.05 ! Climiatological relaxation length: 20 profiles, or ~1/4 orbit
*ASD Section (Optional)
Field Type Description Units Range Default
CORREL= R*4
Correlation length [km] 0:100 50km
TEM= R*4
Temperature SD [K] 0.1:1000 10K R*4
PRE= R*4
Pressure SD [%] 0.1:1000 50%
[gas]= R*4
VMR SD [%] 0.1:1000 100%
*ASD
CORREL=0 ! Uncorrelated A Priori
CH4=10 ! Change CH4 to 10% uncertainty
*CLD Section (Optional)
Field Type Description Default
IDXCLD= R*4
Maximum Cloud Index (-1 or positive) 1.8
RADCLD= R*4
Minimum Cloud Radiance (-1 or positive) 125.0
[nW/(cm2.sr.cm-1)]
TOPCLD= R*4
Maximum Cloud Altitude (-1 or positive) 30.0 [km]
*CLD
RADCLD=-1 ! Remove Cloud Radiance threshold test
CLDIDX=2.2 ! Set more stringent Cloud Index threshold
*CNV Section (Optional)
Field Type Description Default
MAXITR=
I
Maximum number of iterations (GE 1) 10
CHILIM= R*4
Minimum value of ChiSq for convergence 1.0
DELCHI= R*4
Minimum Change in ChiSq for convergence 0.1
*CNV
CHILIM=2.0 MAXITR=5
*FMZ Section (Optional)
Field Type Description
FMZ R*4
Nominal Foward Model Altitudes [km] (increasing monotonically)
*FMZ
3 4 5 6 7 8 9 10 11 12 13 14 15 18 21 24 27 30 33 36 39 42
47 52 60 68 76 85 100 120
*HIT Section (Optional)
Field Type Description
HITFIL C*80
Name of binary spectroscopic
database file
*HIT
../morse_files/hitran_mipas.bin
*L2P Section (Optional)
Field Type Description
L2PFIL C*80
Name of .atm file
generated by MORSE for the same limb scan.
*L2P
morse_pt.atm ! Input retrieved pT profile
morse_h2o.atm ! Input retrieved H2O profile
*LUT Section (Optional)
Field Type Description
LUTDIR C*80
Directory to be searched for
CS*DAT files
*LUT
../morse_files/CS/pt/ ! pT LUTs
../morse_files/CS/o3/ ! O3 LUTs
*NTE Section (Optional)
Field Type Description
NTEFIL C*80
Name of .nte file
*NTE
co2_day_amil2da.nte ! Use mid-lat daytime vib.temps
h2o_day_amil2da.nte
o3_day_amil2da.nte
*OUT Section (Optional)
Field Type Description Default
DIR= C*76
Output directory
[local directory]
ACV= C*76
A Priori Covariance
morse.acv
APR= C*76
A Priori Vector
morse.apr
ATM= C*76
Retrieved Atmosphere
morse.atm
COV= C*76
Retrieval Covariance
morse.cov
DIA= C*76
Diagnostics
morse.dia
RES= C*76
Residual Spectra
morse.res
RTV= C*76
Retrieval Vector
morse.rtv
SWP= C*76
Sweep Diagnostics
morse.swp
*OUT
APR=morse_pt.apr
cov=./cov/morse_pt.cov
dir=../Morse_Outputs/
*SHP Section (Optional)
Field Type Value Description
[gas]= C*3
VOI Voigt line shape
CHI Chi Factor
LOR Lorentz line shape
DOP Doppler line shape
*SHP
co2=chi ! Chi factor for CO2
o3=doppler h2o=dop ! Doppler for O3 and H2O
*=lorentz ! Lorentz for everything else
*XSC Section (Optional)
Field Type Description
XSCDIR C*80
Directories for .xsc files
*XSC
./ ! just aerosol.xsc in current directory
../morse_files/ ! directory for all other .xsc files