*** Results log of t_1dvar_COSMIC_04newop_bangle (1DVAR COSMIC; new interp) *** ./../tools/ropp_1dvar_bangle -y ../data/IT-1DVAR-04_y.nc -b ../data/IT-1DVAR-04_b.nc --bg-corr ../data/IT-1DVAR-04_c.nc -c ../config/metoffice_bangle_1dvar.cf -o IT-1DVAR-04newop_bangle.1.nc -new_op -d --------------------------------------------------------------------- ROPP Bending Angle 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_bangle): Reading configuration file ../config/metoffice_bangle_1dvar.cf. INFO (from ropp_1dvar_bangle): Reading observation data for profile 1 from the file ../data/IT-1DVAR-04_y.nc. WARNING (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 263772223.125 seconds - using yr/../ms timestamp INFO (from ropp_1dvar_bangle): Reading background data for profile 1 from the file ../data/IT-1DVAR-04_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 1.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 50 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 = 149.97 max(relative change in state) = - n_iter = 2 J = 130.08 max(relative change in state) = 0.82267 n_iter = 3 J = 123.54 max(relative change in state) = 0.87336 n_iter = 4 J = 123.22 max(relative change in state) = 0.47011 n_iter = 5 J = 122.24 max(relative change in state) = 0.23468 n_iter = 6 J = 122.07 max(relative change in state) = 0.76949E-01 n_iter = 7 J = 122.04 max(relative change in state) = 0.82616E-01 INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the state vector did not change by more than 0.10000 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 7. INFO (from ropp_1dvar_bangle): Writing 1DVar retrieval for profile 1 to the file IT-1DVAR-04newop_bangle.1.nc. INFO (from ropp_1dvar_bangle): Reading observation data for profile 2 from the file ../data/IT-1DVAR-04_y.nc. WARNING (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 263772222.906 seconds - using yr/../ms timestamp INFO (from ropp_1dvar_bangle): Reading background data for profile 2 from the file ../data/IT-1DVAR-04_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 1.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 50 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 = 87.267 max(relative change in state) = - n_iter = 2 J = 74.868 max(relative change in state) = 0.70135 n_iter = 3 J = 68.956 max(relative change in state) = 0.48658 n_iter = 4 J = 69.756 max(relative change in state) = 0.56033 n_iter = 5 J = 67.602 max(relative change in state) = 0.39501 n_iter = 6 J = 67.503 max(relative change in state) = 0.55271E-01 n_iter = 7 J = 67.424 max(relative change in state) = 0.54618E-01 INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the state vector did not change by more than 0.10000 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 7. INFO (from ropp_1dvar_bangle): Appending 1DVar retrieval for profile 2 to the file IT-1DVAR-04newop_bangle.1.nc. INFO (from ropp_1dvar_bangle): Reading observation data for profile 3 from the file ../data/IT-1DVAR-04_y.nc. WARNING (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 263772223.122 seconds - using yr/../ms timestamp INFO (from ropp_1dvar_bangle): Reading background data for profile 3 from the file ../data/IT-1DVAR-04_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 1.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 50 points rejected for being above maximum impact height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control removes 1 bending angle data points from the observations as their deviation from the background exceeds 10.0 times the expected (1-sigma) error. 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 = 66.070 max(relative change in state) = - n_iter = 2 J = 55.756 max(relative change in state) = 0.58705 n_iter = 3 J = 47.591 max(relative change in state) = 1.2262 n_iter = 4 J = 53.100 max(relative change in state) = 0.81538 n_iter = 5 J = 47.177 max(relative change in state) = 0.69805 n_iter = 6 J = 47.055 max(relative change in state) = 0.22450E-01 n_iter = 7 J = 46.638 max(relative change in state) = 0.20814 n_iter = 8 J = 46.478 max(relative change in state) = 0.11611 n_iter = 9 J = 46.645 max(relative change in state) = 0.15557 n_iter = 10 J = 46.450 max(relative change in state) = 0.11927 n_iter = 11 J = 46.443 max(relative change in state) = 0.76598E-02 n_iter = 12 J = 46.436 max(relative change in state) = 0.61555E-01 INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the state vector did not change by more than 0.10000 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 12. INFO (from ropp_1dvar_bangle): Appending 1DVar retrieval for profile 3 to the file IT-1DVAR-04newop_bangle.1.nc. INFO (from ropp_1dvar_bangle): Reading observation data for profile 4 from the file ../data/IT-1DVAR-04_y.nc. WARNING (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 263772222.896 seconds - using yr/../ms timestamp INFO (from ropp_1dvar_bangle): Reading background data for profile 4 from the file ../data/IT-1DVAR-04_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 1.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 51 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 = 175.20 max(relative change in state) = - n_iter = 2 J = 155.94 max(relative change in state) = 0.92166 n_iter = 3 J = 142.75 max(relative change in state) = 1.4051 n_iter = 4 J = 154.21 max(relative change in state) = 1.5632 n_iter = 5 J = 139.48 max(relative change in state) = 1.1211 n_iter = 6 J = 138.82 max(relative change in state) = 0.10015 n_iter = 7 J = 140.07 max(relative change in state) = 0.53740 n_iter = 8 J = 137.89 max(relative change in state) = 0.34994 n_iter = 9 J = 137.30 max(relative change in state) = 0.15741 n_iter = 10 J = 138.28 max(relative change in state) = 0.32295 n_iter = 11 J = 137.00 max(relative change in state) = 0.29020 n_iter = 12 J = 136.75 max(relative change in state) = 0.57659E-01 n_iter = 13 J = 136.59 max(relative change in state) = 0.29401 n_iter = 14 J = 136.24 max(relative change in state) = 0.24140 n_iter = 15 J = 136.18 max(relative change in state) = 0.11504 n_iter = 16 J = 136.10 max(relative change in state) = 0.13411 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 2 iterations. INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 16. INFO (from ropp_1dvar_bangle): Appending 1DVar retrieval for profile 4 to the file IT-1DVAR-04newop_bangle.1.nc. INFO (from ropp_1dvar_bangle): Reading observation data for profile 5 from the file ../data/IT-1DVAR-04_y.nc. WARNING (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 263772223.118 seconds - using yr/../ms timestamp INFO (from ropp_1dvar_bangle): Reading background data for profile 5 from the file ../data/IT-1DVAR-04_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 1.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 49 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 = 98.895 max(relative change in state) = - n_iter = 2 J = 84.273 max(relative change in state) = 0.61193 n_iter = 3 J = 69.134 max(relative change in state) = 1.0399 n_iter = 4 J = 66.731 max(relative change in state) = 0.41293 n_iter = 5 J = 65.335 max(relative change in state) = 0.40183 n_iter = 6 J = 65.277 max(relative change in state) = 0.26982 n_iter = 7 J = 64.907 max(relative change in state) = 0.13679 n_iter = 8 J = 64.860 max(relative change in state) = 0.51302E-01 n_iter = 9 J = 64.796 max(relative change in state) = 0.65211E-01 INFO (from ropp_1dvar_cost): Convergence assumed to be achieved as the state vector did not change by more than 0.10000 relative to the assumed background errors for the last 2 iterations. INFO (from ropp_1dvar_solve): Minimization of cost function successfully finished (according to additional convergence criteria). Number of required iterations: 9. INFO (from ropp_1dvar_bangle): Appending 1DVar retrieval for profile 5 to the file IT-1DVAR-04newop_bangle.1.nc. *** Running IDL to plot results IDL 8.8.2 (linux x86_64 m64). (c) 2022, L3Harris Geospatial Solutions, Inc. Licensed for use by: University of Maryland (MAIN) License: 94891 A new version is available: IDL 8.8.3 https://harrisgeospatial.flexnetoperations.com % Compiled module: PLOT_1DVAR. PLOT_1DVAR: Plotting increments for file_id: 04newop_bangle % Loaded DLM: NCDF. % Compiled module: NCDF_GETVAR. Processing scenario: OC_20080512000000_CO01_U999_METO % Compiled module: INTERPOL. Processing scenario: OC_20080512000000_CO01_U999_METO Processing scenario: OC_20080512000000_CO04_U999_METO Processing scenario: OC_20080512000000_CO04_U999_METO Processing scenario: OC_20080512000000_CO05_U999_METO Found max difference [%] in observation scenarios: -3.8120445e+11 Found max difference [K] in temperature scenarios: -2.2785997 Found max difference [g/kg] in specific humidity scenarios: 0.79168464 Found max difference [hPa] in pressure scenarios: -1.0891737 Generating output image file ropp_1dvar_04newop_bangle_bangle.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. % Program caused arithmetic error: Floating illegal operand *** Finished IDL *** ropp_1dvar_04newop_bangle.jpg can be compared against ../data/ropp_1dvar_04newop_bangle_reference.jpg