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.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.354576334     
 Ncm (from Ncn) =    10246392.665117813     
 Ncnm (from Nc_in_cloud) =    50000000.000000000     
 Ncnm (from Ncm) =    50000000.000000000     
 
 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           6          29        -300           3          17           6          37        2058          43          66        -263          40          54          43          74        2058          43          66        -263          40          54          43          74        2058          43          66        -263          40          54          43          74        2058
 hmm =    3.6736151266325510E-005
 mixt_frac =   0.64568729952124015     
 precip_frac =   0.88860834007141021     
 precip_frac_1 =   0.92444578204935424     
 precip_frac_2 =   0.82329941634042203     
 hmp2_ip_on_hmm2_ip =    4.6224196234685282     
 hmm_ip =    4.1341218183224935E-005
 hmp2 =    7.1893208293127386E-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           6          29        -300           3          17           6          40        2061          46          69        -260          43          57          46          80        2061          46          69        -260          43          57          46          80        2061          46          69        -260          43          57          46          80        2061
 wm =    1.3990036430270858E-002
 rtm =    6.4684105990897214E-003
 thlm =    292.40606606552501     
 wp2 =   0.46049349775869042     
 rtp2 =    4.0352919802364656E-006
 thlp2 =    1.3463356663843720     
 Skw =  -0.85249615840124526     
 wp3 =  -0.26639599049392881     
 Skrt =   -3.8131512876811073     
 rtp3 =   -3.0909820482975234E-008
 Skthl =    2.1684392401775314     
 thlp3 =    3.3874842468571704     
 wprtp =    9.5158818237846949E-004
 wpthlp =   0.17775656901408005     
 rtpthlp =    2.0977641881265901E-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           6          29        -300           3          17           6          40        2061          46          69        -260          43          57          46          80        2061          46          69        -260          43          57          46          80        2061          46          69        -260          43          57          46          80        2061
 wm =    1.3990036430270858E-002
 rtm =    6.4684105990897214E-003
 thlm =    292.40606606552501     
 wp2 =   0.46049349775869042     
 rtp2 =    4.0352919802364656E-006
 thlp2 =    1.3463356663843720     
 Skw =  -0.85249615840124526     
 wp3 =  -0.26639599049392881     
 Skrt =   -3.8131512876811073     
 rtp3 =   -3.0909820482975234E-008
 Skthl =    2.1684392401775314     
 thlp3 =    3.3874842468571704     
 wprtp =    9.5158818237846949E-004
 wpthlp =   0.17775656901408005     
 rtpthlp =    2.0977641881265901E-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           6          29        -300           3          17           6          46        2067          52          75        -254          49          63          52          92        2067          52          75        -254          49          63          52          92        2067          52          75        -254          49          63          52          92        2067
 wm =    8.3022784799689830E-002
 rtm =    1.7814562682763539E-002
 thlm =    308.53420095539843     
 wp2 =   0.62827218730139489     
 rtp2 =    3.5296840155703257E-005
 thlp2 =    1.4542574526886949     
 Skw =    1.2916081858141348     
 wp3 =   0.64320963393208919     
 Skrt =    1.8454840280338882     
 rtp3 =    3.8700273407515525E-007
 Skthl =   -4.4974834986191778     
 thlp3 =   -7.8873574962213322     
 wprtp =    4.3639869880778099E-003
 wpthlp =   0.19295961691502631     
 rtpthlp =    6.4480897334607849E-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           6          29        -300           3          17           6          59        2080          65          88        -241          62          76          65         118        2080          65          88        -241          62          76          65         118        2080          65          88        -241          62          76          65         118        2080
 wm =   -4.3658041121098501E-003
 rtm =    2.7402320889264435E-003
 thlm =    309.05334368786970     
 wp2 =   0.81148637345525854     
 rtp2 =    6.6377798907301190E-007
 thlp2 =   0.74251886545838031     
 Skw =    3.8822421080326652     
 wp3 =    2.8379481926209320     
 Skrt =   -7.9079904517566524E-002
 rtp3 =   -4.2766174814438903E-011
 Skthl =  -0.49074891970094336     
 thlp3 =  -0.31399345160404668     
 wprtp =    7.0606949324429342E-004
 wpthlp =   0.42344666963401245     
 rtpthlp =    6.3184081878572699E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 Success!
 
 
 Performing PDF parameter values unit test
 =========================================
 
 Performing PDF parameter unit tests for the new hybrid PDF
 
 PDF parameter set 1:
 wm =   -5.0000000000000001E-003
 rtm =    4.2868400000000001E-003
 thlm =    309.13099999999997     
 wp2 =    3.3741899999999998E-002
 rtp2 =    1.1786499999999999E-006
 thlp2 =    5.9920000000000001E-002
 Skw =    8.6882480819056749     
 wp3 =    5.3850099999999998E-002
 Skrt =    1.4880899999999999     
 rtp3 =    1.9041733779221134E-009
 Skthl =   -1.4580900000000001     
 thlp3 =   -2.1386614365822010E-002
 wprtp =    2.3725400000000002E-005
 wpthlp =   -5.5512699999999996E-003
 rtpthlp =   -2.3917778634313013E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 2:
 wm =    1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.10000000000000001     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =   -4.5000000000000000     
 wp3 =  -0.14230249470757708     
 Skrt =    1.0000000000000000     
 rtp3 =    1.1180339887498949E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =   -5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 3:
 wm =   -1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.50000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =    4.5000000000000000     
 wp3 =    1.5909902576697321     
 Skrt =    2.5000000000000000     
 rtp3 =    2.7950849718747374E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 4:
 wm =   0.10000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    5.0000000000000000     
 rtp2 =    6.9999999999999999E-006
 thlp2 =    1.0000000000000000     
 Skw =    2.0000000000000000     
 wp3 =    22.360679774997898     
 Skrt =    1.5000000000000000     
 rtp3 =    2.7780388766178202E-008
 Skthl =   -1.5000000000000000     
 thlp3 =   -1.5000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -1.0000000000000000E-002
 rtpthlp =   -2.3811761799581318E-003
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 5:
 wm =   -1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.29999999999999999     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.20000000000000001     
 Skw =   -2.0000000000000000     
 wp3 =  -0.32863353450309962     
 Skrt =  -0.25000000000000000     
 rtp3 =   -7.0710678118654745E-010
 Skthl =   -3.5000000000000000     
 thlp3 =  -0.31304951684997062     
 wprtp =    5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =    5.6920997883030825E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 6:
 wm =   -1.0000000000000000E-002
 rtm =    5.0000000000000001E-003
 thlm =    305.00000000000000     
 wp2 =    1.0000000000000000     
 rtp2 =    4.9999999999999998E-007
 thlp2 =   0.10000000000000001     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -2.0124611797498108E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 7:
 wm =   0.20000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.25000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.14999999999999999     
 Skw =   0.75000000000000000     
 wp3 =    9.3750000000000000E-002
 Skrt =    3.0000000000000000     
 rtp3 =    3.3541019662496846E-008
 Skthl =  -0.50000000000000000     
 thlp3 =   -2.9047375096555625E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -7.7942286340599491E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 8:
 wm =   -2.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    0.0000000000000000     
 thlp2 =    0.0000000000000000     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    0.0000000000000000     
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 9:
 wm =    1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.69999999999999996     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =   -1.0000000000000000     
 rtp3 =   -2.8284271247461898E-009
 Skthl =   -2.0000000000000000     
 thlp3 =   -1.1713240371477056     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    1.0648943609579309E-003
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 10 (randomly generated):
 Random seed values =         2021           6          29        -300           3          17           6          71        2092          77         100        -229          74          88          77         142        2092          77         100        -229          74          88          77         142        2092          77         100        -229          74          88          77         142        2092
 wm =  -0.33251587004716299     
 rtm =    4.5415620563433968E-003
 thlm =    290.54619859880802     
 wp2 =   0.63392930285821392     
 rtp2 =    2.4310363864626030E-006
 thlp2 =    1.4308702413293375     
 Skw =    1.7281710667349970     
 wp3 =   0.87226404392722656     
 Skrt =    3.4804074572572210     
 rtp3 =    1.3192201594651528E-008
 Skthl =   0.88481384594188917     
 thlp3 =    1.5144407511370710     
 wprtp =    8.3938292300530235E-004
 wpthlp =  -0.58266026037597329     
 rtpthlp =   -1.6785657785717834E-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.454761250379875     
 Vertical integral of <rt> forcing =   -3.9696697370157397E-005
 Flux of <rt> at the surface =    5.0000000000000001E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    1.3037531124138191E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    793918.42642888101     
 Vertical integral of <thl> after advance_xm_wpxp =    793908.03522151220     
 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> =    8.7298224205056840E-014
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            2
 
 Vertical integral of <rt> before advance_xm_wpxp =    28.639554483708842     
 Vertical integral of <rt> after advance_xm_wpxp =    28.690899728050656     
 Vertical integral of <rt> forcing =   -3.6936782447310184E-005
 Flux of <rt> at the surface =    2.0808759691998734E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    3.5731237846975406E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    793385.15959492838     
 Vertical integral of <thl> after advance_xm_wpxp =    793369.92936369812     
 Vertical integral of <thl> forcing =   -5.4929189373047839E-002
 Flux of <thl> at the surface =    4.1617519383997468E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    4.3610254296666540E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            3
 
 Vertical integral of <rt> before advance_xm_wpxp =    37.847868347453371     
 Vertical integral of <rt> after advance_xm_wpxp =    37.971603659147632     
 Vertical integral of <rt> forcing =   -5.7743666783424203E-005
 Flux of <rt> at the surface =    4.7019470576422197E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    7.3698642674702164E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    797701.72306735942     
 Vertical integral of <thl> after advance_xm_wpxp =    797684.35047198320     
 Vertical integral of <thl> forcing =   -6.7312545369553958E-002
 Flux of <thl> at the surface =    9.4038941152844394E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    2.0729251648532454E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            4
 
 Vertical integral of <rt> before advance_xm_wpxp =    10.523778178339809     
 Vertical integral of <rt> after advance_xm_wpxp =    10.605435206353198     
 Vertical integral of <rt> forcing =   -1.5571694241973466E-005
 Flux of <rt> at the surface =    2.8776178761992044E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    1.5731095896406866E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    787551.98586361436     
 Vertical integral of <thl> after advance_xm_wpxp =    787549.97323695174     
 Vertical integral of <thl> forcing =   -1.2463991293956325E-002
 Flux of <thl> at the surface =    5.7552357523984092E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -5.0058221456872332E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            5
 
 Vertical integral of <rt> before advance_xm_wpxp =    15.969537090899822     
 Vertical integral of <rt> after advance_xm_wpxp =    15.961001955828573     
 Vertical integral of <rt> forcing =   -2.9028212874470324E-005
 Flux of <rt> at the surface =    5.7776263697539855E-007
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -3.7608262858090935E-019
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    789839.77897984849     
 Vertical integral of <thl> after advance_xm_wpxp =    789829.30923091248     
 Vertical integral of <thl> forcing =   -3.4910718372894951E-002
 Flux of <thl> at the surface =    1.1555252739507971E-005
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    1.1726036808212825E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Success!
 
 Running tuner tests
 -------------------
 goldstein_price pass
 rastrigin pass
 himmelblau pass
 eggholder pass
 schaffer pass
All tests have succeeded.