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           7          19        -300          11          33          54         407        2428         414         426         107         418         440         461         814        2428         414         426         107         418         440         461         814        2428         414         426         107         418         440         461         814        2428
 hmm =    3.0198256220120073E-005
 mixt_frac =   0.21963623215635672     
 precip_frac =   0.30230163104030966     
 precip_frac_1 =    3.1255735971661085E-002
 precip_frac_2 =   0.37858848800040001     
 hmp2_ip_on_hmm2_ip =   0.64493547592890299     
 hmm_ip =    9.9894453484087754E-005
 hmp2 =    4.0502406810875713E-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           7          19        -300          11          33          54         410        2431         417         429         110         421         443         464         820        2431         417         429         110         421         443         464         820        2431         417         429         110         421         443         464         820        2431
 wm =    6.3770566889241320E-002
 rtm =    2.0805193422237273E-002
 thlm =    301.46993598528809     
 wp2 =    3.3147272113153270E-002
 rtp2 =    8.7582310679591959E-005
 thlp2 =   0.98182714052237863     
 Skw =   -4.3657219450111109     
 wp3 =   -2.6346793169802870E-002
 Skrt =    2.6504715741288369     
 rtp3 =    2.1724397933273062E-006
 Skthl =   -3.4611396177733735     
 thlp3 =   -3.3672213632963306     
 wprtp =   -1.3638704656715507E-003
 wpthlp =    7.7072466998266082E-003
 rtpthlp =   -8.3458048515681199E-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           7          19        -300          11          33          54         410        2431         417         429         110         421         443         464         820        2431         417         429         110         421         443         464         820        2431         417         429         110         421         443         464         820        2431
 wm =    6.3770566889241320E-002
 rtm =    2.0805193422237273E-002
 thlm =    301.46993598528809     
 wp2 =    3.3147272113153270E-002
 rtp2 =    8.7582310679591959E-005
 thlp2 =   0.98182714052237863     
 Skw =   -4.3657219450111109     
 wp3 =   -2.6346793169802870E-002
 Skrt =    2.6504715741288369     
 rtp3 =    2.1724397933273062E-006
 Skthl =   -3.4611396177733735     
 thlp3 =   -3.3672213632963306     
 wprtp =   -1.3638704656715507E-003
 wpthlp =    7.7072466998266082E-003
 rtpthlp =   -8.3458048515681199E-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           7          19        -300          11          33          54         415        2436         422         434         115         426         448         469         830        2436         422         434         115         426         448         469         830        2436         422         434         115         426         448         469         830        2436
 wm =   -7.0846647868444523E-002
 rtm =    2.4184147068069921E-003
 thlm =    292.15659840494726     
 wp2 =   0.21602086563519751     
 rtp2 =    8.1007611312408889E-007
 thlp2 =   0.23020871563038170     
 Skw =    3.9258661132256543     
 wp3 =   0.39416588868673796     
 Skrt =    2.6470525445694921     
 rtp3 =    1.9299733032229847E-009
 Skthl =  -0.21542932047605035     
 thlp3 =   -2.3795095620669553E-002
 wprtp =   -3.2504371214282611E-004
 wpthlp =    7.1598500349429869E-002
 rtpthlp =   -3.8865682943667687E-004
 
 Success!
 
 
 Performing PDF parameter values unit test
 =========================================
 
 Performing PDF parameter unit tests for the new PDF
 
 PDF parameter set 1:
 wm =   -5.0000000000000001E-003
 rtm =    4.2868400000000001E-003
 thlm =    309.13099999999997     
 wp2 =    3.3741899999999998E-002
 rtp2 =    1.1786499999999999E-006
 thlp2 =    5.9920000000000001E-002
 Skw =    8.6882480819056749     
 wp3 =    5.3850099999999998E-002
 Skrt =    1.4880899999999999     
 rtp3 =    1.9041733779221134E-009
 Skthl =   -1.4580900000000001     
 thlp3 =   -2.1386614365822010E-002
 wprtp =    2.3725400000000002E-005
 wpthlp =   -5.5512699999999996E-003
 rtpthlp =   -2.3917778634313013E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 2:
 wm =    1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.10000000000000001     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =   -4.5000000000000000     
 wp3 =  -0.14230249470757708     
 Skrt =    1.0000000000000000     
 rtp3 =    1.1180339887498949E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =   -5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 3:
 wm =   -1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.50000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =    4.5000000000000000     
 wp3 =    1.5909902576697321     
 Skrt =    2.5000000000000000     
 rtp3 =    2.7950849718747374E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 4:
 wm =   0.10000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    5.0000000000000000     
 rtp2 =    6.9999999999999999E-006
 thlp2 =    1.0000000000000000     
 Skw =    2.0000000000000000     
 wp3 =    22.360679774997898     
 Skrt =    1.5000000000000000     
 rtp3 =    2.7780388766178202E-008
 Skthl =   -1.5000000000000000     
 thlp3 =   -1.5000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -1.0000000000000000E-002
 rtpthlp =   -2.3811761799581318E-003
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 5:
 wm =   -1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.29999999999999999     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.20000000000000001     
 Skw =   -2.0000000000000000     
 wp3 =  -0.32863353450309962     
 Skrt =  -0.25000000000000000     
 rtp3 =   -7.0710678118654745E-010
 Skthl =   -3.5000000000000000     
 thlp3 =  -0.31304951684997062     
 wprtp =    5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =    5.6920997883030825E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 6:
 wm =   -1.0000000000000000E-002
 rtm =    5.0000000000000001E-003
 thlm =    305.00000000000000     
 wp2 =    1.0000000000000000     
 rtp2 =    4.9999999999999998E-007
 thlp2 =   0.10000000000000001     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -2.0124611797498108E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 7:
 wm =   0.20000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.25000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.14999999999999999     
 Skw =   0.75000000000000000     
 wp3 =    9.3750000000000000E-002
 Skrt =    3.0000000000000000     
 rtp3 =    3.3541019662496846E-008
 Skthl =  -0.50000000000000000     
 thlp3 =   -2.9047375096555625E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -7.7942286340599491E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 8:
 wm =   -2.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    0.0000000000000000     
 thlp2 =    0.0000000000000000     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    0.0000000000000000     
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 9:
 wm =    1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.69999999999999996     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =   -1.0000000000000000     
 rtp3 =   -2.8284271247461898E-009
 Skthl =   -2.0000000000000000     
 thlp3 =   -1.1713240371477056     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    1.0648943609579309E-003
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 10 (randomly generated):
 Random seed values =         2021           7          19        -300          11          33          54         427        2448         434         446         127         438         460         481         854        2448         434         446         127         438         460         481         854        2448         434         446         127         438         460         481         854        2448
 wm =  -0.14041552108746513     
 rtm =    4.0739459244499678E-003
 thlm =    306.43894547078236     
 wp2 =   0.74889031642100401     
 rtp2 =    1.1147846775937789E-006
 thlp2 =    9.9317975109613488E-002
 Skw =    3.3685160976301822     
 wp3 =    2.1830613942254353     
 Skrt =  -0.84052417687196013     
 rtp3 =   -9.8931978913659027E-010
 Skthl =  -0.33890166396974097     
 thlp3 =   -1.0607559734040558E-002
 wprtp =    8.4033251715212045E-004
 wpthlp =    5.4892621687830601E-002
 rtpthlp =    2.9946904190264532E-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           7          19        -300          11          33          54         438        2459         445         457         138         449         471         492         876        2459         445         457         138         449         471         492         876        2459         445         457         138         449         471         492         876        2459
 wm =   -9.0576422873809848E-002
 rtm =    1.1822386615267260E-002
 thlm =    307.42627587260847     
 wp2 =   0.23304715795898412     
 rtp2 =    8.5183100121334167E-006
 thlp2 =    2.5568856687646599     
 Skw =  -0.12886475551169130     
 wp3 =   -1.4497725436596047E-002
 Skrt =  -0.27384301528798471     
 rtp3 =   -6.8081925582399744E-009
 Skthl =   -1.9745566175715257     
 thlp3 =   -8.0730297790389507     
 wprtp =   -1.8689787067390222E-004
 wpthlp =   0.31947001955674870     
 rtpthlp =   -4.2002475262317516E-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 =    30.753373618421570     
 Vertical integral of <rt> after advance_xm_wpxp =    30.824885422704192     
 Vertical integral of <rt> forcing =   -5.1188909884094147E-005
 Flux of <rt> at the surface =    2.8956159082622315E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -5.5890621991627754E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    794112.28116669890     
 Vertical integral of <thl> after advance_xm_wpxp =    794096.42993908061     
 Vertical integral of <thl> forcing =   -5.8628657211119678E-002
 Flux of <thl> at the surface =    5.7912318165244626E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    3.0105778980882292E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            3
 
 Vertical integral of <rt> before advance_xm_wpxp =    23.069688288914165     
 Vertical integral of <rt> after advance_xm_wpxp =    23.196551535169142     
 Vertical integral of <rt> forcing =   -3.3883219322171509E-005
 Flux of <rt> at the surface =    4.5676070683873225E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    3.0804894225744395E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    793972.08784434409     
 Vertical integral of <thl> after advance_xm_wpxp =    793964.66771038249     
 Vertical integral of <thl> forcing =   -3.3868994008441046E-002
 Flux of <thl> at the surface =    9.1352141367746463E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -3.2512187386757319E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            4
 
 Vertical integral of <rt> before advance_xm_wpxp =    44.568486608245131     
 Vertical integral of <rt> after advance_xm_wpxp =    44.792819830486849     
 Vertical integral of <rt> forcing =   -6.6913364971723607E-005
 Flux of <rt> at the surface =    8.1469077244421544E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -9.8404899680215596E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    799345.63680456916     
 Vertical integral of <thl> after advance_xm_wpxp =    799329.24436809309     
 Vertical integral of <thl> forcing =   -7.0935270370946130E-002
 Flux of <thl> at the surface =    1.6293815448884308E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    1.8115231537052523E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            5
 
 Vertical integral of <rt> before advance_xm_wpxp =    40.907129670589995     
 Vertical integral of <rt> after advance_xm_wpxp =    41.126927664848708     
 Vertical integral of <rt> forcing =   -4.6734766212221454E-005
 Flux of <rt> at the surface =    7.7939474707456691E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    2.9571614254542133E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    797779.95700221346     
 Vertical integral of <thl> after advance_xm_wpxp =    797766.09240741609     
 Vertical integral of <thl> forcing =   -6.1803210932470765E-002
 Flux of <thl> at the surface =    1.5587894941491338E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -2.5708601913976281E-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.