MORSE .dia File

02APR24
Description
MORSE Diagnostics file

Structure
Text file. Formal structure tbd
File Structure
! Diagnostics file for MORSE v.VERSID
! HDRREC
FMT           (= .dia file format version)
NPTMAX     (= Max No.pts in any microwindow)
! -------------------------
IPIX   IMIC   ISWP   ITER   MICLAB   STATUS   = Pixel#, MW#, Sweep#, Iter#, MW label
CHI2X   CHI2Y   GAMMA   NX   NY   = Chi2X, Chi2Y, Gamma, NX, NY
NMEAS   NTRUE   OFFSET   NESR   ALT   = NMeas, Ntrue, Offset, NESR, Alt
Meas
MEAS(1) MEAS(2) ...
... MEAS(NMEAS)
Masks
MASK(1) MASK(2) ...
... MASK(NMEAS)
Y vec
Y(1) Y(2) ...
... Y(NY)
Y-F vec
DY(1) DY(2) ...
... DY(NY)
NRTV     = N.Rtv
For i = 1, NRTV Repeat for each retrieved product
    *RTV[i]
    NLVRTV[i] = N.Lev
! ix ilv iatm z[km] x a x_sd a_sd
    IX ILV[ix] IATM[ix] Z[ix] X[ix] A[ix] X_SD[ix] A_SD[ix]

Dictionary
FieldFormat Description Units Range
FMT F5.2 .dia file format identifier 0.01
HDRREC C79 Copied from *HDR section of the driver table
IX I4 Index in state vector 1:NX
IPIX I4 Pixel# 1:NPix
IMIC I4 Microwindow# 1:NMic
ISWP I4 Sweep# 1:NSwp
ITER I4 Iteration# 1:MAXITR
MICLAB 2X, C8 Microwindow Label
NLVRTV(:) * (NRTV) No. profile levels for species i > 0
NPTMAX I5 Max no. spectral points in any microwindow > 0
CHI2X E10.2 χ2 value for state vector ≥ 0
CHI2Y E10.2 χ2 value for measurement vector ≥ 0
GAMMA E10.2 Levenberg-Marquardt damping parameter ≥ 1
NX I5 No.elements in state vector ≥ 1
NY I5 No.elements in measurement vector ≥ 1
NMEAS I4 Total no. spectral points in microwindow NY
NTRUE I4 Total used spectral points in microwindow (=NY = ∑MASK) 1 ≥ : ≥ NY
OFFSET E10.2 Radiometric offset R.U.
NESR E10.2 Radiometric noise R.U. > 0
ALT E10.2 Sweep tangent altitude km
MEAS(:) E10.2 (NMEAS) Full list of measurements in Microwindow R.U.
MASK(:) I2 (NMEAS) List of logical masks (1=use meas.)
*RTV(:) Cn (NRTV) Retrieved species (preceded by '*')
DY(:) E10.2 (NY) Vector of (meas-forward model) differences R.U.
NRTV * No. different retrieval products > 0
STATUS 2X, C7 Status ∈{'ITERATE', 'CONVRGD', 'STOPPED'}
VERSID C11 MORSE version identifier eg 27APR23
Y(:) E10.2 (NY) Vector of measurements used in microwindow R.U.
Example
morse.dia
At the moment a bit of a mess with all sorts of possibly useful information output each iteration by subroutine wrtdia.for

Notes
  1. IDL procedure inpdia.pro may be useful in reading the MORSE diagnostics file into a structure.