-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathSCIA2D_mlt.conf
340 lines (340 loc) · 8.74 KB
/
SCIA2D_mlt.conf
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
!*****************************************************************************!
! !
! Main configuration file for program SCIA2D !
! !
!*****************************************************************************!
!
!***************************** Directory structure ***************************!
!
! Number of emitting species (must be <=10)
Number of emitters
5
!
! Path to SCIA solar spectra [SolPath]
! Note that currently this setting has no effect because we overwrite it
! with the value from the commend line.
Solar spectra
auto
!
! Relative path to fallback solar spectrum
! that is used when a current spectrum is not available
! File may be zipped, but provide the file name WITHOUT ".gz" suffix here
Fallback solar spectrum
DATA/SCIA_solar_20100218_001751_41668.dat
!
! Correction factors to the reference solar spectrum
! can contain an actual filename or "auto" for automatic
! calculation per orbit.
solar correction factors
DATA/sol_nomfac2002-2012_mfac_2_238-282nm_daily.dat
!
! Line parameters of retrieval species
Line parameters
DATA/LinePars.dat
!
! Atmospheric structure
! (species in columns, altitude grid in column 1,
! first line has to give the number of altitudes and the number of species)
Atmospheric structure
DATA/ACE_air_o3_50_200.dat
!
! Absorption cross sections
! (species in columns, wavelength grid in column 1,
! first line has to give the number of wavelengths and the number of species)
Cross sections
DATA/XSECT_air_o3_new.dat
!
! NO parameters: oscillator strengths, Einstein coefficients,...
NO parameters
DATA/Luqueetal.dat
!
! Ap index file
Ap index file
DATA/spidr_ap_2000-2012.dat
!
! Kp index file
Kp index file
DATA/spidr_kp_2000-2012.dat
!
! f10.7 index file
f10.7 index file
DATA/f107_noontime_flux_obs.txt
!
! f10.7a index file
f10.7a index file
DATA/f107a_noontime_flux_obs.txt
!
! f10.7 adjusted flux (to 1 AU) index file
f10.7 adjusted index file
DATA/spidr_f107_2000-2012.dat
!
! Wavelengths to calculate the absorption along rays
! Note that the order of values below must correspond to the order of species
! in line parameter file (see above)
Absorption wavelengths
280.00 282.00
!
!***************************** Input data ************************************!
!
! Relative path to file that contains the full path of limb scan to be processed
! (e.g. C:\SCIA\NADIR\SCIA_nadir_20060906_115958_23624.dat )
Scan ID data path
Testdaten/orbit12908.txt
!
! Relative path to correction factors (note input format!)
Correction factors
KeyData/corrections.limbrad.irrad.spectralon
!
!***************************** Dark signal background value*******************!
!
! Dark signal estimate, can be a fixed positive value (e.g. empirical avarage
! dark signal) or negative:
! -1: use spectral median within [238.0, 282.0]
! -2: lowess smoothed spectrum (not implemented yet)
! -3: whittaker smoothed spectrum (not implemented yet)
!
dark background
3.9e9
!
!***************************** Altitude/retrieval grid ***********************!
!
! Minimal and maximal retrieval altitude [MinAlt,MaxAlt]
MinAlt and MaxAlt
60. 160.
!
! Altitude layer size [km]
dAlt
2.
!
! Minimal and maximal latitude [MinLat,MaxLat]
MinLat and MaxLat
-90 90
! number of horizontal grid points
NLat
72
!
!
! Altitude grid
! The measurement altitude grid can be extended by a number of altitude bins
! Note that the extension bins have to be located above the top
! measurement altitude to avoid numerical discontinuities
! Lower limit of extension bin no. 1 is upper limit of last measurement bin
!
! Extension format: No. of bins
! Lower/upper limits of bins [km]
!
Altitude grid extensions
0
0. 120.
120. 150.
150. 250.
250. 500.
!
! Altitude above which absorption is neglected
TOA
200.
!
! Altitude below which we have only absorption
BOA
50.
!
! Altitude of the highest tangent point considered
! In general equal to TOA, change only for testing purposes
max TP alt
200.
!
! Altitude of the lowest tangent point considered
! In general equal to BOA, change only for testing purposes
min TP alt
50.
!
!******************************* Selection rules *****************************!
!
! Do a pure nadir retrieval without inclusion of limb measurements
Nadir retrieval
0
!
! Sort out certain measurements (1) or not (0) [iNight,iGEO,iSZA,iNLC]
Night
1
Geolocation
0
Large SZA
1
NLC
0
Skip SAA
1
!
! Maximal SZA to be considered [MaxSZA, deg]
Maximal SZA
95.000
!
! SAA cut-off intensity
SAA cut-off value
8.8e10
!
! Geolocation boundaries to be sorted out (LonMin,LonMax,LatMin,LatMax) [deg]
Geolocation boundaries
0. 0. 0. 0.
!
!***************************** Baseline fit parameters ***********************!
!
! Number of spectral windows for baseline fit (O3, AIR)
! Spectral ranges for background fit [nm,nm]
Baseline windows
3
270. 278.
282. 284.
286. 290.
!
!************************ Retrieval species fit parameters *******************!
!
! Number of retrieval windows
! Spectral intervals for retrieval species [nm,nm]
! WARNING: Specify one window for each retrieval species
Retrieval windows
1
270. 290.
!
! Specify which wavelength windows refers to which species
! Format: [species of window 1] [species of window 2] ...
! Note: Numbering of species follows the order of line specifications
! (e.g. in LinePars.dat)
Assignment of wl windows
1 2
!
!******************************* Regularization ******************************!
!
! Retrieval covariances
! (rows: a priori, latitudinal, altitudinal)
! (columns: species 1, species 2,...)
! original values for MgII (species 1): 1.e-6, 1.e-5, 5.e-7 (3 km) (1.e-5 (1 km))
! original values for MgI (species 2): 5.e-2, 1.e-5, 1.e-5
! original values for unknown (species 3): 1.e-14, 5.e-6, 5.e-6
! original values for FeI (species 4): 1.e-10, 1.e-8, 1.e-8
! original values for NO (species 5): 3.e-17, 3.e-16, 3.e-16
Retrieval covariances
3.D-18 3.D-8 1.D-14 1.D-10 1.D-17
3.D-8 3.D-6 5.D-6 1.D-8 1.D-16
3.D-8 3.D-6 5.D-6 1.D-8 3.D-17
!
!***************************** Miscellaneous *********************************!
!
! Thresholds for rejection of measurements (Nadir, Limb) (0.259 0.879)
! To be given are thresholds for ||yfit-ymeas|| / ||ymeas||
Error thresholds
1.E20 1.E20
!
! FWHM of spectral slit function [FWHM, nm]
Spectral FWHM
.22
!
! Apply correction factors for solar spectrum and limb radiances
Do correction of radiances
0
!
! Operating system ('WINDOWS' or 'LINUX', provide left-justified!)
OS type
LINUX
!
! Maximal number of Levenberg-Marquardt steps
! ( A priori covariance is incresed by scaling factor in each step)
Maximal number of LM steps and scaling factor
20 2.0
!
! Maximal number of iteration steps in improving the solution
! Iterations are aborted if maximal number of steps is reached or
! residual < conv_threshold * (initial residual)
Maximal number of iteration steps and convergence threshold
2000 .01
!
! The temperature to use when calculating the
! NO emission lines [K]
atmosphere temperature
200.
!
number of NO transitions
3
!
! The NO transitions to fit
! v_u v_l v_l_abs
! (default: 0 2 0)
NO transitions
0 2 0
1 4 0
1 5 0
0 3 0
0 1 0
!
! Rough correction of the polarisation effect
! on the NO rotational lines
NO polarisation correction
1
!
! NO Rayleigh background fit method
! 0: None
! 1: per NO gamma band (standard up to v5.0)
! 2: fixed Rayleigh fit window (parameters below)
NO Rayleigh fit method
1
!
! NO Rayleigh background fit window
! for case 2 above (min, max)
NO Rayleigh fit window
238 282
!
! NO apriori
! 0: null apriori
! 1: NOEM apriori
! 2: Regression apriori
NO apriori
0
!
! NO apriori altitude range
! zero below and above
! actual values between
NO apriori altitude range
40.0 160.0
!
! NO apriori scale factor
! use value times scale factor as apriori value
! a negative value switches on "automatic" scaling:
! -1: scale = f10.7_adj(day - 1) / 150
! -2: fit scale factor
! -3: scale according to the slant columns from one forward run
NO apriori scale
1.0
!
! NO apriori altitude smoothing scale
! Smooth the apriori xx km along the vertical direction
NO apriori smoothness
6.0
!
! NO apriori relative covariance scaling
! Scale the NO apriori covariance proportional
! to the NO apriori value (0 = no, 1 = yes)
! (default: 0)
NO apriori covariance relative
0
!
! NO apriori covariance factor
! Divide the NO apriori covariance by the given factor
! (= less weight) where the NO apriori values is not zero.
! (default: 1.0)
NO apriori covariance factor
1.0
!
!
! Number of the retrieval algorithm
! 0: old variant
! 1: new variant
! 2: new variant using Eigen3
Retrieval algorithm
1
!
! Use MLT scans for the retrieval
! 0: false (nominal scans)
! 1: true (MLT scans)
MLT
1