Skipping 319 KB..
Full LogTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTest
[ 97%] [32mBuilding Fortran object src/CMakeFiles/clubb_tuner.dir/clubb_tuner.F90.o[0m
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_tuner.F90:294:6:
294 | use nrtype, only: &
| 1
Warning: Unused parameter 'sp' which has been explicitly imported at (1) [-Wunused-parameter]
[ 97%] [32m[1mLinking Fortran executable clubb_tuner[0m
[ 97%] Built target clubb_tuner
[35m[1mScanning dependencies of target G_unit_test_lib[0m
[ 97%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/KK_integrals_tests.F90.o[0m
[ 97%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/Nc_Ncn_test.F90.o[0m
[ 97%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/corr_cholesky_mtx_tests.F90.o[0m
[ 97%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/fill_holes_tests.F90.o[0m
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90:393:27:
393 | .or. any( field(:,upper_hf_level:nz) /= field(:,upper_hf_level:nz) ) ) then
| 1
Warning: Inequality comparison for REAL(8) at (1) [-Wcompare-reals]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90:392:22:
392 | if ( any( field(:,1:lower_hf_level) /= field(:,1:lower_hf_level) ) &
| 1
Warning: Inequality comparison for REAL(8) at (1) [-Wcompare-reals]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90:444:22:
444 | if ( any( field_rev_rev /= field ) ) then
| 1
Warning: Inequality comparison for REAL(8) at (1) [-Wcompare-reals]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90:215:27:
215 | .or. any( field(:,upper_hf_level:nz) /= field(:,upper_hf_level:nz) ) ) then
| 1
Warning: Inequality comparison for REAL(8) at (1) [-Wcompare-reals]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90:214:22:
214 | if ( any( field(:,1:lower_hf_level) /= field(:,1:lower_hf_level) ) &
| 1
Warning: Inequality comparison for REAL(8) at (1) [-Wcompare-reals]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90:266:22:
266 | if ( any( field_rev_rev /= field ) ) then
| 1
Warning: Inequality comparison for REAL(8) at (1) [-Wcompare-reals]
[ 97%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/hole_filling_tests.F90.o[0m
[ 97%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/mu_sigma_hm_tests.F90.o[0m
[ 97%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/pdf_parameter_tests.F90.o[0m
[ 98%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/read_corr_mtx_test.F90.o[0m
[ 98%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/rev_direction_grid_test.F90.o[0m
[ 98%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/silhs_category_test.F90.o[0m
[ 98%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/smooth_heaviside_tests.F90.o[0m
[ 99%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/smooth_min_max_tests.F90.o[0m
[ 99%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/spurious_source_test.F90.o[0m
[100%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/tuner_tests.F90.o[0m
[100%] [32mBuilding Fortran object src/G_unit_test_types/CMakeFiles/G_unit_test_lib.dir/w_up_in_cloud_tests.F90.o[0m
[100%] [32m[1mLinking Fortran static library libG_unit_test_lib.a[0m
[100%] Built target G_unit_test_lib
[35m[1mScanning dependencies of target G_unit_tests[0m
[100%] [32mBuilding Fortran object src/CMakeFiles/G_unit_tests.dir/G_unit_tests.F90.o[0m
[100%] [32m[1mLinking Fortran executable G_unit_tests[0m
[100%] Built target G_unit_tests
[35m[1mScanning dependencies of target clubb_driver_test[0m
[100%] [32mBuilding Fortran object src/CMakeFiles/clubb_driver_test.dir/clubb_driver_test.F90.o[0m
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90:74:26:
74 | character(len=10) :: arg
| 1
Warning: Unused variable 'arg' declared at (1) [-Wunused-variable]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90:17:6:
17 | use error_code, only: &
| 1
Warning: Unused parameter 'clubb_generalized_grd_test_err' which has been explicitly imported at (1) [-Wunused-parameter]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90:34:6:
34 | use constants_clubb, only: &
| 1
Warning: Unused parameter 'fstdout' which has been explicitly imported at (1) [-Wunused-parameter]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90:21:6:
21 | use model_flags, only: &
| 1
Warning: Unused parameter 'l_test_grid_generalization' which has been explicitly imported at (1) [-Wunused-parameter]
[100%] [32m[1mLinking Fortran executable clubb_driver_test[0m
[100%] Built target clubb_driver_test
[36mInstall the project...[0m
-- Install configuration: "Debug"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/libclubb_driver_lib.a
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_tuner
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_tuner" to "/usr/local/spack/opt/spack/linux-pop22-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/gfortran_DEBUG_PRECdouble/G_unit_tests
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/G_unit_tests" to "/usr/local/spack/opt/spack/linux-pop22-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/gfortran_DEBUG_PRECdouble/clubb_thread_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_thread_test" to "/usr/local/spack/opt/spack/linux-pop22-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/gfortran_DEBUG_PRECdouble/clubb_standalone
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_standalone" to "/usr/local/spack/opt/spack/linux-pop22-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/gfortran_DEBUG_PRECdouble/clubb_driver_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_driver_test" to "/usr/local/spack/opt/spack/linux-pop22-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"
Setting CLUBB installation dir: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gfortran_DEBUG_PRECdouble
Setting CLUBB installation dir: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble
Using inferred toolchain file: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/cmake/toolchains/linux_x86_64_gfortran.cmake
about to cmnake /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gfortran_DEBUG_PRECdouble
Running CMake configure...
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:
rtm_sponge_damp_settings
wp2_sponge_damp_settings
rtm_sponge_damp_profile
wp2_sponge_damp_profile
wp3_sponge_damp_settings
thlm_sponge_damp_settings
uv_sponge_damp_profile
up2_vp2_sponge_damp_settings
thlm_sponge_damp_profile
wp3_sponge_damp_profile
uv_sponge_damp_settings
up2_vp2_sponge_damp_profile
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/gfortran_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/soil_vegetation.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 174 : 103 chars : call stat_update_var_pt( stats_metadata%isfc_soil_T_in_K, 1, sfc_soil_T_in_K(i), stats_sfc(i) )
line 175 : 105 chars : call stat_update_var_pt( stats_metadata%ideep_soil_T_in_K, 1, deep_soil_T_in_K(i), stats_sfc(i) )
line 212 : 101 chars : call stat_update_var_pt( stats_metadata%isoil_heat_flux, 1, soil_heat_flux(i), & ! intent(in)
line 213 : 104 chars : stats_sfc(i) ) ! intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 111 : 107 chars : write(fstderr, *) "This section ensures that everything allocated in init_clubb_case will be deallocated"
line 112 : 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 908 : 106 chars : pdf_params_zm_flip%ice_supersat_frac_1(i,:) = flip( pdf_params_zm%ice_supersat_frac_1(i,:), nzm )
line 909 : 106 chars : pdf_params_zm_flip%ice_supersat_frac_2(i,:) = flip( pdf_params_zm%ice_supersat_frac_2(i,:), nzm )
line 911 : 117 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1080 : 101 chars : rtm_ref_flip, thlm_ref_flip, um_ref_flip, vm_ref_flip, ug_flip, vg_flip, & ! Intent(in)
line 1094 : 111 chars : um_flip, vm_flip, upwp_flip, vpwp_flip, up2_flip, vp2_flip, up3_flip, vp3_flip, & ! Intent(inout)
line 1096 : 111 chars : wp2_flip, wp3_flip, rtp2_flip, rtp3_flip, thlp2_flip, thlp3_flip, rtpthlp_flip, & ! Intent(inout)
line 1102 : 109 chars : wp2rtp_flip, wp2thlp_flip, uprcp_flip, vprcp_flip, rc_coef_zm_flip, wp4_flip, & ! intent(inout)
line 1103 : 105 chars : wpup2_flip, wpvp2_flip, wp2up2_flip, wp2vp2_flip, ice_supersat_frac_flip, & ! intent(inout)
line 1587 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1614 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2up_implicit, nzt, ngrdcol, &
line 1618 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2up_explicit, nzt, ngrdcol, &
line 1622 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2vp_implicit, nzt, ngrdcol, &
line 1626 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2vp_explicit, nzt, ngrdcol, &
line 1660 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpup2_implicit, nzt, ngrdcol, &
line 1664 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpup2_explicit, nzt, ngrdcol, &
line 1668 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpvp2_implicit, nzt, ngrdcol, &
line 1672 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpvp2_explicit, nzt, ngrdcol, &
----------------------------------------------------------------------------------------------------
--------------------- 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.
3115 : 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 198 : 110 chars : l_modify_ic_with_cubic_int, & ! Flag for interpolating the sounding profile with Steffen's monotone cubic
line 199 : 110 chars : ! method to obtain smoother initial condition profile, which is found to be
line 200 : 111 chars : ! beneficial to achive a better numerical solution convergence. If this flag
line 201 : 118 chars : ! is turned off, the initial conditions will be generated with linear interpolation.
line 202 : 111 chars : ! This is done on a case-by-case basis, since using the monotone cubic method
line 203 : 109 chars : ! requires a special sounding.in file with many additional sounding levels.
line 204 : 109 chars : l_modify_bc_for_cnvg_test, & ! Flag to activate modifications on boundary condition for convergence test
line 736 : 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 752 : 102 chars : !$omp rtp2_mc, thlp2_mc, rtpthlp_mc, Ncm_mc, Ncm, Nccnm, em, tau_zm, tau_zt, Kh_zt, Kh_zm, rfrzm, &
line 753 : 105 chars : !$omp rrm, Nrm, Nc_in_cloud, wpNcp, radht, Frad, Frad_SW_up, Frad_LW_up, Frad_SW_down, Frad_LW_down, &
line 754 : 102 chars : !$omp wpsclrp_sfc, wpedsclrp_sfc, sclrm, wpsclrp, sclrp2, sclrp3, sclrprtp, sclrpthlp, sclrpthvp, &
line 799 : 101 chars : !----------------------------------- Optional Input Variables -----------------------------------
line 1017 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 1691 : 102 chars : open(unit=iunit_grid_adaptation, file=fname_grid_adaptation, status='replace', action='write')
line 2004 : 105 chars : allocate( um_ref(ngrdcol, gr%nzt) ) ! Reference u wind for nudging; Michael Falk, 17 Oct 2007
line 2005 : 105 chars : allocate( vm_ref(ngrdcol, gr%nzt) ) ! Reference v wind for nudging; Michael Falk, 17 Oct 2007
line 2495 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 2526 : 106 chars : ! Similar call to above, checking parameter values, but because we've placed this in "init_clubb_case"
line 2534 : 116 chars : write(fstderr, *) "Fatal error calling check_clubb_settings_api and/or check_parameters_api in clubb_driver"
line 2577 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 2589 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 2590 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 2593 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 2595 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 2601 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 2602 : 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 2603 : 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 2604 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 2605 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 2606 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 2607 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 2608 : 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 2609 : 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 2610 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 2725 : 103 chars : invrs_rho_ds_zt(i,k) = invrs_rho_ds_zt_init(k) ! inv. dry, static density: t-levs
line 2731 : 103 chars : um_pert(i,k) = zero ! Variables used to track perturbed version of winds.
line 2733 : 105 chars : wm_zt(i,k) = wm_zt_init(k) ! Imposed large scale w - Thermodynamic levels
line 2926 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 3136 : 101 chars : !----------------------------------- Optional Input Variables -----------------------------------
line 3292 : 103 chars : um(i,:), vm(i,:), rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), &
line 3630 : 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 3632 : 104 chars : !$acc lh_rt_clipped, lh_thl_clipped, lh_rc_clipped, lh_rv_clipped, lh_Nc_clipped, &
line 3636 : 116 chars : !$acc pdf_params%cloud_frac_1, pdf_params%cloud_frac_2, pdf_params%cthl_1, pdf_params%cthl_2, &
line 3637 : 112 chars : !$acc pdf_params%crt_1, pdf_params%crt_2, pdf_params%stdev_chi_1, pdf_params%stdev_chi_2, &
line 3653 : 102 chars : wp2_zt = zm2zt_api( gr%nzm, gr%nzt, ngrdcol, gr, wp2, w_tol_sqd ) ! Positive definite quantity
line 3659 : 109 chars : thlm(i,:), p_in_Pa(i,:), exner(i,:), rho(i,:), rho_zm(i,:), rtm(i,:), &! In
line 3660 : 109 chars : rcm(i,:), cloud_frac(i,:), wm_zt(i,:), wm_zm(i,:), wp2_zt(i,:), & ! In
line 3661 : 106 chars : hydromet(i,:,:), Nc_in_cloud(i,:), & ! In
line 3677 : 103 chars : stats_zt(i), stats_zm(i), & ! Inout
line 3678 : 103 chars : stats_sfc(i), stats_lh_zt(i), & ! Inout
line 3679 : 102 chars : Nccnm(i,:), & ! Inout
line 3680 : 117 chars : hydromet_mc(i,:,:), Ncm_mc(i,:), rcm_mc(i,:), rvm_mc(i,:), & ! Out
line 3681 : 107 chars : thlm_mc(i,:), hydromet_vel_zt(i,:,:), & ! Out
line 3682 : 102 chars : hydromet_vel_covar_zt_impc(i,:,:), & ! Out
line 3683 : 102 chars : hydromet_vel_covar_zt_expc(i,:,:), & ! Out
line 3684 : 110 chars : wprtp_mc(i,:), wpthlp_mc(i,:), rtp2_mc(i,:), & ! Out
line 3685 : 105 chars : thlp2_mc(i,:), rtpthlp_mc(i,:) ) ! Out
line 3710 : 103 chars : stats_zt(i), stats_zm(i), stats_sfc(i), & ! Inout
line 3711 : 103 chars : hydromet(i,:,:), hydromet_vel_zt(i,:,:), hydrometp2(i,:,:), & ! Inout
line 3712 : 103 chars : K_hm(i,:,:), Ncm(i,:), Nc_in_cloud(i,:), rvm_mc(i,:), & ! Inout
line 3713 : 103 chars : thlm_mc(i,:), err_info, & ! Inout
line 3714 : 101 chars : wphydrometp(i,:,:), wpNcp(i,:) ) ! Out
line 3723 : 106 chars : !$acc update device( rcm_mc, rvm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3781 : 104 chars : !$acc update host( rvm_mc, rcm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3787 : 105 chars : call stat_update_var( stats_metadata%irtm_mc, rvm_mc(i,:)+rcm_mc(i,:), stats_zt(i) ) ! kg/kg/s
line 3789 : 106 chars : call stat_update_var( stats_metadata%iwprtp_mc, wprtp_mc(i,:), stats_zm(i) ) ! m*(kg/kg)/s^2
line 3791 : 104 chars : call stat_update_var( stats_metadata%irtp2_mc, rtp2_mc(i,:), stats_zm(i) ) ! (kg/kg)^2/s
line 3793 : 104 chars : call stat_update_var( stats_metadata%irtpthlp_mc, rtpthlp_mc(i,:), stats_zm(i) ) ! K*(kg/kg)/s
line 3810 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, X_nl_all_levs, &
line 3836 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, thlm, rtm, rcm )
line 3882 : 106 chars : call update_radiation_variables( gr%nzm, gr%nzt, radht(i,:), Frad_SW_up(i,:), Frad_LW_up(i,:), &
line 3981 : 104 chars : clubb_config_flags%grid_remap_method, & ! intent(in)
line 4049 : 101 chars : .and. (stats_metadata%l_stats_last) ) then ! only adapt grid when the average of the last
line 4154 : 103 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_scheme = ', time_microphys_scheme
line 4155 : 104 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_advance = ', time_microphys_advance
line 4157 : 103 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_output_multi_col = ', time_output_multi_col
line 4284 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 4296 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 4297 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 4300 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 4302 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 4308 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 4309 : 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 4310 : 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 4311 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 4312 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 4313 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 4314 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 4315 : 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 4316 : 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 4317 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 4346 : 119 chars : if ( thlm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( thlm_sponge_damp_profile )
line 4347 : 118 chars : if ( rtm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( rtm_sponge_damp_profile )
line 4348 : 117 chars : if ( uv_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( uv_sponge_damp_profile )
line 4349 : 118 chars : if ( wp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp2_sponge_damp_profile )
line 4350 : 118 chars : if ( wp3_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp3_sponge_damp_profile )
line 4351 : 122 chars : if ( up2_vp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( up2_vp2_sponge_damp_profile )
line 5027 : 104 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 5035 : 106 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 5049 : 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/sounding.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 377 : 102 chars : um(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 379 : 101 chars : vm(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 381 : 101 chars : ugm(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 383 : 101 chars : vgm(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 385 : 101 chars : thlm(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 387 : 101 chars : rtm(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 389 : 101 chars : press(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 392 : 101 chars : wm(i,k) = mono_cubic_interp( gr%zt(i,k), km1, k00, kp1, kp2, z(km1), z(k00), z(kp1), &
line 409 : 108 chars : edsclr_snd(km1,edsclr), edsclr_snd(k00,edsclr), &
line 410 : 107 chars : edsclr_snd(kp1,edsclr), edsclr_snd(kp2,edsclr) )
line 434 : 118 chars : um(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), u(k_above), u(k_above-1) )
line 435 : 118 chars : vm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), v(k_above), v(k_above-1) )
line 436 : 120 chars : ugm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), ug(k_above), ug(k_above-1) )
line 437 : 120 chars : vgm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), vg(k_above), vg(k_above-1) )
line 438 : 126 chars : thlm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), theta(k_above), theta(k_above-1) )
line 439 : 120 chars : rtm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), rt(k_above), rt(k_above-1) )
line 440 : 112 chars : press(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), p_in_Pa(k_above), &
line 442 : 122 chars : wm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), subs(k_above), subs(k_above-1) )
line 446 : 105 chars : sclrm(i,k,sclr) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), &
line 452 : 109 chars : edsclrm(i,k,edsclr) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), &
line 481 : 124 chars : wm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), subs(k_above), subs(k_above-1) )
line 507 : 124 chars : wm(i,k) = lin_interpolate_two_points( gr%zt(i,k), z(k_above), z(k_above-1), subs(k_above), subs(k_above-1) )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_tuner.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 650 : 102 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 8 check(s) failed.
CLUBBStandardsCheck.py has begun.
--------------------- 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 776 : 105 chars : call sfc_varnce_check( sclr_dim, wp2(i,gr%k_lb_zm), up2(i,gr%k_lb_zm), & ! intent(in)
line 777 : 105 chars : vp2(i,gr%k_lb_zm), thlp2(i,gr%k_lb_zm), rtp2(i,gr%k_lb_zm), & ! intent(in)
line 778 : 105 chars : rtpthlp(i,gr%k_lb_zm), sclrp2(i,gr%k_lb_zm,:), & ! intent(in)
line 779 : 105 chars : sclrprtp(i,gr%k_lb_zm,:), sclrpthlp(i,gr%k_lb_zm,:), & ! intent(in)
line 780 : 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 2523 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
line 2590 : 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 233 : 113 chars : write(*,'(A6,I5,1X,I5,A4,E30.20,A3,E30.20)') "field(", i, k, ") = ", field(i,k), " < ", threshold
line 241 : 105 chars : write(fstderr, *) "column", i, " field average = ", sum(field(i,:) * rho_ds(i,:) * dz(i,:)) &
line 583 : 102 chars : ! rho_ds_dz(i,k-num_hf_draw_points:k+num_hf_draw_points))
line 822 : 110 chars : * min( threshold, field(i,k_start:k_end:grid_dir_indx) ) ) * invrs_denom
line 827 : 106 chars : + mass_fraction * ( min( threshold, field(i,k_start:k_end:grid_dir_indx) ) - threshold )
line 832 : 110 chars : * max( threshold, field(i,k_start:k_end:grid_dir_indx) ) ) * invrs_denom
line 839 : 108 chars : + mass_fraction * ( max( threshold, field(i,k_start:k_end:grid_dir_indx) ) - threshold )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/output_netcdf.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1746 : 105 chars : status = nf90_def_var( ncid_zm, "columns", NF90_INT, (/column_id_zm/), column_var_id_zm )
line 1747 : 104 chars : status = nf90_def_var( ncid_zm, "altitude", NF90_PREC, (/vertical_id_zm/), vertical_var_id_zm )
line 1748 : 104 chars : status = nf90_def_var( ncid_zm, "time", NF90_PREC, (/time_id_zm/), time_var_id_zm )
line 1800 : 105 chars : status = nf90_def_var( ncid_zt, "columns", NF90_INT, (/column_id_zt/), column_var_id_zt )
line 1801 : 104 chars : status = nf90_def_var( ncid_zt, "altitude", NF90_PREC, (/vertical_id_zt/), vertical_var_id_zt )
line 1802 : 104 chars : status = nf90_def_var( ncid_zt, "time", NF90_PREC, (/time_id_zt/), time_var_id_zt )
line 1848 : 101 chars : status = nf90_put_var( ncid_zm, column_var_id_zm, column_list, start=(/1/), count=(/ngrdcol/) )
line 1849 : 101 chars : status = nf90_put_var( ncid_zt, column_var_id_zt, column_list, start=(/1/), count=(/ngrdcol/) )
line 1851 : 101 chars : !=================================== Allocate Output Arrays ===================================
line 1978 : 116 chars : !$acc upwp_out, vpwp_out, up2_out, vp2_out, wp3_out, rcm_out, cloud_frac_out, rtm_out, thlm_out )
----------------------------------------------------------------------------------------------------
--------------------- 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 688 : 104 chars : !$acc parallel loop gang vector collapse(2) default(present) reduction(.or.:l_any_adjustment_needed)
line 839 : 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/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 1002 : 103 chars : call windm_edsclrm_rhs( nzm, nzt, ngrdcol, gr, windm_edsclrm_scalar, dt, & ! intent(in)
line 1003 : 103 chars : lhs_diff, edsclrm(:,:,edsclr), & ! intent(in)
line 1004 : 103 chars : edsclrm_forcing(:,:,edsclr), & ! intent(in)
line 1005 : 103 chars : rho_ds_zm, invrs_rho_ds_zt, & ! intent(in)
line 1006 : 103 chars : l_imp_sfc_momentum_flux, wpedsclrp(:,gr%k_lb_zm,edsclr), & ! intent(in)
line 1007 : 103 chars : stats_metadata, & ! intent(in)
line 1008 : 106 chars : stats_zt, & ! intent(inout)
line 1009 : 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."
----------------------------------------------------------------------------------------------------
--------------------- 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 1330 : 101 chars : sclrprtp(:,:,sclr) ) ! Intent(inout)
line 1365 : 102 chars : sclrpthlp(:,:,sclr) ) ! Intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- 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 307 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 527 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 823 : 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 784 : 101 chars : call wpxp_term_pr1_lhs( nzm, ngrdcol, gr, C6rt_Skw_fnc, C6thl_Skw_fnc, C7_Skw_fnc, & ! Intent(in)
line 785 : 101 chars : invrs_tau_C6_zm, l_scalar_calc, & ! Intent(in)
line 786 : 102 chars : lhs_pr1_wprtp, lhs_pr1_wpthlp, lhs_pr1_wpsclrp ) ! Intent(out)
line 865 : 101 chars : stats_zt, stats_zm, stats_sfc, & ! InOut
line 866 : 101 chars : rtm, wprtp, thlm, wpthlp, sclrm, wpsclrp, & ! InOut
line 867 : 101 chars : err_info ) ! InOut
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/stats_zm_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 2449 : 104 chars : call stat_assign( var_index=stats_metadata%inorm_min_grid_dens, var_name="norm_min_grid_dens", &
----------------------------------------------------------------------------------------------------
--------------------- 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 3024 : 103 chars : call stat_update_var( stats_metadata%icovar_chi_eta_1, pdf_params%covar_chi_eta_1(icol,:), & ! In
line 3026 : 103 chars : call stat_update_var( stats_metadata%icovar_chi_eta_2, pdf_params%covar_chi_eta_2(icol,:), & ! In
line 3036 : 101 chars : call stat_update_var( stats_metadata%icorr_chi_eta_1, pdf_params%corr_chi_eta_1(icol,:), & ! In
line 3038 : 101 chars : call stat_update_var( stats_metadata%icorr_chi_eta_2, pdf_params%corr_chi_eta_2(icol,:), & ! In
line 3099 : 101 chars : chip2 = compute_variance_binormal( chi, pdf_params%chi_1(icol,:), pdf_params%chi_2(icol,:), &
line 3100 : 106 chars : pdf_params%stdev_chi_1(icol,:), pdf_params%stdev_chi_2(icol,:), &
line 4119 : 108 chars : 'floor stuff= ', floor( stats_metadata%stats_tout / stats_metadata%stats_tsamp )
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 16 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 )
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 1 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 299 : 104 chars : !$acc enter data create( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 465 : 118 chars : !$acc update device( wm_zt, wm_zm, thlm_forcing, rtm_forcing, um_ref, vm_ref, sclrm_forcing, edsclrm_forcing )
line 867 : 105 chars : call stat_update_var_pt( stats_metadata%ish, 1, wpthlp_sfc(i)*rho_zm(i,1)*Cp,& ! intent(in)
line 870 : 105 chars : call stat_update_var_pt( stats_metadata%ilh, 1, wprtp_sfc(i)*rho_zm(i,1)*Lv, & ! intent(in)
line 873 : 103 chars : call stat_update_var_pt( stats_metadata%iwpthlp_sfc, 1, wpthlp_sfc(i), & ! intent(in)
line 876 : 103 chars : call stat_update_var_pt( stats_metadata%iwprtp_sfc, 1, wprtp_sfc(i), & ! intent(in)
line 879 : 103 chars : call stat_update_var_pt( stats_metadata%iupwp_sfc, 1, upwp_sfc(i), & ! intent(in)
line 882 : 103 chars : call stat_update_var_pt( stats_metadata%ivpwp_sfc, 1, vpwp_sfc(i), & ! intent(in)
line 885 : 102 chars : call stat_update_var_pt( stats_metadata%iustar, 1, ustar(i), & ! intent(in)
line 888 : 103 chars : call stat_update_var_pt( stats_metadata%iT_sfc, 1, T_sfc(i), & ! intent(in)
line 894 : 103 chars : !$acc exit data delete( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 1097 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1101 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1105 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1106 : 112 chars : exner_zm(i,km1), exner_zm(i,k00), exner_zm(i,kp1), exner_zm(i,kp2) )
line 1109 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1110 : 108 chars : thlm_zm(i,km1), thlm_zm(i,k00), thlm_zm(i,kp1), thlm_zm(i,kp2) )
line 1113 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1114 : 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 538 : 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 501 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 6 check(s) failed.
[91m===============================================================
[91mCLUBBStandardsCheck FAILED
[91m THIS IS PRINTED IN ALL RED, CAPITAL LETTERS, AND USES
[91m AN EXCLAMATION MARK TO ENSURE THE DEVELOPERS FEEL SHAME!
[91m IF YOU ARE ONE OF THESE "DEVELOPERS" CHECK THE
[91m LOG FILE FOR DETAILS: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gfortran_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
[93mBuild completed successfully, but some source code checks have failed.[0m
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Run CLUBB's G Unit Tests)
[Pipeline] sh
+ source /etc/profile.d/larson-group.sh
++ export GIT_EDITOR=vi
++ GIT_EDITOR=vi
++ export SVN_EDITOR=vi
++ SVN_EDITOR=vi
++ export OMP_STACKSIZE=1048579
++ OMP_STACKSIZE=1048579
++ export LMOD_ROOT=/opt/lmod/
++ LMOD_ROOT=/opt/lmod/
++ source /opt/lmod//lmod/lmod/init/bash
+++ '[' -z '' ']'
+++ case "$-" in
+++ __lmod_vx=x
+++ '[' -n x ']'
+++ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for this output (/usr/local/spack/opt/spack/linux-pop22-skylake_avx512/gcc-12.2.0/lmod-8.7.37-fq24mybyn2li6got2bxzk62ejh5d3p4z/lmod/lmod/init/bash)
Shell debugging restarted
+++ unset __lmod_vx
+++ find /usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core -print -quit
++ export MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
++ MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
+ module load gcc netcdf-fortran
+ '[' -z '' ']'
+ case "$-" in
+ __lmod_sh_dbg=x
+ '[' -n x ']'
+ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for Lmod's output
Lmod Warning:
-------------------------------------------------------------------------------
The following dependent module(s) are not currently loaded:
gcc-runtime/13.1.0-wiuv7wc (required by: mpich/4.2.1, hdf5/1.14.3-nhc23hl,
netcdf-c/4.9.2-xduzp77, netcdf-fortran/4.5.3), glibc/2.35-tdr5czg (required by:
mpich/4.2.1, hdf5/1.14.3-nhc23hl, netcdf-c/4.9.2-xduzp77, netcdf-fortran/4.5.3)
-------------------------------------------------------------------------------
Shell debugging restarted
+ unset __lmod_sh_dbg
+ return 0
+ python run_scripts/run_G_unit_tests.py
Running G_unit_tests
Performing KK integrals tests
Percent Difference Tolerance (maximum acceptable percent difference between the CLUBB results and the MATLAB results using the same equations): 1.50000E-12
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 1
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: 3.940566852245362E-04
MATLAB result: 3.940566852245361E-04
Percent difference: 1.37569E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: 3.927115990964176E-04
MATLAB result: 3.927115990964175E-04
Percent difference: 2.76081E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: 3.959157231937519E-04
MATLAB result: 3.959157231937519E-04
Percent difference: 1.36923E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: 3.961224408774508E-04
MATLAB result: 3.961224408774508E-04
Percent difference: 1.36852E-14%
Agreement with regards to integral: Trivariate NNL
Obtained result: 5.957899876478518E-10
MATLAB result: 5.957899876478559E-10
Percent difference: 6.94188E-13%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 5.299646183560617E-10
MATLAB result: 5.299646183560650E-10
Percent difference: 6.24329E-13%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.543555152546510E-10
MATLAB result: 5.543555152546526E-10
Percent difference: 2.98430E-13%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.286452253001267E-10
MATLAB result: 5.286452253001284E-10
Percent difference: 3.12944E-13%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -9.854231981928352E-04
MATLAB result: -9.854231981928350E-04
Percent difference: 2.20048E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549015E-04
MATLAB result: -9.922448817549015E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 8.940070763287877E-09
MATLAB result: 8.940070763287870E-09
Percent difference: 7.40201E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 8.942709549399747E-09
MATLAB result: 8.942709549399743E-09
Percent difference: 3.69991E-14%
Agreement with regards to integral: Bivariate LL
Obtained result: 1.007487885941809E-02
MATLAB result: 1.007487885941809E-02
Percent difference: 1.72183E-14%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 2
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -4.441972023148827E-05
MATLAB result: -4.441972023148829E-05
Percent difference: 4.57652E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -4.559373570842309E-05
MATLAB result: -4.559373570842311E-05
Percent difference: 4.45868E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -5.040842768062482E-05
MATLAB result: -5.040842768062482E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -5.038775591225493E-05
MATLAB result: -5.038775591225493E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.629614112194526E-09
MATLAB result: 4.629614112194527E-09
Percent difference: 1.78672E-14%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.604609365066508E-09
MATLAB result: 4.604609365066509E-09
Percent difference: 1.79642E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 4.685272472815394E-09
MATLAB result: 4.685272472815394E-09
Percent difference: 1.76549E-14%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 4.683452323676534E-09
MATLAB result: 4.683452323676534E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -1.088125285831538E-04
MATLAB result: -1.088125285831538E-04
Percent difference: 2.49099E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549016E-05
MATLAB result: -9.922448817549016E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 7.907812698669825E-10
MATLAB result: 7.907812698669821E-10
Percent difference: 5.23015E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 6.330953526469324E-10
MATLAB result: 6.330953526469324E-10
Percent difference: 0.00000E+00%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 3
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -7.622010070476262E-05
MATLAB result: -7.622010070476262E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -7.695119321547280E-05
MATLAB result: -7.695119321547280E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.851589269710427E-09
MATLAB result: 4.851589269710427E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.837259568295437E-09
MATLAB result: 4.837259568295437E-09
Percent difference: 1.71002E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -4.609761356905442E-05
MATLAB result: -4.609761356905441E-05
Percent difference: 1.46998E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 3.254808634091251E-10
MATLAB result: 3.254808634091251E-10
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
==================================================
Success!
==================================================
Performing correlation cholesky matrix setup tests
correlation matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
--
correlation matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 1.0000000000000000 0.86602540378443860
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
--
==================================================
Performing correlation matrix approximation tests
correlation cholesky matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.86602540378443860 1.0000000000000000 0.96650635094610959
0.86602540378443860 0.96650635094610959 1.0000000000000000
--
==================================================
Success!
==================================================
Performing hole_filling_hm_one_lev_tests
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
==================================================
Performing fill_holes_hydromet_tests
testset_in(:, 1 ) = 1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 6.0000000000000000
testset_in(:, 2 ) = 6.0000000000000000 5.0000000000000000 4.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000
testset_in(:, 3 ) = 3.0000000000000000 3.0000000000000000 3.0000000000000000 4.0000000000000000 4.0000000000000000 4.0000000000000000
testset_in(:, 4 ) = 1.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 5 ) = 1.0000000000000000 -2.0000000000000000 -4.0000000000000000 -3.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 6 ) = 2.0000000000000000 -1.0000000000000000 -1.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 7 ) = 1.0000000000000000 4.0000000000000000 2.0000000000000000 3.0000000000000000 -1.0000000000000000 0.0000000000000000
Hole-filling was conservative
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
Nc-Ncn equations unit test
==========================
Iteration = 1
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.15865525393145707
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 4.7896918412674905E-002
Nc_in_cloud (from Ncn) = 66978771.649335034
Ncm (from Ncn) = 3208076.7610693830
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 1 is a success!
======================================================
Iteration = 2
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.0000000000000000
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 0.0000000000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 2 is a success!
======================================================
Iteration = 3
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 0.99957093966680322
cloud fraction = 0.95210308158732504
Nc_in_cloud (from Ncn) = 51650441.529586732
Ncm (from Ncn) = 49176544.545665480
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 3 is a success!
======================================================
Iteration = 4
mu_chi_1 = 1.0000000000000000E-008
mu_chi_2 = -1.0000000000000000E-008
sigma_chi_1 = 5.0000000000000002E-005
sigma_chi_2 = 5.0000000000000002E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.10000000000000001
mixture fraction = 0.50000000000000000
cloud fraction (1st PDF component) = 0.50007978845554835
cloud fraction (2nd PDF component) = 0.49992021154445165
cloud fraction = 0.50000000000000000
Nc_in_cloud (from Ncn) = 53982783.648312390
Ncm (from Ncn) = 26991391.824156195
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 4 is a success!
======================================================
Iteration = 5
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 9.9999999999999995E-007
sigma_chi_2 = 1.5000000000000000E-006
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.50000000000000000
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 1.0000000000000000
cloud fraction = 1.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 50000000.000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 5 is a success!
======================================================
Iteration = 6
mu_chi_1 = -9.9999999999999995E-008
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.0000000000000000
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.49601064368536840
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 9.9545377003631155E-002
Nc_in_cloud (from Ncn) = 49999999.999999993
Ncm (from Ncn) = 4977268.8501815572
Ncnm (from Nc_in_cloud) = 49999999.999999993
Ncnm (from Ncm) = 49999999.999999993
Test iteration 6 is a success!
======================================================
Iteration = 7
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 0.0000000000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.48670438618290790
cloud fraction = 0.49732907640173213
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 24866453.820086606
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 7 is a success!
======================================================
Iteration = 8
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.10796556745540581
Nc_in_cloud (from Ncn) = 81037694.237492040
Ncm (from Ncn) = 8749280.6436284967
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 8 is a success!
======================================================
Iteration = 9
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.20034324826655747
Nc_in_cloud (from Ncn) = 52151569.416871324
Ncm (from Ncn) = 10448214.819174858
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 9 is a success!
======================================================
Iteration = 10
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.16861219748026607
Nc_in_cloud (from Ncn) = 60768988.354576342
Ncm (from Ncn) = 10246392.665117815
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 10 is a success!
======================================================
Success!
Performing corr_mtx_test
==================================================
Expected:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Actual:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Success!
==================================================
==================================================
Performing SILHS category test
==================================================
Success!
Performing hydrometeor mu and sigma values unit test
====================================================
PDF parameter set 1:
hmm = 2.0000000000000002E-005
mixt_frac = 2.7000000000000000E-002
precip_frac = 8.0000000000000002E-002
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 5.4470709146968145E-002
hmp2_ip_on_hmm2_ip = 3.6000000000000001
hmm_ip = 2.5000000000000001E-004
hmp2 = 2.2599999999999997E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 2:
hmm = 5.0000000000000004E-006
mixt_frac = 0.50000000000000000
precip_frac = 0.50000000000000000
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 0.50000000000000000
hmp2_ip_on_hmm2_ip = 1.0000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 7.5000000000000012E-011
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 3:
hmm = 1.0000000000000000E-010
mixt_frac = 0.25000000000000000
precip_frac = 0.40000000000000002
precip_frac_1 = 0.75000000000000000
precip_frac_2 = 0.28333333333333338
hmp2_ip_on_hmm2_ip = 0.50000000000000000
hmm_ip = 2.5000000000000002E-010
hmp2 = 2.7500000000000010E-020
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 4:
hmm = 4.0000000000000003E-005
mixt_frac = 1.0000000000000000E-002
precip_frac = 0.10000000000000001
precip_frac_1 = 1.0000000000000000E-002
precip_frac_2 = 0.10090909090909091
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 5:
hmm = 4.0000000000000003E-005
mixt_frac = 0.98999999999999999
precip_frac = 0.10000000000000001
precip_frac_1 = 0.10090909090909091
precip_frac_2 = 1.0000000000000278E-002
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 6:
hmm = 1.0000000000000001E-005
mixt_frac = 0.12500000000000000
precip_frac = 1.0000000000000000
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 1.0000000000000000
hmp2_ip_on_hmm2_ip = 1.5000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 1.5000000000000005E-010
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 7:
hmm = 6.0000000000000002E-005
mixt_frac = 5.0000000000000003E-002
precip_frac = 0.11000000000000000
precip_frac_1 = 0.90000000000000002
precip_frac_2 = 6.8421052631578952E-002
hmp2_ip_on_hmm2_ip = 5.0000000000000000
hmm_ip = 5.4545454545454548E-004
hmp2 = 1.9276363636363637E-007
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 8:
hmm = 7.9999999999999996E-007
mixt_frac = 0.80000000000000004
precip_frac = 0.25000000000000000
precip_frac_1 = 0.10000000000000001
precip_frac_2 = 0.85000000000000009
hmp2_ip_on_hmm2_ip = 0.75000000000000000
hmm_ip = 3.1999999999999999E-006
hmp2 = 3.8399999999999998E-012
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 9:
hmm = 1.3499999999999999E-005
mixt_frac = 0.10000000000000001
precip_frac = 0.10000000000000001
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 5.5555555555555559E-002
hmp2_ip_on_hmm2_ip = 2.0000000000000000
hmm_ip = 1.3499999999999997E-004
hmp2 = 5.2852499999999985E-009
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 10 (randomly generated):
Random seed values = 2025 11 14 -360 3 12 54 949
hmm = 2.4553947299393144E-005
mixt_frac = 0.12496832104129976
precip_frac = 0.74387703780731296
precip_frac_1 = 0.32747407509859106
precip_frac_2 = 0.80334594662242798
hmp2_ip_on_hmm2_ip = 2.8116078133423610
hmm_ip = 3.3008072640297538E-005
hmp2 = 2.4863297639315251E-009
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Checking values with precipitation in 1st PDF component only.
Checking values with precipitation in 2nd PDF component only.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for ADG1
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 10 (randomly generated):
Random seed values = 2025 11 14 -360 3 12 54 950
wm = -2.8921054012137182E-002
rtm = 3.3462922661489929E-003
thlm = 304.14597280643653
wp2 = 0.32068088381519000
rtp2 = 1.4923659000415826E-006
thlp2 = 2.4872132860353489
Skw = -0.19313842710226847
wp3 = -3.5073434710539711E-002
Skrt = -1.6284154284224979
rtp3 = -2.9687811510654807E-009
Skthl = -1.1763964435456669
thlp3 = -4.6144850816697005
wprtp = 6.8984394383493674E-004
wpthlp = -0.37762816898298485
rtpthlp = -1.7339504486315964E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for LY93
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 10 (randomly generated):
Random seed values = 2025 11 14 -360 3 12 54 950
wm = -2.8921054012137182E-002
rtm = 3.3462922661489929E-003
thlm = 304.14597280643653
wp2 = 0.32068088381519000
rtp2 = 1.4923659000415826E-006
thlp2 = 2.4872132860353489
Skw = -0.19313842710226847
wp3 = -3.5073434710539711E-002
Skrt = -1.6284154284224979
rtp3 = -2.9687811510654807E-009
Skthl = -1.1763964435456669
thlp3 = -4.6144850816697005
wprtp = 6.8984394383493674E-004
wpthlp = -0.37762816898298485
rtpthlp = -1.7339504486315964E-003
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for TSDADG
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
PDF parameter set 10 (randomly generated):
Random seed values = 2025 11 14 -360 3 12 54 950
wm = -2.8921054012137182E-002
rtm = 3.3462922661489929E-003
thlm = 304.14597280643653
wp2 = 0.32068088381519000
rtp2 = 1.4923659000415826E-006
thlp2 = 2.4872132860353489
Skw = -0.19313842710226847
wp3 = -3.5073434710539711E-002
Skrt = -1.6284154284224979
rtp3 = -2.9687811510654807E-009
Skthl = -1.1763964435456669
thlp3 = -4.6144850816697005
wprtp = 6.8984394383493674E-004
wpthlp = -0.37762816898298485
rtpthlp = -1.7339504486315964E-003
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2025 11 14 -360 3 12 54 952
wm = -2.8921054012137182E-002
rtm = 3.3462922661490059E-003
thlm = 304.14597280643648
wp2 = 0.32068088340591305
rtp2 = 1.9229404434963161E-006
thlp2 = 2.4869876035380205
Skw = 0.79534142777916284
wp3 = 0.14443192897923857
Skrt = 1.4652380083695125
rtp3 = 3.9071215858202215E-009
Skthl = 4.3652541888394865
thlp3 = 17.120638943171834
wprtp = 7.2527082044207699E-004
wpthlp = 0.23877171425225843
rtpthlp = 1.9681683176768611E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new hybrid PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2025 11 14 -360 3 12 54 953
wm = -2.8921054012137182E-002
rtm = 3.3462922661489994E-003
thlm = 304.14597280643642
wp2 = 0.32068088344683998
rtp2 = 1.8921851188656233E-006
thlp2 = 2.4870036338050552
Skw = 0.69651275404844526
wp3 = 0.12648489960597131
Skrt = -1.2065887122365648
rtp3 = -3.1405426239245006E-009
Skthl = 1.7922321911369625
thlp3 = 7.0292486048200287
wprtp = 2.6122396558150884E-004
wpthlp = -0.29760194219259800
rtpthlp = -1.9523718205877259E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing spurious source unit test
====================================
Test 1
Vertical integral of <rt> before advance_xm_wpxp = 28.316670259590918
Vertical integral of <rt> after advance_xm_wpxp = 28.454761250379871
Vertical integral of <rt> forcing = -3.9696697370157397E-005
Flux of <rt> at the surface = 5.0000000000000001E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 1.2197274440461925E-018
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 793918.42642888101
Vertical integral of <thl> after advance_xm_wpxp = 793908.03522151255
Vertical integral of <thl> forcing = -4.4637357896125940E-002
Flux of <thl> at the surface = 1.0000000000000000E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 1.2514503322513804E-012
The spurious source of <thl> is within acceptable limits.
Test 2
Vertical integral of <rt> before advance_xm_wpxp = 30.538572189940858
Vertical integral of <rt> after advance_xm_wpxp = 30.675947141863642
Vertical integral of <rt> forcing = -4.3435804138101911E-005
Flux of <rt> at the surface = 5.0135231054742615E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -4.5577148825859393E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 793526.20879985753
Vertical integral of <thl> after advance_xm_wpxp = 793514.18036282028
Vertical integral of <thl> forcing = -5.0121836335417891E-002
Flux of <thl> at the surface = 1.0027046210948522E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = 2.9019842084920811E-013
The spurious source of <thl> is within acceptable limits.
Test 3
Vertical integral of <rt> before advance_xm_wpxp = 24.686709013987230
Vertical integral of <rt> after advance_xm_wpxp = 24.705459874640248
Vertical integral of <rt> forcing = -3.8022584664444255E-005
Flux of <rt> at the surface = 1.0052545350780095E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 3.6673138484322187E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 792147.05323497509
Vertical integral of <thl> after advance_xm_wpxp = 792132.43931249762
Vertical integral of <thl> forcing = -5.0723583993502311E-002
Flux of <thl> at the surface = 2.0105090701560192E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -1.5436610323327216E-012
The spurious source of <thl> is within acceptable limits.
Test 4
Vertical integral of <rt> before advance_xm_wpxp = 34.401263410406429
Vertical integral of <rt> after advance_xm_wpxp = 34.406239496184604
Vertical integral of <rt> forcing = -6.8429344850707552E-005
Flux of <rt> at the surface = 8.5016297444666662E-005
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -4.2853090867489563E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 798397.67855735740
Vertical integral of <thl> after advance_xm_wpxp = 798375.94280580617
Vertical integral of <thl> forcing = -7.4152831119599485E-002
Flux of <thl> at the surface = 1.7003259488933332E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -4.9557580261705425E-014
The spurious source of <thl> is within acceptable limits.
Test 5
Vertical integral of <rt> before advance_xm_wpxp = 20.800632473975561
Vertical integral of <rt> after advance_xm_wpxp = 20.845095542640561
Vertical integral of <rt> forcing = -3.4395721906668523E-005
Flux of <rt> at the surface = 1.8260595078999807E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 1.7279472123987727E-018
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 791086.92245564435
Vertical integral of <thl> after advance_xm_wpxp = 791076.16835060238
Vertical integral of <thl> forcing = -3.9499135821001272E-002
Flux of <thl> at the surface = 3.6521190157999617E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -1.3577056146019117E-012
The spurious source of <thl> is within acceptable limits.
Success!
Running tuner tests
-------------------
goldstein_price pass
rastrigin pass
himmelblau pass
eggholder pass
schaffer pass
==================================================
Performing 'mean vertical velocity in clouds' tests
Success!
==================================================
==================================================
Performing 'Smooth Heaviside Peskin' tests
Case 1: heaviside_smth_range = 0, smooth Heaviside collapses to Heaviside
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
Case 2: heaviside_smth_range = 2.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Expected outcome: 0.0000000000000000 9.0845056908100005E-002 0.90915494309000000 1.0000000000000000
True outcome: 0.0000000000000000 9.0845056909999999E-002 0.90915494309000000 1.0000000000000000
Success!
==================================================
==================================================
Performing 'Smooth Min / Max' tests
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
(Simple) min:
Expected outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
True outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
(Simple) max:
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Smooth min:
Expected outcome: -2.2247448714824141 -1.3660254037844388 -0.36602540378443871 -0.22474487138241406
True outcome: -2.2247448714824136 -1.3660254037844388 -0.36602540378443871 -0.22474487138241384
Smooth max:
Expected outcome: 0.22474487138241406 0.36602540378443871 1.3660254037844388 2.2247448714824141
True outcome: 0.22474487138241384 0.36602540378443871 1.3660254037844388 2.2247448714824136
Testing smooth min and max with min_max_smth_mag=1e-7, input_var1 = equidistant grid on [-1, 1), input_var2 = 0
Input: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
True outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
Smooth max:
Expected outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
True outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
Input_var1: -1000000.0000000000 -10000000.000000000 -100000000.00000000 -1000000000.0000000
Input_var2: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
True outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
Smooth max:
Expected outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
True outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
Input_var1: -1000000.0000000000
Input_var2: 1.0000000000000000
Smooth min:
Expected outcome: -1000000.0000002501
True outcome: -1000000.0000002501
Smooth max:
Expected outcome: 1.0000002500019036
True outcome: 1.0000002500019036
Success!
==================================================
====================== Running fill holes tests ======================
-------------- easy_fill_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: easy_fill_test
-----------------------------------------------
-------------- below_thresh_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: below_thresh_test
-----------------------------------------------
====================== End fill holes tests ======================
Running /home/jenkins/workspace/clubb_G_unit_gfortran_tests/run_scripts/../install/latest/G_unit_tests ...
All tests have succeeded.
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Declarative: Post Actions)
[Pipeline] script
[Pipeline] {
[Pipeline] cleanWs
[WS-CLEANUP] Deleting project workspace...
[WS-CLEANUP] Deferred wipeout is used...
[WS-CLEANUP] done
[Pipeline] }
[Pipeline] // script
[Pipeline] }
[Pipeline] // stage
[Pipeline] }
[Pipeline] // withEnv
[Pipeline] }
[Pipeline] // node
[Pipeline] End of Pipeline
Finished: SUCCESS
![[Jenkins]](/jenkins/static/9dada5b8/images/svgs/logo.svg){kind=logo}
