70 likes | 221 Views
sfit. Refmod maker Binput to sfit4.ctl converter update Output levels Comment line / version / run date / # reference 2 – every thing currently except very redundant data Better read checks add tag name Check units Check normalization – write out scale factors
E N D
sfit • Refmod maker • Binput to sfit4.ctl converter update • Output levels • Comment line / version / run date / # reference • 2 – every thing currently except very redundant data • Better read checks add tag name • Check units • Check normalization – write out scale factors • 1 – basic output, detail file, statevec, spectra, AK • 0 – just for operational runs – binary format enough to build hdf (default) • Ability to recalculate spectrum from previous retrieved state • Error analysis mode • New input block – a posteriori • Entry for each kb ?? Or all we can? • Separate snr ad-hoc versus real? • Put snr in t15asc for each spectrum • Allow different resolutions for a given µw
Possible Components Already in K • BckGrdSlp - background slope • BckGrdCur - background curvature • DWNumShft - differential wavenumber shift • SWNumShft - single wavenumber shift • IWNumShft - independent wavenumber shift • EmpApdFcn - empirical apodization function • EmpPhsFnc - empirical phase function • SolLnShft - solar line shift • Sollnstrn - solar line strength • SPhsErr - simple phase error • ZeroLev - zero level offset • Interfering species - n other fit species • Channel - 4 parameters / beam IRWG/NORS Error Analysis Workshop, NCAR, January 2013
Components Requiring new Kb • not fitted: • Target gas absorption line(s) intensity • Target gas absorption line(s) air broadened line half-width • Retrieved gas absorption line(s) intensity • Retrieved gas absorption line(s) air broadened line half-width • Other possible line parameters: • Self broadening • Lower state energy • Temperature dependence exponent • pressure shift • Galatry Beta • Line mixing , SDV params: GAMMA2, ETA2, LMTK1, LMTK2, YLM • Solar pointing • Temperature • Kb from linefit output IRWG/NORS Error Analysis Workshop, NCAR, January 2013
Sb prescriptions defined by physical error source 1/ • Pointing • Time or ephemeris calculation (at racking and/or at synthetic spec. calc.) • Bias on (quad) sensors • SZA change during measurement integration time • ? • ILS both phase & modulation • Error along scan path • How many points or order of polynomial? (4th order) • Errors from linfit output: cell p, T • instrument variability 120/M > 125HR… • Solar intensity variation during course of scan • Temperature • Error by layer • Some inter-layer correlation? If so what? From ensemble / NCEP • ? • Off axis values for Sb should be considered
Sb prescriptions defined by physical error source 2/ • background slope • background curvature • differential wavenumber shift • single wavenumber shift • independent wavenumber shift • solar line shift • solar line strength • simple phase error • Separate from ILS phase ? • zero level offset • Fitted only
Documentation • SFIT to hdf work flow - Bavo • Ncar batch outline • Sfit4 input file key • Sfit4.ctl – mostly done • Isotope.input – done • Hbin.input - done • Ils.dat - done but old • Wvmr.input - Reference.prf / waccm / zpt • Sfit4 output keys • Nothing yet – start next week • Sfit 4 theory to implementation • How do options work • Sb definition / outline
Actions 1 Feb 2013 • Kb from sfit – Jim, Mathias, Nicholas • Sb outline – Martine to start • Expand and pass to group - Jim • error calculation code – Stephanie • refmaker.f90 - Jim, Mathias • Sfit4 workflow - Bavo • Move codes to IRWG ws – Jim • Sfit4 input/output key - Mathias • Sfit4 implementation doc – Jim, • New batching code - Hideaki, Nicholas, Jim • Input definition (control, spec list) - everyone • Output definition (hdf, list) - everyone • Task manager - Bavo • Specmaker- runlog or ckopus/prepspsec • Provide testcases (3+ NDACC) - everyone • Speclistmaker - after runlog