Console Output
+ source /etc/profile.d/larson-group.sh
++ export GIT_EDITOR=vi
++ GIT_EDITOR=vi
++ export SVN_EDITOR=vi
++ SVN_EDITOR=vi
++ export OMP_STACKSIZE=1048579
++ OMP_STACKSIZE=1048579
++ export LMOD_ROOT=/opt/lmod/
++ LMOD_ROOT=/opt/lmod/
++ source /opt/lmod//lmod/lmod/init/bash
+++ '[' -z '' ']'
+++ case "$-" in
+++ __lmod_vx=x
+++ '[' -n x ']'
+++ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for this output (/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-12.2.0/lmod-8.7.37-qorv6nkkjzyzntx6wbjupmoxrot26t6r/lmod/lmod/init/bash)
Shell debugging restarted
+++ unset __lmod_vx
+++ find /usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core -print -quit
++ export MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
++ MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
+ module load gcc netcdf-fortran
+ '[' -z '' ']'
+ case "$-" in
+ __lmod_sh_dbg=x
+ '[' -n x ']'
+ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for Lmod's output
Lmod Warning:
-------------------------------------------------------------------------------
The following dependent module(s) are not currently loaded:
gcc-runtime/13.1.0-wiuv7wc (required by: mpich/4.2.1, hdf5/1.14.3-cb5lzz7,
netcdf-c/4.9.2, netcdf-fortran/4.5.3), glibc/2.35-tdr5czg (required by:
mpich/4.2.1, hdf5/1.14.3-cb5lzz7, netcdf-c/4.9.2, netcdf-fortran/4.5.3)
-------------------------------------------------------------------------------
Shell debugging restarted
+ unset __lmod_sh_dbg
+ return 0
+ python run_scripts/run_G_unit_tests.py
Running G_unit_tests
Performing KK integrals tests
Percent Difference Tolerance (maximum acceptable percent difference between the CLUBB results and the MATLAB results using the same equations): 1.50000E-12
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 1
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: 3.940566852245362E-04
MATLAB result: 3.940566852245361E-04
Percent difference: 1.37569E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: 3.927115990964176E-04
MATLAB result: 3.927115990964175E-04
Percent difference: 2.76081E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: 3.959157231937519E-04
MATLAB result: 3.959157231937519E-04
Percent difference: 1.36923E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: 3.961224408774508E-04
MATLAB result: 3.961224408774508E-04
Percent difference: 1.36852E-14%
Agreement with regards to integral: Trivariate NNL
Obtained result: 5.957899876478518E-10
MATLAB result: 5.957899876478559E-10
Percent difference: 6.94188E-13%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 5.299646183560617E-10
MATLAB result: 5.299646183560650E-10
Percent difference: 6.24329E-13%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.543555152546510E-10
MATLAB result: 5.543555152546526E-10
Percent difference: 2.98430E-13%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.286452253001267E-10
MATLAB result: 5.286452253001284E-10
Percent difference: 3.12944E-13%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -9.854231981928352E-04
MATLAB result: -9.854231981928350E-04
Percent difference: 2.20048E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549015E-04
MATLAB result: -9.922448817549015E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 8.940070763287877E-09
MATLAB result: 8.940070763287870E-09
Percent difference: 7.40201E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 8.942709549399747E-09
MATLAB result: 8.942709549399743E-09
Percent difference: 3.69991E-14%
Agreement with regards to integral: Bivariate LL
Obtained result: 1.007487885941809E-02
MATLAB result: 1.007487885941809E-02
Percent difference: 1.72183E-14%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 2
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -4.441972023148827E-05
MATLAB result: -4.441972023148829E-05
Percent difference: 4.57652E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -4.559373570842309E-05
MATLAB result: -4.559373570842311E-05
Percent difference: 4.45868E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -5.040842768062482E-05
MATLAB result: -5.040842768062482E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -5.038775591225493E-05
MATLAB result: -5.038775591225493E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.629614112194526E-09
MATLAB result: 4.629614112194527E-09
Percent difference: 1.78672E-14%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.604609365066508E-09
MATLAB result: 4.604609365066509E-09
Percent difference: 1.79642E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 4.685272472815394E-09
MATLAB result: 4.685272472815394E-09
Percent difference: 1.76549E-14%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 4.683452323676534E-09
MATLAB result: 4.683452323676534E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -1.088125285831538E-04
MATLAB result: -1.088125285831538E-04
Percent difference: 2.49099E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549016E-05
MATLAB result: -9.922448817549016E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 7.907812698669825E-10
MATLAB result: 7.907812698669821E-10
Percent difference: 5.23015E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 6.330953526469324E-10
MATLAB result: 6.330953526469324E-10
Percent difference: 0.00000E+00%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 3
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -7.622010070476262E-05
MATLAB result: -7.622010070476262E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -7.695119321547280E-05
MATLAB result: -7.695119321547280E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.851589269710427E-09
MATLAB result: 4.851589269710427E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.837259568295437E-09
MATLAB result: 4.837259568295437E-09
Percent difference: 1.71002E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -4.609761356905442E-05
MATLAB result: -4.609761356905441E-05
Percent difference: 1.46998E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 3.254808634091251E-10
MATLAB result: 3.254808634091251E-10
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
==================================================
Success!
==================================================
Performing correlation cholesky matrix setup tests
correlation matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
--
correlation matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 1.0000000000000000 0.86602540378443860
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
--
==================================================
Performing correlation matrix approximation tests
correlation cholesky matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.86602540378443860 1.0000000000000000 0.96650635094610959
0.86602540378443860 0.96650635094610959 1.0000000000000000
--
==================================================
Success!
==================================================
Performing hole_filling_hm_one_lev_tests
Hole-filling was conservative
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: Hole filling failed! A hole could not be filled.
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
==================================================
Performing fill_holes_hydromet_tests
testset_in(:, 1 ) = 1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 6.0000000000000000
testset_in(:, 2 ) = 6.0000000000000000 5.0000000000000000 4.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000
testset_in(:, 3 ) = 3.0000000000000000 3.0000000000000000 3.0000000000000000 4.0000000000000000 4.0000000000000000 4.0000000000000000
testset_in(:, 4 ) = 1.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 5 ) = 1.0000000000000000 -2.0000000000000000 -4.0000000000000000 -3.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 6 ) = 2.0000000000000000 -1.0000000000000000 -1.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 7 ) = 1.0000000000000000 4.0000000000000000 2.0000000000000000 3.0000000000000000 -1.0000000000000000 0.0000000000000000
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: Hole filling failed! A hole could not be filled.
Warning: One-level hole filling was not successful! total_mass ~= 0
Hole-filling was conservative
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
Nc-Ncn equations unit test
==========================
Iteration = 1
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.15865525393145707
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 4.7896918412674905E-002
Nc_in_cloud (from Ncn) = 66978771.649335034
Ncm (from Ncn) = 3208076.7610693830
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 1 is a success!
======================================================
Iteration = 2
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.0000000000000000
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 0.0000000000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 2 is a success!
======================================================
Iteration = 3
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 0.99957093966680322
cloud fraction = 0.95210308158732504
Nc_in_cloud (from Ncn) = 51650441.529586732
Ncm (from Ncn) = 49176544.545665480
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 3 is a success!
======================================================
Iteration = 4
mu_chi_1 = 1.0000000000000000E-008
mu_chi_2 = -1.0000000000000000E-008
sigma_chi_1 = 5.0000000000000002E-005
sigma_chi_2 = 5.0000000000000002E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.10000000000000001
mixture fraction = 0.50000000000000000
cloud fraction (1st PDF component) = 0.50007978845554835
cloud fraction (2nd PDF component) = 0.49992021154445165
cloud fraction = 0.50000000000000000
Nc_in_cloud (from Ncn) = 53982783.648312390
Ncm (from Ncn) = 26991391.824156195
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 4 is a success!
======================================================
Iteration = 5
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 9.9999999999999995E-007
sigma_chi_2 = 1.5000000000000000E-006
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.50000000000000000
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 1.0000000000000000
cloud fraction = 1.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 50000000.000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 5 is a success!
======================================================
Iteration = 6
mu_chi_1 = -9.9999999999999995E-008
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.0000000000000000
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.49601064368536840
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 9.9545377003631155E-002
Nc_in_cloud (from Ncn) = 49999999.999999993
Ncm (from Ncn) = 4977268.8501815572
Ncnm (from Nc_in_cloud) = 49999999.999999993
Ncnm (from Ncm) = 49999999.999999993
Test iteration 6 is a success!
======================================================
Iteration = 7
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 0.0000000000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.48670438618290790
cloud fraction = 0.49732907640173213
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 24866453.820086606
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 7 is a success!
======================================================
Iteration = 8
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.10796556745540581
Nc_in_cloud (from Ncn) = 81037694.237492040
Ncm (from Ncn) = 8749280.6436284967
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 8 is a success!
======================================================
Iteration = 9
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.20034324826655747
Nc_in_cloud (from Ncn) = 52151569.416871324
Ncm (from Ncn) = 10448214.819174858
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 9 is a success!
======================================================
Iteration = 10
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.16861219748026607
Nc_in_cloud (from Ncn) = 60768988.354576342
Ncm (from Ncn) = 10246392.665117815
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 10 is a success!
======================================================
Success!
Performing corr_mtx_test
==================================================
Expected:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Actual:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Success!
==================================================
==================================================
Performing SILHS category test
==================================================
Success!
Performing hydrometeor mu and sigma values unit test
====================================================
PDF parameter set 1:
hmm = 2.0000000000000002E-005
mixt_frac = 2.7000000000000000E-002
precip_frac = 8.0000000000000002E-002
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 5.4470709146968145E-002
hmp2_ip_on_hmm2_ip = 3.6000000000000001
hmm_ip = 2.5000000000000001E-004
hmp2 = 2.2599999999999997E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 2:
hmm = 5.0000000000000004E-006
mixt_frac = 0.50000000000000000
precip_frac = 0.50000000000000000
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 0.50000000000000000
hmp2_ip_on_hmm2_ip = 1.0000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 7.5000000000000012E-011
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 3:
hmm = 1.0000000000000000E-010
mixt_frac = 0.25000000000000000
precip_frac = 0.40000000000000002
precip_frac_1 = 0.75000000000000000
precip_frac_2 = 0.28333333333333338
hmp2_ip_on_hmm2_ip = 0.50000000000000000
hmm_ip = 2.5000000000000002E-010
hmp2 = 2.7500000000000010E-020
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 4:
hmm = 4.0000000000000003E-005
mixt_frac = 1.0000000000000000E-002
precip_frac = 0.10000000000000001
precip_frac_1 = 1.0000000000000000E-002
precip_frac_2 = 0.10090909090909091
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 5:
hmm = 4.0000000000000003E-005
mixt_frac = 0.98999999999999999
precip_frac = 0.10000000000000001
precip_frac_1 = 0.10090909090909091
precip_frac_2 = 1.0000000000000278E-002
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 6:
hmm = 1.0000000000000001E-005
mixt_frac = 0.12500000000000000
precip_frac = 1.0000000000000000
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 1.0000000000000000
hmp2_ip_on_hmm2_ip = 1.5000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 1.5000000000000005E-010
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 7:
hmm = 6.0000000000000002E-005
mixt_frac = 5.0000000000000003E-002
precip_frac = 0.11000000000000000
precip_frac_1 = 0.90000000000000002
precip_frac_2 = 6.8421052631578952E-002
hmp2_ip_on_hmm2_ip = 5.0000000000000000
hmm_ip = 5.4545454545454548E-004
hmp2 = 1.9276363636363637E-007
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 8:
hmm = 7.9999999999999996E-007
mixt_frac = 0.80000000000000004
precip_frac = 0.25000000000000000
precip_frac_1 = 0.10000000000000001
precip_frac_2 = 0.85000000000000009
hmp2_ip_on_hmm2_ip = 0.75000000000000000
hmm_ip = 3.1999999999999999E-006
hmp2 = 3.8399999999999998E-012
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 9:
hmm = 1.3499999999999999E-005
mixt_frac = 0.10000000000000001
precip_frac = 0.10000000000000001
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 5.5555555555555559E-002
hmp2_ip_on_hmm2_ip = 2.0000000000000000
hmm_ip = 1.3499999999999997E-004
hmp2 = 5.2852499999999985E-009
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 10 (randomly generated):
Random seed values = 2026 3 18 -300 20 15 53 785
hmm = 2.4554684906872281E-005
mixt_frac = 0.12497883530071006
precip_frac = 0.99339786930577845
precip_frac_1 = 0.98963591167507392
precip_frac_2 = 0.99393518783257961
hmp2_ip_on_hmm2_ip = 0.59921612332917074
hmm_ip = 2.4717875551748427E-005
hmp2 = 3.6769511648006417E-010
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Checking values with precipitation in 1st PDF component only.
Checking values with precipitation in 2nd PDF component only.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for ADG1
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 9.0000000000000012E-011
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 3 18 -300 20 15 53 786
wm = -2.8906301862554384E-002
rtm = 3.3465068428716140E-003
thlm = 309.84909780561730
wp2 = 0.80380702116872071
rtp2 = 2.0209825723327141E-007
thlp2 = 0.23911341863092139
Skw = -2.3778906195294334
wp3 = -1.7136399044038948
Skrt = -2.0719127568698159
rtp3 = -1.8824147495421997E-010
Skthl = -0.85808862414233111
thlp3 = -0.10033167610136987
wprtp = 2.2802959790051166E-004
wpthlp = -0.16942099767961608
rtpthlp = -1.9784531382113419E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
Running ADG1 expected-failure tests.
ADG1_w_closure correctly failed when given wp2 = 0.
ADG1_pdf_driver correctly failed when given thlp2 = 0.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for LY93
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 9.0000000000000012E-011
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 10 (randomly generated):
Random seed values = 2026 3 18 -300 20 15 53 786
wm = -2.8906301862554384E-002
rtm = 3.3465068428716140E-003
thlm = 309.84909780561730
wp2 = 0.80380702116872071
rtp2 = 2.0209825723327141E-007
thlp2 = 0.23911341863092139
Skw = -2.3778906195294334
wp3 = -1.7136399044038948
Skrt = -2.0719127568698159
rtp3 = -1.8824147495421997E-010
Skthl = -0.85808862414233111
thlp3 = -0.10033167610136987
wprtp = 2.2802959790051166E-004
wpthlp = -0.16942099767961608
rtpthlp = -1.9784531382113419E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for TSDADG
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 9.0000000000000012E-011
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
PDF parameter set 10 (randomly generated):
Random seed values = 2026 3 18 -300 20 15 53 787
wm = -2.8906301862554384E-002
rtm = 3.3465068428716210E-003
thlm = 309.84909780561730
wp2 = 0.80380702120964831
rtp2 = 1.7133898841338700E-007
thlp2 = 0.23912944914351975
Skw = -2.4767192925234602
wp3 = -1.7848613292562781
Skrt = 4.2561639637362614
rtp3 = 3.0185797066487199E-010
Skthl = -3.4310627201898121
thlp3 = -0.40121600773780713
wprtp = -3.1394891300777310E-004
wpthlp = 9.3899179137538175E-002
rtpthlp = -1.8217431298633166E-004
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 9.0000000000000012E-011
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 3 18 -300 20 15 53 792
wm = -2.8906301862554384E-002
rtm = 3.3465068428715763E-003
thlm = 309.84909780561736
wp2 = 0.80380702141428739
rtp2 = 5.0969094543211094E-007
thlp2 = 0.23920934828489249
Skw = -2.9706695428704157
wp3 = -2.1408292845699743
Skrt = -2.3528730746127540
rtp3 = -8.5616777509903251E-010
Skthl = -3.8157518655128317
thlp3 = -0.44642380728519265
wprtp = 3.9503306360634213E-004
wpthlp = 0.18563237301295521
rtpthlp = 3.1425705957966328E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new hybrid PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 9.0000000000000012E-011
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 3 18 -300 20 15 53 795
wm = -2.8906301862554384E-002
rtm = 3.3465068428715720E-003
thlm = 309.84909780561736
wp2 = 0.80380702153706884
rtp2 = 4.1741313897245653E-007
thlp2 = 0.23925769127978236
Skw = -3.2673486823749500
wp3 = -2.3546327327183381
Skrt = -1.3684497952188313
rtp3 = -3.6904390908143722E-010
Skthl = -0.84726995844017061
thlp3 = -9.9156381441434308E-002
wprtp = 2.3398100831989834E-004
wpthlp = -0.23015283845908227
rtpthlp = -2.8441894485375176E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing spurious source unit test
====================================
Test 1
Vertical integral of <rt> before advance_xm_wpxp = 28.316670259590918
Vertical integral of <rt> after advance_xm_wpxp = 28.454761250379875
Vertical integral of <rt> forcing = -3.9696697370157397E-005
Flux of <rt> at the surface = 5.0000000000000001E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 1.3037531124138191E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 793918.42642888101
Vertical integral of <thl> after advance_xm_wpxp = 793908.03522151255
Vertical integral of <thl> forcing = -4.4637357896125940E-002
Flux of <thl> at the surface = 1.0000000000000000E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.2514503322513804E-012
The spurious source of <thl> is within acceptable limits.
Test 2
Vertical integral of <rt> before advance_xm_wpxp = 40.649314645548607
Vertical integral of <rt> after advance_xm_wpxp = 40.698949526677758
Vertical integral of <rt> forcing = -7.4112518778961763E-005
Flux of <rt> at the surface = 2.3956212254280064E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -1.6940658945086007E-018
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 798859.44169500307
Vertical integral of <thl> after advance_xm_wpxp = 798835.94851462403
Vertical integral of <thl> forcing = -8.3101843714126294E-002
Flux of <thl> at the surface = 4.7912424508560125E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -1.9197143874549738E-013
The spurious source of <thl> is within acceptable limits.
Test 3
Vertical integral of <rt> before advance_xm_wpxp = 31.009437736985099
Vertical integral of <rt> after advance_xm_wpxp = 31.279950090205219
Vertical integral of <rt> forcing = -3.3005391628401000E-005
Flux of <rt> at the surface = 9.3471323569541890E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 4.8409627001477773E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 795646.21794054087
Vertical integral of <thl> after advance_xm_wpxp = 795641.60150694009
Vertical integral of <thl> forcing = -3.4082376715235507E-002
Flux of <thl> at the surface = 1.8694264713908379E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -1.2685130723610882E-012
The spurious source of <thl> is within acceptable limits.
Test 4
Vertical integral of <rt> before advance_xm_wpxp = 24.372996085230319
Vertical integral of <rt> after advance_xm_wpxp = 24.392615086757392
Vertical integral of <rt> forcing = -4.6350805137021549E-005
Flux of <rt> at the surface = 1.1174747689390508E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 2.4137050864958542E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 791668.49656725035
Vertical integral of <thl> after advance_xm_wpxp = 791653.81938223296
Vertical integral of <thl> forcing = -5.1158899595623988E-002
Flux of <thl> at the surface = 2.2349495378781014E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -1.9653723093426834E-013
The spurious source of <thl> is within acceptable limits.
Test 5
Vertical integral of <rt> before advance_xm_wpxp = 32.146406578019139
Vertical integral of <rt> after advance_xm_wpxp = 32.288881886416320
Vertical integral of <rt> forcing = -4.0609056079218530E-005
Flux of <rt> at the surface = 5.1552675073658096E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -9.3065203980724487E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 794007.22380042728
Vertical integral of <thl> after advance_xm_wpxp = 793994.57675988507
Vertical integral of <thl> forcing = -5.2467336821761558E-002
Flux of <thl> at the surface = 1.0310535014731620E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -3.5375868900899832E-013
The spurious source of <thl> is within acceptable limits.
Success!
Running tuner tests
-------------------
goldstein_price pass
rastrigin pass
himmelblau pass
eggholder pass
schaffer pass
==================================================
Performing 'mean vertical velocity in clouds' tests
Success!
==================================================
==================================================
Performing 'Smooth Heaviside Peskin' tests
Case 1: heaviside_smth_range = 0, smooth Heaviside collapses to Heaviside
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
Case 2: heaviside_smth_range = 2.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Expected outcome: 0.0000000000000000 9.0845056908100005E-002 0.90915494309000000 1.0000000000000000
True outcome: 0.0000000000000000 9.0845056909999999E-002 0.90915494309000000 1.0000000000000000
Success!
==================================================
==================================================
Performing 'Smooth Min / Max' tests
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
(Simple) min:
Expected outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
True outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
(Simple) max:
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Smooth min:
Expected outcome: -2.2247448714824141 -1.3660254037844388 -0.36602540378443871 -0.22474487138241406
True outcome: -2.2247448714824136 -1.3660254037844388 -0.36602540378443871 -0.22474487138241384
Smooth max:
Expected outcome: 0.22474487138241406 0.36602540378443871 1.3660254037844388 2.2247448714824141
True outcome: 0.22474487138241384 0.36602540378443871 1.3660254037844388 2.2247448714824136
Testing smooth min and max with min_max_smth_mag=1e-7, input_var1 = equidistant grid on [-1, 1), input_var2 = 0
Input: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
True outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
Smooth max:
Expected outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
True outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
Fatal error in CLUBB - MPI Process 1 / Chunk index 1
Column index not available. Latitude/Longitude for column 1: ******/******
Latitude range from -90 to 90 degrees, Longitude range from -180 to 180 degrees
Error in ADG1_w_closure: wp2 must not be less than the w tolerance value.
Fatal error in CLUBB - MPI Process 1 / Chunk index 1
Column index not available. Latitude/Longitude for column 1: ******/******
Latitude range from -90 to 90 degrees, Longitude range from -180 to 180 degrees
Error in ADG1_ADG2_responder_params: xp2 must not be less than its tolerance value.
Input_var1: -1000000.0000000000 -10000000.000000000 -100000000.00000000 -1000000000.0000000
Input_var2: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
True outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
Smooth max:
Expected outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
True outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
Input_var1: -1000000.0000000000
Input_var2: 1.0000000000000000
Smooth min:
Expected outcome: -1000000.0000002501
True outcome: -1000000.0000002501
Smooth max:
Expected outcome: 1.0000002500019036
True outcome: 1.0000002500019036
Success!
==================================================
====================== Running fill holes tests ======================
-------------- easy_fill_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: easy_fill_test
-----------------------------------------------
-------------- below_thresh_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: below_thresh_test
-----------------------------------------------
====================== End fill holes tests ======================
Running /home/jenkins/workspace/clubb_G_unit_gfortran_tests/run_scripts/../install/latest/G_unit_tests ...
All tests have succeeded.
