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
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 16 -300 11 48 42 987
hmm = 2.4553410857589477E-005
mixt_frac = 0.12497620673580050
precip_frac = 0.38840210129683761
precip_frac_1 = 0.51210439675456254
precip_frac_2 = 0.37073418898425220
hmp2_ip_on_hmm2_ip = 0.85989614437528905
hmm_ip = 6.3216472762654937E-005
hmp2 = 2.2840236660307533E-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).
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
#0 0x746cef44251f in ???
at ./signal/../sysdeps/unix/sysv/linux/x86_64/libc_sigaction.c:0
#1 0x4b5f58 in adg1_adg2_responder_params
at /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/adg1_adg2_3d_luhar_pdf.F90:1190
#2 0x4caf6e in __adg1_adg2_3d_luhar_pdf_MOD_adg1_pdf_driver
at /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/adg1_adg2_3d_luhar_pdf.F90:176
#3 0x4261db in __pdf_parameter_tests_MOD_pdf_parameter_unit_tests
at /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/pdf_parameter_tests.F90:1616
#4 0x403f2e in g_unit_tests
at /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_tests.F90:267
#5 0x4040da in main
at /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_tests.F90:118
Running /home/jenkins/workspace/clubb_G_unit_gfortran_tests/run_scripts/../install/latest/G_unit_tests ...
Error: A G unit test has failed.
