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-fq24mybyn2li6got2bxzk62ejh5d3p4z/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-nhc23hl,
netcdf-c/4.9.2-xduzp77, netcdf-fortran/4.5.3), glibc/2.35-tdr5czg (required by:
mpich/4.2.1, hdf5/1.14.3-nhc23hl, netcdf-c/4.9.2-xduzp77, 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
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
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 14 -300 3 11 53 219
hmm = 2.4553947299393127E-005
mixt_frac = 0.12496700675912630
precip_frac = 0.74387665293612981
precip_frac_1 = 0.53412687709832540
precip_frac_2 = 0.77383186817637706
hmp2_ip_on_hmm2_ip = 2.3406478066918504
hmm_ip = 3.3008089718204071E-005
hmp2 = 2.1046282162880373E-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.
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 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
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 = 0.0000000000000000
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 = 0.0000000000000000
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 14 -300 3 11 53 219
wm = -2.8921054012137515E-002
rtm = 3.3462654440638020E-003
thlm = 304.14597280234517
wp2 = 0.52942108797810650
rtp2 = 1.2484316455152604E-006
thlp2 = 0.24287476779206951
Skw = -4.2485016832771221
wp3 = -1.6365815428249368
Skrt = 0.32205648018380373
rtp3 = 4.4924080337369112E-010
Skthl = -3.5939232569687620
thlp3 = -0.43017225563491818
wprtp = -7.6798764794273032E-004
wpthlp = -0.22298303005648598
rtpthlp = 4.9558264927973409E-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).
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 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
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 = 0.0000000000000000
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 = 0.0000000000000000
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 14 -300 3 11 53 220
wm = -2.8921054012137515E-002
rtm = 3.3462654440638332E-003
thlm = 304.14597280234517
wp2 = 0.52942108810088739
rtp2 = 1.4637154662308678E-006
thlp2 = 0.24289092427877335
Skw = -3.5567009671613761
wp3 = -1.3700903500028432
Skrt = 1.0250366434754170
rtp3 = 1.8151995855423685E-009
Skthl = -3.1874062976278612
thlp3 = -0.38155254759247115
wprtp = -6.4379619930305580E-004
wpthlp = 6.9158488416989541E-003
rtpthlp = -5.3663189804921256E-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 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
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 = 0.0000000000000000
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
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 14 -300 3 11 53 221
wm = -2.8921054012137515E-002
rtm = 3.3462654440638397E-003
thlm = 304.14597280234517
wp2 = 0.52942108805996047
rtp2 = 1.4329606346346066E-006
thlp2 = 0.24290708027435370
Skw = -3.6555296408928131
wp3 = -1.4081605203242815
Skrt = 4.0709922869339721E-002
rtp3 = 6.9831547996455927E-011
Skthl = 3.2396199844318829
thlp3 = 0.38784154671023996
wprtp = -1.2461620441741854E-004
wpthlp = -0.20748714036149254
rtpthlp = 5.3098191114652252E-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 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
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 = 0.0000000000000000
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 = 0.0000000000000000
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 14 -300 3 11 53 224
wm = -2.8921054012137515E-002
rtm = 3.3462654440635600E-003
thlm = 304.14597280234500
wp2 = 0.52942108695491785
rtp2 = 1.5867347907825447E-006
thlp2 = 0.24346952165227220
Skw = 1.0917474897781334
wp3 = 0.42055621468847426
Skrt = -0.66033905729845754
rtp3 = -1.3198467327954894E-009
Skthl = 2.3589862692871710
thlp3 = 0.28339508687848841
wprtp = 3.6131602342128659E-004
wpthlp = -0.14946071280867557
rtpthlp = -5.5939294239888649E-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 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
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 = 0.0000000000000000
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 = 0.0000000000000000
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 14 -300 3 11 53 226
wm = -2.8921054012137515E-002
rtm = 3.3462654440635465E-003
thlm = 304.14597280234500
wp2 = 0.52942108703677171
rtp2 = 1.5252251280482918E-006
thlp2 = 0.24350183413456405
Skw = 0.89409014231669826
wp3 = 0.34441587404666979
Skrt = -2.6291856175596844
rtp3 = -4.9524736710926646E-009
Skthl = -3.1824689809350675
thlp3 = -0.38239969265944890
wprtp = -7.0299380307869183E-004
wpthlp = 0.27393576654501817
rtpthlp = -5.4847975721970484E-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 = 22.009464516326343
Vertical integral of <rt> after advance_xm_wpxp = 22.281504292772677
Vertical integral of <rt> forcing = -1.8095933255454536E-005
Flux of <rt> at the surface = 9.2489518807652583E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 4.2910689107902855E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 790530.31915425370
Vertical integral of <thl> after advance_xm_wpxp = 790531.38738241000
Vertical integral of <thl> forcing = -1.4937143241047555E-002
Flux of <thl> at the surface = 1.8497903761530515E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 5.3653782805529460E-013
The spurious source of <thl> is within acceptable limits.
Test 3
Vertical integral of <rt> before advance_xm_wpxp = 35.907285164465996
Vertical integral of <rt> after advance_xm_wpxp = 36.136937074056398
Vertical integral of <rt> forcing = -4.8851534487917105E-005
Flux of <rt> at the surface = 8.1435789978928511E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -2.8914316687472796E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 795357.54135273700
Vertical integral of <thl> after advance_xm_wpxp = 795347.05827640486
Vertical integral of <thl> forcing = -5.1230745769885816E-002
Flux of <thl> at the surface = 1.6287157995785700E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 3.1015467971684529E-013
The spurious source of <thl> is within acceptable limits.
Test 4
Vertical integral of <rt> before advance_xm_wpxp = 39.757100997474922
Vertical integral of <rt> after advance_xm_wpxp = 40.040342265118163
Vertical integral of <rt> forcing = -3.7527074329577332E-005
Flux of <rt> at the surface = 9.8166463314032721E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 5.4128793461338809E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 796039.53751379240
Vertical integral of <thl> after advance_xm_wpxp = 796029.86503898213
Vertical integral of <thl> forcing = -5.1874875365658622E-002
Flux of <thl> at the surface = 1.9633292662806546E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.9469911793912331E-012
The spurious source of <thl> is within acceptable limits.
Test 5
Vertical integral of <rt> before advance_xm_wpxp = 30.089412703879969
Vertical integral of <rt> after advance_xm_wpxp = 30.216996946567004
Vertical integral of <rt> forcing = -4.5943301484973959E-005
Flux of <rt> at the surface = 4.7122411044173428E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 2.4699480741935398E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 793474.95864257228
Vertical integral of <thl> after advance_xm_wpxp = 793462.61414223991
Vertical integral of <thl> forcing = -5.0572816651033574E-002
Flux of <thl> at the surface = 9.4244822088346853E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 9.5801144794904758E-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
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.
