*** Results log of t_1dvar_GRAS_05_refrac (1DVAR GRAS; default) *** ./../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-05_refrac.1.nc -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.261 max(relative change in state) = - n_iter = 2 J = 33.255 max(relative change in state) = 0.45574 n_iter = 3 J = 30.546 max(relative change in state) = 0.32552 n_iter = 4 J = 29.733 max(relative change in state) = 0.30806 n_iter = 5 J = 27.661 max(relative change in state) = 0.43586 n_iter = 6 J = 24.994 max(relative change in state) = 0.26013 n_iter = 7 J = 24.612 max(relative change in state) = 0.15771 n_iter = 8 J = 23.993 max(relative change in state) = 0.14904 n_iter = 9 J = 23.682 max(relative change in state) = 0.26229 n_iter = 10 J = 23.438 max(relative change in state) = 0.22042 n_iter = 11 J = 23.163 max(relative change in state) = 0.88407E-01 n_iter = 12 J = 23.060 max(relative change in state) = 0.16648 n_iter = 13 J = 22.904 max(relative change in state) = 0.43429E-01 n_iter = 14 J = 22.664 max(relative change in state) = 0.16437 n_iter = 15 J = 22.521 max(relative change in state) = 0.15812 n_iter = 16 J = 22.600 max(relative change in state) = 0.31883 n_iter = 17 J = 22.457 max(relative change in state) = 0.21146 n_iter = 18 J = 22.412 max(relative change in state) = 0.37856E-01 n_iter = 19 J = 22.397 max(relative change in state) = 0.36141E-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): Writing 1DVar retrieval for profile 1 to the file IT-1DVAR-05_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.250 max(relative change in state) = - n_iter = 2 J = 33.395 max(relative change in state) = 0.38006 n_iter = 3 J = 31.951 max(relative change in state) = 0.42285 n_iter = 4 J = 33.624 max(relative change in state) = 0.77864 n_iter = 5 J = 31.884 max(relative change in state) = 0.74049 n_iter = 6 J = 31.831 max(relative change in state) = 0.19729E-01 n_iter = 7 J = 31.516 max(relative change in state) = 0.40648 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-05_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.538 max(relative change in state) = - n_iter = 2 J = 13.820 max(relative change in state) = 0.35951 n_iter = 3 J = 12.318 max(relative change in state) = 0.56305 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-05_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.264 max(relative change in state) = - n_iter = 2 J = 29.713 max(relative change in state) = 0.84063 n_iter = 3 J = 29.285 max(relative change in state) = 0.55684 n_iter = 4 J = 29.113 max(relative change in state) = 0.60296 n_iter = 5 J = 28.972 max(relative change in state) = 0.14102 n_iter = 6 J = 28.967 max(relative change in state) = 0.18731E-01 n_iter = 7 J = 28.954 max(relative change in state) = 0.10166 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-05_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.913 max(relative change in state) = - n_iter = 2 J = 34.188 max(relative change in state) = 0.23822 n_iter = 3 J = 32.534 max(relative change in state) = 0.16404 n_iter = 4 J = 31.294 max(relative change in state) = 0.13750 n_iter = 5 J = 30.321 max(relative change in state) = 0.28835 n_iter = 6 J = 29.785 max(relative change in state) = 0.19816 n_iter = 7 J = 28.770 max(relative change in state) = 0.11817 n_iter = 8 J = 28.179 max(relative change in state) = 0.10136 n_iter = 9 J = 27.956 max(relative change in state) = 0.14192 n_iter = 10 J = 27.799 max(relative change in state) = 0.54229E-01 n_iter = 11 J = 27.735 max(relative change in state) = 0.16386E-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: 11. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 5 to the file IT-1DVAR-05_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.25 max(relative change in state) = - n_iter = 2 J = 134.04 max(relative change in state) = 0.15261 n_iter = 3 J = 118.20 max(relative change in state) = 1.1976 n_iter = 4 J = 122.03 max(relative change in state) = 1.0678 n_iter = 5 J = 115.53 max(relative change in state) = 1.2998 n_iter = 6 J = 112.59 max(relative change in state) = 0.33902 n_iter = 7 J = 112.27 max(relative change in state) = 1.0733 n_iter = 8 J = 110.58 max(relative change in state) = 0.30820 n_iter = 9 J = 109.73 max(relative change in state) = 0.33449 n_iter = 10 J = 108.81 max(relative change in state) = 0.11501 n_iter = 11 J = 108.27 max(relative change in state) = 0.38067 n_iter = 12 J = 107.22 max(relative change in state) = 0.24395 n_iter = 13 J = 106.03 max(relative change in state) = 0.88194E-01 n_iter = 14 J = 105.04 max(relative change in state) = 0.11162 n_iter = 15 J = 105.31 max(relative change in state) = 0.21578 n_iter = 16 J = 104.19 max(relative change in state) = 0.11964 n_iter = 17 J = 103.38 max(relative change in state) = 0.18118 n_iter = 18 J = 103.22 max(relative change in state) = 0.46415E-01 n_iter = 19 J = 103.04 max(relative change in state) = 0.78855E-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-05_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.331 max(relative change in state) = - n_iter = 2 J = 21.447 max(relative change in state) = 0.80049 n_iter = 3 J = 21.200 max(relative change in state) = 0.26455 n_iter = 4 J = 21.758 max(relative change in state) = 0.26259 n_iter = 5 J = 21.170 max(relative change in state) = 0.20283 n_iter = 6 J = 21.168 max(relative change in state) = 0.53533E-02 n_iter = 7 J = 21.160 max(relative change in state) = 0.57471E-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-05_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.166 max(relative change in state) = - n_iter = 2 J = 31.097 max(relative change in state) = 0.27418 n_iter = 3 J = 30.259 max(relative change in state) = 0.31122 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-05_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.296 max(relative change in state) = - n_iter = 2 J = 8.9038 max(relative change in state) = 0.60779 n_iter = 3 J = 8.6610 max(relative change in state) = 0.28785 n_iter = 4 J = 9.8471 max(relative change in state) = 0.35283 n_iter = 5 J = 8.6586 max(relative change in state) = 0.33731 n_iter = 6 J = 8.6585 max(relative change in state) = 0.70519E-03 n_iter = 7 J = 8.6567 max(relative change in state) = 0.25903E-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-05_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.406 max(relative change in state) = - n_iter = 2 J = 24.558 max(relative change in state) = 1.0160 n_iter = 3 J = 24.526 max(relative change in state) = 0.20712 n_iter = 4 J = 24.463 max(relative change in state) = 0.20330 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-05_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: 05_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.7358206 Found max difference [g/kg] in specific humidity scenarios: -1.1631501 Found max difference [hPa] in pressure scenarios: 1.3071630 Generating output image file ropp_1dvar_05_refrac_bangle.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. *** Finished IDL *** ropp_1dvar_05_refrac.jpg can be compared against ../data/ropp_1dvar_05_refrac_reference.jpg