Ticket #646: fwarrick_results.txt

File fwarrick_results.txt, 33.2 KB (added by warrick, 5 years ago)

Reproducing Chris' results at the Met Office

Line 
1###########################################################################
2run_1dvar-00
3
4---------------------------------------------------------------------
5 ROPP Bending Angle 1D-Var
6---------------------------------------------------------------------
7
8INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-00.cf.
9
10INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
11 ./data/metop-c/first.nc.
12
13INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
14 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
15
16INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
17
18INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
19
20INFO (from ropp_1dvar_solve): Using background error covariance matrix for preconditioning.
21
22INFO (from ropp_1dvar_cost): Using absolute decrease of cost function and relative change
23 of state vector as additional convergence criteria.
24
25 n_iter = 1 J = 97.039 max(relative change in state) = -
26 n_iter = 2 J = 79.446 max(relative change in state) = 0.42320
27INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the cost function is increasing
28
29INFO (from ropp_1dvar_cost): Minimization of cost function successfully finished (according to additional
30 convergence criteria). Number of required iterations: 3.
31
32INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
33 first-1dvar-default.nc.
34
35###########################################################################
36run_1dvar-01
37
38---------------------------------------------------------------------
39 ROPP Bending Angle 1D-Var
40---------------------------------------------------------------------
41
42INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-01.cf.
43
44INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
45 ./data/metop-c/first.nc.
46
47INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
48 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
49
50INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
51
52INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
53
54 n_iter = 1 J = 97.039 max(relative change in state) = -
55 n_iter = 2 J = 32.400 max(relative change in state) = 0.13905
56 n_iter = 3 J = 32.382 max(relative change in state) = 0.16073
57 n_iter = 4 J = 32.378 max(relative change in state) = 0.12803
58
59INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the cost function did not change by more
60 than 0.10000 for the last 2 iterations.
61
62INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
63 first-1dvar-default.nc.
64
65##############################################################################
66run_1dvar-02
67
68---------------------------------------------------------------------
69 ROPP Bending Angle 1D-Var
70---------------------------------------------------------------------
71
72INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-02.cf.
73
74INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
75 ./data/metop-c/first.nc.
76
77INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
78 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
79
80INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
81
82INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
83
84INFO (from ropp_1dvar_solve): Using background error covariance matrix for preconditioning.
85
86INFO (from ropp_1dvar_cost): Using absolute decrease of cost function and relative change
87 of state vector as additional convergence criteria.
88
89 n_iter = 1 J = 97.039 max(relative change in state) = -
90 n_iter = 2 J = 79.446 max(relative change in state) = 0.42320
91 n_iter = 3 J = 16747. max(relative change in state) = 71.630
92 n_iter = 4 J = 76.892 max(relative change in state) = 71.347
93 n_iter = 5 J = 76.486 max(relative change in state) = 0.11244E-01
94 n_iter = 6 J = 373.25 max(relative change in state) = 9.9224
95 n_iter = 7 J = 53.400 max(relative change in state) = 8.7370
96 n_iter = 8 J = 49.459 max(relative change in state) = 0.30230
97 n_iter = 9 J = 43.141 max(relative change in state) = 0.74605
98 n_iter = 10 J = 43.110 max(relative change in state) = 0.25562
99 n_iter = 11 J = 43.023 max(relative change in state) = 0.12781
100 n_iter = 12 J = 43.016 max(relative change in state) = 0.75320E-02
101 n_iter = 13 J = 42.914 max(relative change in state) = 0.76796E-01
102 n_iter = 14 J = 43.489 max(relative change in state) = 0.46498
103 n_iter = 15 J = 42.885 max(relative change in state) = 0.39930
104 n_iter = 16 J = 42.878 max(relative change in state) = 0.60095E-02
105 n_iter = 17 J = 42.769 max(relative change in state) = 0.97121E-01
106 n_iter = 18 J = 43.071 max(relative change in state) = 0.35938
107 n_iter = 19 J = 42.741 max(relative change in state) = 0.30358
108 n_iter = 20 J = 42.730 max(relative change in state) = 0.87931E-02
109 n_iter = 21 J = 42.565 max(relative change in state) = 0.11115
110 n_iter = 22 J = 42.476 max(relative change in state) = 0.10533
111 n_iter = 23 J = 42.462 max(relative change in state) = 0.40446E-01
112 n_iter = 24 J = 42.445 max(relative change in state) = 0.29540E-01
113 n_iter = 25 J = 42.358 max(relative change in state) = 0.52629E-01
114 n_iter = 26 J = 42.076 max(relative change in state) = 0.16068
115 n_iter = 27 J = 42.040 max(relative change in state) = 0.28554E-01
116 n_iter = 28 J = 41.725 max(relative change in state) = 0.16654
117 n_iter = 29 J = 40.846 max(relative change in state) = 0.12660
118 n_iter = 30 J = 39.286 max(relative change in state) = 0.75896
119 n_iter = 31 J = 38.406 max(relative change in state) = 0.59504
120 n_iter = 32 J = 38.350 max(relative change in state) = 0.10336
121 n_iter = 33 J = 38.342 max(relative change in state) = 0.28838E-01
122 n_iter = 34 J = 38.334 max(relative change in state) = 0.29887E-01
123 n_iter = 35 J = 38.322 max(relative change in state) = 0.22172E-01
124 n_iter = 36 J = 38.245 max(relative change in state) = 0.61646E-01
125
126INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the cost function did not change by more
127 than 0.10000 for the last 5 iterations.
128
129INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional
130 convergence criteria). Number of required iterations: 36.
131
132... (from ropp_1dvar_diagnostics): Warning: Too many iterations, niter = 36 2J/m = 0.361
133INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
134 first-1dvar-default.nc.
135
136##############################################################################
137run_1dvar-03
138
139---------------------------------------------------------------------
140 ROPP Bending Angle 1D-Var
141---------------------------------------------------------------------
142
143INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-03.cf.
144
145INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
146 ./data/metop-c/first.nc.
147
148INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
149 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
150
151INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
152
153INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
154
155INFO (from ropp_1dvar_solve): Using background error covariance matrix for preconditioning.
156
157INFO (from ropp_1dvar_cost): Using absolute decrease of cost function and relative change
158 of state vector as additional convergence criteria.
159
160 n_iter = 1 J = 97.039 max(relative change in state) = -
161 n_iter = 2 J = 79.446 max(relative change in state) = 0.42320
162 n_iter = 3 J = 16747. max(relative change in state) = 71.630
163 n_iter = 4 J = 76.892 max(relative change in state) = 71.347
164 n_iter = 5 J = 76.486 max(relative change in state) = 0.11244E-01
165 n_iter = 6 J = 373.25 max(relative change in state) = 9.9224
166 n_iter = 7 J = 53.400 max(relative change in state) = 8.7370
167 n_iter = 8 J = 49.459 max(relative change in state) = 0.30230
168 n_iter = 9 J = 43.141 max(relative change in state) = 0.74605
169 n_iter = 10 J = 43.110 max(relative change in state) = 0.25562
170 n_iter = 11 J = 43.023 max(relative change in state) = 0.12781
171 n_iter = 12 J = 43.016 max(relative change in state) = 0.75320E-02
172 n_iter = 13 J = 42.914 max(relative change in state) = 0.76796E-01
173 n_iter = 14 J = 43.489 max(relative change in state) = 0.46498
174 n_iter = 15 J = 42.885 max(relative change in state) = 0.39930
175 n_iter = 16 J = 42.878 max(relative change in state) = 0.60095E-02
176 n_iter = 17 J = 42.769 max(relative change in state) = 0.97121E-01
177 n_iter = 18 J = 43.071 max(relative change in state) = 0.35938
178 n_iter = 19 J = 42.741 max(relative change in state) = 0.30358
179 n_iter = 20 J = 42.730 max(relative change in state) = 0.87931E-02
180 n_iter = 21 J = 42.565 max(relative change in state) = 0.11115
181 n_iter = 22 J = 42.476 max(relative change in state) = 0.10533
182 n_iter = 23 J = 42.462 max(relative change in state) = 0.40446E-01
183 n_iter = 24 J = 42.445 max(relative change in state) = 0.29540E-01
184 n_iter = 25 J = 42.358 max(relative change in state) = 0.52629E-01
185 n_iter = 26 J = 42.076 max(relative change in state) = 0.16068
186 n_iter = 27 J = 42.040 max(relative change in state) = 0.28554E-01
187 n_iter = 28 J = 41.725 max(relative change in state) = 0.16654
188 n_iter = 29 J = 40.846 max(relative change in state) = 0.12660
189 n_iter = 30 J = 39.286 max(relative change in state) = 0.75896
190 n_iter = 31 J = 38.406 max(relative change in state) = 0.59504
191 n_iter = 32 J = 38.350 max(relative change in state) = 0.10336
192 n_iter = 33 J = 38.342 max(relative change in state) = 0.28838E-01
193 n_iter = 34 J = 38.334 max(relative change in state) = 0.29887E-01
194 n_iter = 35 J = 38.322 max(relative change in state) = 0.22172E-01
195 n_iter = 36 J = 38.245 max(relative change in state) = 0.61646E-01
196 n_iter = 37 J = 37.917 max(relative change in state) = 0.74321E-01
197 n_iter = 38 J = 50.348 max(relative change in state) = 2.0487
198 n_iter = 39 J = 37.722 max(relative change in state) = 1.8466
199 n_iter = 40 J = 37.700 max(relative change in state) = 0.86831E-02
200 n_iter = 41 J = 37.600 max(relative change in state) = 0.89404E-01
201 n_iter = 42 J = 317.92 max(relative change in state) = 9.3296
202 n_iter = 43 J = 37.579 max(relative change in state) = 9.2572
203 n_iter = 44 J = 37.576 max(relative change in state) = 0.34352E-03
204 n_iter = 45 J = 91.314 max(relative change in state) = 4.0968
205 n_iter = 46 J = 36.909 max(relative change in state) = 3.8664
206 n_iter = 47 J = 36.614 max(relative change in state) = 0.35645E-01
207 n_iter = 48 J = 35.339 max(relative change in state) = 0.35465
208 n_iter = 49 J = 3636.7 max(relative change in state) = 33.363
209 n_iter = 50 J = 35.247 max(relative change in state) = 33.198
210INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional
211 convergence criteria). Number of required iterations: 50.
212
213... (from ropp_1dvar_diagnostics): Warning: Too many iterations, niter = 50 2J/m = 0.333
214INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
215 first-1dvar-default.nc.
216
217##############################################################################
218run_1dvar-04
219
220---------------------------------------------------------------------
221 ROPP Bending Angle 1D-Var
222---------------------------------------------------------------------
223
224INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-04.cf.
225
226INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
227 ./data/metop-c/first.nc.
228
229INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
230 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
231
232INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
233
234INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
235
236 n_iter = 1 J = 97.039 max(relative change in state) = -
237 n_iter = 2 J = 32.400 max(relative change in state) = 0.13905
238 n_iter = 3 J = 32.382 max(relative change in state) = 0.16073
239 n_iter = 4 J = 32.378 max(relative change in state) = 0.12803
240 n_iter = 5 J = 32.376 max(relative change in state) = 0.12886
241 n_iter = 6 J = 32.378 max(relative change in state) = 0.12885
242 n_iter = 7 J = 32.416 max(relative change in state) = 0.12841
243 n_iter = 8 J = 32.482 max(relative change in state) = 0.98499E-01
244 n_iter = 9 J = 32.559 max(relative change in state) = 0.16935E-01
245 n_iter = 10 J = 32.593 max(relative change in state) = 0.80768E-03
246 n_iter = 11 J = 32.596 max(relative change in state) = 0.88798E-05
247 n_iter = 12 J = 32.596 max(relative change in state) = 0.88887E-07
248
249INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more
250 than 0.10000E-03 relative to the assumed background errors for the last 2 iterations.
251
252INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
253 first-1dvar-default.nc.
254
255##############################################################################
256run_1dvar-06 NumRec vn1
257
258---------------------------------------------------------------------
259 ROPP Bending Angle 1D-Var
260---------------------------------------------------------------------
261
262INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf.
263
264INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
265 ./data/metop-c/first.nc.
266
267INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
268 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
269
270INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
271
272INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
273
274 n_iter = 0 J = 97.039 max(relative change in state) = -
275 n_iter = 1 J = 32.394 max(relative change in state) = 7876.7
276 n_iter = 1 J = 32.560 lambda -> 0.10000E-03
277 n_iter = 1 J = 32.404 lambda -> 0.10000E-02
278 n_iter = 1 J = 32.469 lambda -> 0.10000E-01
279 n_iter = 1 J = 32.397 lambda -> 0.10000
280 n_iter = 2 J = 32.383 max(relative change in state) = 0.41798E-01
281 n_iter = 3 J = 32.366 max(relative change in state) = 0.34476E-01
282 n_iter = 3 J = 32.375 lambda -> 0.10000E-01
283 n_iter = 3 J = 32.368 lambda -> 0.10000
284 n_iter = 3 J = 32.367 lambda -> 1.0000
285 n_iter = 4 J = 32.366 max(relative change in state) = 0.32681E-03
286 n_iter = 4 J = 32.368 lambda -> 1.0000
287 n_iter = 5 J = 32.366 max(relative change in state) = 0.24915E-03
288 n_iter = 5 J = 32.368 lambda -> 1.0000
289 n_iter = 6 J = 32.366 max(relative change in state) = 0.24482E-03
290 n_iter = 6 J = 32.368 lambda -> 1.0000
291 n_iter = 7 J = 32.366 max(relative change in state) = 0.24052E-03
292 n_iter = 7 J = 32.368 lambda -> 1.0000
293 n_iter = 8 J = 32.365 max(relative change in state) = 0.23628E-03
294 n_iter = 8 J = 32.368 lambda -> 1.0000
295 n_iter = 9 J = 32.365 max(relative change in state) = 0.23211E-03
296 n_iter = 9 J = 32.368 lambda -> 1.0000
297 n_iter = 9 J = 32.366 lambda -> 10.000
298 n_iter = 10 J = 32.365 max(relative change in state) = 0.43513E-04
299 n_iter = 10 J = 32.366 lambda -> 10.000
300 n_iter = 11 J = 32.365 max(relative change in state) = 0.35401E-04
301
302INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more
303 than 0.10000E-03 relative to the assumed background errors for the last 2 iterations.
304
305INFO (from ropp_1dvar_levmarq): Finished after 11 iterations (26 forward model / gradient evaluations).
306
307INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
308 first-1dvar-default.nc.
309
310
311##############################################################################
312run_1dvar-06 TR1 vn1
313
314---------------------------------------------------------------------
315 ROPP Bending Angle 1D-Var
316---------------------------------------------------------------------
317
318INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf.
319
320INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
321 ./data/metop-c/first.nc.
322
323INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
324 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
325
326INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
327
328INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
329
330 n_iter = 0 J = 97.039 max(relative change in state) = -
331 n_iter = 1 J = 35.030 max(relative change in state) = 7876.7
332 n_iter = 2 J = 32.532 max(relative change in state) = 0.41053
333 n_iter = 3 J = 32.382 max(relative change in state) = 0.17524
334 n_iter = 4 J = 32.370 max(relative change in state) = 0.88772E-01
335 n_iter = 4 J = 32.379 mu -> 0.74074E-03
336 n_iter = 5 J = 32.367 max(relative change in state) = 0.16165E-01
337 n_iter = 5 J = 32.404 mu -> 0.49383E-03
338 n_iter = 6 J = 32.361 max(relative change in state) = 0.26131E-01
339 n_iter = 6 J = 32.372 mu -> 0.32922E-03
340 n_iter = 6 J = 32.364 mu -> 0.65844E-03
341 n_iter = 6 J = 32.369 mu -> 0.13169E-02
342 n_iter = 6 J = 32.366 mu -> 0.26337E-02
343 n_iter = 6 J = 32.366 mu -> 0.52675E-02
344 n_iter = 6 J = 32.365 mu -> 0.10535E-01
345 n_iter = 7 J = 32.360 max(relative change in state) = 0.12883E-01
346 n_iter = 7 J = 32.367 mu -> 0.70233E-02
347 n_iter = 7 J = 32.364 mu -> 0.14047E-01
348 n_iter = 7 J = 32.363 mu -> 0.28093E-01
349 n_iter = 7 J = 32.362 mu -> 0.56187E-01
350 n_iter = 8 J = 32.360 max(relative change in state) = 0.18178E-02
351 n_iter = 8 J = 32.362 mu -> 0.11237
352 n_iter = 8 J = 32.361 mu -> 0.22475
353 n_iter = 8 J = 32.361 mu -> 0.44949
354 n_iter = 8 J = 32.360 mu -> 0.89898
355 n_iter = 9 J = 32.360 max(relative change in state) = 0.18119E-03
356 n_iter = 9 J = 32.360 mu -> 3.5959
357 n_iter = 9 J = 32.360 mu -> 7.1919
358 n_iter = 9 J = 32.360 mu -> 14.384
359 n_iter = 9 J = 32.360 mu -> 28.768
360 n_iter = 10 J = 32.360 max(relative change in state) = 0.60608E-05
361 n_iter = 11 J = 32.360 max(relative change in state) = 0.31354E-05
362
363INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more
364 than 0.10000E-03 relative to the assumed background errors for the last 2 iterations.
365
366INFO (from ropp_1dvar_levmarq): Finished after 11 iterations (31 forward model / gradient evaluations).
367
368INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
369 first-1dvar-default.nc.
370
371##############################################################################
372run_1dvar-06 TR2 vn1
373
374---------------------------------------------------------------------
375 ROPP Bending Angle 1D-Var
376---------------------------------------------------------------------
377
378INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf.
379
380INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
381 ./data/metop-c/first.nc.
382
383INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
384 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
385
386INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
387
388INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
389
390 n_iter = 0 J = 97.039 max(relative change in state) = -
391 n_iter = 1 J = 35.030 max(relative change in state) = 7876.7
392 n_iter = 2 J = 32.719 max(relative change in state) = 0.36452
393 n_iter = 3 J = 32.433 max(relative change in state) = 0.18787
394 n_iter = 4 J = 32.370 max(relative change in state) = 0.12908
395 n_iter = 4 J = 32.371 mu -> 4.0000
396 n_iter = 5 J = 32.370 max(relative change in state) = 0.12539E-03
397 n_iter = 6 J = 32.370 max(relative change in state) = 0.21140E-02
398 n_iter = 7 J = 32.370 max(relative change in state) = 0.91153E-03
399 n_iter = 8 J = 32.370 max(relative change in state) = 0.42446E-03
400 n_iter = 9 J = 32.370 max(relative change in state) = 0.20538E-03
401 n_iter = 10 J = 32.370 max(relative change in state) = 0.10131E-03
402 n_iter = 11 J = 32.370 max(relative change in state) = 0.50460E-04
403 n_iter = 12 J = 32.370 max(relative change in state) = 0.25256E-04
404
405INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more
406 than 0.10000E-03 relative to the assumed background errors for the last 2 iterations.
407
408INFO (from ropp_1dvar_levmarq): Finished after 12 iterations (13 forward model / gradient evaluations).
409
410INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
411 first-1dvar-default.nc.
412
413##############################################################################
414run_1dvar-06 TR1 vn2
415
416---------------------------------------------------------------------
417 ROPP Bending Angle 1D-Var
418---------------------------------------------------------------------
419
420INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf.
421
422INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
423 ./data/metop-c/first.nc.
424
425INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
426 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
427
428INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
429
430INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
431
432 n_iter = 0 J = 97.039 max(relative change in state) = -
433 n_iter = 1 J = 35.030 max(relative change in state) = 7876.7
434 n_iter = 2 J = 32.532 max(relative change in state) = 0.41053
435 n_iter = 3 J = 32.382 max(relative change in state) = 0.17524
436 n_iter = 4 J = 32.370 max(relative change in state) = 0.88772E-01
437 n_iter = 4 J = 32.379 mu -> 0.44444E-02
438 n_iter = 5 J = 32.368 max(relative change in state) = 0.94089E-02
439 n_iter = 5 J = 32.377 mu -> 0.17778E-01
440 n_iter = 6 J = 32.366 max(relative change in state) = 0.91381E-02
441 n_iter = 6 J = 32.375 mu -> 0.71111E-01
442 n_iter = 7 J = 32.363 max(relative change in state) = 0.76370E-02
443 n_iter = 7 J = 32.377 mu -> 0.28444
444 n_iter = 8 J = 32.360 max(relative change in state) = 0.68296E-02
445 n_iter = 9 J = 32.359 max(relative change in state) = 0.59932E-02
446 n_iter = 10 J = 32.359 max(relative change in state) = 0.56609E-03
447 n_iter = 11 J = 32.359 max(relative change in state) = 0.32166E-04
448 n_iter = 12 J = 32.359 max(relative change in state) = 0.26036E-05
449
450INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more
451 than 0.10000E-03 relative to the assumed background errors for the last 2 iterations.
452
453INFO (from ropp_1dvar_levmarq): Finished after 12 iterations (16 forward model / gradient evaluations).
454
455INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
456 first-1dvar-default.nc.
457
458##############################################################################
459run_1dvar-06 TR2 vn2
460
461---------------------------------------------------------------------
462 ROPP Bending Angle 1D-Var
463---------------------------------------------------------------------
464
465INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf.
466
467INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
468 ./data/metop-c/first.nc.
469
470INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
471 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
472
473INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
474
475INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
476
477 n_iter = 0 J = 97.039 max(relative change in state) = -
478 n_iter = 1 J = 35.030 max(relative change in state) = 7876.7
479 n_iter = 2 J = 32.646 max(relative change in state) = 0.38236
480 n_iter = 3 J = 32.411 max(relative change in state) = 0.18946
481 n_iter = 4 J = 32.365 max(relative change in state) = 0.11897
482 n_iter = 4 J = 32.371 mu -> 0.43383E-03
483 n_iter = 5 J = 32.364 max(relative change in state) = 0.17311E-01
484 n_iter = 5 J = 32.672 mu -> 0.28922E-03
485 n_iter = 5 J = 32.494 mu -> 0.11569E-02
486 n_iter = 5 J = 32.547 mu -> 0.92550E-02
487 n_iter = 5 J = 32.416 mu -> 0.14808
488 n_iter = 5 J = 32.377 mu -> 4.7386
489 n_iter = 6 J = 32.362 max(relative change in state) = 0.11154E-01
490 n_iter = 7 J = 32.362 max(relative change in state) = 0.25784E-02
491 n_iter = 8 J = 32.362 max(relative change in state) = 0.11699E-02
492 n_iter = 9 J = 32.362 max(relative change in state) = 0.55834E-03
493 n_iter = 10 J = 32.362 max(relative change in state) = 0.27307E-03
494 n_iter = 11 J = 32.362 max(relative change in state) = 0.13517E-03
495 n_iter = 12 J = 32.362 max(relative change in state) = 0.67315E-04
496 n_iter = 13 J = 32.362 max(relative change in state) = 0.33623E-04
497
498INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more
499 than 0.10000E-03 relative to the assumed background errors for the last 2 iterations.
500
501INFO (from ropp_1dvar_levmarq): Finished after 13 iterations (19 forward model / gradient evaluations).
502
503INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
504 first-1dvar-default.nc.
505
506##############################################################################
507run_1dvar-06 NumRec vn2
508
509---------------------------------------------------------------------
510 ROPP Bending Angle 1D-Var
511---------------------------------------------------------------------
512
513INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf.
514
515INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file
516 ./data/metop-c/first.nc.
517
518INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file
519 ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc.
520
521INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km.
522
523INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass.
524
525 n_iter = 0 J = 97.039 max(relative change in state) = -
526 n_iter = 1 J = 32.394 max(relative change in state) = 7876.7
527 n_iter = 1 J = 32.560 lambda -> 0.10000E-03
528 n_iter = 1 J = 32.404 lambda -> 0.10000E-02
529 n_iter = 1 J = 32.469 lambda -> 0.10000E-01
530 n_iter = 1 J = 32.397 lambda -> 0.10000
531 n_iter = 2 J = 32.383 max(relative change in state) = 0.41798E-01
532 n_iter = 3 J = 32.366 max(relative change in state) = 0.34476E-01
533 n_iter = 3 J = 32.375 lambda -> 0.10000E-01
534 n_iter = 3 J = 32.368 lambda -> 0.10000
535 n_iter = 3 J = 32.367 lambda -> 1.0000
536 n_iter = 4 J = 32.366 max(relative change in state) = 0.32681E-03
537 n_iter = 4 J = 32.368 lambda -> 1.0000
538 n_iter = 5 J = 32.366 max(relative change in state) = 0.24915E-03
539 n_iter = 5 J = 32.368 lambda -> 1.0000
540 n_iter = 6 J = 32.366 max(relative change in state) = 0.24482E-03
541 n_iter = 6 J = 32.368 lambda -> 1.0000
542 n_iter = 7 J = 32.366 max(relative change in state) = 0.24052E-03
543 n_iter = 7 J = 32.368 lambda -> 1.0000
544 n_iter = 8 J = 32.365 max(relative change in state) = 0.23628E-03
545 n_iter = 8 J = 32.368 lambda -> 1.0000
546 n_iter = 9 J = 32.365 max(relative change in state) = 0.23211E-03
547 n_iter = 9 J = 32.368 lambda -> 1.0000
548 n_iter = 9 J = 32.366 lambda -> 10.000
549 n_iter = 10 J = 32.365 max(relative change in state) = 0.43513E-04
550 n_iter = 10 J = 32.366 lambda -> 10.000
551 n_iter = 11 J = 32.365 max(relative change in state) = 0.35401E-04
552
553INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more
554 than 0.10000E-03 relative to the assumed background errors for the last 2 iterations.
555
556INFO (from ropp_1dvar_levmarq): Finished after 11 iterations (26 forward model / gradient evaluations).
557
558INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file
559 first-1dvar-default.nc.
560
561##############################################################################
562