*** Results log of t_1dvar_01_refrac (1DVAR basic; default) *** ./../tools/ropp_1dvar_refrac -y ../data/IT-1DVAR-01_y.nc -b ../data/IT-1DVAR-01_b.nc --bg-corr ../data/IT-1DVAR-01_c.nc -o IT-1DVAR-01_refrac.1.nc -d --------------------------------------------------------------------- ROPP Refractivity 1D-Var --------------------------------------------------------------------- INFO (from ropp_1dvar_refrac): Reading observation data for profile 1 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 1 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.22872E-20 max(relative change in state) = - n_iter = 2 J = 0.19324E-20 max(relative change in state) = 0.11653E-11 n_iter = 3 J = 0.13371E-20 max(relative change in state) = 0.39790E-11 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: 3. INFO (from ropp_1dvar_refrac): Writing 1DVar retrieval for profile 1 to the file IT-1DVAR-01_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 2 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 2 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.24509E-20 max(relative change in state) = - n_iter = 2 J = 0.18491E-20 max(relative change in state) = 0.11653E-11 n_iter = 3 J = 0.13755E-20 max(relative change in state) = 0.28706E-11 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: 3. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 2 to the file IT-1DVAR-01_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 3 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 3 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.23332E-20 max(relative change in state) = - n_iter = 2 J = 0.18772E-20 max(relative change in state) = 0.13074E-11 n_iter = 3 J = 0.11751E-20 max(relative change in state) = 0.29416E-11 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: 3. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 3 to the file IT-1DVAR-01_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 4 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 4 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.23553E-20 max(relative change in state) = - n_iter = 2 J = 0.19606E-20 max(relative change in state) = 0.12506E-11 n_iter = 3 J = 0.13810E-20 max(relative change in state) = 0.31832E-11 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: 3. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 4 to the file IT-1DVAR-01_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 5 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 5 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.46577E-20 max(relative change in state) = - n_iter = 2 J = 0.37583E-20 max(relative change in state) = 0.21032E-11 n_iter = 3 J = 0.26742E-20 max(relative change in state) = 0.34674E-11 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: 3. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 5 to the file IT-1DVAR-01_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 6 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 6 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.23260E-20 max(relative change in state) = - n_iter = 2 J = 0.19549E-20 max(relative change in state) = 0.13074E-11 n_iter = 3 J = 0.12434E-20 max(relative change in state) = 0.34106E-11 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: 3. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 6 to the file IT-1DVAR-01_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 7 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 7 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.22167E-20 max(relative change in state) = - n_iter = 2 J = 0.18528E-20 max(relative change in state) = 0.11653E-11 n_iter = 3 J = 0.11485E-20 max(relative change in state) = 0.42348E-11 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: 3. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 7 to the file IT-1DVAR-01_refrac.1.nc. INFO (from ropp_1dvar_refrac): Reading observation data for profile 8 from the file ../data/IT-1DVAR-01_y.nc. INFO (from ropp_1dvar_refrac): Reading background data for profile 8 from the file ../data/IT-1DVAR-01_b.nc. 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 = 0.22525E-20 max(relative change in state) = - n_iter = 2 J = 0.17474E-20 max(relative change in state) = 0.11084E-11 n_iter = 3 J = 0.11693E-20 max(relative change in state) = 0.38085E-11 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: 3. INFO (from ropp_1dvar_refrac): Appending 1DVar retrieval for profile 8 to the file IT-1DVAR-01_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: 01_refrac % Loaded DLM: NCDF. % Compiled module: NCDF_GETVAR. Processing scenario: OC_20000101000000_9999_U999_METO % Compiled module: INTERPOL. Processing scenario: OC_20000101000000_9999_U999_METO Processing scenario: OC_20000101000000_9999_U999_METO Processing scenario: OC_20000101000000_9999_U999_METO Processing scenario: OC_20000101000000_9999_U999_METO Processing scenario: OC_20000101000000_9999_U999_METO Processing scenario: OC_20000101000000_9999_U999_METO Processing scenario: OC_20000101000000_9999_U999_METO Found max difference [%] in observation scenarios: -2.5495935e+11 Found max difference [K] in temperature scenarios: -1.1539214e-11 Found max difference [g/kg] in specific humidity scenarios: -1.6591173e-12 Found max difference [hPa] in pressure scenarios: 3.4106051e-12 Generating output image file ropp_1dvar_01_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_01_refrac.jpg can be compared against ../data/ropp_1dvar_01_refrac_reference.jpg