*** Results log of t_1dvar_COSMIC_04_refrac (1DVAR COSMIC; default) *** ./../tools/ropp_1dvar_refrac -y ../data/IT-1DVAR-04_y.nc -b ../data/IT-1DVAR-04_b.nc --bg-corr ../data/IT-1DVAR-04_c.nc --obs-corr ../data/IT-1DVAR-04_o.nc -c ../config/metoffice_refrac_1dvar.cf -o IT-1DVAR-04_refrac.1.nc -d --------------------------------------------------------------------- ROPP Refractivity 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_refrac): Reading configuration file ../config/metoffice_refrac_1dvar.cf... INFO (from ropp_1dvar_refrac): 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_refrac): 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): 48 points rejected for being above maximum geopotential height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all refractivity 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 criterium. n_iter = 1 J = 55.619 max(relative change in state) = - n_iter = 2 J = 53.171 max(relative change in state) = 0.91119 n_iter = 3 J = 52.717 max(relative change in state) = 0.19503 n_iter = 4 J = 52.666 max(relative change in state) = 0.89036E-01 n_iter = 5 J = 52.206 max(relative change in state) = 0.22761 n_iter = 6 J = 52.232 max(relative change in state) = 0.11037 n_iter = 7 J = 52.075 max(relative change in state) = 0.11088 n_iter = 8 J = 52.079 max(relative change in state) = 0.30579E-01 n_iter = 9 J = 52.027 max(relative change in state) = 0.59843E-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_refrac): Writing 1DVar retrieval for profile 1 to the file IT-1DVAR-04_refrac.1.nc. INFO (from ropp_1dvar_refrac): 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_refrac): 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): 48 points rejected for being above maximum geopotential height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all refractivity 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 criterium. n_iter = 1 J = 26.132 max(relative change in state) = - n_iter = 2 J = 25.665 max(relative change in state) = 0.55513 n_iter = 3 J = 24.905 max(relative change in state) = 0.39578 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: 4. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 2 to the file IT-1DVAR-04_refrac.1.nc. INFO (from ropp_1dvar_refrac): 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_refrac): 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): 48 points rejected for being above maximum geopotential height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all refractivity 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 criterium. n_iter = 1 J = 37.143 max(relative change in state) = - n_iter = 2 J = 32.769 max(relative change in state) = 0.69294 n_iter = 3 J = 31.863 max(relative change in state) = 0.40386 n_iter = 4 J = 36.141 max(relative change in state) = 0.68427 n_iter = 5 J = 31.783 max(relative change in state) = 0.61766 n_iter = 6 J = 31.737 max(relative change in state) = 0.16517E-01 n_iter = 7 J = 32.853 max(relative change in state) = 0.32041 n_iter = 8 J = 31.695 max(relative change in state) = 0.31733 n_iter = 9 J = 31.656 max(relative change in state) = 0.16311E-01 n_iter = 10 J = 31.780 max(relative change in state) = 0.11181 n_iter = 11 J = 31.611 max(relative change in state) = 0.10013 n_iter = 12 J = 31.569 max(relative change in state) = 0.17193E-01 n_iter = 13 J = 31.637 max(relative change in state) = 0.10037 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: 13. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 3 to the file IT-1DVAR-04_refrac.1.nc. INFO (from ropp_1dvar_refrac): 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_refrac): 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): 49 points rejected for being above maximum geopotential height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all refractivity 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 criterium. n_iter = 1 J = 33.327 max(relative change in state) = - n_iter = 2 J = 32.277 max(relative change in state) = 0.41943 n_iter = 3 J = 31.864 max(relative change in state) = 0.27890 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: 4. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 4 to the file IT-1DVAR-04_refrac.1.nc. INFO (from ropp_1dvar_refrac): 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_refrac): 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): 47 points rejected for being above maximum geopotential height of 50.0 km. INFO (from ropp_qc_bgqc): Background quality control lets all refractivity 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 criterium. n_iter = 1 J = 28.623 max(relative change in state) = - n_iter = 2 J = 27.087 max(relative change in state) = 0.48814 n_iter = 3 J = 26.672 max(relative change in state) = 0.25898 n_iter = 4 J = 26.910 max(relative change in state) = 0.22307 n_iter = 5 J = 26.698 max(relative change in state) = 0.83727E-01 n_iter = 6 J = 26.459 max(relative change in state) = 0.13445 n_iter = 7 J = 26.412 max(relative change in state) = 0.15589 n_iter = 8 J = 26.322 max(relative change in state) = 0.15502 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: 8. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 5 to the file IT-1DVAR-04_refrac.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: 04_refrac % 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: -3.5585007 Found max difference [g/kg] in specific humidity scenarios: -0.77814325 Found max difference [hPa] in pressure scenarios: -0.45877572 Generating output image file ropp_1dvar_04_refrac_bangle.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. % Program caused arithmetic error: Floating illegal operand *** Finished IDL *** ropp_1dvar_04_refrac.jpg can be compared against ../data/ropp_1dvar_04_refrac_reference.jpg