Microsoft word - sfit4_init.doc

Keyword1.keyword2.keyword3 = value
First keyword:file defines in and output files. All output files are also set to default values by the code.gas defines parameters for the retrieval gasesfw all parameters for the forward model are defined in this sectionkb is set to true if Kb matrizes are calculated at the end of the talk.rt in this section the retrieval parameters are defined, except for the SA matrices for the gases which are defined in the gas section. If rt is set to F, no retrieval is performed but only a forward calculation.
In the section band the parameters for al MW bands are defined. sp contains additional noise information for the spectra, i.e. for deweighting
In out more details of output are definedRecommendation: Do not rely on default values for in and output to make the program forward compatible Dependency Default Description
file.in.modulation_fcn fw.modulation_fcn = 4
file containing spectral data, created with hbin
file.out.solarspectrum fw.solar_spectrum=T
Nr of layers the gas is retrieved on. Must match the number of layers in file.statlayers, in file.refprofile and file.isotope
names of the gases for which profiles are retrieved
names of the gases for which columns are retrieved
gas.profile.x.correlati gas.profile.x.correlation
correlation in the sa-matrix 1 – Gaussian shape 2 – Exponential shape 4 - the sa matrix is read in from file.sa_matrix 5 – the inverse Sa matrix is read from file.sa_matrix
gas.profile.x.correlati gas.profile.x.correlation
= T gas.profile.x.correlation.type =1,2
gas.profile.x.correlati gas.profile.x.correlation
gas.profile.x.correlati gas.profile.x.correlation
diagonals of Sa matrix in fractions of the a priori Nr of entries must correspond to the number of layers defined in the statlayers
factor the entries in file.sa_matrix are multiplied
Half width of integration interval for crosssection calculation
Lineshape model 0 – depends on the spectroscopic values given 2 - Galatry if BETA_T is given, if not Voigt 3 - SDV, if parameters are given
if T and parameters found, linemixing is included
currently only CO2 possible, gas for which linemixing is calculated
if T inclusion of solar lines (files.solarlines)
Pressure induced line shift 0 - read from linelist 1 - no shift
0 - no empirical apodiziation 1 - tabular function is read in 2 - polynomial 3 - fourier series 4 - linfit output is read in
Empirical phase error 0 - no empirical phase 1 - tabular function is read in 2 - polynomial 4 - linefit output
if T emitted radiation from the atmosphere is calculated
Temperatur (in K) of the radiating object outside the atmosphere Moon = 370.0 Sun = 6000.0 None = 2.7
Reflexion of solar light off object .e. only emission is calculated, no reflection .m. reflection of solar ligh of the moon (pre-alpha)
spectra are normalized to one (T) or not normalized (F)
T if Kb matrix entries are calculated, if the respective statevector entries are not retrieved, i.e. given kb.slope = T
the Kb row for the slopes are only calculated it slope is not retrieved.
Calculates AB matrix for a wrong assumed profile if the retrieved gas is a column
For which gas an error for the retreived profile is calculated?
T if Kb calculation for line intensities
for which gases lineparameters are calculated, default: all gases which are retrieved predefined values: target – calculation only for the target gas retreival – kb are calculated for each gas which is retrieved.
1 if all lineparameters of a gas are perturnbed together with the same perturbation (this is the only parameter supported so far) Kb are calculated for:
Temperature depency of pressure braodening
Switch on (T) or off (F) Retrieval, if F only a forward model calculation is performed
Increase gamma by value if step was succesful
convergence is reached if change in cost function is smaller than value
Convergence is reached if the proposed change in the spectrum is smaller than the noise * rt.tolerance
0 - no shift for any bandpass 1 - single shift for each bandpass 2 - independent shift for each bandpass 3 - idependent shift for each fit
apriori of all types of wavenumber shift in [cm^-1]
Empirical phase function retrieved if T
Empirical phase function retrieved if T
if T retrieval of line shifts for each retrieved gas
diagonals of sa matrix for temperature for each layer in state vector
= 1.2 MWs that are included in the calculation
Imposed apodization code 0 – Boxcar 1 -3 Norton Beer 4 – Triangle 5 - Happ – Genzel 6 – KPNO Atmospheric Spectra 7 – Hamming function
T if an offset is retrieved in this band
0 - use the a priori as given 1 - allow to retrieve for each bad 2 - use zero level from first band in list
Its standart deviation, if set to o.o not retrieved but fixed
Channel model PS - phase model IP - interferogram pertubation model
gases which are retrieved from this band, must be contained in gas (see above)
T if temperature is retrieved in this band
which additional snr windows are taken into account e.g. = 1, the lines containing an 1 are read in, all other lines are ignored
Output level, a predefined set of putput files
Type of GASFILE 1 - only the final spectrum, the spectrum of each gas and the solar spectrum will be printed out.
The files are named like gas1.1.1, allgases.1.1 and solar 1.1 2 - spectra will be printed out for each iteration. The names are like above,but an extra number is appended denoting the iteration number the numbers appended to the files are "nr of window", "nr of scan" and "nr of iteration" The information of for band nr, gas and iteration number are also contained in the file header
Additionally to the predefined output acc to the level given in <output> the following quantities can be written out. If such a key is given the resp quantity is written out to the file defined by the parameter string (e.g. output.k-matrix = kk.out – the Kmatrix is written out to kk.out)
0,1,2 for the level of the output (see below)
<filename> matrices written in file (now only: K.out)
<filename> write out SA-matrix (now only: SA.out)
<filename> write out SM-matrix (error on profile due to the measurement noise (now only: SM.out)
T calculated spectra for each retrieved gas and each band are printed out
Output description
Files which may appear but are not described here are a legacy and are subject to modification or removal in the future, so dont relay on them, but notify the maintainers of sfit4 if you need the information contained in them.
contains al profiles of the retrieval gases together with the alitude grid, pressure , temperature and airmass (vertical). For each gas there are five columns:
contains the retrieved, measured spectra and the difference thereof
Averaging kernel in units of the internal statevector, i.e. x/xa
Contains the GyKb matrix (see formula 3.16 page 48 in Rodgers (2000)
contains the Kb matrix for al parameters which are not retrieved (and contained in the K-­‐matrix)
contains the ful Sє matrix of the retrieval noise
contains the spectra calculated for each retrieval gas, each iteration (if output.gas_spectra.type = 2) and each microwindow