########################################################################### run_1dvar-00 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-00.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. INFO (from ropp_1dvar_solve): Using background error covariance matrix for preconditioning. INFO (from ropp_1dvar_cost): Using absolute decrease of cost function and relative change of state vector as additional convergence criteria. n_iter = 1 J = 97.039 max(relative change in state) = - n_iter = 2 J = 79.446 max(relative change in state) = 0.42320 INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the cost function is increasing INFO (from ropp_1dvar_cost): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 3. INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ########################################################################### run_1dvar-01 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-01.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 1 J = 97.039 max(relative change in state) = - n_iter = 2 J = 32.400 max(relative change in state) = 0.13905 n_iter = 3 J = 32.382 max(relative change in state) = 0.16073 n_iter = 4 J = 32.378 max(relative change in state) = 0.12803 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the cost function did not change by more than 0.10000 for the last 2 iterations. INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-02 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-02.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. INFO (from ropp_1dvar_solve): Using background error covariance matrix for preconditioning. INFO (from ropp_1dvar_cost): Using absolute decrease of cost function and relative change of state vector as additional convergence criteria. n_iter = 1 J = 97.039 max(relative change in state) = - n_iter = 2 J = 79.446 max(relative change in state) = 0.42320 n_iter = 3 J = 16747. max(relative change in state) = 71.630 n_iter = 4 J = 76.892 max(relative change in state) = 71.347 n_iter = 5 J = 76.486 max(relative change in state) = 0.11244E-01 n_iter = 6 J = 373.25 max(relative change in state) = 9.9224 n_iter = 7 J = 53.400 max(relative change in state) = 8.7370 n_iter = 8 J = 49.459 max(relative change in state) = 0.30230 n_iter = 9 J = 43.141 max(relative change in state) = 0.74605 n_iter = 10 J = 43.110 max(relative change in state) = 0.25562 n_iter = 11 J = 43.023 max(relative change in state) = 0.12781 n_iter = 12 J = 43.016 max(relative change in state) = 0.75320E-02 n_iter = 13 J = 42.914 max(relative change in state) = 0.76796E-01 n_iter = 14 J = 43.489 max(relative change in state) = 0.46498 n_iter = 15 J = 42.885 max(relative change in state) = 0.39930 n_iter = 16 J = 42.878 max(relative change in state) = 0.60095E-02 n_iter = 17 J = 42.769 max(relative change in state) = 0.97121E-01 n_iter = 18 J = 43.071 max(relative change in state) = 0.35938 n_iter = 19 J = 42.741 max(relative change in state) = 0.30358 n_iter = 20 J = 42.730 max(relative change in state) = 0.87931E-02 n_iter = 21 J = 42.565 max(relative change in state) = 0.11115 n_iter = 22 J = 42.476 max(relative change in state) = 0.10533 n_iter = 23 J = 42.462 max(relative change in state) = 0.40446E-01 n_iter = 24 J = 42.445 max(relative change in state) = 0.29540E-01 n_iter = 25 J = 42.358 max(relative change in state) = 0.52629E-01 n_iter = 26 J = 42.076 max(relative change in state) = 0.16068 n_iter = 27 J = 42.040 max(relative change in state) = 0.28554E-01 n_iter = 28 J = 41.725 max(relative change in state) = 0.16654 n_iter = 29 J = 40.846 max(relative change in state) = 0.12660 n_iter = 30 J = 39.286 max(relative change in state) = 0.75896 n_iter = 31 J = 38.406 max(relative change in state) = 0.59504 n_iter = 32 J = 38.350 max(relative change in state) = 0.10336 n_iter = 33 J = 38.342 max(relative change in state) = 0.28838E-01 n_iter = 34 J = 38.334 max(relative change in state) = 0.29887E-01 n_iter = 35 J = 38.322 max(relative change in state) = 0.22172E-01 n_iter = 36 J = 38.245 max(relative change in state) = 0.61646E-01 INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the cost function did not change by more than 0.10000 for the last 5 iterations. INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 36. ... (from ropp_1dvar_diagnostics): Warning: Too many iterations, niter = 36 2J/m = 0.361 INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-03 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-03.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. INFO (from ropp_1dvar_solve): Using background error covariance matrix for preconditioning. INFO (from ropp_1dvar_cost): Using absolute decrease of cost function and relative change of state vector as additional convergence criteria. n_iter = 1 J = 97.039 max(relative change in state) = - n_iter = 2 J = 79.446 max(relative change in state) = 0.42320 n_iter = 3 J = 16747. max(relative change in state) = 71.630 n_iter = 4 J = 76.892 max(relative change in state) = 71.347 n_iter = 5 J = 76.486 max(relative change in state) = 0.11244E-01 n_iter = 6 J = 373.25 max(relative change in state) = 9.9224 n_iter = 7 J = 53.400 max(relative change in state) = 8.7370 n_iter = 8 J = 49.459 max(relative change in state) = 0.30230 n_iter = 9 J = 43.141 max(relative change in state) = 0.74605 n_iter = 10 J = 43.110 max(relative change in state) = 0.25562 n_iter = 11 J = 43.023 max(relative change in state) = 0.12781 n_iter = 12 J = 43.016 max(relative change in state) = 0.75320E-02 n_iter = 13 J = 42.914 max(relative change in state) = 0.76796E-01 n_iter = 14 J = 43.489 max(relative change in state) = 0.46498 n_iter = 15 J = 42.885 max(relative change in state) = 0.39930 n_iter = 16 J = 42.878 max(relative change in state) = 0.60095E-02 n_iter = 17 J = 42.769 max(relative change in state) = 0.97121E-01 n_iter = 18 J = 43.071 max(relative change in state) = 0.35938 n_iter = 19 J = 42.741 max(relative change in state) = 0.30358 n_iter = 20 J = 42.730 max(relative change in state) = 0.87931E-02 n_iter = 21 J = 42.565 max(relative change in state) = 0.11115 n_iter = 22 J = 42.476 max(relative change in state) = 0.10533 n_iter = 23 J = 42.462 max(relative change in state) = 0.40446E-01 n_iter = 24 J = 42.445 max(relative change in state) = 0.29540E-01 n_iter = 25 J = 42.358 max(relative change in state) = 0.52629E-01 n_iter = 26 J = 42.076 max(relative change in state) = 0.16068 n_iter = 27 J = 42.040 max(relative change in state) = 0.28554E-01 n_iter = 28 J = 41.725 max(relative change in state) = 0.16654 n_iter = 29 J = 40.846 max(relative change in state) = 0.12660 n_iter = 30 J = 39.286 max(relative change in state) = 0.75896 n_iter = 31 J = 38.406 max(relative change in state) = 0.59504 n_iter = 32 J = 38.350 max(relative change in state) = 0.10336 n_iter = 33 J = 38.342 max(relative change in state) = 0.28838E-01 n_iter = 34 J = 38.334 max(relative change in state) = 0.29887E-01 n_iter = 35 J = 38.322 max(relative change in state) = 0.22172E-01 n_iter = 36 J = 38.245 max(relative change in state) = 0.61646E-01 n_iter = 37 J = 37.917 max(relative change in state) = 0.74321E-01 n_iter = 38 J = 50.348 max(relative change in state) = 2.0487 n_iter = 39 J = 37.722 max(relative change in state) = 1.8466 n_iter = 40 J = 37.700 max(relative change in state) = 0.86831E-02 n_iter = 41 J = 37.600 max(relative change in state) = 0.89404E-01 n_iter = 42 J = 317.92 max(relative change in state) = 9.3296 n_iter = 43 J = 37.579 max(relative change in state) = 9.2572 n_iter = 44 J = 37.576 max(relative change in state) = 0.34352E-03 n_iter = 45 J = 91.314 max(relative change in state) = 4.0968 n_iter = 46 J = 36.909 max(relative change in state) = 3.8664 n_iter = 47 J = 36.614 max(relative change in state) = 0.35645E-01 n_iter = 48 J = 35.339 max(relative change in state) = 0.35465 n_iter = 49 J = 3636.7 max(relative change in state) = 33.363 n_iter = 50 J = 35.247 max(relative change in state) = 33.198 INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 50. ... (from ropp_1dvar_diagnostics): Warning: Too many iterations, niter = 50 2J/m = 0.333 INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-04 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-04.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 1 J = 97.039 max(relative change in state) = - n_iter = 2 J = 32.400 max(relative change in state) = 0.13905 n_iter = 3 J = 32.382 max(relative change in state) = 0.16073 n_iter = 4 J = 32.378 max(relative change in state) = 0.12803 n_iter = 5 J = 32.376 max(relative change in state) = 0.12886 n_iter = 6 J = 32.378 max(relative change in state) = 0.12885 n_iter = 7 J = 32.416 max(relative change in state) = 0.12841 n_iter = 8 J = 32.482 max(relative change in state) = 0.98499E-01 n_iter = 9 J = 32.559 max(relative change in state) = 0.16935E-01 n_iter = 10 J = 32.593 max(relative change in state) = 0.80768E-03 n_iter = 11 J = 32.596 max(relative change in state) = 0.88798E-05 n_iter = 12 J = 32.596 max(relative change in state) = 0.88887E-07 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more than 0.10000E-03 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-06 NumRec vn1 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 0 J = 97.039 max(relative change in state) = - n_iter = 1 J = 32.394 max(relative change in state) = 7876.7 n_iter = 1 J = 32.560 lambda -> 0.10000E-03 n_iter = 1 J = 32.404 lambda -> 0.10000E-02 n_iter = 1 J = 32.469 lambda -> 0.10000E-01 n_iter = 1 J = 32.397 lambda -> 0.10000 n_iter = 2 J = 32.383 max(relative change in state) = 0.41798E-01 n_iter = 3 J = 32.366 max(relative change in state) = 0.34476E-01 n_iter = 3 J = 32.375 lambda -> 0.10000E-01 n_iter = 3 J = 32.368 lambda -> 0.10000 n_iter = 3 J = 32.367 lambda -> 1.0000 n_iter = 4 J = 32.366 max(relative change in state) = 0.32681E-03 n_iter = 4 J = 32.368 lambda -> 1.0000 n_iter = 5 J = 32.366 max(relative change in state) = 0.24915E-03 n_iter = 5 J = 32.368 lambda -> 1.0000 n_iter = 6 J = 32.366 max(relative change in state) = 0.24482E-03 n_iter = 6 J = 32.368 lambda -> 1.0000 n_iter = 7 J = 32.366 max(relative change in state) = 0.24052E-03 n_iter = 7 J = 32.368 lambda -> 1.0000 n_iter = 8 J = 32.365 max(relative change in state) = 0.23628E-03 n_iter = 8 J = 32.368 lambda -> 1.0000 n_iter = 9 J = 32.365 max(relative change in state) = 0.23211E-03 n_iter = 9 J = 32.368 lambda -> 1.0000 n_iter = 9 J = 32.366 lambda -> 10.000 n_iter = 10 J = 32.365 max(relative change in state) = 0.43513E-04 n_iter = 10 J = 32.366 lambda -> 10.000 n_iter = 11 J = 32.365 max(relative change in state) = 0.35401E-04 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more than 0.10000E-03 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_levmarq): Finished after 11 iterations (26 forward model / gradient evaluations). INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-06 TR1 vn1 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 0 J = 97.039 max(relative change in state) = - n_iter = 1 J = 35.030 max(relative change in state) = 7876.7 n_iter = 2 J = 32.532 max(relative change in state) = 0.41053 n_iter = 3 J = 32.382 max(relative change in state) = 0.17524 n_iter = 4 J = 32.370 max(relative change in state) = 0.88772E-01 n_iter = 4 J = 32.379 mu -> 0.74074E-03 n_iter = 5 J = 32.367 max(relative change in state) = 0.16165E-01 n_iter = 5 J = 32.404 mu -> 0.49383E-03 n_iter = 6 J = 32.361 max(relative change in state) = 0.26131E-01 n_iter = 6 J = 32.372 mu -> 0.32922E-03 n_iter = 6 J = 32.364 mu -> 0.65844E-03 n_iter = 6 J = 32.369 mu -> 0.13169E-02 n_iter = 6 J = 32.366 mu -> 0.26337E-02 n_iter = 6 J = 32.366 mu -> 0.52675E-02 n_iter = 6 J = 32.365 mu -> 0.10535E-01 n_iter = 7 J = 32.360 max(relative change in state) = 0.12883E-01 n_iter = 7 J = 32.367 mu -> 0.70233E-02 n_iter = 7 J = 32.364 mu -> 0.14047E-01 n_iter = 7 J = 32.363 mu -> 0.28093E-01 n_iter = 7 J = 32.362 mu -> 0.56187E-01 n_iter = 8 J = 32.360 max(relative change in state) = 0.18178E-02 n_iter = 8 J = 32.362 mu -> 0.11237 n_iter = 8 J = 32.361 mu -> 0.22475 n_iter = 8 J = 32.361 mu -> 0.44949 n_iter = 8 J = 32.360 mu -> 0.89898 n_iter = 9 J = 32.360 max(relative change in state) = 0.18119E-03 n_iter = 9 J = 32.360 mu -> 3.5959 n_iter = 9 J = 32.360 mu -> 7.1919 n_iter = 9 J = 32.360 mu -> 14.384 n_iter = 9 J = 32.360 mu -> 28.768 n_iter = 10 J = 32.360 max(relative change in state) = 0.60608E-05 n_iter = 11 J = 32.360 max(relative change in state) = 0.31354E-05 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more than 0.10000E-03 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_levmarq): Finished after 11 iterations (31 forward model / gradient evaluations). INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-06 TR2 vn1 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 0 J = 97.039 max(relative change in state) = - n_iter = 1 J = 35.030 max(relative change in state) = 7876.7 n_iter = 2 J = 32.719 max(relative change in state) = 0.36452 n_iter = 3 J = 32.433 max(relative change in state) = 0.18787 n_iter = 4 J = 32.370 max(relative change in state) = 0.12908 n_iter = 4 J = 32.371 mu -> 4.0000 n_iter = 5 J = 32.370 max(relative change in state) = 0.12539E-03 n_iter = 6 J = 32.370 max(relative change in state) = 0.21140E-02 n_iter = 7 J = 32.370 max(relative change in state) = 0.91153E-03 n_iter = 8 J = 32.370 max(relative change in state) = 0.42446E-03 n_iter = 9 J = 32.370 max(relative change in state) = 0.20538E-03 n_iter = 10 J = 32.370 max(relative change in state) = 0.10131E-03 n_iter = 11 J = 32.370 max(relative change in state) = 0.50460E-04 n_iter = 12 J = 32.370 max(relative change in state) = 0.25256E-04 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more than 0.10000E-03 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_levmarq): Finished after 12 iterations (13 forward model / gradient evaluations). INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-06 TR1 vn2 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 0 J = 97.039 max(relative change in state) = - n_iter = 1 J = 35.030 max(relative change in state) = 7876.7 n_iter = 2 J = 32.532 max(relative change in state) = 0.41053 n_iter = 3 J = 32.382 max(relative change in state) = 0.17524 n_iter = 4 J = 32.370 max(relative change in state) = 0.88772E-01 n_iter = 4 J = 32.379 mu -> 0.44444E-02 n_iter = 5 J = 32.368 max(relative change in state) = 0.94089E-02 n_iter = 5 J = 32.377 mu -> 0.17778E-01 n_iter = 6 J = 32.366 max(relative change in state) = 0.91381E-02 n_iter = 6 J = 32.375 mu -> 0.71111E-01 n_iter = 7 J = 32.363 max(relative change in state) = 0.76370E-02 n_iter = 7 J = 32.377 mu -> 0.28444 n_iter = 8 J = 32.360 max(relative change in state) = 0.68296E-02 n_iter = 9 J = 32.359 max(relative change in state) = 0.59932E-02 n_iter = 10 J = 32.359 max(relative change in state) = 0.56609E-03 n_iter = 11 J = 32.359 max(relative change in state) = 0.32166E-04 n_iter = 12 J = 32.359 max(relative change in state) = 0.26036E-05 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more than 0.10000E-03 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_levmarq): Finished after 12 iterations (16 forward model / gradient evaluations). INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-06 TR2 vn2 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 0 J = 97.039 max(relative change in state) = - n_iter = 1 J = 35.030 max(relative change in state) = 7876.7 n_iter = 2 J = 32.646 max(relative change in state) = 0.38236 n_iter = 3 J = 32.411 max(relative change in state) = 0.18946 n_iter = 4 J = 32.365 max(relative change in state) = 0.11897 n_iter = 4 J = 32.371 mu -> 0.43383E-03 n_iter = 5 J = 32.364 max(relative change in state) = 0.17311E-01 n_iter = 5 J = 32.672 mu -> 0.28922E-03 n_iter = 5 J = 32.494 mu -> 0.11569E-02 n_iter = 5 J = 32.547 mu -> 0.92550E-02 n_iter = 5 J = 32.416 mu -> 0.14808 n_iter = 5 J = 32.377 mu -> 4.7386 n_iter = 6 J = 32.362 max(relative change in state) = 0.11154E-01 n_iter = 7 J = 32.362 max(relative change in state) = 0.25784E-02 n_iter = 8 J = 32.362 max(relative change in state) = 0.11699E-02 n_iter = 9 J = 32.362 max(relative change in state) = 0.55834E-03 n_iter = 10 J = 32.362 max(relative change in state) = 0.27307E-03 n_iter = 11 J = 32.362 max(relative change in state) = 0.13517E-03 n_iter = 12 J = 32.362 max(relative change in state) = 0.67315E-04 n_iter = 13 J = 32.362 max(relative change in state) = 0.33623E-04 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more than 0.10000E-03 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_levmarq): Finished after 13 iterations (19 forward model / gradient evaluations). INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ############################################################################## run_1dvar-06 NumRec vn2 --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ecmwf_bangle_1dvar-06.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ./data/metop-c/first.nc. INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ./data/ecmwf/GRAS_1B_M03_20181113100752Z_20181113100945Z_N_C_20181113121522Z_G22_NN.fw.nc. INFO (from ropp_qc_cutoff): 33 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all bending angle values pass. n_iter = 0 J = 97.039 max(relative change in state) = - n_iter = 1 J = 32.394 max(relative change in state) = 7876.7 n_iter = 1 J = 32.560 lambda -> 0.10000E-03 n_iter = 1 J = 32.404 lambda -> 0.10000E-02 n_iter = 1 J = 32.469 lambda -> 0.10000E-01 n_iter = 1 J = 32.397 lambda -> 0.10000 n_iter = 2 J = 32.383 max(relative change in state) = 0.41798E-01 n_iter = 3 J = 32.366 max(relative change in state) = 0.34476E-01 n_iter = 3 J = 32.375 lambda -> 0.10000E-01 n_iter = 3 J = 32.368 lambda -> 0.10000 n_iter = 3 J = 32.367 lambda -> 1.0000 n_iter = 4 J = 32.366 max(relative change in state) = 0.32681E-03 n_iter = 4 J = 32.368 lambda -> 1.0000 n_iter = 5 J = 32.366 max(relative change in state) = 0.24915E-03 n_iter = 5 J = 32.368 lambda -> 1.0000 n_iter = 6 J = 32.366 max(relative change in state) = 0.24482E-03 n_iter = 6 J = 32.368 lambda -> 1.0000 n_iter = 7 J = 32.366 max(relative change in state) = 0.24052E-03 n_iter = 7 J = 32.368 lambda -> 1.0000 n_iter = 8 J = 32.365 max(relative change in state) = 0.23628E-03 n_iter = 8 J = 32.368 lambda -> 1.0000 n_iter = 9 J = 32.365 max(relative change in state) = 0.23211E-03 n_iter = 9 J = 32.368 lambda -> 1.0000 n_iter = 9 J = 32.366 lambda -> 10.000 n_iter = 10 J = 32.365 max(relative change in state) = 0.43513E-04 n_iter = 10 J = 32.366 lambda -> 10.000 n_iter = 11 J = 32.365 max(relative change in state) = 0.35401E-04 INFO (from ropp_1dvar_levmarq): Convergence assumed to be achieved as the state vector did not change by more than 0.10000E-03 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_levmarq): Finished after 11 iterations (26 forward model / gradient evaluations). INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file first-1dvar-default.nc. ##############################################################################