Skip to content

Console Output

+ run_scripts/run_G_unit_tests.bash
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.940566852245361E-04
MATLAB result:    3.940566852245361E-04
Percent difference:   1.37569E-14%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x1)
Obtained result:    3.927115990964175E-04
MATLAB result:    3.927115990964175E-04
Percent difference:   0.00000E+00%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x2)
Obtained result:    3.959157231937518E-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.961224408774507E-04
MATLAB result:    3.961224408774508E-04
Percent difference:   2.73704E-14%
 
Agreement with regards to integral:  Trivariate NNL
Obtained result:    5.957899876478482E-10
MATLAB result:    5.957899876478559E-10
Percent difference:   1.30160E-12%
 
Agreement with regards to integral:  Trivariate NNL (const. x1)
Obtained result:    5.299646183560589E-10
MATLAB result:    5.299646183560650E-10
Percent difference:   1.15111E-12%
 
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.286452253001269E-10
MATLAB result:    5.286452253001284E-10
Percent difference:   2.73826E-13%
 
--------------------------------------------------
 
Agreement with regards to integral:  Trivariate NLL
Obtained result:   -9.854231981928350E-04
MATLAB result:   -9.854231981928350E-04
Percent difference:   0.00000E+00%
 
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.940070763287882E-09
MATLAB result:    8.940070763287870E-09
Percent difference:   1.29535E-13%
 
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.007487885941808E-02
MATLAB result:    1.007487885941809E-02
Percent difference:   1.03310E-13%
 
==================================================
 
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals):  2
 
Agreement with regards to integral:  Quadrivariate NNLL
Obtained result:   -4.441972023148828E-05
MATLAB result:   -4.441972023148829E-05
Percent difference:   3.05102E-14%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x1)
Obtained result:   -4.559373570842311E-05
MATLAB result:   -4.559373570842311E-05
Percent difference:   1.48623E-14%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x2)
Obtained result:   -5.040842768062482E-05
MATLAB result:   -5.040842768062482E-05
Percent difference:   1.34427E-14%
 
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:   0.00000E+00%
 
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.622010070476263E-05
MATLAB result:   -7.622010070476262E-05
Percent difference:   1.77808E-14%
 
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.609761356905441E-05
MATLAB result:   -4.609761356905441E-05
Percent difference:   0.00000E+00%
 
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.95210308158732515     
 Nc_in_cloud (from Ncn) =    51650441.529586725     
 Ncm (from Ncn) =    49176544.545665473     
 Ncnm (from Nc_in_cloud) =    50000000.000000000     
 Ncnm (from Ncm) =    49999999.999999993     
 
 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.354576334     
 Ncm (from Ncn) =    10246392.665117815     
 Ncnm (from Nc_in_cloud) =    50000000.000000000     
 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 =         2021          11          16        -360           3          17          16         153        2174         164         169        -207         156         170         169         306        2174         164         169        -207         156         170         169         306        2174         164         169        -207         156         170         169         306        2174
 hmm =    1.7957099143356954E-005
 mixt_frac =    7.1320791205783704E-002
 precip_frac =   0.80525810072937776     
 precip_frac_1 =   0.77557266785875290     
 precip_frac_2 =   0.80753788533052984     
 hmp2_ip_on_hmm2_ip =    3.9768537013877854     
 hmm_ip =    2.2299805648762765E-005
 hmp2 =    1.6704729952467098E-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 =         2021          11          16        -360           3          17          16         156        2177         167         172        -204         159         173         172         312        2177         167         172        -204         159         173         172         312        2177         167         172        -204         159         173         172         312        2177
 wm =  -0.16738295691290073     
 rtm =    1.5381998584300560E-002
 thlm =    304.94033072526753     
 wp2 =   0.62617409560942450     
 rtp2 =    3.2828698400685833E-005
 thlp2 =    1.8809606056952481     
 Skw =  -0.88636461607051620     
 wp3 =  -0.43919264279886033     
 Skrt =   0.13000307838110137     
 rtp3 =    2.4453111735622018E-008
 Skthl =   0.57985419335336330     
 thlp3 =    1.4958510549816599     
 wprtp =    1.6615336650220777E-003
 wpthlp =  -0.59996543004683001     
 rtpthlp =   -7.0722758449195088E-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).
 
 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 =         2021          11          16        -360           3          17          16         157        2178         168         173        -203         160         174         173         314        2178         168         173        -203         160         174         173         314        2178         168         173        -203         160         174         173         314        2178
 wm =   0.35476000262786955     
 rtm =    1.0441122727653919E-002
 thlm =    306.83444787208504     
 wp2 =    6.1253349821658531E-002
 rtp2 =    5.5910304055744853E-006
 thlp2 =    1.5501269987171238     
 Skw =    3.5575365820191323     
 wp3 =    5.3931703864419651E-002
 Skrt =   0.27434600766776696     
 rtp3 =    3.6269070508070095E-009
 Skthl =   0.52521048611358490     
 thlp3 =    1.0136413278610448     
 wprtp =   -8.2303581901907859E-005
 wpthlp =   0.26090832491845911     
 rtpthlp =   -2.6495497386700975E-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.
 
 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 =         2021          11          16        -360           3          17          16         162        2183         173         178        -198         165         179         178         324        2183         173         178        -198         165         179         178         324        2183         173         178        -198         165         179         178         324        2183
 wm =   0.28547614242094532     
 rtm =    9.6467676814046433E-003
 thlm =    307.19105638569096     
 wp2 =   0.61242954493370061     
 rtp2 =    1.4089467226513835E-005
 thlp2 =    2.2364917082403588     
 Skw =    4.2848701672114968     
 wp3 =    2.0536284299473824     
 Skrt =   -4.1821649601739708     
 rtp3 =   -2.2117855468228660E-007
 Skthl =    1.8805215714470980     
 thlp3 =    6.2896902496084133     
 wprtp =   -2.4736864148689423E-003
 wpthlp =   -9.2124355104445790E-002
 rtpthlp =    5.0521174834671258E-003
 
 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 =         2021          11          16        -360           3          17          16         185        2206         196         201        -175         188         202         201         370        2206         196         201        -175         188         202         201         370        2206         196         201        -175         188         202         201         370        2206
 wm =  -0.16085801713286096     
 rtm =    4.5736389836652817E-003
 thlm =    304.57489201994298     
 wp2 =   0.77330572510985141     
 rtp2 =    2.8678217261653444E-006
 thlp2 =    2.2579807134792849     
 Skw =    1.8234135251790538     
 wp3 =    1.2399723758481926     
 Skrt =    2.7143327989875901     
 rtp3 =    1.3182303928755170E-008
 Skthl =    2.7584188191137704     
 thlp3 =    9.3592392487998985     
 wprtp =    1.3323293920402732E-003
 wpthlp =   0.36496138300689579     
 rtpthlp =    2.2902278880880737E-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).
 
 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 =         2021          11          16        -360           3          17          16         198        2219         209         214        -162         201         215         214         396        2219         209         214        -162         201         215         214         396        2219         209         214        -162         201         215         214         396        2219
 wm =   0.18091963321812921     
 rtm =    1.0637584244912080E-002
 thlm =    301.23167229264897     
 wp2 =   0.11087143863066085     
 rtp2 =    2.0756676454667372E-005
 thlp2 =   0.43664702559629054     
 Skw =   -3.6741127285273802     
 wp3 =  -0.13563819486773873     
 Skrt =   -4.3889700076793412     
 rtp3 =   -4.1504896649551598E-007
 Skthl =    1.9059615986393723     
 thlp3 =   0.54993315713943614     
 wprtp =   -5.0429970962015131E-004
 wpthlp =   0.13483773317345746     
 rtpthlp =   -2.7094844968200138E-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).
 
 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.454761250379878     
 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> =    2.4909544912854464E-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 =    18.504037052486161     
 Vertical integral of <rt> after advance_xm_wpxp =    18.528355085615985     
 Vertical integral of <rt> forcing =   -3.2008735805343402E-005
 Flux of <rt> at the surface =    1.1306884623810664E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -1.8600843521704435E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    790366.29120395682     
 Vertical integral of <thl> after advance_xm_wpxp =    790355.92115061113     
 Vertical integral of <thl> forcing =   -3.6828221409155773E-002
 Flux of <thl> at the surface =    2.2613769247621330E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -1.2288781103819701E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            3
 
 Vertical integral of <rt> before advance_xm_wpxp =    29.095372494829498     
 Vertical integral of <rt> after advance_xm_wpxp =    29.260485314512174     
 Vertical integral of <rt> forcing =   -3.8738304535650652E-005
 Flux of <rt> at the surface =    5.8911437014459766E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -2.6698478497455547E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    793014.33036481973     
 Vertical integral of <thl> after advance_xm_wpxp =    793004.79479640105     
 Vertical integral of <thl> forcing =   -4.3567515467102601E-002
 Flux of <thl> at the surface =    1.1782287402891953E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    1.9447846111297906E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            4
 
 Vertical integral of <rt> before advance_xm_wpxp =    17.373843454211215     
 Vertical integral of <rt> after advance_xm_wpxp =    17.547949797083060     
 Vertical integral of <rt> forcing =   -1.1407528393578277E-005
 Flux of <rt> at the surface =    5.9176200463306207E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -1.2180333781516839E-018
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    788983.43860784650     
 Vertical integral of <thl> after advance_xm_wpxp =    788982.25998345262     
 Vertical integral of <thl> forcing =   -1.5763988073085423E-002
 Flux of <thl> at the surface =    1.1835240092661242E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    8.1053913580930725E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            5
 
 Vertical integral of <rt> before advance_xm_wpxp =    23.347581712617110     
 Vertical integral of <rt> after advance_xm_wpxp =    23.336963999484500     
 Vertical integral of <rt> forcing =   -4.5295703294175358E-005
 Flux of <rt> at the surface =    9.9033261854417634E-006
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    3.1536730692172110E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    793122.96396363096     
 Vertical integral of <thl> after advance_xm_wpxp =    793107.55013637478     
 Vertical integral of <thl> forcing =   -5.1577490710995939E-002
 Flux of <thl> at the surface =    1.9806652370883527E-004
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    9.9711905399146872E-015
 
 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!
  
 ==================================================
  
All tests have succeeded.