*** Results log of t_1dvar_GRAS_05comp_refrac (1DVAR GRAS; comp factors) *** ./../tools/ropp_1dvar_refrac -y ../data/IT-1DVAR-05_y.nc -b ../data/IT-1DVAR-05_b.nc --bg-corr ../data/IT-1DVAR-05_c.nc --obs-corr ../errors/ropp_ob_refrac_error_corr_247L.nc -c ../config/ecmwf_refrac_1dvar.cf -o IT-1DVAR-05comp_refrac.1.nc -comp -d --------------------------------------------------------------------- ROPP Refractivity 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_refrac): Reading configuration file ../config/ecmwf_refrac_1dvar.cf... INFO (from ropp_1dvar_refrac): Reading observation data for profile 1 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 1 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 32 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 = 40.880 max(relative change in state) = - n_iter = 2 J = 33.737 max(relative change in state) = 0.45255 n_iter = 3 J = 30.923 max(relative change in state) = 0.33484 n_iter = 4 J = 30.130 max(relative change in state) = 0.30964 n_iter = 5 J = 27.805 max(relative change in state) = 0.44510 n_iter = 6 J = 25.393 max(relative change in state) = 0.19707 n_iter = 7 J = 24.887 max(relative change in state) = 0.19441 n_iter = 8 J = 24.102 max(relative change in state) = 0.16191 n_iter = 9 J = 23.869 max(relative change in state) = 0.18752 n_iter = 10 J = 23.488 max(relative change in state) = 0.20692 n_iter = 11 J = 23.271 max(relative change in state) = 0.62004E-01 n_iter = 12 J = 23.056 max(relative change in state) = 0.67122E-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_refrac): Writing 1DVar retrieval for profile 1 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 2 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 2 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 31 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 = 36.200 max(relative change in state) = - n_iter = 2 J = 34.350 max(relative change in state) = 0.37425 n_iter = 3 J = 32.721 max(relative change in state) = 0.42703 n_iter = 4 J = 34.522 max(relative change in state) = 0.81557 n_iter = 5 J = 32.642 max(relative change in state) = 0.77496 n_iter = 6 J = 32.581 max(relative change in state) = 0.21335E-01 n_iter = 7 J = 32.191 max(relative change in state) = 0.44828 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: 8. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 2 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 3 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 3 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 32 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 = 15.968 max(relative change in state) = - n_iter = 2 J = 14.220 max(relative change in state) = 0.36273 n_iter = 3 J = 12.564 max(relative change in state) = 0.61804 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 3 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 4 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 4 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 31 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 = 34.332 max(relative change in state) = - n_iter = 2 J = 29.747 max(relative change in state) = 0.82233 n_iter = 3 J = 29.274 max(relative change in state) = 0.54791 n_iter = 4 J = 29.154 max(relative change in state) = 0.62942 n_iter = 5 J = 28.972 max(relative change in state) = 0.17845 n_iter = 6 J = 28.966 max(relative change in state) = 0.18660E-01 n_iter = 7 J = 28.953 max(relative change in state) = 0.11075 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: 7. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 4 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 5 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 5 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 32 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 = 38.441 max(relative change in state) = - n_iter = 2 J = 34.668 max(relative change in state) = 0.24222 n_iter = 3 J = 32.997 max(relative change in state) = 0.16645 n_iter = 4 J = 31.707 max(relative change in state) = 0.14879 n_iter = 5 J = 30.717 max(relative change in state) = 0.28818 n_iter = 6 J = 30.177 max(relative change in state) = 0.19237 n_iter = 7 J = 29.286 max(relative change in state) = 0.12162 n_iter = 8 J = 28.814 max(relative change in state) = 0.83455E-01 n_iter = 9 J = 28.618 max(relative change in state) = 0.99576E-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): Appending 1DVar retrieval for profile 5 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 6 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 6 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 32 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 = 136.91 max(relative change in state) = - n_iter = 2 J = 133.94 max(relative change in state) = 0.14235 n_iter = 3 J = 116.24 max(relative change in state) = 1.4096 n_iter = 4 J = 132.66 max(relative change in state) = 3.2969 n_iter = 5 J = 114.18 max(relative change in state) = 2.9332 n_iter = 6 J = 113.14 max(relative change in state) = 0.10929 n_iter = 7 J = 111.36 max(relative change in state) = 0.42423 n_iter = 8 J = 108.36 max(relative change in state) = 0.46125 n_iter = 9 J = 107.30 max(relative change in state) = 0.40757 n_iter = 10 J = 107.57 max(relative change in state) = 0.24655 n_iter = 11 J = 105.82 max(relative change in state) = 0.19267 n_iter = 12 J = 105.34 max(relative change in state) = 0.99332E-01 n_iter = 13 J = 104.67 max(relative change in state) = 0.18175 n_iter = 14 J = 104.41 max(relative change in state) = 0.60459E-01 n_iter = 15 J = 104.01 max(relative change in state) = 0.15869 n_iter = 16 J = 103.71 max(relative change in state) = 0.51141 n_iter = 17 J = 103.39 max(relative change in state) = 0.99943E-01 n_iter = 18 J = 103.22 max(relative change in state) = 0.67829E-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: 18. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 6 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 7 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 7 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 31 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 = 24.173 max(relative change in state) = - n_iter = 2 J = 21.332 max(relative change in state) = 0.80255 n_iter = 3 J = 21.103 max(relative change in state) = 0.26143 n_iter = 4 J = 21.540 max(relative change in state) = 0.24266 n_iter = 5 J = 21.070 max(relative change in state) = 0.17969 n_iter = 6 J = 21.068 max(relative change in state) = 0.58651E-02 n_iter = 7 J = 21.061 max(relative change in state) = 0.56725E-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_refrac): Appending 1DVar retrieval for profile 7 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 8 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 8 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 32 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.463 max(relative change in state) = - n_iter = 2 J = 31.469 max(relative change in state) = 0.25102 n_iter = 3 J = 30.539 max(relative change in state) = 0.31394 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 8 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 9 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 9 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 31 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 = 10.229 max(relative change in state) = - n_iter = 2 J = 8.8384 max(relative change in state) = 0.60323 n_iter = 3 J = 8.5874 max(relative change in state) = 0.29299 n_iter = 4 J = 9.7492 max(relative change in state) = 0.32936 n_iter = 5 J = 8.5851 max(relative change in state) = 0.31559 n_iter = 6 J = 8.5850 max(relative change in state) = 0.68829E-03 n_iter = 7 J = 8.5834 max(relative change in state) = 0.26410E-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_refrac): Appending 1DVar retrieval for profile 9 to the file IT-1DVAR-05comp_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 10 from the file ../data/IT-1DVAR-05_y.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_1dvar_refrac): Reading background data for profile 10 from the file ../data/IT-1DVAR-05_b.nc. ... (from ropp_io_read_ncdf_get): 'start_time' and yr/mo/dy/hr/mn/sc/ms timestamps differ by 2.000 seconds (probably a leap-second issue) - using yr/../ms timestamp INFO (from ropp_qc_cutoff): 31 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.553 max(relative change in state) = - n_iter = 2 J = 24.741 max(relative change in state) = 1.0592 n_iter = 3 J = 24.668 max(relative change in state) = 0.19690 n_iter = 4 J = 24.672 max(relative change in state) = 0.23288 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: 4. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 10 to the file IT-1DVAR-05comp_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: 05comp_refrac % Loaded DLM: NCDF. % Compiled module: NCDF_GETVAR. Processing scenario: BG_20090501000049_UNKN_U999_DMI_ % Compiled module: INTERPOL. Processing scenario: BG_20090501000102_UNKN_U999_DMI_ Processing scenario: BG_20090501000152_UNKN_U999_DMI_ Processing scenario: BG_20090501000623_UNKN_U999_DMI_ Processing scenario: BG_20090501000755_UNKN_U999_DMI_ Processing scenario: BG_20090501001552_UNKN_U999_DMI_ Processing scenario: BG_20090501001716_UNKN_U999_DMI_ Processing scenario: BG_20090501002017_UNKN_U999_DMI_ Processing scenario: BG_20090501002135_UNKN_U999_DMI_ Processing scenario: BG_20090501002245_UNKN_U999_DMI_ Found max difference [%] in observation scenarios: 999900.00 Found max difference [K] in temperature scenarios: -1.8317260 Found max difference [g/kg] in specific humidity scenarios: -1.2436550 Found max difference [hPa] in pressure scenarios: -1.0066380 Generating output image file ropp_1dvar_05comp_refrac_bangle.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. *** Finished IDL *** ropp_1dvar_05comp_refrac.jpg can be compared against ../data/ropp_1dvar_05comp_refrac_reference.jpg