Skipping 316 KB..
Full LogTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTest
/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/clubb_driver.F90:4514:33:
4514 | l_add_dycore_grid, &
| 1
Warning: Unused dummy argument 'l_add_dycore_grid' at (1) [-Wunused-dummy-argument]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/clubb_driver.F90:3122:18:
3122 | l_stats_last
| 1
Warning: Unused variable 'l_stats_last' declared at (1) [-Wunused-variable]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/clubb_driver.F90:3142:14:
3142 | ret_code, &
| 1
Warning: Unused variable 'ret_code' declared at (1) [-Wunused-variable]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/clubb_driver.F90:2637:49:
2637 | subroutine set_case_initial_conditions(err_info)
| 1
Warning: Unused dummy argument 'err_info' at (1) [-Wunused-dummy-argument]
[ 97%] Linking Fortran static library libclubb_driver_lib.a
[ 97%] Built target clubb_driver_lib
Scanning dependencies of target clubb_tuner
Scanning dependencies of target clubb_standalone
Scanning dependencies of target clubb_thread_test
Scanning dependencies of target clubb_driver_test
[ 98%] Building Fortran object src/CMakeFiles/clubb_tuner.dir/error.F90.o
[ 99%] Building Fortran object src/CMakeFiles/clubb_standalone.dir/clubb_standalone.F90.o
[ 99%] Building Fortran object src/CMakeFiles/clubb_thread_test.dir/clubb_thread_test.F90.o
[ 99%] Building Fortran object src/CMakeFiles/clubb_driver_test.dir/clubb_driver_test.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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]
[ 99%] Linking Fortran executable clubb_thread_test
[ 99%] Linking Fortran executable clubb_standalone
[100%] Linking Fortran executable clubb_driver_test
[100%] Built target clubb_thread_test
[100%] Built target clubb_standalone
[100%] Built target clubb_driver_test
[100%] Building Fortran object src/CMakeFiles/clubb_tuner.dir/clubb_tuner.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/clubb_tuner.F90:294:6:
294 | use nrtype, only: &
| 1
Warning: Unused parameter 'sp' which has been explicitly imported at (1) [-Wunused-parameter]
[100%] Linking Fortran executable clubb_tuner
[100%] Built target clubb_tuner
Install the project...
-- Install configuration: "Debug"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/libclubb_driver_lib.a
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_tuner
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_tuner" to "/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i:/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/G_unit_tests
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/G_unit_tests" to "/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i:/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_thread_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_thread_test" to "/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i:/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_standalone
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_standalone" to "/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i:/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_driver_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/install/gfortran_DEBUG_PRECdouble/clubb_driver_test" to "/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i:/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-13.1.0/netcdf-fortran-4.5.3-rr4amopxiouymfpd5tw5ducrmg4u7a6i/lib"
check_for_missing_threadprivate.py has begun.
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/error_code.F90 failed:
Missing threadprivate for:
clubb_debug_level
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/fill_holes.F90 failed:
Missing threadprivate for:
ret_code
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/sponge_layer_damping.F90 failed:
Missing threadprivate for:
thlm_sponge_damp_settings
rtm_sponge_damp_profile
wp2_sponge_damp_profile
wp3_sponge_damp_profile
uv_sponge_damp_settings
wp2_sponge_damp_settings
up2_vp2_sponge_damp_settings
wp3_sponge_damp_settings
thlm_sponge_damp_profile
up2_vp2_sponge_damp_profile
rtm_sponge_damp_settings
uv_sponge_damp_profile
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/corr_varnce_module.F90 failed:
Missing threadprivate for:
corr_array_n_cloud_def
corr_array_n_below_def
[91mcheck_for_missing_threadprivate.py did not pass.[0m
[91mfailedFiles = ['/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/error_code.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/fill_holes.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/sponge_layer_damping.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/corr_varnce_module.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_branch/build/gfortran_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/clubb_driver_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 70 : 107 chars : write(fstderr, *) "This section ensures that everything allocated in init_clubb_case will be deallocated"
line 71 : 102 chars : write(fstderr, *) "in clean_up_clubb. This may cause a runtime error if there is a mismatch between"
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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 1070 : 101 chars : rtm_ref_flip, thlm_ref_flip, um_ref_flip, vm_ref_flip, ug_flip, vg_flip, & ! Intent(in)
line 1083 : 111 chars : um_flip, vm_flip, upwp_flip, vpwp_flip, up2_flip, vp2_flip, up3_flip, vp3_flip, & ! Intent(inout)
line 1085 : 111 chars : wp2_flip, wp3_flip, rtp2_flip, rtp3_flip, thlp2_flip, thlp3_flip, rtpthlp_flip, & ! Intent(inout)
line 1091 : 109 chars : wp2rtp_flip, wp2thlp_flip, uprcp_flip, vprcp_flip, rc_coef_zm_flip, wp4_flip, & ! intent(inout)
line 1092 : 105 chars : wpup2_flip, wpvp2_flip, wp2up2_flip, wp2vp2_flip, ice_supersat_frac_flip, & ! intent(inout)
line 1581 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1608 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2up_implicit, nzt, ngrdcol, &
line 1612 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2up_explicit, nzt, ngrdcol, &
line 1616 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2vp_implicit, nzt, ngrdcol, &
line 1620 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2vp_explicit, nzt, ngrdcol, &
line 1654 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpup2_implicit, nzt, ngrdcol, &
line 1658 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpup2_explicit, nzt, ngrdcol, &
line 1662 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpvp2_implicit, nzt, ngrdcol, &
line 1666 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpvp2_explicit, nzt, ngrdcol, &
line 2470 : 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_branch/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.
3103 : use gptl
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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 1745 : 102 chars : open(unit=iunit_grid_adaptation, file=fname_grid_adaptation, status='replace', action='write')
line 2060 : 105 chars : allocate( um_ref(ngrdcol, gr%nzt) ) ! Reference u wind for nudging; Michael Falk, 17 Oct 2007
line 2061 : 105 chars : allocate( vm_ref(ngrdcol, gr%nzt) ) ! Reference v wind for nudging; Michael Falk, 17 Oct 2007
line 2498 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 2529 : 106 chars : ! Similar call to above, checking parameter values, but because we've placed this in "init_clubb_case"
line 2537 : 116 chars : write(fstderr, *) "Fatal error calling check_clubb_settings_api and/or check_parameters_api in clubb_driver"
line 2580 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 2592 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 2593 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 2596 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 2598 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 2604 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 2605 : 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 2606 : 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 2607 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 2608 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 2609 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 2610 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 2611 : 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 2612 : 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 2613 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 2722 : 103 chars : invrs_rho_ds_zt(i,k) = invrs_rho_ds_zt_init(k) ! inv. dry, static density: t-levs
line 2728 : 103 chars : um_pert(i,k) = zero ! Variables used to track perturbed version of winds.
line 2730 : 105 chars : wm_zt(i,k) = wm_zt_init(k) ! Imposed large scale w - Thermodynamic levels
line 2924 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 3115 : 101 chars : !----------------------------------- Optional Input Variables -----------------------------------
line 3270 : 103 chars : um(i,:), vm(i,:), rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), &
line 3603 : 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 3605 : 104 chars : !$acc lh_rt_clipped, lh_thl_clipped, lh_rc_clipped, lh_rv_clipped, lh_Nc_clipped, &
line 3609 : 116 chars : !$acc pdf_params%cloud_frac_1, pdf_params%cloud_frac_2, pdf_params%cthl_1, pdf_params%cthl_2, &
line 3610 : 112 chars : !$acc pdf_params%crt_1, pdf_params%crt_2, pdf_params%stdev_chi_1, pdf_params%stdev_chi_2, &
line 3626 : 102 chars : wp2_zt = zm2zt_api( gr%nzm, gr%nzt, ngrdcol, gr, wp2, w_tol_sqd ) ! Positive definite quantity
line 3632 : 109 chars : thlm(i,:), p_in_Pa(i,:), exner(i,:), rho(i,:), rho_zm(i,:), rtm(i,:), &! In
line 3633 : 109 chars : rcm(i,:), cloud_frac(i,:), wm_zt(i,:), wm_zm(i,:), wp2_zt(i,:), & ! In
line 3634 : 106 chars : hydromet(i,:,:), Nc_in_cloud(i,:), & ! In
line 3651 : 102 chars : Nccnm(i,:), & ! Inout
line 3652 : 117 chars : hydromet_mc(i,:,:), Ncm_mc(i,:), rcm_mc(i,:), rvm_mc(i,:), & ! Out
line 3653 : 107 chars : thlm_mc(i,:), hydromet_vel_zt(i,:,:), & ! Out
line 3654 : 102 chars : hydromet_vel_covar_zt_impc(i,:,:), & ! Out
line 3655 : 102 chars : hydromet_vel_covar_zt_expc(i,:,:), & ! Out
line 3656 : 110 chars : wprtp_mc(i,:), wpthlp_mc(i,:), rtp2_mc(i,:), & ! Out
line 3657 : 105 chars : thlp2_mc(i,:), rtpthlp_mc(i,:) ) ! Out
line 3681 : 103 chars : stats, i, & ! Inout
line 3682 : 103 chars : hydromet(i,:,:), hydromet_vel_zt(i,:,:), hydrometp2(i,:,:), & ! Inout
line 3683 : 103 chars : K_hm(i,:,:), Ncm(i,:), Nc_in_cloud(i,:), rvm_mc(i,:), & ! Inout
line 3684 : 103 chars : thlm_mc(i,:), err_info, & ! Inout
line 3685 : 101 chars : wphydrometp(i,:,:), wpNcp(i,:) ) ! Out
line 3694 : 106 chars : !$acc update device( rcm_mc, rvm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3747 : 104 chars : !$acc update host( rvm_mc, rcm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3772 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, X_nl_all_levs, &
line 3798 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, thlm, rtm, rcm )
line 3899 : 112 chars : call stats_update( "brunt_term_time_avg", richardson_num_term_weight*richardson_num_term_time_avg, &
line 3925 : 101 chars : .and. stats%l_last_sample ) then ! only adapt grid when the average of the last
line 4053 : 103 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_scheme = ', time_microphys_scheme
line 4054 : 104 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_advance = ', time_microphys_advance
line 4173 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 4185 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 4186 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 4189 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 4191 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 4197 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 4198 : 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 4199 : 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 4200 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 4201 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 4202 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 4203 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 4204 : 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 4205 : 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 4206 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 4228 : 119 chars : if ( thlm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( thlm_sponge_damp_profile )
line 4229 : 118 chars : if ( rtm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( rtm_sponge_damp_profile )
line 4230 : 117 chars : if ( uv_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( uv_sponge_damp_profile )
line 4231 : 118 chars : if ( wp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp2_sponge_damp_profile )
line 4232 : 118 chars : if ( wp3_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp3_sponge_damp_profile )
line 4233 : 122 chars : if ( up2_vp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( up2_vp2_sponge_damp_profile )
line 4920 : 104 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 4928 : 106 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 4942 : 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_branch/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
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 6 check(s) failed.
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/advance_wp2_wp3_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 925 : 101 chars : stats, & ! intent(inout)
line 1712 : 105 chars : stats_tmp_zt(i,k) = - gamma_over_implicit_ts * wp3_term_ta_lhs_result(5,i,k) * wp3(i,k-1) &
line 1713 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(4,i,k) * wp2(i,k) &
line 1714 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(3,i,k) * wp3(i,k) &
line 1715 : 105 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(2,i,k) * wp2(i,k+1) &
line 1716 : 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_branch/src/CLUBB_core/sfc_varnce_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 722 : 105 chars : call sfc_varnce_check( sclr_dim, wp2(i,gr%k_lb_zm), up2(i,gr%k_lb_zm), & ! intent(in)
line 723 : 105 chars : vp2(i,gr%k_lb_zm), thlp2(i,gr%k_lb_zm), rtp2(i,gr%k_lb_zm), & ! intent(in)
line 724 : 105 chars : rtpthlp(i,gr%k_lb_zm), sclrp2(i,gr%k_lb_zm,:), & ! intent(in)
line 725 : 105 chars : sclrprtp(i,gr%k_lb_zm,:), sclrpthlp(i,gr%k_lb_zm,:), & ! intent(in)
line 726 : 108 chars : err_info ) ! intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/clubb_api_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 2410 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
line 2477 : 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_branch/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 )
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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_branch/src/CLUBB_core/mono_flux_limiter.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 698 : 104 chars : !$acc parallel loop gang vector collapse(2) default(present) reduction(.or.:l_any_adjustment_needed)
line 849 : 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_branch/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_branch/src/CLUBB_core/advance_clubb_core_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1868 : 101 chars : stats, & ! intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/src/CLUBB_core/advance_windm_edsclrm_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 492 : 102 chars : stats ) ! intent(inout)
line 500 : 102 chars : stats ) ! intent(inout)
line 952 : 103 chars : call windm_edsclrm_rhs( nzm, nzt, ngrdcol, gr, windm_edsclrm_scalar, dt, & ! intent(in)
line 953 : 103 chars : lhs_diff, edsclrm(:,:,edsclr), & ! intent(in)
line 954 : 103 chars : edsclrm_forcing(:,:,edsclr), & ! intent(in)
line 955 : 103 chars : rho_ds_zm, invrs_rho_ds_zt, & ! intent(in)
line 956 : 103 chars : l_imp_sfc_momentum_flux, wpedsclrp(:,gr%k_lb_zm,edsclr), & ! intent(in)
line 957 : 106 chars : stats, & ! intent(inout)
line 958 : 104 chars : rhs(:,:,edsclr) ) ! intent(out)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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."
----------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/stats_netcdf.F90 ---------------------
WARNING: Missing 'implicit none' after these blocks:
line 219: subroutine stats_init( registry_path, output_path, ncol, stats_tsamp, stats_tout, dt_main, &
line 402: subroutine stats_finalize( stats, err_info )
line 469: subroutine stats_open_netcdf( output_path, ncol, defs, ndefs, time_units, &
line 767: subroutine stats_begin_timestep( itime, stats )
line 803: subroutine stats_end_timestep( time_value, stats, err_info )
line 853: subroutine stats_update_scalar( name, values, stats, icol, level )
line 912: subroutine stats_update_1d( name, values, stats, icol, level )
line 971: subroutine stats_update_2d( name, values, stats, icol, level )
line 1033: subroutine stats_begin_budget_scalar( name, values, stats, icol )
line 1078: subroutine stats_begin_budget_1d( name, values, stats, icol )
line 1130: subroutine stats_begin_budget_2d( name, values, stats, icol )
line 1174: subroutine stats_update_budget_scalar( name, values, stats, icol, level )
line 1223: subroutine stats_update_budget_1d( name, values, stats, icol, level )
line 1272: subroutine stats_update_budget_2d( name, values, stats, icol, level )
line 1315: subroutine stats_finalize_budget_scalar( name, values, stats, icol, l_count_sample )
line 1379: subroutine stats_finalize_budget_1d( name, values, stats, icol, l_count_sample )
line 1450: subroutine stats_finalize_budget_2d( name, values, stats, icol, l_count_sample )
line 1514: subroutine stats_write( time_value, stats, ierr )
line 1635: subroutine stats_grid_init( ncol, zt_src, zm_src, zt_tgt, zm_tgt, &
line 1696: subroutine stats_update_grid( zt_src, zm_src, rho_vals, rho_levels, p_sfc, stats )
line 1773: subroutine stats_lh_samples_init( num_samples, nzt, nl_var_names, u_var_names, zt_vals, &
line 1887: subroutine stats_lh_samples_write_lognormal( samples, stats, err_info )
line 1935: subroutine stats_lh_samples_write_uniform( uniform_vals, mixture_comp, sample_weights, &
line 2143: subroutine stats_lookup_append( stats, id )
line 2174: subroutine stats_lookup_reject_append( stats, name )
line 2277: subroutine stats_read_registry_namelist( path, defs, ndefs, ierr )
line 2343: subroutine stats_expand_registry( sclr_dim, edsclr_dim, defs, ndefs, hydromet_list )
line 2633: subroutine add_expanded_def( base_def, name, out_defs, nout )
line 2674: subroutine parse_registry_line( line, def, l_valid )
line 2736: subroutine split_registry_fields( line, fields, nfields )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/stats_netcdf.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 134 : 105 chars : real(kind=stat_rknd), allocatable, dimension(:,:) :: buffer ! (col, height); height=1 for rank-1 vars
line 179 : 103 chars : logical :: l_last_sample = .false. ! True if we want to output at the end of this sample
line 180 : 104 chars : logical :: l_different_output_grid = .false. ! True if we need to remap the variables before output
line 181 : 102 chars : logical :: l_output_rad_files = .false. ! True if we want to output radiation grid variables
line 435 : 101 chars : if ( allocated( stats%grid%rho_lin_spline_vals ) ) deallocate( stats%grid%rho_lin_spline_vals )
line 436 : 105 chars : if ( allocated( stats%grid%rho_lin_spline_levels ) ) deallocate( stats%grid%rho_lin_spline_levels )
line 581 : 104 chars : ret_code = nf90_def_var( stats%ncid, "time", NF90_DOUBLE, (/ stats%time_dimid /), stats%time_varid )
line 586 : 101 chars : ret_code = nf90_def_var( stats%ncid, "col", NF90_DOUBLE, (/ stats%col_dimid /), stats%col_varid )
line 606 : 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 613 : 112 chars : ret_code = nf90_def_var( stats%ncid, "rad_zt", NF90_DOUBLE, (/ stats%rad_zt_dimid /), stats%rad_zt_varid )
line 617 : 112 chars : ret_code = nf90_def_var( stats%ncid, "rad_zm", NF90_DOUBLE, (/ stats%rad_zm_dimid /), stats%rad_zm_varid )
line 624 : 109 chars : ret_code = nf90_def_var( stats%ncid, "param", NF90_DOUBLE, (/ stats%param_dimid /), stats%param_varid )
line 632 : 101 chars : (/ stats%param_strlen_dimid, stats%param_dimid /), stats%param_name_varid )
line 722 : 110 chars : ret_code = nf90_put_att( stats%ncid, stats%vars(i)%varid, "long_name", trim( stats%vars(i)%long_name ) )
line 723 : 102 chars : ret_code = nf90_put_att( stats%ncid, stats%vars(i)%varid, "units", trim( stats%vars(i)%units ) )
line 736 : 108 chars : ret_code = nf90_put_var( stats%ncid, stats%col_varid, [( real( i, kind=core_rknd ), i=1, stats%ncol )] )
line 756 : 109 chars : ret_code = nf90_put_var( stats%ncid, stats%param_varid, [( real( i, kind=core_rknd ), i=1, nparams )] )
line 1594 : 101 chars : stats%grid%gr_source, stats%grid%gr_target, stats%vars(i)%size, &
line 1595 : 107 chars : real( stats%vars(i)%buffer, kind=core_rknd ), stats%vars(i)%out_size, &
line 1596 : 112 chars : size( stats%grid%rho_lin_spline_vals, 2 ), stats%grid%rho_lin_spline_vals, &
line 1751 : 101 chars : if ( allocated( stats%grid%rho_lin_spline_vals ) ) deallocate( stats%grid%rho_lin_spline_vals )
line 1752 : 105 chars : if ( allocated( stats%grid%rho_lin_spline_levels ) ) deallocate( stats%grid%rho_lin_spline_levels )
line 2166 : 116 chars : if ( stats%lookup%cache_len > 0 ) tmp(1:stats%lookup%cache_len) = stats%lookup%cache(1:stats%lookup%cache_len)
line 2200 : 105 chars : tmp(1:stats%lookup%reject_cache_len) = stats%lookup%reject_cache(1:stats%lookup%reject_cache_len)
line 2266 : 101 chars : write(date(32:33),'(i2.2)') nint( ( ( real( mod( nint( st_time ),3600 ),kind=time_precision ) / &
line 2334 : 104 chars : write( fstderr,* ) "Invalid stats_registry entry (expected: name | grid | long_name [| units]):"
line 2424 : 101 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2425 : 101 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2432 : 101 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2433 : 101 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2440 : 101 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2441 : 101 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2448 : 113 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_1", &
line 2452 : 113 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_2", &
line 2458 : 101 chars : call add_expanded_def(base_def,"corr_w_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2459 : 101 chars : call add_expanded_def(base_def,"corr_w_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2466 : 103 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2467 : 103 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2474 : 103 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2475 : 103 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2482 : 103 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2483 : 103 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2490 : 115 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_1_n", &
line 2494 : 115 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_2_n", &
line 2531 : 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_branch/src/CLUBB_core/advance_xp2_xpyp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1310 : 101 chars : sclrprtp(:,:,sclr) ) ! Intent(inout)
line 1343 : 102 chars : sclrpthlp(:,:,sclr) ) ! Intent(inout)
line 3011 : 101 chars : stats_tmp(i,k) = (-gamma_over_implicit_ts * lhs_ta(2+gr%grid_dir_indx,i,k)) * xapxbp(i,km1) &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/CLUBB_core/advance_xm_wpxp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 777 : 101 chars : call wpxp_term_pr1_lhs( nzm, ngrdcol, gr, C6rt_Skw_fnc, C6thl_Skw_fnc, C7_Skw_fnc, & ! Intent(in)
line 778 : 101 chars : invrs_tau_C6_zm, l_scalar_calc, & ! Intent(in)
line 779 : 102 chars : lhs_pr1_wprtp, lhs_pr1_wpthlp, lhs_pr1_wpsclrp ) ! Intent(out)
line 852 : 102 chars : stats, & ! InOut
line 853 : 101 chars : rtm, wprtp, thlm, wpthlp, sclrm, wpsclrp, & ! InOut
line 854 : 101 chars : err_info ) ! InOut
line 894 : 101 chars : stats, & ! InOut
line 4833 : 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 4835 : 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_branch/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_branch/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)
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 19 check(s) failed.
CLUBBStandardsCheck.py has begun.
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/src/SILHS/latin_hypercube_driver_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 407 : 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_branch/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_branch/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)
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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 909 : 103 chars : !$acc exit data delete( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 1112 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1116 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1120 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1121 : 112 chars : exner_zm(i,km1), exner_zm(i,k00), exner_zm(i,kp1), exner_zm(i,kp2) )
line 1124 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1125 : 108 chars : thlm_zm(i,km1), thlm_zm(i,k00), thlm_zm(i,kp1), thlm_zm(i,kp2) )
line 1128 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1129 : 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_branch/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_branch/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_branch/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)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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) )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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_branch/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)
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 11 check(s) failed.
No matches for /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/src/KK_microphys/
CLUBBStandardsCheck.py has begun.
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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"
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests_branch/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_branch/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-skylake_avx512/gcc-12.2.0/lmod-8.7.37-fq24mybyn2li6got2bxzk62ejh5d3p4z/lmod/lmod/init/bash)
Shell debugging restarted
+++ unset __lmod_vx
+++ find /usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core -print -quit
++ export MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
++ MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
+ module load gcc netcdf-fortran
+ '[' -z '' ']'
+ case "$-" in
+ __lmod_sh_dbg=x
+ '[' -n x ']'
+ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for Lmod's output
Lmod Warning:
-------------------------------------------------------------------------------
The following dependent module(s) are not currently loaded:
gcc-runtime/13.1.0-wiuv7wc (required by: mpich/4.2.1, hdf5/1.14.3-nhc23hl,
netcdf-c/4.9.2-xduzp77, netcdf-fortran/4.5.3), glibc/2.35-tdr5czg (required by:
mpich/4.2.1, hdf5/1.14.3-nhc23hl, netcdf-c/4.9.2-xduzp77, 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 13 -360 21 51 57 300
hmm = 2.4554617851646228E-005
mixt_frac = 0.12495517821733178
precip_frac = 0.27215525713949634
precip_frac_1 = 0.30471981394725883
precip_frac_2 = 0.26750508391614913
hmp2_ip_on_hmm2_ip = 0.39359898640856472
hmm_ip = 9.0222831297579796E-005
hmp2 = 2.4844320205484118E-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 13 -360 21 51 57 301
wm = -2.8907642967075486E-002
rtm = 3.3460240452516756E-003
thlm = 293.36440089142599
wp2 = 0.29769678180583481
rtp2 = 1.4364156227613172E-007
thlp2 = 2.1590458347221055
Skw = -3.0585641135702049
wp3 = -0.49679680978429036
Skrt = 2.8861442891419795
rtp3 = 1.5712244161450199E-010
Skthl = 3.6837177161702801
thlp3 = 11.686356643718725
wprtp = -1.8122254540922970E-004
wpthlp = -0.70036711247698080
rtpthlp = 5.0120281378645420E-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).
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 13 -360 21 51 57 302
wm = -2.8907642967075486E-002
rtm = 3.3460240452516578E-003
thlm = 293.36440089142593
wp2 = 0.29769678184676185
rtp2 = 5.1390379343851672E-008
thlp2 = 2.1589977436754397
Skw = -3.1573927880376131
wp3 = -0.51284936539241144
Skrt = 1.9016244487501393
rtp3 = 2.2153762096053673E-011
Skthl = -2.9410138142692461
thlp3 = -9.3298647301186097
wprtp = 6.4862353628165884E-005
wpthlp = 0.74198124028836909
rtpthlp = 2.9978473538607754E-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.
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 13 -360 21 51 57 303
wm = -2.8907642967075486E-002
rtm = 3.3460240452516622E-003
thlm = 293.36440089142599
wp2 = 0.29769678188768878
rtp2 = 8.2140773807350252E-008
thlp2 = 2.1590137741880362
Skw = -3.2562214610316400
wp3 = -0.52890192076783427
Skrt = -0.77020227111924733
rtp3 = -1.8131885989550379E-011
Skthl = -0.17023753221022186
thlp3 = -0.54005555138904349
wprtp = -1.1701507396471073E-005
wpthlp = 0.26266425785289466
rtpthlp = -3.7900907677620565E-004
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 2 13 -360 21 51 57 307
wm = -2.8907642967075486E-002
rtm = 3.3460240452516379E-003
thlm = 293.36440089142576
wp2 = 0.29769678237882458
rtp2 = 9.4315181420312722E-007
thlp2 = 2.1592056384072946
Skw = -2.0710498554118222
wp3 = -0.33639672926206204
Skrt = 2.0413321410022984
rtp3 = 1.8697609287037594E-009
Skthl = 0.45310470949674553
thlp3 = 1.4376050333110049
wprtp = -3.7084954025056144E-004
wpthlp = 0.77705457012042389
rtpthlp = -1.2843409046013384E-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 13 -360 21 51 57 309
wm = -2.8907642967075486E-002
rtm = 3.3460240452516245E-003
thlm = 293.36440089142582
wp2 = 0.29769678213326045
rtp2 = 8.8165102573434069E-007
thlp2 = 2.1591730749590350
Skw = -2.2687071999264914
wp3 = -0.36850183918556639
Skrt = 1.7599855836878344
rtp3 = 1.4569819991127269E-009
Skthl = -3.4008505392549702
thlp3 = -10.789932549385666
wprtp = 4.4753006020209463E-004
wpthlp = -0.47347743748064886
rtpthlp = -1.2417512216958916E-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 = 22.142180434274024
Vertical integral of <rt> after advance_xm_wpxp = 22.238970503607820
Vertical integral of <rt> forcing = -3.3939610803605483E-005
Flux of <rt> at the surface = 3.5657317524961799E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -2.5959865767449797E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 791632.59694389009
Vertical integral of <thl> after advance_xm_wpxp = 791623.86191870086
Vertical integral of <thl> forcing = -3.6248214135958183E-002
Flux of <thl> at the surface = 7.1314635049923589E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.8290924330699454E-013
The spurious source of <thl> is within acceptable limits.
Test 3
Vertical integral of <rt> before advance_xm_wpxp = 10.359280447986295
Vertical integral of <rt> after advance_xm_wpxp = 10.479609584312751
Vertical integral of <rt> forcing = -5.5682660694087893E-006
Flux of <rt> at the surface = 4.0666538715762834E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -3.3525564052325207E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 785117.61519881606
Vertical integral of <thl> after advance_xm_wpxp = 785117.76585894160
Vertical integral of <thl> forcing = -7.6311073251430469E-003
Flux of <thl> at the surface = 8.1333077431525666E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 4.5283655297767567E-013
The spurious source of <thl> is within acceptable limits.
Test 4
Vertical integral of <rt> before advance_xm_wpxp = 26.734224989283177
Vertical integral of <rt> after advance_xm_wpxp = 26.889697158268412
Vertical integral of <rt> forcing = -2.9133900610661250E-005
Flux of <rt> at the surface = 5.4737446389475886E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 1.6723818510588906E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 791538.53223353240
Vertical integral of <thl> after advance_xm_wpxp = 791530.10444203240
Vertical integral of <thl> forcing = -3.9040127610256645E-002
Flux of <thl> at the surface = 1.0947489277895177E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -9.6660179860208473E-013
The spurious source of <thl> is within acceptable limits.
Test 5
Vertical integral of <rt> before advance_xm_wpxp = 17.971720063150372
Vertical integral of <rt> after advance_xm_wpxp = 18.047590627208393
Vertical integral of <rt> forcing = -2.5817201659435672E-005
Flux of <rt> at the surface = 2.7871908185278136E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 6.0085129146431049E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 788951.44936560316
Vertical integral of <thl> after advance_xm_wpxp = 788944.18695791590
Vertical integral of <thl> forcing = -2.9782407260605888E-002
Flux of <thl> at the surface = 5.5743816370556276E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -6.4919944420260833E-013
The spurious source of <thl> is within acceptable limits.
Success!
Running tuner tests
-------------------
goldstein_price pass
rastrigin pass
himmelblau pass
eggholder pass
schaffer pass
==================================================
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_branch/run_scripts/../install/latest/G_unit_tests ...
All tests have succeeded.
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Declarative: Post Actions)
[Pipeline] script
[Pipeline] {
[Pipeline] }
[Pipeline] // script
[Pipeline] }
[Pipeline] // stage
[Pipeline] }
[Pipeline] // withEnv
[Pipeline] }
[Pipeline] // node
[Pipeline] End of Pipeline
Finished: SUCCESS
