Skipping 358 KB..
Full LogTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTest
2641 | subroutine set_case_initial_conditions(err_info)
| 1
Warning: Unused dummy argument 'err_info' at (1) [-Wunused-dummy-argument]
[901/915] Linking Fortran static library src/libclubb_driver_lib.a
[902/915] Generating Fortran dyndep file src/CMakeFiles/clubb_tuner.dir/Fortran.dd
[903/915] Generating Fortran dyndep file src/CMakeFiles/clubb_standalone.dir/Fortran.dd
[904/915] Generating Fortran dyndep file src/CMakeFiles/clubb_driver_test.dir/Fortran.dd
[905/915] Generating Fortran dyndep file src/CMakeFiles/clubb_thread_test.dir/Fortran.dd
[906/915] Building Fortran object src/CMakeFiles/clubb_thread_test.dir/clubb_thread_test.F90.o
[907/915] Building Fortran object src/CMakeFiles/clubb_standalone.dir/clubb_standalone.F90.o
[908/915] Building Fortran object src/CMakeFiles/clubb_driver_test.dir/clubb_driver_test.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90:17:6:
17 | use error_code, only: &
| 1
Warning: Unused parameter 'clubb_generalized_grd_test_err' which has been explicitly imported at (1) [-Wunused-parameter]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90:21:6:
21 | use model_flags, only: &
| 1
Warning: Unused parameter 'l_test_grid_generalization' which has been explicitly imported at (1) [-Wunused-parameter]
[909/915] Linking Fortran executable src/clubb_thread_test
[910/915] Linking Fortran executable src/clubb_driver_test
[911/915] Building Fortran object src/CMakeFiles/clubb_tuner.dir/error.F90.o
[912/915] Linking Fortran executable src/clubb_standalone
[913/915] Building Fortran object src/CMakeFiles/clubb_tuner.dir/clubb_tuner.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_tuner.F90:294:6:
294 | use nrtype, only: &
| 1
Warning: Unused parameter 'sp' which has been explicitly imported at (1) [-Wunused-parameter]
[914/915] Linking Fortran executable src/clubb_tuner
[914/915] Install the project...
-- Install configuration: "Debug"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/libclubb_driver_lib.a
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_tuner
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_tuner" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/G_unit_tests
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/G_unit_tests" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_thread_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_thread_test" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_standalone
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_standalone" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_driver_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_driver_test" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
check_for_missing_threadprivate.py has begun.
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/corr_varnce_module.F90 failed:
Missing threadprivate for:
corr_array_n_below_def
corr_array_n_cloud_def
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/error_code.F90 failed:
Missing threadprivate for:
clubb_debug_level
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sponge_layer_damping.F90 failed:
Missing threadprivate for:
wp2_sponge_damp_settings
rtm_sponge_damp_settings
up2_vp2_sponge_damp_settings
thlm_sponge_damp_settings
wp3_sponge_damp_settings
thlm_sponge_damp_profile
wp2_sponge_damp_profile
wp3_sponge_damp_profile
up2_vp2_sponge_damp_profile
rtm_sponge_damp_profile
uv_sponge_damp_settings
uv_sponge_damp_profile
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90 failed:
Missing threadprivate for:
ret_code
[91mcheck_for_missing_threadprivate.py did not pass.[0m
[91mfailedFiles = ['/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/corr_varnce_module.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/error_code.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sponge_layer_damping.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90'][0m
[91m===============================================================
[91mTHREADPRIVATE CHECK FAILED
[91m THIS IS PRINTED IN ALL RED, CAPITAL LETTERS, AND USES
[91m AN EXCLAMATION MARK TO ENSURE THE DEVELOPERS FEEL SHAME!
[91m IF YOU ARE ONE OF THESE "DEVELOPERS" CHECK THE
[91m LOG FILE FOR DETAILS: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gfortran_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
CLUBBStandardsCheck.py has begun.
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/generalized_grid_test.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 3: module generalized_grid_test
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/generalized_grid_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 903 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1071 : 101 chars : rtm_ref_flip, thlm_ref_flip, um_ref_flip, vm_ref_flip, ug_flip, vg_flip, & ! Intent(in)
line 1084 : 111 chars : um_flip, vm_flip, upwp_flip, vpwp_flip, up2_flip, vp2_flip, up3_flip, vp3_flip, & ! Intent(inout)
line 1086 : 111 chars : wp2_flip, wp3_flip, rtp2_flip, rtp3_flip, thlp2_flip, thlp3_flip, rtpthlp_flip, & ! Intent(inout)
line 1092 : 109 chars : wp2rtp_flip, wp2thlp_flip, uprcp_flip, vprcp_flip, rc_coef_zm_flip, wp4_flip, & ! intent(inout)
line 1093 : 105 chars : wpup2_flip, wpvp2_flip, wp2up2_flip, wp2vp2_flip, ice_supersat_frac_flip, & ! intent(inout)
line 1582 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1609 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2up_implicit, nzt, ngrdcol, &
line 1613 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2up_explicit, nzt, ngrdcol, &
line 1617 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2vp_implicit, nzt, ngrdcol, &
line 1621 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2vp_explicit, nzt, ngrdcol, &
line 1655 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpup2_implicit, nzt, ngrdcol, &
line 1659 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpup2_explicit, nzt, ngrdcol, &
line 1663 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpvp2_implicit, nzt, ngrdcol, &
line 1667 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpvp2_explicit, nzt, ngrdcol, &
line 2471 : 105 chars : tolerance = max( 1.e-8_core_rknd, epsilon(tolerance) ) ! max statement for single precision runs
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90 ---------------------
WARNING: 'use' statements without 'only' found in the following lines:
WARNING: Check that comma is on same line as 'use', as CLUBB requires.
3107 : use gptl
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 200 : 110 chars : l_modify_ic_with_cubic_int, & ! Flag for interpolating the sounding profile with Steffen's monotone cubic
line 201 : 110 chars : ! method to obtain smoother initial condition profile, which is found to be
line 202 : 111 chars : ! beneficial to achive a better numerical solution convergence. If this flag
line 203 : 118 chars : ! is turned off, the initial conditions will be generated with linear interpolation.
line 204 : 111 chars : ! This is done on a case-by-case basis, since using the monotone cubic method
line 205 : 109 chars : ! requires a special sounding.in file with many additional sounding levels.
line 206 : 109 chars : l_modify_bc_for_cnvg_test, & ! Flag to activate modifications on boundary condition for convergence test
line 738 : 101 chars : !$omp sens_ht, latent_ht, lmin, vert_decorr_coef, deep_soil_T_in_K, sfc_soil_T_in_K, veg_T_in_K, &
line 745 : 101 chars : !$omp exner, rtm, thlm, rcm, wp3, wp3_zm, delta_zm, p_in_Pa, cloud_frac, wpthvp, wp2thvp, wp2up, &
line 754 : 102 chars : !$omp rtp2_mc, thlp2_mc, rtpthlp_mc, Ncm_mc, Ncm, Nccnm, em, tau_zm, tau_zt, Kh_zt, Kh_zm, rfrzm, &
line 755 : 105 chars : !$omp rrm, Nrm, Nc_in_cloud, wpNcp, radht, Frad, Frad_SW_up, Frad_LW_up, Frad_SW_down, Frad_LW_down, &
line 756 : 102 chars : !$omp wpsclrp_sfc, wpedsclrp_sfc, sclrm, wpsclrp, sclrp2, sclrp3, sclrprtp, sclrpthlp, sclrpthvp, &
line 1016 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 1746 : 102 chars : open(unit=iunit_grid_adaptation, file=fname_grid_adaptation, status='replace', action='write')
line 2062 : 105 chars : allocate( um_ref(ngrdcol, gr%nzt) ) ! Reference u wind for nudging; Michael Falk, 17 Oct 2007
line 2063 : 105 chars : allocate( vm_ref(ngrdcol, gr%nzt) ) ! Reference v wind for nudging; Michael Falk, 17 Oct 2007
line 2502 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 2533 : 106 chars : ! Similar call to above, checking parameter values, but because we've placed this in "init_clubb_case"
line 2541 : 116 chars : write(fstderr, *) "Fatal error calling check_clubb_settings_api and/or check_parameters_api in clubb_driver"
line 2584 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 2596 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 2597 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 2600 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 2602 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 2608 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 2609 : 106 chars : !$acc pdf_params_zm%corr_w_rt_1, pdf_params_zm%corr_w_rt_2, pdf_params_zm%corr_w_thl_1, &
line 2610 : 110 chars : !$acc pdf_params_zm%corr_w_thl_2, pdf_params_zm%corr_rt_thl_1, pdf_params_zm%corr_rt_thl_2, &
line 2611 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 2612 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 2613 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 2614 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 2615 : 108 chars : !$acc pdf_params_zm%corr_w_chi_1, pdf_params_zm%corr_w_chi_2, pdf_params_zm%corr_w_eta_1, &
line 2616 : 112 chars : !$acc pdf_params_zm%corr_w_eta_2, pdf_params_zm%corr_chi_eta_1, pdf_params_zm%corr_chi_eta_2, &
line 2617 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 2726 : 103 chars : invrs_rho_ds_zt(i,k) = invrs_rho_ds_zt_init(k) ! inv. dry, static density: t-levs
line 2732 : 103 chars : um_pert(i,k) = zero ! Variables used to track perturbed version of winds.
line 2734 : 105 chars : wm_zt(i,k) = wm_zt_init(k) ! Imposed large scale w - Thermodynamic levels
line 2928 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 3119 : 101 chars : !----------------------------------- Optional Input Variables -----------------------------------
line 3276 : 103 chars : um(i,:), vm(i,:), rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), &
line 3609 : 111 chars : !$acc X_nl_all_levs, X_mixt_comp_all_levs, lh_sample_point_weights, mu_x_1_n, mu_x_2_n, &
line 3611 : 104 chars : !$acc lh_rt_clipped, lh_thl_clipped, lh_rc_clipped, lh_rv_clipped, lh_Nc_clipped, &
line 3615 : 116 chars : !$acc pdf_params%cloud_frac_1, pdf_params%cloud_frac_2, pdf_params%cthl_1, pdf_params%cthl_2, &
line 3616 : 112 chars : !$acc pdf_params%crt_1, pdf_params%crt_2, pdf_params%stdev_chi_1, pdf_params%stdev_chi_2, &
line 3632 : 102 chars : wp2_zt = zm2zt_api( gr%nzm, gr%nzt, ngrdcol, gr, wp2, w_tol_sqd ) ! Positive definite quantity
line 3638 : 109 chars : thlm(i,:), p_in_Pa(i,:), exner(i,:), rho(i,:), rho_zm(i,:), rtm(i,:), &! In
line 3639 : 109 chars : rcm(i,:), cloud_frac(i,:), wm_zt(i,:), wm_zm(i,:), wp2_zt(i,:), & ! In
line 3640 : 106 chars : hydromet(i,:,:), Nc_in_cloud(i,:), & ! In
line 3657 : 102 chars : Nccnm(i,:), & ! Inout
line 3658 : 117 chars : hydromet_mc(i,:,:), Ncm_mc(i,:), rcm_mc(i,:), rvm_mc(i,:), & ! Out
line 3659 : 107 chars : thlm_mc(i,:), hydromet_vel_zt(i,:,:), & ! Out
line 3660 : 102 chars : hydromet_vel_covar_zt_impc(i,:,:), & ! Out
line 3661 : 102 chars : hydromet_vel_covar_zt_expc(i,:,:), & ! Out
line 3662 : 110 chars : wprtp_mc(i,:), wpthlp_mc(i,:), rtp2_mc(i,:), & ! Out
line 3663 : 105 chars : thlp2_mc(i,:), rtpthlp_mc(i,:) ) ! Out
line 3687 : 103 chars : stats, i, & ! Inout
line 3688 : 103 chars : hydromet(i,:,:), hydromet_vel_zt(i,:,:), hydrometp2(i,:,:), & ! Inout
line 3689 : 103 chars : K_hm(i,:,:), Ncm(i,:), Nc_in_cloud(i,:), rvm_mc(i,:), & ! Inout
line 3690 : 103 chars : thlm_mc(i,:), err_info, & ! Inout
line 3691 : 101 chars : wphydrometp(i,:,:), wpNcp(i,:) ) ! Out
line 3700 : 106 chars : !$acc update device( rcm_mc, rvm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3753 : 104 chars : !$acc update host( rvm_mc, rcm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3780 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, X_nl_all_levs, &
line 3806 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, thlm, rtm, rcm )
line 3907 : 112 chars : call stats_update( "brunt_term_time_avg", richardson_num_term_weight*richardson_num_term_time_avg, &
line 3933 : 101 chars : .and. stats%l_last_sample ) then ! only adapt grid when the average of the last
line 4062 : 103 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_scheme = ', time_microphys_scheme
line 4063 : 104 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_advance = ', time_microphys_advance
line 4182 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 4194 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 4195 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 4198 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 4200 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 4206 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 4207 : 106 chars : !$acc pdf_params_zm%corr_w_rt_1, pdf_params_zm%corr_w_rt_2, pdf_params_zm%corr_w_thl_1, &
line 4208 : 110 chars : !$acc pdf_params_zm%corr_w_thl_2, pdf_params_zm%corr_rt_thl_1, pdf_params_zm%corr_rt_thl_2, &
line 4209 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 4210 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 4211 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 4212 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 4213 : 108 chars : !$acc pdf_params_zm%corr_w_chi_1, pdf_params_zm%corr_w_chi_2, pdf_params_zm%corr_w_eta_1, &
line 4214 : 112 chars : !$acc pdf_params_zm%corr_w_eta_2, pdf_params_zm%corr_chi_eta_1, pdf_params_zm%corr_chi_eta_2, &
line 4215 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 4237 : 119 chars : if ( thlm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( thlm_sponge_damp_profile )
line 4238 : 118 chars : if ( rtm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( rtm_sponge_damp_profile )
line 4239 : 117 chars : if ( uv_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( uv_sponge_damp_profile )
line 4240 : 118 chars : if ( wp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp2_sponge_damp_profile )
line 4241 : 118 chars : if ( wp3_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp3_sponge_damp_profile )
line 4242 : 122 chars : if ( up2_vp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( up2_vp2_sponge_damp_profile )
line 4929 : 104 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 4937 : 106 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 4951 : 101 chars : 'be used for cases with forcings from an input file and for the atex ', &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_tuner.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 654 : 102 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 71 : 107 chars : write(fstderr, *) "This section ensures that everything allocated in init_clubb_case will be deallocated"
line 72 : 102 chars : write(fstderr, *) "in clean_up_clubb. This may cause a runtime error if there is a mismatch between"
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 6 check(s) failed.
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_clubb_core_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1876 : 101 chars : stats, & ! intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_xm_wpxp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 779 : 101 chars : call wpxp_term_pr1_lhs( nzm, ngrdcol, gr, C6rt_Skw_fnc, C6thl_Skw_fnc, C7_Skw_fnc, & ! Intent(in)
line 780 : 101 chars : invrs_tau_C6_zm, l_scalar_calc, & ! Intent(in)
line 781 : 102 chars : lhs_pr1_wprtp, lhs_pr1_wpthlp, lhs_pr1_wpsclrp ) ! Intent(out)
line 854 : 102 chars : stats, & ! InOut
line 855 : 101 chars : rtm, wprtp, thlm, wpthlp, sclrm, wpsclrp, & ! InOut
line 856 : 101 chars : err_info ) ! InOut
line 896 : 101 chars : stats, & ! InOut
line 4881 : 108 chars : tmp_zm_stats(i,k) = (-gamma_over_implicit_ts*lhs_ta_wpxp(2+gr%grid_dir_indx,i,k)) * wpxp(i,km1) &
line 4883 : 105 chars : + (-gamma_over_implicit_ts*lhs_ta_wpxp(2-gr%grid_dir_indx,i,k)) * wpxp(i,kp1)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sfc_varnce_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 723 : 105 chars : call sfc_varnce_check( sclr_dim, wp2(i,gr%k_lb_zm), up2(i,gr%k_lb_zm), & ! intent(in)
line 724 : 105 chars : vp2(i,gr%k_lb_zm), thlp2(i,gr%k_lb_zm), rtp2(i,gr%k_lb_zm), & ! intent(in)
line 725 : 105 chars : rtpthlp(i,gr%k_lb_zm), sclrp2(i,gr%k_lb_zm,:), & ! intent(in)
line 726 : 105 chars : sclrprtp(i,gr%k_lb_zm,:), sclrpthlp(i,gr%k_lb_zm,:), & ! intent(in)
line 727 : 108 chars : err_info ) ! intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/grid_adaptation_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 295 : 128 chars : write(fstderr, *) "The fixed minimum grid density momentum heights need to cover at least the domain of the physics grid."
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/pdf_parameter_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 295 : 101 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 298 : 102 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 300 : 101 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/clubb_api_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 2430 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
line 2497 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/stats_clubb_utilities.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 129 : 101 chars : wp2_zt, thlp2_zt, wpthlp_zt, wprtp_zt, rtp2_zt, rtpthlp_zt, up2_zt, vp2_zt, upwp_zt, vpwp_zt, &
line 130 : 101 chars : wpup2, wpvp2, a3_coef_zt, wp3_on_wp2_zt, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_windm_edsclrm_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 493 : 102 chars : stats ) ! intent(inout)
line 501 : 102 chars : stats ) ! intent(inout)
line 953 : 103 chars : call windm_edsclrm_rhs( nzm, nzt, ngrdcol, gr, windm_edsclrm_scalar, dt, & ! intent(in)
line 954 : 103 chars : lhs_diff, edsclrm(:,:,edsclr), & ! intent(in)
line 955 : 103 chars : edsclrm_forcing(:,:,edsclr), & ! intent(in)
line 956 : 103 chars : rho_ds_zm, invrs_rho_ds_zt, & ! intent(in)
line 957 : 103 chars : l_imp_sfc_momentum_flux, wpedsclrp(:,gr%k_lb_zm,edsclr), & ! intent(in)
line 958 : 106 chars : stats, & ! intent(inout)
line 959 : 104 chars : rhs(:,:,edsclr) ) ! intent(out)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/setup_clubb_pdf_params.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 55 : 104 chars : subroutine setup_pdf_parameters_api( gr, nzm, nzt, ngrdcol, pdf_dim, hydromet_dim, dt, & ! Intent(in)
line 3823 : 104 chars : ! Correlation (in-precip) of chi (old s) and the precipitating hydrometeor in PDF component 1.
line 3827 : 104 chars : ! Correlation (in-precip) of chi (old s) and the precipitating hydrometeor in PDF component 2.
line 3841 : 104 chars : ! Correlation (in-precip) of eta (old t) and the precipitating hydrometeor in PDF component 1.
line 3845 : 104 chars : ! Correlation (in-precip) of eta (old t) and the precipitating hydrometeor in PDF component 2.
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/parameters_tunable.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 839 : 106 chars : / sqrt( 4.0_core_rknd * ( 1.0_core_rknd - 0.4_core_rknd )**3 &
line 840 : 118 chars : + clubb_params(ngrdcol,iSkw_max_mag)**2 ) ) ) ! Known magic number
line 902 : 117 chars : mult_factor_zt(i) = 1.0_core_rknd + clubb_params(i,imult_coef) * log( avg_deltaz(i) / grid_spacing_thresh )
line 917 : 104 chars : nu_vert_res_dep%nu10(i) = clubb_params(i,inu10) * mult_factor_zt(i) !We're unsure of the grid
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/mono_flux_limiter.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 700 : 104 chars : !$acc parallel loop gang vector collapse(2) default(present) reduction(.or.:l_any_adjustment_needed)
line 851 : 121 chars : ! MONOFLUX TEST COMMENT DO NOT REMOVE if( any(wpxp_net_adjust /= 0.0) ) write(fstderr,*) "MONOFLUX: wpxp adjusted "
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/model_flags.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 315 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 541 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 845 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/err_info_type_module.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 4: module err_info_type_module
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/err_info_type_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 265 : 102 chars : write(err_info%err_header_global,'(A38, I3, A15, I3, A, A28, A33, F6.2, A1, F6.2, A, A39, A40)') &
line 276 : 105 chars : write(err_info%err_header(i),'(A38, I3, A15, I3, A, A22, I2, A11, F6.2, A13, F6.2, A, A39, A40)') &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/penta_bicgstab_solver.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 419 : 110 chars : ! It's interesting to see how close E_inv is to the lhs diagonal, it's usually (visually) very close, but
line 420 : 104 chars : ! using it instead of E_inv requires so many more iterations (use it would be jacobi preconditioner)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_wp2_wp3_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 928 : 101 chars : stats, & ! intent(inout)
line 1742 : 105 chars : stats_tmp_zt(i,k) = - gamma_over_implicit_ts * wp3_term_ta_lhs_result(5,i,k) * wp3(i,k-1) &
line 1743 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(4,i,k) * wp2(i,k) &
line 1744 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(3,i,k) * wp3(i,k) &
line 1745 : 105 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(2,i,k) * wp2(i,k+1) &
line 1746 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(1,i,k) * wp3(i,k+1)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 234 : 113 chars : write(*,'(A6,I5,1X,I5,A4,E30.20,A3,E30.20)') "field(", i, k, ") = ", field(i,k), " < ", threshold
line 240 : 105 chars : write(fstderr, *) "column", i, " field average = ", sum(field(i,:) * rho_ds(i,:) * dz(i,:)) &
line 584 : 102 chars : ! rho_ds_dz(i,k-num_hf_draw_points:k+num_hf_draw_points))
line 823 : 110 chars : * min( threshold, field(i,k_start:k_end:grid_dir_indx) ) ) * invrs_denom
line 828 : 106 chars : + mass_fraction * ( min( threshold, field(i,k_start:k_end:grid_dir_indx) ) - threshold )
line 833 : 110 chars : * max( threshold, field(i,k_start:k_end:grid_dir_indx) ) ) * invrs_denom
line 840 : 108 chars : + mass_fraction * ( max( threshold, field(i,k_start:k_end:grid_dir_indx) ) - threshold )
----------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/stats_netcdf.F90 ---------------------
WARNING: Missing 'implicit none' after these blocks:
line 281: subroutine stats_init_api( registry_path, output_path, ncol, stats_tsamp, stats_tout, dt_main, &
line 489: subroutine stats_finalize_api( stats, err_info )
line 565: subroutine stats_type_initialize( ncol, defs, ndefs, zt, zm, stats, ierr, rad_zt, rad_zm )
line 716: subroutine stats_open_netcdf( output_path, time_units, &
line 951: subroutine stats_begin_timestep_api( itime, stats )
line 1001: subroutine stats_end_timestep_api( time_value, stats, err_info )
line 1052: subroutine stats_update_scalar( name, values, stats, icol, level )
line 1112: subroutine stats_update_1d( name, values, stats, icol, level )
line 1172: subroutine stats_update_2d( name, values, stats, icol, level )
line 1235: subroutine stats_begin_budget_scalar( name, values, stats, icol )
line 1281: subroutine stats_begin_budget_1d( name, values, stats, icol )
line 1334: subroutine stats_begin_budget_2d( name, values, stats, icol )
line 1379: subroutine stats_update_budget_scalar( name, values, stats, icol, level )
line 1429: subroutine stats_update_budget_1d( name, values, stats, icol, level )
line 1479: subroutine stats_update_budget_2d( name, values, stats, icol, level )
line 1524: subroutine stats_finalize_budget_scalar( name, values, stats, icol, l_count_sample )
line 1590: subroutine stats_finalize_budget_1d( name, values, stats, icol, l_count_sample )
line 1663: subroutine stats_finalize_budget_2d( name, values, stats, icol, l_count_sample )
line 1729: subroutine stats_avg_and_write( time_value, stats, ierr )
line 1855: subroutine stats_grid_init( ncol, zt_src, zm_src, zt_tgt, zm_tgt, &
line 1917: subroutine stats_update_grid( zt_src, zm_src, rho_vals, rho_levels, p_sfc, stats )
line 1995: subroutine stats_lh_samples_init( num_samples, nzt, nl_var_names, u_var_names, zt_vals, &
line 2114: subroutine stats_lh_samples_write_lognormal( samples, stats, err_info )
line 2167: subroutine stats_lh_samples_write_uniform( uniform_vals, mixture_comp, sample_weights, &
line 2383: subroutine stats_lookup_append( stats, id )
line 2425: subroutine stats_lookup_reject_append( stats, name )
line 2539: subroutine stats_read_registry_namelist( path, defs, ndefs, ierr )
line 2614: subroutine stats_expand_registry( sclr_dim, edsclr_dim, defs, ndefs, hydromet_list )
line 2915: subroutine add_expanded_def( base_def, name, out_defs, nout )
line 2965: subroutine parse_registry_line( line, def, l_valid )
line 3039: subroutine split_registry_fields( line, fields, nfields )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/stats_netcdf.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 235 : 103 chars : logical :: l_last_sample = .false. ! True if we want to output at the end of this sample
line 237 : 104 chars : logical :: l_different_output_grid = .false. ! True if we need to remap the variables before output
line 238 : 102 chars : logical :: l_output_rad_files = .false. ! True if we want to output radiation grid variables
line 526 : 101 chars : if ( allocated( stats%grid%rho_lin_spline_vals ) ) deallocate( stats%grid%rho_lin_spline_vals )
line 527 : 105 chars : if ( allocated( stats%grid%rho_lin_spline_levels ) ) deallocate( stats%grid%rho_lin_spline_levels )
line 803 : 104 chars : ret_code = nf90_def_var( stats%ncid, "time", NF90_DOUBLE, (/ stats%time_dimid /), stats%time_varid )
line 808 : 101 chars : ret_code = nf90_def_var( stats%ncid, "col", NF90_DOUBLE, (/ stats%col_dimid /), stats%col_varid )
line 828 : 121 chars : ret_code = nf90_def_var( stats%ncid, "lh_zt", NF90_DOUBLE, (/ stats%lh_2d%lh_zt_dimid /), stats%lh_2d%lh_zt_varid )
line 835 : 112 chars : ret_code = nf90_def_var( stats%ncid, "rad_zt", NF90_DOUBLE, (/ stats%rad_zt_dimid /), stats%rad_zt_varid )
line 839 : 112 chars : ret_code = nf90_def_var( stats%ncid, "rad_zm", NF90_DOUBLE, (/ stats%rad_zm_dimid /), stats%rad_zm_varid )
line 849 : 101 chars : (/ stats%param_strlen_dimid, stats%param_dimid /), stats%param_name_varid )
line 852 : 109 chars : ret_code = nf90_def_var( stats%ncid, "param", NF90_DOUBLE, (/ stats%param_dimid /), stats%param_varid )
line 912 : 110 chars : ret_code = nf90_put_att( stats%ncid, stats%vars(i)%varid, "long_name", trim( stats%vars(i)%long_name ) )
line 913 : 102 chars : ret_code = nf90_put_att( stats%ncid, stats%vars(i)%varid, "units", trim( stats%vars(i)%units ) )
line 920 : 108 chars : ret_code = nf90_put_var( stats%ncid, stats%col_varid, [( real( i, kind=core_rknd ), i=1, stats%ncol )] )
line 940 : 109 chars : ret_code = nf90_put_var( stats%ncid, stats%param_varid, [( real( i, kind=core_rknd ), i=1, nparams )] )
line 1815 : 101 chars : stats%grid%gr_source, stats%grid%gr_output, stats%vars(i)%nz, &
line 1816 : 107 chars : real( stats%vars(i)%buffer, kind=core_rknd ), stats%vars(i)%out_nz, &
line 1817 : 114 chars : size( stats%grid%rho_lin_spline_vals, 2 ), stats%grid%rho_lin_spline_vals, &
line 1973 : 101 chars : if ( allocated( stats%grid%rho_lin_spline_vals ) ) deallocate( stats%grid%rho_lin_spline_vals )
line 1974 : 105 chars : if ( allocated( stats%grid%rho_lin_spline_levels ) ) deallocate( stats%grid%rho_lin_spline_levels )
line 2416 : 116 chars : if ( stats%lookup%cache_len > 0 ) tmp(1:stats%lookup%cache_len) = stats%lookup%cache(1:stats%lookup%cache_len)
line 2461 : 105 chars : tmp(1:stats%lookup%reject_cache_len) = stats%lookup%reject_cache(1:stats%lookup%reject_cache_len)
line 2528 : 101 chars : write(date(32:33),'(i2.2)') nint( ( ( real( mod( nint( st_time ),3600 ),kind=time_precision ) / &
line 2605 : 101 chars : write( fstderr,* ) "Invalid stats registry entry (expected: name | grid | units | long_name)"
line 2706 : 101 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2707 : 101 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2714 : 101 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2715 : 101 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2722 : 101 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2723 : 101 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2730 : 113 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_1", &
line 2734 : 113 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_2", &
line 2740 : 101 chars : call add_expanded_def(base_def,"corr_w_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2741 : 101 chars : call add_expanded_def(base_def,"corr_w_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2748 : 103 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2749 : 103 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2756 : 103 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2757 : 103 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2764 : 103 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2765 : 103 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2772 : 115 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_1_n", &
line 2776 : 115 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_2_n", &
line 2813 : 112 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"p", &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_xp2_xpyp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1320 : 101 chars : sclrprtp(:,:,sclr) ) ! Intent(inout)
line 1353 : 102 chars : sclrpthlp(:,:,sclr) ) ! Intent(inout)
line 3042 : 101 chars : stats_tmp(i,k) = (-gamma_over_implicit_ts * lhs_ta(2+gr%grid_dir_indx,i,k)) * xapxbp(i,km1) &
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 20 check(s) failed.
CLUBBStandardsCheck.py has begun.
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/SILHS/mt95.F90 ---------------------
WARNING: Missing 'implicit none' after these blocks:
line 362: subroutine init_by_type( put, get )
line 382: subroutine init_by_scalar( put )
line 413: subroutine init_by_array( put )
line 461: subroutine next_state( )
line 593: subroutine genrand_int32_0d( y )
line 614: subroutine genrand_int32_1d( y )
line 629: subroutine genrand_int32_2d( y )
line 644: subroutine genrand_int32_3d( y )
line 659: subroutine genrand_int32_4d( y )
line 674: subroutine genrand_int32_5d( y )
line 689: subroutine genrand_int32_6d( y )
line 704: subroutine genrand_int32_7d( y )
line 720: subroutine genrand_int31_0d( y )
line 732: subroutine genrand_int31_1d( y )
line 747: subroutine genrand_int31_2d( y )
line 762: subroutine genrand_int31_3d( y )
line 777: subroutine genrand_int31_4d( y )
line 792: subroutine genrand_int31_5d( y )
line 807: subroutine genrand_int31_6d( y )
line 822: subroutine genrand_int31_7d( y )
line 838: subroutine genrand_real1_0d( r )
line 853: subroutine genrand_real1_1d( r )
line 868: subroutine genrand_real1_2d( r )
line 883: subroutine genrand_real1_3d( r )
line 898: subroutine genrand_real1_4d( r )
line 913: subroutine genrand_real1_5d( r )
line 928: subroutine genrand_real1_6d( r )
line 943: subroutine genrand_real1_7d( r )
line 959: subroutine genrand_real2_0d( r )
line 974: subroutine genrand_real2_1d( r )
line 989: subroutine genrand_real2_2d( r )
line 1004: subroutine genrand_real2_3d( r )
line 1019: subroutine genrand_real2_4d( r )
line 1034: subroutine genrand_real2_5d( r )
line 1049: subroutine genrand_real2_6d( r )
line 1064: subroutine genrand_real2_7d( r )
line 1080: subroutine genrand_real3_0d( r )
line 1095: subroutine genrand_real3_1d( r )
line 1110: subroutine genrand_real3_2d( r )
line 1125: subroutine genrand_real3_3d( r )
line 1140: subroutine genrand_real3_4d( r )
line 1155: subroutine genrand_real3_5d( r )
line 1170: subroutine genrand_real3_6d( r )
line 1185: subroutine genrand_real3_7d( r )
line 1201: subroutine genrand_res53_0d( r )
line 1218: subroutine genrand_res53_1d( r )
line 1233: subroutine genrand_res53_2d( r )
line 1248: subroutine genrand_res53_3d( r )
line 1263: subroutine genrand_res53_4d( r )
line 1278: subroutine genrand_res53_5d( r )
line 1293: subroutine genrand_res53_6d( r )
line 1308: subroutine genrand_res53_7d( r )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/SILHS/latin_hypercube_driver_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 113 : 101 chars : ! The uniform distribution corresponds to all the same variables as X_nl,
line 410 : 105 chars : !$acc update host(X_u_all_levs,l_in_precip,lh_sample_point_weights,X_mixt_comp_all_levs,k_lh_start)
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 2 check(s) failed.
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/wangara.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 104 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 105 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/mpace_a.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 272 : 103 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 273 : 102 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/prescribe_forcings.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 1: module prescribe_forcings_module
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/prescribe_forcings.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 297 : 104 chars : !$acc enter data create( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 471 : 118 chars : !$acc update device( wm_zt, wm_zm, thlm_forcing, rtm_forcing, um_ref, vm_ref, sclrm_forcing, edsclrm_forcing )
line 911 : 103 chars : !$acc exit data delete( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 1114 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1118 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1122 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1123 : 112 chars : exner_zm(i,km1), exner_zm(i,k00), exner_zm(i,kp1), exner_zm(i,kp2) )
line 1126 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1127 : 108 chars : thlm_zm(i,km1), thlm_zm(i,k00), thlm_zm(i,kp1), thlm_zm(i,kp2) )
line 1130 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1131 : 104 chars : rtm_zm(i,km1), rtm_zm(i,k00), rtm_zm(i,kp1), rtm_zm(i,kp2) )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/rico.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 159 : 103 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 160 : 102 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/dycoms2_rf02.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 121 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 122 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/bomex.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 127 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 128 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/dycoms2_rf01.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 91 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 92 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/gabls2.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 152 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 153 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/mpace_b.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 157 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 158 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/lba.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 96 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 97 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
line 176 : 101 chars : wpthlp_sfc(i) = convert_sens_ht_to_km_s( ( 270._core_rknd * ft**1.5_core_rknd ), rho_sfc(i) )
line 177 : 101 chars : wprtp_sfc(i) = convert_latent_ht_to_m_s( ( 554._core_rknd * ft**1.3_core_rknd ), rho_sfc(i) )
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 11 check(s) failed.
No matches for /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/KK_microphys/
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/pdf_parameter_tests.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 542 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/rev_direction_grid_test.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 1: module rev_direction_grid_test
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/rev_direction_grid_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 233 : 113 chars : !seed used = -1102123834 1480475276 -248259584 525496705 1614213535 -1267052688 6290033 -78466652
line 234 : 127 chars : !seed_output = (/ -1102123834, 1480475276, -248259584, 525496705, 1614213535, -1267052688, 6290033, -78466652 /)
line 750 : 101 chars : call xpyp_term_ta_pdf_lhs_godunov( gr_ascending%nzm, gr_ascending%nzt, 1, gr_ascending, & ! In
line 754 : 104 chars : call xpyp_term_ta_pdf_lhs_godunov( gr_descending%nzm, gr_descending%nzt, 1, gr_descending, & ! In
line 795 : 101 chars : call xpyp_term_ta_pdf_rhs_godunov( gr_ascending%nzm, gr_ascending%nzt, 1, gr_ascending, & ! In
line 801 : 104 chars : call xpyp_term_ta_pdf_rhs_godunov( gr_descending%nzm, gr_descending%nzt, 1, gr_descending, & ! In
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 19: module fill_holes_tests
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 162 : 102 chars : field_initial(i,k) = 10_core_rknd * ( sin( 8.0 * pi * ( k + i ) / nz ) + one ) &
line 180 : 110 chars : if ( l_print_name ) print *, "filling holes of: easy_fill_test_"//trim(fill_type_names(fill_type))
line 195 : 121 chars : write(*,'(A6,I5,1X,I5,A4,E30.20,A3,E30.20)') "field(", i, k, ") = ", field(i,k), " < ", threshold
line 241 : 118 chars : if ( l_print_name ) print *, "filling holes of: reverse_easy_fill_test_"//trim(fill_type_names(fill_type))
line 344 : 102 chars : field_initial(i,k) = 10_core_rknd * ( sin( 8.0 * pi * ( k + i ) / nz ) + one ) &
line 363 : 113 chars : if ( l_print_name ) print *, "filling holes of: below_thresh_test_"//trim(fill_type_names(fill_type))
line 385 : 103 chars : print *, "-- error", two * abs( initial_mass - new_mass ) / ( initial_mass + new_mass )
line 419 : 121 chars : if ( l_print_name ) print *, "filling holes of: reverse_below_thresh_test_"//trim(fill_type_names(fill_type))
line 430 : 116 chars : print *, "ERROR in below_thresh_test REVERSE MODE: method was not conservative in grid reverse mode"
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/spurious_source_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 500 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 6 check(s) failed.
[91m===============================================================
[91mCLUBBStandardsCheck FAILED
[91m THIS IS PRINTED IN ALL RED, CAPITAL LETTERS, AND USES
[91m AN EXCLAMATION MARK TO ENSURE THE DEVELOPERS FEEL SHAME!
[91m IF YOU ARE ONE OF THESE "DEVELOPERS" CHECK THE
[91m LOG FILE FOR DETAILS: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gfortran_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
[93mBuild completed successfully, but some source code checks have failed.[0m
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Run CLUBB's G Unit Tests)
[Pipeline] sh
+ source /etc/profile.d/larson-group.sh
++ export GIT_EDITOR=vi
++ GIT_EDITOR=vi
++ export SVN_EDITOR=vi
++ SVN_EDITOR=vi
++ export OMP_STACKSIZE=1048579
++ OMP_STACKSIZE=1048579
++ export LMOD_ROOT=/opt/lmod/
++ LMOD_ROOT=/opt/lmod/
++ source /opt/lmod//lmod/lmod/init/bash
+++ '[' -z '' ']'
+++ case "$-" in
+++ __lmod_vx=x
+++ '[' -n x ']'
+++ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for this output (/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-12.2.0/lmod-8.7.37-bi3kyxcdrfgw3y7vv2k7c5rjxg75qzju/lmod/lmod/init/bash)
Shell debugging restarted
+++ unset __lmod_vx
+++ find /usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core -print -quit
++ export MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
++ MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
+ module load gcc netcdf-fortran
+ '[' -z '' ']'
+ case "$-" in
+ __lmod_sh_dbg=x
+ '[' -n x ']'
+ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for Lmod's output
Lmod Warning:
-------------------------------------------------------------------------------
The following dependent module(s) are not currently loaded:
gcc-runtime/13.1.0-gzedt3f (required by: mpich/4.2.1, hdf5/1.14.3-42x3zrv,
netcdf-c/4.9.2, netcdf-fortran/4.5.3)
-------------------------------------------------------------------------------
Shell debugging restarted
+ unset __lmod_sh_dbg
+ return 0
+ python run_scripts/run_G_unit_tests.py
Running G_unit_tests
Performing KK integrals tests
Percent Difference Tolerance (maximum acceptable percent difference between the CLUBB results and the MATLAB results using the same equations): 1.50000E-12
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 1
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: 3.940566852245362E-04
MATLAB result: 3.940566852245361E-04
Percent difference: 1.37569E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: 3.927115990964176E-04
MATLAB result: 3.927115990964175E-04
Percent difference: 2.76081E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: 3.959157231937519E-04
MATLAB result: 3.959157231937519E-04
Percent difference: 1.36923E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: 3.961224408774508E-04
MATLAB result: 3.961224408774508E-04
Percent difference: 1.36852E-14%
Agreement with regards to integral: Trivariate NNL
Obtained result: 5.957899876478518E-10
MATLAB result: 5.957899876478559E-10
Percent difference: 6.94188E-13%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 5.299646183560617E-10
MATLAB result: 5.299646183560650E-10
Percent difference: 6.24329E-13%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.543555152546510E-10
MATLAB result: 5.543555152546526E-10
Percent difference: 2.98430E-13%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.286452253001267E-10
MATLAB result: 5.286452253001284E-10
Percent difference: 3.12944E-13%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -9.854231981928352E-04
MATLAB result: -9.854231981928350E-04
Percent difference: 2.20048E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549015E-04
MATLAB result: -9.922448817549015E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 8.940070763287877E-09
MATLAB result: 8.940070763287870E-09
Percent difference: 7.40201E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 8.942709549399747E-09
MATLAB result: 8.942709549399743E-09
Percent difference: 3.69991E-14%
Agreement with regards to integral: Bivariate LL
Obtained result: 1.007487885941809E-02
MATLAB result: 1.007487885941809E-02
Percent difference: 1.72183E-14%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 2
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -4.441972023148827E-05
MATLAB result: -4.441972023148829E-05
Percent difference: 4.57652E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -4.559373570842309E-05
MATLAB result: -4.559373570842311E-05
Percent difference: 4.45868E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -5.040842768062482E-05
MATLAB result: -5.040842768062482E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -5.038775591225493E-05
MATLAB result: -5.038775591225493E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.629614112194526E-09
MATLAB result: 4.629614112194527E-09
Percent difference: 1.78672E-14%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.604609365066508E-09
MATLAB result: 4.604609365066509E-09
Percent difference: 1.79642E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 4.685272472815394E-09
MATLAB result: 4.685272472815394E-09
Percent difference: 1.76549E-14%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 4.683452323676534E-09
MATLAB result: 4.683452323676534E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -1.088125285831538E-04
MATLAB result: -1.088125285831538E-04
Percent difference: 2.49099E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549016E-05
MATLAB result: -9.922448817549016E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 7.907812698669825E-10
MATLAB result: 7.907812698669821E-10
Percent difference: 5.23015E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 6.330953526469324E-10
MATLAB result: 6.330953526469324E-10
Percent difference: 0.00000E+00%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 3
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -7.622010070476262E-05
MATLAB result: -7.622010070476262E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -7.695119321547280E-05
MATLAB result: -7.695119321547280E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.851589269710427E-09
MATLAB result: 4.851589269710427E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.837259568295437E-09
MATLAB result: 4.837259568295437E-09
Percent difference: 1.71002E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -4.609761356905442E-05
MATLAB result: -4.609761356905441E-05
Percent difference: 1.46998E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 3.254808634091251E-10
MATLAB result: 3.254808634091251E-10
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
==================================================
Success!
==================================================
Performing correlation cholesky matrix setup tests
correlation matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
--
correlation matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 1.0000000000000000 0.86602540378443860
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
--
==================================================
Performing correlation matrix approximation tests
correlation cholesky matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.86602540378443860 1.0000000000000000 0.96650635094610959
0.86602540378443860 0.96650635094610959 1.0000000000000000
--
==================================================
Success!
==================================================
Performing hole_filling_hm_one_lev_tests
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
==================================================
Performing fill_holes_hydromet_tests
testset_in(:, 1 ) = 1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 6.0000000000000000
testset_in(:, 2 ) = 6.0000000000000000 5.0000000000000000 4.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000
testset_in(:, 3 ) = 3.0000000000000000 3.0000000000000000 3.0000000000000000 4.0000000000000000 4.0000000000000000 4.0000000000000000
testset_in(:, 4 ) = 1.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 5 ) = 1.0000000000000000 -2.0000000000000000 -4.0000000000000000 -3.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 6 ) = 2.0000000000000000 -1.0000000000000000 -1.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 7 ) = 1.0000000000000000 4.0000000000000000 2.0000000000000000 3.0000000000000000 -1.0000000000000000 0.0000000000000000
Hole-filling was conservative
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
Nc-Ncn equations unit test
==========================
Iteration = 1
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.15865525393145707
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 4.7896918412674905E-002
Nc_in_cloud (from Ncn) = 66978771.649335034
Ncm (from Ncn) = 3208076.7610693830
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 1 is a success!
======================================================
Iteration = 2
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.0000000000000000
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 0.0000000000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 2 is a success!
======================================================
Iteration = 3
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 0.99957093966680322
cloud fraction = 0.95210308158732504
Nc_in_cloud (from Ncn) = 51650441.529586732
Ncm (from Ncn) = 49176544.545665480
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 3 is a success!
======================================================
Iteration = 4
mu_chi_1 = 1.0000000000000000E-008
mu_chi_2 = -1.0000000000000000E-008
sigma_chi_1 = 5.0000000000000002E-005
sigma_chi_2 = 5.0000000000000002E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.10000000000000001
mixture fraction = 0.50000000000000000
cloud fraction (1st PDF component) = 0.50007978845554835
cloud fraction (2nd PDF component) = 0.49992021154445165
cloud fraction = 0.50000000000000000
Nc_in_cloud (from Ncn) = 53982783.648312390
Ncm (from Ncn) = 26991391.824156195
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 4 is a success!
======================================================
Iteration = 5
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 9.9999999999999995E-007
sigma_chi_2 = 1.5000000000000000E-006
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.50000000000000000
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 1.0000000000000000
cloud fraction = 1.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 50000000.000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 5 is a success!
======================================================
Iteration = 6
mu_chi_1 = -9.9999999999999995E-008
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.0000000000000000
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.49601064368536840
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 9.9545377003631155E-002
Nc_in_cloud (from Ncn) = 49999999.999999993
Ncm (from Ncn) = 4977268.8501815572
Ncnm (from Nc_in_cloud) = 49999999.999999993
Ncnm (from Ncm) = 49999999.999999993
Test iteration 6 is a success!
======================================================
Iteration = 7
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 0.0000000000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.48670438618290790
cloud fraction = 0.49732907640173213
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 24866453.820086606
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 7 is a success!
======================================================
Iteration = 8
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.10796556745540581
Nc_in_cloud (from Ncn) = 81037694.237492040
Ncm (from Ncn) = 8749280.6436284967
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 8 is a success!
======================================================
Iteration = 9
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.20034324826655747
Nc_in_cloud (from Ncn) = 52151569.416871324
Ncm (from Ncn) = 10448214.819174858
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 9 is a success!
======================================================
Iteration = 10
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.16861219748026607
Nc_in_cloud (from Ncn) = 60768988.354576342
Ncm (from Ncn) = 10246392.665117815
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 10 is a success!
======================================================
Success!
Performing corr_mtx_test
==================================================
Expected:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Actual:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Success!
==================================================
==================================================
Performing SILHS category test
==================================================
Success!
Performing hydrometeor mu and sigma values unit test
====================================================
PDF parameter set 1:
hmm = 2.0000000000000002E-005
mixt_frac = 2.7000000000000000E-002
precip_frac = 8.0000000000000002E-002
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 5.4470709146968145E-002
hmp2_ip_on_hmm2_ip = 3.6000000000000001
hmm_ip = 2.5000000000000001E-004
hmp2 = 2.2599999999999997E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 2:
hmm = 5.0000000000000004E-006
mixt_frac = 0.50000000000000000
precip_frac = 0.50000000000000000
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 0.50000000000000000
hmp2_ip_on_hmm2_ip = 1.0000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 7.5000000000000012E-011
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 3:
hmm = 1.0000000000000000E-010
mixt_frac = 0.25000000000000000
precip_frac = 0.40000000000000002
precip_frac_1 = 0.75000000000000000
precip_frac_2 = 0.28333333333333338
hmp2_ip_on_hmm2_ip = 0.50000000000000000
hmm_ip = 2.5000000000000002E-010
hmp2 = 2.7500000000000010E-020
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 4:
hmm = 4.0000000000000003E-005
mixt_frac = 1.0000000000000000E-002
precip_frac = 0.10000000000000001
precip_frac_1 = 1.0000000000000000E-002
precip_frac_2 = 0.10090909090909091
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 5:
hmm = 4.0000000000000003E-005
mixt_frac = 0.98999999999999999
precip_frac = 0.10000000000000001
precip_frac_1 = 0.10090909090909091
precip_frac_2 = 1.0000000000000278E-002
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 6:
hmm = 1.0000000000000001E-005
mixt_frac = 0.12500000000000000
precip_frac = 1.0000000000000000
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 1.0000000000000000
hmp2_ip_on_hmm2_ip = 1.5000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 1.5000000000000005E-010
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 7:
hmm = 6.0000000000000002E-005
mixt_frac = 5.0000000000000003E-002
precip_frac = 0.11000000000000000
precip_frac_1 = 0.90000000000000002
precip_frac_2 = 6.8421052631578952E-002
hmp2_ip_on_hmm2_ip = 5.0000000000000000
hmm_ip = 5.4545454545454548E-004
hmp2 = 1.9276363636363637E-007
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 8:
hmm = 7.9999999999999996E-007
mixt_frac = 0.80000000000000004
precip_frac = 0.25000000000000000
precip_frac_1 = 0.10000000000000001
precip_frac_2 = 0.85000000000000009
hmp2_ip_on_hmm2_ip = 0.75000000000000000
hmm_ip = 3.1999999999999999E-006
hmp2 = 3.8399999999999998E-012
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 9:
hmm = 1.3499999999999999E-005
mixt_frac = 0.10000000000000001
precip_frac = 0.10000000000000001
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 5.5555555555555559E-002
hmp2_ip_on_hmm2_ip = 2.0000000000000000
hmm_ip = 1.3499999999999997E-004
hmp2 = 5.2852499999999985E-009
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 10 (randomly generated):
Random seed values = 2026 2 17 -360 22 24 36 520
hmm = 2.4554819017322670E-005
mixt_frac = 0.12495780678209140
precip_frac = 0.68578329175505837
precip_frac_1 = 0.42439014681447595
precip_frac_2 = 0.72311076505041116
hmp2_ip_on_hmm2_ip = 4.1846560916713900
hmm_ip = 3.5805507822277082E-005
hmp2 = 3.9553989080034591E-009
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Checking values with precipitation in 1st PDF component only.
Checking values with precipitation in 2nd PDF component only.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for ADG1
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 2 17 -360 22 24 36 521
wm = -2.8903619653546597E-002
rtm = 3.3460776894304305E-003
thlm = 302.81825012141576
wp2 = 0.41857590583228110
rtp2 = 2.3543529220009550E-006
thlp2 = 2.5196941751792372
Skw = -2.1480522837360678
wp3 = -0.58170954038938694
Skrt = 4.4607573241340308
rtp3 = 1.6114466972066737E-008
Skthl = -2.5941033195279495
thlp3 = -10.375499587309903
wprtp = 4.3598940731977016E-004
wpthlp = -0.89158513056116373
rtpthlp = -2.1920588424000935E-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 = 2026 2 17 -360 22 24 36 521
wm = -2.8903619653546597E-002
rtm = 3.3460776894304305E-003
thlm = 302.81825012141576
wp2 = 0.41857590583228110
rtp2 = 2.3543529220009550E-006
thlp2 = 2.5196941751792372
Skw = -2.1480522837360678
wp3 = -0.58170954038938694
Skrt = 4.4607573241340308
rtp3 = 1.6114466972066737E-008
Skthl = -2.5941033195279495
thlp3 = -10.375499587309903
wprtp = 4.3598940731977016E-004
wpthlp = -0.89158513056116373
rtpthlp = -2.1920588424000935E-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 = 2026 2 17 -360 22 24 36 522
wm = -2.8903619653546597E-002
rtm = 3.3460776894304505E-003
thlm = 302.81825012141570
wp2 = 0.41857590579135417
rtp2 = 2.2620987810632908E-006
thlp2 = 2.5196459586495910
Skw = -2.2469775177272928
wp3 = -0.60849927574264184
Skrt = 3.4764306035279535
rtp3 = 1.1827716682231791E-008
Skthl = -3.4796260732593014
thlp3 = -13.916879228012627
wprtp = 6.2415358048165321E-004
wpthlp = 0.14916630673499065
rtpthlp = 2.1486617514014511E-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 = 2026 2 17 -360 22 24 36 524
wm = -2.8903619653546597E-002
rtm = 3.3460776894304128E-003
thlm = 302.81825012141576
wp2 = 0.41857590570949776
rtp2 = 2.2005960205908564E-006
thlp2 = 2.5196133952013273
Skw = -2.4448279812910405
wp3 = -0.66207874522895138
Skrt = 1.5077771652625618
rtp3 = 4.9220684236038275E-009
Skthl = -1.3131708245090308
thlp3 = -5.2519682023284018
wprtp = 2.0632641727959527E-004
wpthlp = -0.28036352965967376
rtpthlp = -2.1192374762718405E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new hybrid PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 2 17 -360 22 24 36 526
wm = -2.8903619653546597E-002
rtm = 3.3460776894304236E-003
thlm = 302.81825012141576
wp2 = 0.41857590562764391
rtp2 = 2.1390932601184869E-006
thlp2 = 2.5195813346672620
Skw = -2.6426784478015479
wp3 = -0.71565821548183950
Skrt = -3.8358762759495928
rtp3 = -1.2000771586864277E-008
Skthl = 4.2282844238729211
thlp3 = 16.910512386755048
wprtp = -9.3065491994814809E-004
wpthlp = -0.74376862610322836
rtpthlp = 2.0893998553552970E-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.454761250379871
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.2197274440461925E-018
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 793918.42642888101
Vertical integral of <thl> after advance_xm_wpxp = 793908.03522151255
Vertical integral of <thl> forcing = -4.4637357896125940E-002
Flux of <thl> at the surface = 1.0000000000000000E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.2514503322513804E-012
The spurious source of <thl> is within acceptable limits.
Test 2
Vertical integral of <rt> before advance_xm_wpxp = 21.986606976570357
Vertical integral of <rt> after advance_xm_wpxp = 22.122120497121372
Vertical integral of <rt> forcing = -3.8098058894043593E-005
Flux of <rt> at the surface = 4.8980979406412425E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -3.0079834022894714E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 794468.69724861323
Vertical integral of <thl> after advance_xm_wpxp = 794462.37745053775
Vertical integral of <thl> forcing = -3.0862189465910104E-002
Flux of <thl> at the surface = 9.7961958812824851E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -2.9768895681847596E-013
The spurious source of <thl> is within acceptable limits.
Test 3
Vertical integral of <rt> before advance_xm_wpxp = 15.815978654659368
Vertical integral of <rt> after advance_xm_wpxp = 16.075632368806346
Vertical integral of <rt> forcing = -1.9853192243070051E-006
Flux of <rt> at the surface = 8.6749769971425795E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -2.3804166916687602E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 788594.81624482269
Vertical integral of <thl> after advance_xm_wpxp = 788599.45790801046
Vertical integral of <thl> forcing = -1.8777433689123157E-003
Flux of <thl> at the surface = 1.7349953994285159E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 5.1221071654949935E-013
The spurious source of <thl> is within acceptable limits.
Test 4
Vertical integral of <rt> before advance_xm_wpxp = 10.457232112436433
Vertical integral of <rt> after advance_xm_wpxp = 10.619120600283068
Vertical integral of <rt> forcing = -3.2311709705559692E-006
Flux of <rt> at the surface = 5.4285946379266890E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 4.8606985678188025E-018
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 785026.07074451179
Vertical integral of <thl> after advance_xm_wpxp = 785028.44177234545
Vertical integral of <thl> forcing = -2.9537631640834975E-003
Flux of <thl> at the surface = 1.0857189275853378E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 4.4919016423117242E-013
The spurious source of <thl> is within acceptable limits.
Test 5
Vertical integral of <rt> before advance_xm_wpxp = 42.650300837573837
Vertical integral of <rt> after advance_xm_wpxp = 42.833199098120026
Vertical integral of <rt> forcing = -6.4296346423834682E-005
Flux of <rt> at the surface = 6.7395721491111265E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 1.8146833861976130E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 801908.68049078609
Vertical integral of <thl> after advance_xm_wpxp = 801891.71896738349
Vertical integral of <thl> forcing = -7.0017555641658499E-002
Flux of <thl> at the surface = 1.3479144298222253E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.4406253967536031E-012
The spurious source of <thl> is within acceptable limits.
Success!
Running tuner tests
-------------------
goldstein_price pass
rastrigin pass
himmelblau pass
eggholder pass
schaffer pass
==================================================
Performing 'mean vertical velocity in clouds' tests
Success!
==================================================
==================================================
Performing 'Smooth Heaviside Peskin' tests
Case 1: heaviside_smth_range = 0, smooth Heaviside collapses to Heaviside
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
Case 2: heaviside_smth_range = 2.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Expected outcome: 0.0000000000000000 9.0845056908100005E-002 0.90915494309000000 1.0000000000000000
True outcome: 0.0000000000000000 9.0845056909999999E-002 0.90915494309000000 1.0000000000000000
Success!
==================================================
==================================================
Performing 'Smooth Min / Max' tests
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
(Simple) min:
Expected outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
True outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
(Simple) max:
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Smooth min:
Expected outcome: -2.2247448714824141 -1.3660254037844388 -0.36602540378443871 -0.22474487138241406
True outcome: -2.2247448714824136 -1.3660254037844388 -0.36602540378443871 -0.22474487138241384
Smooth max:
Expected outcome: 0.22474487138241406 0.36602540378443871 1.3660254037844388 2.2247448714824141
True outcome: 0.22474487138241384 0.36602540378443871 1.3660254037844388 2.2247448714824136
Testing smooth min and max with min_max_smth_mag=1e-7, input_var1 = equidistant grid on [-1, 1), input_var2 = 0
Input: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
True outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
Smooth max:
Expected outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
True outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
Input_var1: -1000000.0000000000 -10000000.000000000 -100000000.00000000 -1000000000.0000000
Input_var2: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
True outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
Smooth max:
Expected outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
True outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
Input_var1: -1000000.0000000000
Input_var2: 1.0000000000000000
Smooth min:
Expected outcome: -1000000.0000002501
True outcome: -1000000.0000002501
Smooth max:
Expected outcome: 1.0000002500019036
True outcome: 1.0000002500019036
Success!
==================================================
====================== Running fill holes tests ======================
-------------- easy_fill_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: easy_fill_test
-----------------------------------------------
-------------- below_thresh_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: below_thresh_test
-----------------------------------------------
====================== End fill holes tests ======================
Running /home/jenkins/workspace/clubb_G_unit_gfortran_tests/run_scripts/../install/latest/G_unit_tests ...
All tests have succeeded.
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Declarative: Post Actions)
[Pipeline] script
[Pipeline] {
[Pipeline] cleanWs
[WS-CLEANUP] Deleting project workspace...
[WS-CLEANUP] Deferred wipeout is used...
[WS-CLEANUP] done
[Pipeline] }
[Pipeline] // script
[Pipeline] }
[Pipeline] // stage
[Pipeline] }
[Pipeline] // withEnv
[Pipeline] }
[Pipeline] // node
[Pipeline] End of Pipeline
Finished: SUCCESS
