*** Results log of t_1dvar_GRAS_05newop_refrac (1DVAR GRAS; new interp) *** ./../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-05newop_refrac.1.nc -new_op -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.208 max(relative change in state) = - n_iter = 2 J = 33.198 max(relative change in state) = 0.45412 n_iter = 3 J = 30.509 max(relative change in state) = 0.32779 n_iter = 4 J = 29.765 max(relative change in state) = 0.30541 n_iter = 5 J = 27.534 max(relative change in state) = 0.40030 n_iter = 6 J = 25.426 max(relative change in state) = 0.15641 n_iter = 7 J = 24.951 max(relative change in state) = 0.22667 n_iter = 8 J = 24.092 max(relative change in state) = 0.19808 n_iter = 9 J = 23.795 max(relative change in state) = 0.23373 n_iter = 10 J = 23.420 max(relative change in state) = 0.15851 n_iter = 11 J = 23.152 max(relative change in state) = 0.17089 n_iter = 12 J = 23.441 max(relative change in state) = 0.14194 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: 13. INFO (from ropp_1dvar_refrac): Writing 1DVar retrieval for profile 1 to the file IT-1DVAR-05newop_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 = 35.854 max(relative change in state) = - n_iter = 2 J = 33.958 max(relative change in state) = 0.38461 n_iter = 3 J = 32.486 max(relative change in state) = 0.45354 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-05newop_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.580 max(relative change in state) = - n_iter = 2 J = 13.859 max(relative change in state) = 0.36004 n_iter = 3 J = 12.351 max(relative change in state) = 0.56532 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-05newop_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.209 max(relative change in state) = - n_iter = 2 J = 29.663 max(relative change in state) = 0.83821 n_iter = 3 J = 29.234 max(relative change in state) = 0.55875 n_iter = 4 J = 29.059 max(relative change in state) = 0.60380 n_iter = 5 J = 28.919 max(relative change in state) = 0.14013 n_iter = 6 J = 28.913 max(relative change in state) = 0.18758E-01 n_iter = 7 J = 28.901 max(relative change in state) = 0.10178 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-05newop_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 = 37.842 max(relative change in state) = - n_iter = 2 J = 34.124 max(relative change in state) = 0.23699 n_iter = 3 J = 32.421 max(relative change in state) = 0.16966 n_iter = 4 J = 31.203 max(relative change in state) = 0.13018 n_iter = 5 J = 30.217 max(relative change in state) = 0.29302 n_iter = 6 J = 29.685 max(relative change in state) = 0.19595 n_iter = 7 J = 28.771 max(relative change in state) = 0.12584 n_iter = 8 J = 28.282 max(relative change in state) = 0.88643E-01 n_iter = 9 J = 28.083 max(relative change in state) = 0.97118E-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-05newop_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 = 137.82 max(relative change in state) = - n_iter = 2 J = 134.51 max(relative change in state) = 0.15655 n_iter = 3 J = 118.61 max(relative change in state) = 1.1945 n_iter = 4 J = 122.70 max(relative change in state) = 1.0784 n_iter = 5 J = 115.62 max(relative change in state) = 1.3212 n_iter = 6 J = 112.37 max(relative change in state) = 0.38700 n_iter = 7 J = 112.81 max(relative change in state) = 1.0324 n_iter = 8 J = 110.80 max(relative change in state) = 0.35631 n_iter = 9 J = 110.36 max(relative change in state) = 0.22384 n_iter = 10 J = 109.44 max(relative change in state) = 0.15704 n_iter = 11 J = 108.94 max(relative change in state) = 0.81437 n_iter = 12 J = 108.00 max(relative change in state) = 0.38431 n_iter = 13 J = 106.70 max(relative change in state) = 0.37572 n_iter = 14 J = 106.33 max(relative change in state) = 0.11696 n_iter = 15 J = 105.43 max(relative change in state) = 0.16335 n_iter = 16 J = 105.68 max(relative change in state) = 0.24524 n_iter = 17 J = 104.70 max(relative change in state) = 0.12529 n_iter = 18 J = 104.56 max(relative change in state) = 0.52418E-01 n_iter = 19 J = 104.08 max(relative change in state) = 0.91167E-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: 19. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 6 to the file IT-1DVAR-05newop_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.211 max(relative change in state) = - n_iter = 2 J = 21.350 max(relative change in state) = 0.80288 n_iter = 3 J = 21.107 max(relative change in state) = 0.26317 n_iter = 4 J = 21.637 max(relative change in state) = 0.25751 n_iter = 5 J = 21.076 max(relative change in state) = 0.19727 n_iter = 6 J = 21.075 max(relative change in state) = 0.54227E-02 n_iter = 7 J = 21.066 max(relative change in state) = 0.57388E-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-05newop_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.492 max(relative change in state) = - n_iter = 2 J = 31.380 max(relative change in state) = 0.29271 n_iter = 3 J = 30.586 max(relative change in state) = 0.33261 n_iter = 4 J = 32.690 max(relative change in state) = 1.4159 n_iter = 5 J = 30.142 max(relative change in state) = 1.0874 n_iter = 6 J = 29.998 max(relative change in state) = 0.90451E-01 n_iter = 7 J = 29.753 max(relative change in state) = 0.17467 n_iter = 8 J = 29.578 max(relative change in state) = 0.30618 n_iter = 9 J = 29.356 max(relative change in state) = 0.12113 n_iter = 10 J = 29.335 max(relative change in state) = 0.25803 n_iter = 11 J = 29.216 max(relative change in state) = 0.11052 n_iter = 12 J = 29.188 max(relative change in state) = 0.39617E-01 n_iter = 13 J = 29.116 max(relative change in state) = 0.10466 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 8 to the file IT-1DVAR-05newop_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.277 max(relative change in state) = - n_iter = 2 J = 8.8818 max(relative change in state) = 0.60416 n_iter = 3 J = 8.6320 max(relative change in state) = 0.29159 n_iter = 4 J = 9.7278 max(relative change in state) = 0.28016 n_iter = 5 J = 8.6300 max(relative change in state) = 0.26938 n_iter = 6 J = 8.6298 max(relative change in state) = 0.82231E-03 n_iter = 7 J = 8.6281 max(relative change in state) = 0.26111E-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-05newop_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.399 max(relative change in state) = - n_iter = 2 J = 24.547 max(relative change in state) = 0.98697 n_iter = 3 J = 24.502 max(relative change in state) = 0.20235 n_iter = 4 J = 24.466 max(relative change in state) = 0.20539 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-05newop_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: 05newop_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.7276414 Found max difference [g/kg] in specific humidity scenarios: -1.3344102 Found max difference [hPa] in pressure scenarios: 1.4187368 Generating output image file ropp_1dvar_05newop_refrac_bangle.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. *** Finished IDL *** ropp_1dvar_05newop_refrac.jpg can be compared against ../data/ropp_1dvar_05newop_refrac_reference.jpg