Skipping 309 KB..
Full LogTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTest
[100%] Built target clubb_standalone
[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/install/gcc_DEBUG_PRECdouble/libclubb_driver_lib.a
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gcc_DEBUG_PRECdouble/clubb_tuner
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gcc_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/install/gcc_DEBUG_PRECdouble/G_unit_tests
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gcc_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/install/gcc_DEBUG_PRECdouble/clubb_thread_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gcc_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/install/gcc_DEBUG_PRECdouble/clubb_standalone
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gcc_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/install/gcc_DEBUG_PRECdouble/clubb_driver_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gcc_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/src/CLUBB_core/error_code.F90 failed:
Missing threadprivate for:
clubb_debug_level
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90 failed:
Missing threadprivate for:
ret_code
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sponge_layer_damping.F90 failed:
Missing threadprivate for:
thlm_sponge_damp_settings
thlm_sponge_damp_profile
wp3_sponge_damp_profile
uv_sponge_damp_profile
wp2_sponge_damp_settings
rtm_sponge_damp_settings
up2_vp2_sponge_damp_profile
wp3_sponge_damp_settings
up2_vp2_sponge_damp_settings
rtm_sponge_damp_profile
wp2_sponge_damp_profile
uv_sponge_damp_settings
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/corr_varnce_module.F90 failed:
Missing threadprivate for:
corr_array_n_below_def
corr_array_n_cloud_def
[91mcheck_for_missing_threadprivate.py did not pass.[0m
[91mfailedFiles = ['/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/error_code.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sponge_layer_damping.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/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/build/gcc_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 71 : 107 chars : write(fstderr, *) "This section ensures that everything allocated in init_clubb_case will be deallocated"
line 72 : 102 chars : write(fstderr, *) "in clean_up_clubb. This may cause a runtime error if there is a mismatch between"
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/generalized_grid_test.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 3: module generalized_grid_test
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/generalized_grid_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 903 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1071 : 101 chars : rtm_ref_flip, thlm_ref_flip, um_ref_flip, vm_ref_flip, ug_flip, vg_flip, & ! Intent(in)
line 1084 : 111 chars : um_flip, vm_flip, upwp_flip, vpwp_flip, up2_flip, vp2_flip, up3_flip, vp3_flip, & ! Intent(inout)
line 1086 : 111 chars : wp2_flip, wp3_flip, rtp2_flip, rtp3_flip, thlp2_flip, thlp3_flip, rtpthlp_flip, & ! Intent(inout)
line 1092 : 109 chars : wp2rtp_flip, wp2thlp_flip, uprcp_flip, vprcp_flip, rc_coef_zm_flip, wp4_flip, & ! intent(inout)
line 1093 : 105 chars : wpup2_flip, wpvp2_flip, wp2up2_flip, wp2vp2_flip, ice_supersat_frac_flip, & ! intent(inout)
line 1582 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1609 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2up_implicit, nzt, ngrdcol, &
line 1613 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2up_explicit, nzt, ngrdcol, &
line 1617 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2vp_implicit, nzt, ngrdcol, &
line 1621 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2vp_explicit, nzt, ngrdcol, &
line 1655 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpup2_implicit, nzt, ngrdcol, &
line 1659 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpup2_explicit, nzt, ngrdcol, &
line 1663 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpvp2_implicit, nzt, ngrdcol, &
line 1667 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpvp2_explicit, nzt, ngrdcol, &
line 2471 : 105 chars : tolerance = max( 5.e-8_core_rknd, epsilon(tolerance) ) ! max statement for single precision runs
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90 ---------------------
WARNING: 'use' statements without 'only' found in the following lines:
WARNING: Check that comma is on same line as 'use', as CLUBB requires.
3107 : use gptl
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 197 : 110 chars : l_modify_ic_with_cubic_int, & ! Flag for interpolating the sounding profile with Steffen's monotone cubic
line 198 : 110 chars : ! method to obtain smoother initial condition profile, which is found to be
line 199 : 111 chars : ! beneficial to achive a better numerical solution convergence. If this flag
line 200 : 118 chars : ! is turned off, the initial conditions will be generated with linear interpolation.
line 201 : 111 chars : ! This is done on a case-by-case basis, since using the monotone cubic method
line 202 : 109 chars : ! requires a special sounding.in file with many additional sounding levels.
line 203 : 109 chars : l_modify_bc_for_cnvg_test, & ! Flag to activate modifications on boundary condition for convergence test
line 714 : 108 chars : !$omp wm_zm_init, em_init, exner_init, thvm_init, p_in_Pa_init, p_in_Pa_zm_init, rho_init, rho_zm_init, &
line 722 : 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 729 : 101 chars : !$omp exner, rtm, thlm, rcm, wp3, wp3_zm, delta_zm, p_in_Pa, cloud_frac, wpthvp, wp2thvp, wp2up, &
line 738 : 102 chars : !$omp rtp2_mc, thlp2_mc, rtpthlp_mc, Ncm_mc, Ncm, Nccnm, em, tau_zm, tau_zt, Kh_zt, Kh_zm, rfrzm, &
line 740 : 102 chars : !$omp wpsclrp_sfc, wpedsclrp_sfc, sclrm, wpsclrp, sclrp2, sclrp3, sclrprtp, sclrpthlp, sclrpthvp, &
line 1002 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 1733 : 102 chars : open(unit=iunit_grid_adaptation, file=fname_grid_adaptation, status='replace', action='write')
line 2049 : 105 chars : allocate( um_ref(ngrdcol, gr%nzt) ) ! Reference u wind for nudging; Michael Falk, 17 Oct 2007
line 2050 : 105 chars : allocate( vm_ref(ngrdcol, gr%nzt) ) ! Reference v wind for nudging; Michael Falk, 17 Oct 2007
line 2500 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 2531 : 106 chars : ! Similar call to above, checking parameter values, but because we've placed this in "init_clubb_case"
line 2539 : 116 chars : write(fstderr, *) "Fatal error calling check_clubb_settings_api and/or check_parameters_api in clubb_driver"
line 2582 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 2594 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 2595 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 2598 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 2600 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 2606 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 2607 : 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 2608 : 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 2609 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 2610 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 2611 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 2612 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 2613 : 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 2614 : 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 2615 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 2731 : 103 chars : invrs_rho_ds_zt(i,k) = invrs_rho_ds_zt_init(k) ! inv. dry, static density: t-levs
line 2737 : 103 chars : um_pert(i,k) = zero ! Variables used to track perturbed version of winds.
line 2739 : 105 chars : wm_zt(i,k) = wm_zt_init(k) ! Imposed large scale w - Thermodynamic levels
line 2921 : 103 chars : open( unit=iunit_grid_adaptation+10, file=trim( output_dir ) // trim( runtype ) // '_grid.txt', &
line 2933 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 3119 : 101 chars : !----------------------------------- Optional Input Variables -----------------------------------
line 3276 : 103 chars : um(i,:), vm(i,:), rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), &
line 3591 : 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 3593 : 104 chars : !$acc lh_rt_clipped, lh_thl_clipped, lh_rc_clipped, lh_rv_clipped, lh_Nc_clipped, &
line 3597 : 116 chars : !$acc pdf_params%cloud_frac_1, pdf_params%cloud_frac_2, pdf_params%cthl_1, pdf_params%cthl_2, &
line 3598 : 112 chars : !$acc pdf_params%crt_1, pdf_params%crt_2, pdf_params%stdev_chi_1, pdf_params%stdev_chi_2, &
line 3614 : 102 chars : wp2_zt = zm2zt_api( gr%nzm, gr%nzt, ngrdcol, gr, wp2, w_tol_sqd ) ! Positive definite quantity
line 3620 : 109 chars : thlm(i,:), p_in_Pa(i,:), exner(i,:), rho(i,:), rho_zm(i,:), rtm(i,:), &! In
line 3621 : 109 chars : rcm(i,:), cloud_frac(i,:), wm_zt(i,:), wm_zm(i,:), wp2_zt(i,:), & ! In
line 3622 : 106 chars : hydromet(i,:,:), Nc_in_cloud(i,:), & ! In
line 3639 : 102 chars : Nccnm(i,:), & ! Inout
line 3640 : 117 chars : hydromet_mc(i,:,:), Ncm_mc(i,:), rcm_mc(i,:), rvm_mc(i,:), & ! Out
line 3641 : 107 chars : thlm_mc(i,:), hydromet_vel_zt(i,:,:), & ! Out
line 3642 : 102 chars : hydromet_vel_covar_zt_impc(i,:,:), & ! Out
line 3643 : 102 chars : hydromet_vel_covar_zt_expc(i,:,:), & ! Out
line 3644 : 110 chars : wprtp_mc(i,:), wpthlp_mc(i,:), rtp2_mc(i,:), & ! Out
line 3645 : 105 chars : thlp2_mc(i,:), rtpthlp_mc(i,:) ) ! Out
line 3669 : 103 chars : stats, i, & ! Inout
line 3670 : 103 chars : hydromet(i,:,:), hydromet_vel_zt(i,:,:), hydrometp2(i,:,:), & ! Inout
line 3671 : 103 chars : K_hm(i,:,:), Ncm(i,:), Nc_in_cloud(i,:), rvm_mc(i,:), & ! Inout
line 3672 : 103 chars : thlm_mc(i,:), err_info, & ! Inout
line 3673 : 101 chars : wphydrometp(i,:,:), wpNcp(i,:) ) ! Out
line 3682 : 106 chars : !$acc update device( rcm_mc, rvm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3735 : 104 chars : !$acc update host( rvm_mc, rcm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3824 : 112 chars : call stats_update( "brunt_term_time_avg", richardson_num_term_weight*richardson_num_term_time_avg, &
line 3850 : 101 chars : .and. stats%l_last_sample ) then ! only adapt grid when the average of the last
line 3979 : 103 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_scheme = ', time_microphys_scheme
line 3980 : 104 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_advance = ', time_microphys_advance
line 4100 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 4112 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 4113 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 4116 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 4118 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 4124 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 4125 : 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 4126 : 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 4127 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 4128 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 4129 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 4130 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 4131 : 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 4132 : 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 4133 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 4155 : 119 chars : if ( thlm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( thlm_sponge_damp_profile )
line 4156 : 118 chars : if ( rtm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( rtm_sponge_damp_profile )
line 4157 : 117 chars : if ( uv_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( uv_sponge_damp_profile )
line 4158 : 118 chars : if ( wp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp2_sponge_damp_profile )
line 4159 : 118 chars : if ( wp3_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp3_sponge_damp_profile )
line 4160 : 122 chars : if ( up2_vp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( up2_vp2_sponge_damp_profile )
line 4835 : 104 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 4843 : 106 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 4857 : 101 chars : 'be used for cases with forcings from an input file and for the atex ', &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_tuner.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 654 : 102 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 6 check(s) failed.
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_wp2_wp3_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 928 : 101 chars : stats, & ! intent(inout)
line 1742 : 105 chars : stats_tmp_zt(i,k) = - gamma_over_implicit_ts * wp3_term_ta_lhs_result(5,i,k) * wp3(i,k-1) &
line 1743 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(4,i,k) * wp2(i,k) &
line 1744 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(3,i,k) * wp3(i,k) &
line 1745 : 105 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(2,i,k) * wp2(i,k+1) &
line 1746 : 103 chars : - gamma_over_implicit_ts * wp3_term_ta_lhs_result(1,i,k) * wp3(i,k+1)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sfc_varnce_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 723 : 105 chars : call sfc_varnce_check( sclr_dim, wp2(i,gr%k_lb_zm), up2(i,gr%k_lb_zm), & ! intent(in)
line 724 : 105 chars : vp2(i,gr%k_lb_zm), thlp2(i,gr%k_lb_zm), rtp2(i,gr%k_lb_zm), & ! intent(in)
line 725 : 105 chars : rtpthlp(i,gr%k_lb_zm), sclrp2(i,gr%k_lb_zm,:), & ! intent(in)
line 726 : 105 chars : sclrprtp(i,gr%k_lb_zm,:), sclrpthlp(i,gr%k_lb_zm,:), & ! intent(in)
line 727 : 108 chars : err_info ) ! intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/clubb_api_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 2430 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
line 2497 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/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/src/CLUBB_core/err_info_type_module.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 4: module err_info_type_module
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/err_info_type_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 265 : 102 chars : write(err_info%err_header_global,'(A38, I3, A15, I3, A, A28, A33, F6.2, A1, F6.2, A, A39, A40)') &
line 276 : 105 chars : write(err_info%err_header(i),'(A38, I3, A15, I3, A, A22, I2, A11, F6.2, A13, F6.2, A, A39, A40)') &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/mono_flux_limiter.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 700 : 104 chars : !$acc parallel loop gang vector collapse(2) default(present) reduction(.or.:l_any_adjustment_needed)
line 851 : 121 chars : ! MONOFLUX TEST COMMENT DO NOT REMOVE if( any(wpxp_net_adjust /= 0.0) ) write(fstderr,*) "MONOFLUX: wpxp adjusted "
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/parameters_tunable.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 839 : 106 chars : / sqrt( 4.0_core_rknd * ( 1.0_core_rknd - 0.4_core_rknd )**3 &
line 840 : 118 chars : + clubb_params(ngrdcol,iSkw_max_mag)**2 ) ) ) ! Known magic number
line 902 : 117 chars : mult_factor_zt(i) = 1.0_core_rknd + clubb_params(i,imult_coef) * log( avg_deltaz(i) / grid_spacing_thresh )
line 917 : 104 chars : nu_vert_res_dep%nu10(i) = clubb_params(i,inu10) * mult_factor_zt(i) !We're unsure of the grid
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_clubb_core_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1876 : 101 chars : stats, & ! intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/pdf_parameter_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 295 : 101 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 298 : 102 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 300 : 101 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_windm_edsclrm_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 493 : 102 chars : stats ) ! intent(inout)
line 501 : 102 chars : stats ) ! intent(inout)
line 953 : 103 chars : call windm_edsclrm_rhs( nzm, nzt, ngrdcol, gr, windm_edsclrm_scalar, dt, & ! intent(in)
line 954 : 103 chars : lhs_diff, edsclrm(:,:,edsclr), & ! intent(in)
line 955 : 103 chars : edsclrm_forcing(:,:,edsclr), & ! intent(in)
line 956 : 103 chars : rho_ds_zm, invrs_rho_ds_zt, & ! intent(in)
line 957 : 103 chars : l_imp_sfc_momentum_flux, wpedsclrp(:,gr%k_lb_zm,edsclr), & ! intent(in)
line 958 : 106 chars : stats, & ! intent(inout)
line 959 : 104 chars : rhs(:,:,edsclr) ) ! intent(out)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/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/src/CLUBB_core/stats_netcdf.F90 ---------------------
WARNING: Missing 'implicit none' after these blocks:
line 281: subroutine stats_init_api( registry_path, output_path, ncol, stats_tsamp, stats_tout, dt_main, &
line 489: subroutine stats_finalize_api( stats, err_info )
line 565: subroutine stats_type_initialize( ncol, defs, ndefs, zt, zm, stats, ierr, rad_zt, rad_zm )
line 716: subroutine stats_open_netcdf( output_path, time_units, &
line 951: subroutine stats_begin_timestep_api( itime, stats )
line 1001: subroutine stats_end_timestep_api( time_value, stats, err_info )
line 1052: subroutine stats_update_scalar( name, values, stats, icol, level )
line 1112: subroutine stats_update_1d( name, values, stats, icol )
line 1171: subroutine stats_update_2d( name, values, stats )
line 1229: subroutine stats_begin_budget_scalar( name, values, stats, icol )
line 1275: subroutine stats_begin_budget_1d( name, values, stats, icol )
line 1328: subroutine stats_begin_budget_2d( name, values, stats )
line 1369: subroutine stats_update_budget_scalar( name, values, stats, icol, level )
line 1419: subroutine stats_update_budget_1d( name, values, stats, icol )
line 1468: subroutine stats_update_budget_2d( name, values, stats )
line 1508: subroutine stats_finalize_budget_scalar( name, values, stats, icol, l_count_sample )
line 1574: subroutine stats_finalize_budget_1d( name, values, stats, icol, l_count_sample )
line 1647: subroutine stats_finalize_budget_2d( name, values, stats, l_count_sample )
line 1710: subroutine stats_avg_and_write( time_value, stats, ierr )
line 1836: subroutine stats_grid_init( ncol, zt_src, zm_src, zt_tgt, zm_tgt, &
line 1898: subroutine stats_update_grid( zt_src, zm_src, rho_vals, rho_levels, p_sfc, stats )
line 1976: subroutine stats_lh_samples_init( num_samples, nzt, nl_var_names, u_var_names, zt_vals, &
line 2095: subroutine stats_lh_samples_write_lognormal( samples, stats, err_info )
line 2148: subroutine stats_lh_samples_write_uniform( uniform_vals, mixture_comp, sample_weights, &
line 2364: subroutine stats_lookup_append( stats, id )
line 2406: subroutine stats_lookup_reject_append( stats, name )
line 2520: subroutine stats_read_registry_namelist( path, defs, ndefs, ierr )
line 2595: subroutine stats_expand_registry( sclr_dim, edsclr_dim, defs, ndefs, hydromet_list )
line 2896: subroutine add_expanded_def( base_def, name, out_defs, nout )
line 2946: subroutine parse_registry_line( line, def, l_valid )
line 3020: subroutine split_registry_fields( line, fields, nfields )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/stats_netcdf.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 235 : 103 chars : logical :: l_last_sample = .false. ! True if we want to output at the end of this sample
line 237 : 104 chars : logical :: l_different_output_grid = .false. ! True if we need to remap the variables before output
line 238 : 102 chars : logical :: l_output_rad_files = .false. ! True if we want to output radiation grid variables
line 526 : 101 chars : if ( allocated( stats%grid%rho_lin_spline_vals ) ) deallocate( stats%grid%rho_lin_spline_vals )
line 527 : 105 chars : if ( allocated( stats%grid%rho_lin_spline_levels ) ) deallocate( stats%grid%rho_lin_spline_levels )
line 803 : 104 chars : ret_code = nf90_def_var( stats%ncid, "time", NF90_DOUBLE, (/ stats%time_dimid /), stats%time_varid )
line 808 : 101 chars : ret_code = nf90_def_var( stats%ncid, "col", NF90_DOUBLE, (/ stats%col_dimid /), stats%col_varid )
line 828 : 121 chars : ret_code = nf90_def_var( stats%ncid, "lh_zt", NF90_DOUBLE, (/ stats%lh_2d%lh_zt_dimid /), stats%lh_2d%lh_zt_varid )
line 835 : 112 chars : ret_code = nf90_def_var( stats%ncid, "rad_zt", NF90_DOUBLE, (/ stats%rad_zt_dimid /), stats%rad_zt_varid )
line 839 : 112 chars : ret_code = nf90_def_var( stats%ncid, "rad_zm", NF90_DOUBLE, (/ stats%rad_zm_dimid /), stats%rad_zm_varid )
line 849 : 109 chars : ret_code = nf90_def_var( stats%ncid, "param", NF90_DOUBLE, (/ stats%param_dimid /), stats%param_varid )
line 853 : 101 chars : (/ stats%param_strlen_dimid, stats%param_dimid /), stats%param_name_varid )
line 912 : 110 chars : ret_code = nf90_put_att( stats%ncid, stats%vars(i)%varid, "long_name", trim( stats%vars(i)%long_name ) )
line 913 : 102 chars : ret_code = nf90_put_att( stats%ncid, stats%vars(i)%varid, "units", trim( stats%vars(i)%units ) )
line 920 : 108 chars : ret_code = nf90_put_var( stats%ncid, stats%col_varid, [( real( i, kind=core_rknd ), i=1, stats%ncol )] )
line 940 : 109 chars : ret_code = nf90_put_var( stats%ncid, stats%param_varid, [( real( i, kind=core_rknd ), i=1, nparams )] )
line 1796 : 101 chars : stats%grid%gr_source, stats%grid%gr_output, stats%vars(i)%nz, &
line 1797 : 107 chars : real( stats%vars(i)%buffer, kind=core_rknd ), stats%vars(i)%out_nz, &
line 1798 : 114 chars : size( stats%grid%rho_lin_spline_vals, 2 ), stats%grid%rho_lin_spline_vals, &
line 1954 : 101 chars : if ( allocated( stats%grid%rho_lin_spline_vals ) ) deallocate( stats%grid%rho_lin_spline_vals )
line 1955 : 105 chars : if ( allocated( stats%grid%rho_lin_spline_levels ) ) deallocate( stats%grid%rho_lin_spline_levels )
line 2397 : 116 chars : if ( stats%lookup%cache_len > 0 ) tmp(1:stats%lookup%cache_len) = stats%lookup%cache(1:stats%lookup%cache_len)
line 2442 : 105 chars : tmp(1:stats%lookup%reject_cache_len) = stats%lookup%reject_cache(1:stats%lookup%reject_cache_len)
line 2509 : 101 chars : write(date(32:33),'(i2.2)') nint( ( ( real( mod( nint( st_time ),3600 ),kind=time_precision ) / &
line 2586 : 101 chars : write( fstderr,* ) "Invalid stats registry entry (expected: name | grid | units | long_name)"
line 2687 : 101 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2688 : 101 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2695 : 101 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2696 : 101 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2703 : 101 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_1",out_defs,k)
line 2704 : 101 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_2",out_defs,k)
line 2711 : 113 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_1", &
line 2715 : 113 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_2", &
line 2721 : 101 chars : call add_expanded_def(base_def,"corr_w_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2722 : 101 chars : call add_expanded_def(base_def,"corr_w_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2729 : 103 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2730 : 103 chars : call add_expanded_def(base_def,"corr_chi_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2737 : 103 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2738 : 103 chars : call add_expanded_def(base_def,"corr_eta_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2745 : 103 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_1_n",out_defs,k)
line 2746 : 103 chars : call add_expanded_def(base_def,"corr_Ncn_"//trim(hydromet_list(j)( 1:2 ))//"_2_n",out_defs,k)
line 2753 : 115 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_1_n", &
line 2757 : 115 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"_2_n", &
line 2794 : 112 chars : trim( hydromet_list(j2)( 1:min( 2,len_trim( hydromet_list(j2) ) ) ) )//"p", &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/advance_xp2_xpyp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1320 : 101 chars : sclrprtp(:,:,sclr) ) ! Intent(inout)
line 1353 : 102 chars : sclrpthlp(:,:,sclr) ) ! Intent(inout)
line 3042 : 101 chars : stats_tmp(i,k) = (-gamma_over_implicit_ts * lhs_ta(2+gr%grid_dir_indx,i,k)) * xapxbp(i,km1) &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/setup_clubb_pdf_params.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 55 : 104 chars : subroutine setup_pdf_parameters_api( gr, nzm, nzt, ngrdcol, pdf_dim, hydromet_dim, dt, & ! Intent(in)
line 3823 : 104 chars : ! Correlation (in-precip) of chi (old s) and the precipitating hydrometeor in PDF component 1.
line 3827 : 104 chars : ! Correlation (in-precip) of chi (old s) and the precipitating hydrometeor in PDF component 2.
line 3841 : 104 chars : ! Correlation (in-precip) of eta (old t) and the precipitating hydrometeor in PDF component 1.
line 3845 : 104 chars : ! Correlation (in-precip) of eta (old t) and the precipitating hydrometeor in PDF component 2.
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/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/src/CLUBB_core/advance_xm_wpxp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 779 : 101 chars : call wpxp_term_pr1_lhs( nzm, ngrdcol, gr, C6rt_Skw_fnc, C6thl_Skw_fnc, C7_Skw_fnc, & ! Intent(in)
line 780 : 101 chars : invrs_tau_C6_zm, l_scalar_calc, & ! Intent(in)
line 781 : 102 chars : lhs_pr1_wprtp, lhs_pr1_wpthlp, lhs_pr1_wpsclrp ) ! Intent(out)
line 854 : 102 chars : stats, & ! InOut
line 855 : 101 chars : rtm, wprtp, thlm, wpthlp, sclrm, wpsclrp, & ! InOut
line 856 : 101 chars : err_info ) ! InOut
line 896 : 101 chars : stats, & ! InOut
line 4881 : 108 chars : tmp_zm_stats(i,k) = (-gamma_over_implicit_ts*lhs_ta_wpxp(2+gr%grid_dir_indx,i,k)) * wpxp(i,km1) &
line 4883 : 105 chars : + (-gamma_over_implicit_ts*lhs_ta_wpxp(2-gr%grid_dir_indx,i,k)) * wpxp(i,kp1)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/stats_clubb_utilities.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 129 : 101 chars : wp2_zt, thlp2_zt, wpthlp_zt, wprtp_zt, rtp2_zt, rtpthlp_zt, up2_zt, vp2_zt, upwp_zt, vpwp_zt, &
line 130 : 101 chars : wpup2, wpvp2, a3_coef_zt, wp3_on_wp2_zt, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/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: 20 check(s) failed.
CLUBBStandardsCheck.py has begun.
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/SILHS/mt95.F90 ---------------------
WARNING: Missing 'implicit none' after these blocks:
line 362: subroutine init_by_type( put, get )
line 382: subroutine init_by_scalar( put )
line 413: subroutine init_by_array( put )
line 461: subroutine next_state( )
line 593: subroutine genrand_int32_0d( y )
line 614: subroutine genrand_int32_1d( y )
line 629: subroutine genrand_int32_2d( y )
line 644: subroutine genrand_int32_3d( y )
line 659: subroutine genrand_int32_4d( y )
line 674: subroutine genrand_int32_5d( y )
line 689: subroutine genrand_int32_6d( y )
line 704: subroutine genrand_int32_7d( y )
line 720: subroutine genrand_int31_0d( y )
line 732: subroutine genrand_int31_1d( y )
line 747: subroutine genrand_int31_2d( y )
line 762: subroutine genrand_int31_3d( y )
line 777: subroutine genrand_int31_4d( y )
line 792: subroutine genrand_int31_5d( y )
line 807: subroutine genrand_int31_6d( y )
line 822: subroutine genrand_int31_7d( y )
line 838: subroutine genrand_real1_0d( r )
line 853: subroutine genrand_real1_1d( r )
line 868: subroutine genrand_real1_2d( r )
line 883: subroutine genrand_real1_3d( r )
line 898: subroutine genrand_real1_4d( r )
line 913: subroutine genrand_real1_5d( r )
line 928: subroutine genrand_real1_6d( r )
line 943: subroutine genrand_real1_7d( r )
line 959: subroutine genrand_real2_0d( r )
line 974: subroutine genrand_real2_1d( r )
line 989: subroutine genrand_real2_2d( r )
line 1004: subroutine genrand_real2_3d( r )
line 1019: subroutine genrand_real2_4d( r )
line 1034: subroutine genrand_real2_5d( r )
line 1049: subroutine genrand_real2_6d( r )
line 1064: subroutine genrand_real2_7d( r )
line 1080: subroutine genrand_real3_0d( r )
line 1095: subroutine genrand_real3_1d( r )
line 1110: subroutine genrand_real3_2d( r )
line 1125: subroutine genrand_real3_3d( r )
line 1140: subroutine genrand_real3_4d( r )
line 1155: subroutine genrand_real3_5d( r )
line 1170: subroutine genrand_real3_6d( r )
line 1185: subroutine genrand_real3_7d( r )
line 1201: subroutine genrand_res53_0d( r )
line 1218: subroutine genrand_res53_1d( r )
line 1233: subroutine genrand_res53_2d( r )
line 1248: subroutine genrand_res53_3d( r )
line 1263: subroutine genrand_res53_4d( r )
line 1278: subroutine genrand_res53_5d( r )
line 1293: subroutine genrand_res53_6d( r )
line 1308: subroutine genrand_res53_7d( r )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/SILHS/latin_hypercube_driver_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 113 : 101 chars : ! The uniform distribution corresponds to all the same variables as X_nl,
line 410 : 105 chars : !$acc update host(X_u_all_levs,l_in_precip,lh_sample_point_weights,X_mixt_comp_all_levs,k_lh_start)
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 2 check(s) failed.
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/wangara.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 104 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 105 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/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/src/Benchmark_cases/prescribe_forcings.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 1: module prescribe_forcings_module
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/prescribe_forcings.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 297 : 104 chars : !$acc enter data create( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 471 : 118 chars : !$acc update device( wm_zt, wm_zm, thlm_forcing, rtm_forcing, um_ref, vm_ref, sclrm_forcing, edsclrm_forcing )
line 911 : 103 chars : !$acc exit data delete( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 1114 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1118 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1122 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1123 : 112 chars : exner_zm(i,km1), exner_zm(i,k00), exner_zm(i,kp1), exner_zm(i,kp2) )
line 1126 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1127 : 108 chars : thlm_zm(i,km1), thlm_zm(i,k00), thlm_zm(i,kp1), thlm_zm(i,kp2) )
line 1130 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1131 : 104 chars : rtm_zm(i,km1), rtm_zm(i,k00), rtm_zm(i,kp1), rtm_zm(i,kp2) )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/dycoms2_rf02.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 121 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 122 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/mpace_b.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 157 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 158 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/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/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/src/Benchmark_cases/dycoms2_rf01.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 91 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 92 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/gabls2.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 152 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 153 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/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/src/KK_microphys/
CLUBBStandardsCheck.py has begun.
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 19: module fill_holes_tests
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 162 : 102 chars : field_initial(i,k) = 10_core_rknd * ( sin( 8.0 * pi * ( k + i ) / nz ) + one ) &
line 180 : 110 chars : if ( l_print_name ) print *, "filling holes of: easy_fill_test_"//trim(fill_type_names(fill_type))
line 195 : 121 chars : write(*,'(A6,I5,1X,I5,A4,E30.20,A3,E30.20)') "field(", i, k, ") = ", field(i,k), " < ", threshold
line 241 : 118 chars : if ( l_print_name ) print *, "filling holes of: reverse_easy_fill_test_"//trim(fill_type_names(fill_type))
line 344 : 102 chars : field_initial(i,k) = 10_core_rknd * ( sin( 8.0 * pi * ( k + i ) / nz ) + one ) &
line 363 : 113 chars : if ( l_print_name ) print *, "filling holes of: below_thresh_test_"//trim(fill_type_names(fill_type))
line 385 : 103 chars : print *, "-- error", two * abs( initial_mass - new_mass ) / ( initial_mass + new_mass )
line 419 : 121 chars : if ( l_print_name ) print *, "filling holes of: reverse_below_thresh_test_"//trim(fill_type_names(fill_type))
line 430 : 116 chars : print *, "ERROR in below_thresh_test REVERSE MODE: method was not conservative in grid reverse mode"
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/rev_direction_grid_test.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 1: module rev_direction_grid_test
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/rev_direction_grid_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 233 : 113 chars : !seed used = -1102123834 1480475276 -248259584 525496705 1614213535 -1267052688 6290033 -78466652
line 234 : 127 chars : !seed_output = (/ -1102123834, 1480475276, -248259584, 525496705, 1614213535, -1267052688, 6290033, -78466652 /)
line 750 : 101 chars : call xpyp_term_ta_pdf_lhs_godunov( gr_ascending%nzm, gr_ascending%nzt, 1, gr_ascending, & ! In
line 754 : 104 chars : call xpyp_term_ta_pdf_lhs_godunov( gr_descending%nzm, gr_descending%nzt, 1, gr_descending, & ! In
line 795 : 101 chars : call xpyp_term_ta_pdf_rhs_godunov( gr_ascending%nzm, gr_ascending%nzt, 1, gr_ascending, & ! In
line 801 : 104 chars : call xpyp_term_ta_pdf_rhs_godunov( gr_descending%nzm, gr_descending%nzt, 1, gr_descending, & ! In
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/pdf_parameter_tests.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 558 : 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/src/G_unit_test_types/spurious_source_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 502 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 6 check(s) failed.
[91m===============================================================
[91mCLUBBStandardsCheck FAILED
[91m THIS IS PRINTED IN ALL RED, CAPITAL LETTERS, AND USES
[91m AN EXCLAMATION MARK TO ENSURE THE DEVELOPERS FEEL SHAME!
[91m IF YOU ARE ONE OF THESE "DEVELOPERS" CHECK THE
[91m LOG FILE FOR DETAILS: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gcc_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
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: Hole filling failed! A hole could not be filled.
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
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: One-level hole filling was not able to fill holes completely! The holes were filled partially. |total_hole| > total_mass
Warning: Hole filling failed! A hole could not be filled.
Warning: One-level hole filling was not successful! total_mass ~= 0
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 3 23 -300 3 12 8 557
hmm = 2.4553947299393144E-005
mixt_frac = 0.12500380666669780
precip_frac = 0.74392849564487618
precip_frac_1 = 0.26358714013352053
precip_frac_2 = 0.81255107754391431
hmp2_ip_on_hmm2_ip = 0.44458680667019534
hmm_ip = 3.3005789458445861E-005
hmp2 = 5.6782920046720086E-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.
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 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
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 = 9.0000000000000012E-011
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 = 4.0000000000000002E-004
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 3 23 -300 3 12 8 558
wm = -2.8921054012137182E-002
rtm = 3.3470164625856872E-003
thlm = 304.14691157958941
wp2 = 0.25614862627471469
rtp2 = 1.1193909100565365E-007
thlp2 = 4.9420642790630742E-002
Skw = 3.4175851112978481
wp3 = 0.44305473843142329
Skrt = -3.8700589941932488
rtp3 = -1.4494066488761685E-010
Skthl = -3.1812255166904597
thlp3 = -3.4950793644405273E-002
wprtp = -1.6048311941582372E-004
wpthlp = 4.1791058220055938E-002
rtpthlp = -6.6940290431306289E-005
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).
Running ADG1 expected-failure tests.
ADG1_w_closure correctly failed when given wp2 = 0.
ADG1_pdf_driver correctly failed when given thlp2 = 0.
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 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
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 = 9.0000000000000012E-011
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 = 4.0000000000000002E-004
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 3 23 -300 3 12 8 558
wm = -2.8921054012137182E-002
rtm = 3.3470164625856872E-003
thlm = 304.14691157958941
wp2 = 0.25614862627471469
rtp2 = 1.1193909100565365E-007
thlp2 = 4.9420642790630742E-002
Skw = 3.4175851112978481
wp3 = 0.44305473843142329
Skrt = -3.8700589941932488
rtp3 = -1.4494066488761685E-010
Skthl = -3.1812255166904597
thlp3 = -3.4950793644405273E-002
wprtp = -1.6048311941582372E-004
wpthlp = 4.1791058220055938E-002
rtpthlp = -6.6940290431306289E-005
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 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
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 = 9.0000000000000012E-011
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
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 3 23 -300 3 12 8 559
wm = -2.8921054012137182E-002
rtm = 3.3470164625856807E-003
thlm = 304.14691157958941
wp2 = 0.25614862631564161
rtp2 = 8.1170452949327002E-008
thlp2 = 4.9436673303228762E-002
Skw = 3.3187564375671306
wp3 = 0.43024261805667868
Skrt = 2.4581142859373646
rtp3 = 5.6845884060852122E-011
Skthl = 3.2457524846854326
thlp3 = 3.5677076199781298E-002
wprtp = 6.5322512898017713E-005
wpthlp = 1.7953854129037251E-002
rtpthlp = 5.7011978593170175E-005
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 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
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 = 9.0000000000000012E-011
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 = 4.0000000000000002E-004
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 3 23 -300 3 12 8 561
wm = -2.8921054012137182E-002
rtm = 3.3470164625858424E-003
thlm = 304.14691157958924
wp2 = 0.25614862557894158
rtp2 = 2.0818156655238908E-006
thlp2 = 4.9339235888956057E-002
Skw = 2.3287316214456606
wp3 = 0.30189608791365558
Skrt = 3.3006090025414156
rtp3 = 9.9141957207821914E-009
Skthl = -4.4403328407536300
thlp3 = -4.8663591115292747E-002
wprtp = 2.0177599500040985E-004
wpthlp = -1.5889255424072005E-003
rtpthlp = -2.8844290128355026E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new hybrid PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 4.0000000000000002E-004
rtp2 = 1.0000000000000001E-016
thlp2 = 1.0000000000000000E-004
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 = 9.0000000000000012E-011
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 = 4.0000000000000002E-004
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 3 23 -300 3 12 8 562
wm = -2.8921054012137182E-002
rtm = 3.3470164625858624E-003
thlm = 304.14691157958924
wp2 = 0.25614862545616068
rtp2 = 2.0510470274675863E-006
thlp2 = 4.9355517858649556E-002
Skw = 2.2297098301392522
wp3 = 0.28905893157400042
Skrt = 2.3160891621495754
rtp3 = 6.8032855806279320E-009
Skthl = 3.2786365937524220
thlp3 = 3.5949830862756521E-002
wprtp = -2.3409115825821059E-004
wpthlp = -2.4175903677867684E-002
rtpthlp = 2.8635065118909071E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing spurious source unit test
====================================
Test 1
Vertical integral of <rt> before advance_xm_wpxp = 28.316670259590918
Vertical integral of <rt> after advance_xm_wpxp = 28.454761250379875
Vertical integral of <rt> forcing = -3.9696697370157397E-005
Flux of <rt> at the surface = 5.0000000000000001E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 1.3037531124138191E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 793918.42642888101
Vertical integral of <thl> after advance_xm_wpxp = 793908.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 = 25.419730600395713
Vertical integral of <rt> after advance_xm_wpxp = 25.647863423019910
Vertical integral of <rt> forcing = -2.0168479036889614E-005
Flux of <rt> at the surface = 7.8061122111760758E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -6.2080738770162180E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 789943.66256586206
Vertical integral of <thl> after advance_xm_wpxp = 789940.03819918819
Vertical integral of <thl> forcing = -2.7693446670342244E-002
Flux of <thl> at the surface = 1.5612224422352152E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.7442991495641991E-012
The spurious source of <thl> is within acceptable limits.
Test 3
Vertical integral of <rt> before advance_xm_wpxp = 37.824094928450137
Vertical integral of <rt> after advance_xm_wpxp = 37.849428403625879
Vertical integral of <rt> forcing = -7.0035013892268117E-005
Flux of <rt> at the surface = 1.5447993114477078E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -2.8880435369582624E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 797774.97980713518
Vertical integral of <thl> after advance_xm_wpxp = 797751.99356319103
Vertical integral of <thl> forcing = -7.9710411770085268E-002
Flux of <thl> at the surface = 3.0895986228954152E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 2.4966140266258208E-014
The spurious source of <thl> is within acceptable limits.
Test 4
Vertical integral of <rt> before advance_xm_wpxp = 43.790599941911118
Vertical integral of <rt> after advance_xm_wpxp = 44.029887449185836
Vertical integral of <rt> forcing = -6.0925650285205875E-005
Flux of <rt> at the surface = 8.5855067453415642E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 1.0734279133964297E-016
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 801963.59585088491
Vertical integral of <thl> after advance_xm_wpxp = 801949.02522760606
Vertical integral of <thl> forcing = -6.5739757754635156E-002
Flux of <thl> at the surface = 1.7171013490683128E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.1163986401996340E-012
The spurious source of <thl> is within acceptable limits.
Test 5
Vertical integral of <rt> before advance_xm_wpxp = 13.493058120739910
Vertical integral of <rt> after advance_xm_wpxp = 13.503540960263788
Vertical integral of <rt> forcing = -2.5703997402426697E-005
Flux of <rt> at the surface = 6.0646795815372133E-005
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -2.0003530082357557E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 787962.99554359692
Vertical integral of <thl> after advance_xm_wpxp = 787955.00173000421
Vertical integral of <thl> forcing = -2.7858981225715566E-002
Flux of <thl> at the surface = 1.2129359163074427E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 3.8232958465833633E-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
Fatal error in CLUBB - MPI Process 1 / Chunk index 1
Column index not available. Latitude/Longitude for column 1: ******/******
Latitude range from -90 to 90 degrees, Longitude range from -180 to 180 degrees
Error in ADG1_w_closure: wp2 must not be less than the w tolerance value.
Fatal error in CLUBB - MPI Process 1 / Chunk index 1
Column index not available. Latitude/Longitude for column 1: ******/******
Latitude range from -90 to 90 degrees, Longitude range from -180 to 180 degrees
Error in ADG1_ADG2_responder_params: xp2 must not be less than its tolerance value.
Input_var1: -1000000.0000000000 -10000000.000000000 -100000000.00000000 -1000000000.0000000
Input_var2: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
True outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
Smooth max:
Expected outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
True outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
Input_var1: -1000000.0000000000
Input_var2: 1.0000000000000000
Smooth min:
Expected outcome: -1000000.0000002501
True outcome: -1000000.0000002501
Smooth max:
Expected outcome: 1.0000002500019036
True outcome: 1.0000002500019036
Success!
==================================================
====================== Running fill holes tests ======================
-------------- easy_fill_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: easy_fill_test
-----------------------------------------------
-------------- below_thresh_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: below_thresh_test
-----------------------------------------------
====================== End fill holes tests ======================
Running /home/jenkins/workspace/clubb_G_unit_gfortran_tests/run_scripts/../install/latest/G_unit_tests ...
All tests have succeeded.
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Declarative: Post Actions)
[Pipeline] script
[Pipeline] {
[Pipeline] cleanWs
[WS-CLEANUP] Deleting project workspace...
[WS-CLEANUP] Deferred wipeout is used...
[WS-CLEANUP] done
[Pipeline] }
[Pipeline] // script
[Pipeline] }
[Pipeline] // stage
[Pipeline] }
[Pipeline] // withEnv
[Pipeline] }
[Pipeline] // node
[Pipeline] End of Pipeline
Finished: SUCCESS
