Skipping 366 KB..
Full LogTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTestTest
[918/935] Building Fortran object src/CMakeFiles/clubb_driver_lib.dir/input_fields.F90.o
[919/935] Building Fortran object src/CMakeFiles/clubb_driver_lib.dir/generalized_grid_test.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/generalized_grid_test.F90:2457:8:
2457 | use constants_clubb, only: &
| 1
Warning: Unused parameter 'eps' which has been explicitly imported at (1) [-Wunused-parameter]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/generalized_grid_test.F90:711:29:
711 | i, sclr, edsclr, hm_idx
| 1
Warning: Unused variable 'hm_idx' declared at (1) [-Wunused-variable]
[920/935] Building Fortran object src/CMakeFiles/clubb_driver_lib.dir/clubb_driver.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90:400:25:
400 | rho_ds_zm_dycore_init
| 1
Warning: Unused PRIVATE module variable 'rho_ds_zm_dycore_init' declared at (1) [-Wunused-value]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90:4666:41:
4666 | grid_adapt_in_time_method, &
| 1
Warning: Unused dummy argument 'grid_adapt_in_time_method' at (1) [-Wunused-dummy-argument]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90:4665:33:
4665 | l_add_dycore_grid, &
| 1
Warning: Unused dummy argument 'l_add_dycore_grid' at (1) [-Wunused-dummy-argument]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90:4208:49:
4208 | subroutine clean_up_clubb(ngrdcol, clubb_params, err_info)
| 1
Warning: Unused dummy argument 'clubb_params' at (1) [-Wunused-dummy-argument]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90:4208:59:
4208 | subroutine clean_up_clubb(ngrdcol, clubb_params, err_info)
| 1
Warning: Unused dummy argument 'err_info' at (1) [-Wunused-dummy-argument]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90:3198:14:
3198 | ret_code, &
| 1
Warning: Unused variable 'ret_code' declared at (1) [-Wunused-variable]
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90:2668:72:
2668 | subroutine set_case_initial_conditions(ngrdcol, clubb_params, err_info)
| 1
Warning: Unused dummy argument 'err_info' at (1) [-Wunused-dummy-argument]
[921/935] Linking Fortran static library src/libclubb_driver_lib.a
[922/935] Generating Fortran dyndep file src/CMakeFiles/clubb_tuner.dir/Fortran.dd
[923/935] Generating Fortran dyndep file src/CMakeFiles/clubb_standalone.dir/Fortran.dd
[924/935] Generating Fortran dyndep file src/CMakeFiles/clubb_thread_test.dir/Fortran.dd
[925/935] Generating Fortran dyndep file src/CMakeFiles/clubb_driver_test.dir/Fortran.dd
[926/935] Building Fortran object src/CMakeFiles/clubb_thread_test.dir/clubb_thread_test.F90.o
[927/935] Building Fortran object src/CMakeFiles/clubb_driver_test.dir/clubb_driver_test.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90: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]
[928/935] Building Fortran object src/CMakeFiles/clubb_standalone.dir/clubb_standalone.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_standalone.F90:66:26:
66 | character(len=10) :: arg
| 1
Warning: Unused variable 'arg' declared at (1) [-Wunused-variable]
[929/935] Linking Fortran executable src/clubb_thread_test
[930/935] Building Fortran object src/CMakeFiles/clubb_tuner.dir/error.F90.o
[931/935] Linking Fortran executable src/clubb_standalone
[932/935] Linking Fortran executable src/clubb_driver_test
[933/935] Building Fortran object src/CMakeFiles/clubb_tuner.dir/clubb_tuner.F90.o
/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_tuner.F90:294:6:
294 | use nrtype, only: &
| 1
Warning: Unused parameter 'sp' which has been explicitly imported at (1) [-Wunused-parameter]
[934/935] Linking Fortran executable src/clubb_tuner
[934/935] Install the project...
-- Install configuration: "Debug"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/libclubb_driver_lib.a
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_tuner
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_tuner" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/G_unit_tests
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/G_unit_tests" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_thread_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_thread_test" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_standalone
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_standalone" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
-- Installing: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_driver_test
-- Set runtime path of "/home/jenkins/workspace/clubb_G_unit_gfortran_tests/install/gfortran_DEBUG_PRECdouble/clubb_driver_test" to "/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir:/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-13.1.0/netcdf-fortran-4.5.3-ivafocv4fav4rj3ckswo6vj33utvtdir/lib"
check_for_missing_threadprivate.py has begun.
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/corr_varnce_module.F90 failed:
Missing threadprivate for:
corr_array_n_cloud_def
corr_array_n_below_def
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/error_code.F90 failed:
Missing threadprivate for:
clubb_debug_level
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sponge_layer_damping.F90 failed:
Missing threadprivate for:
uv_sponge_damp_settings
thlm_sponge_damp_profile
rtm_sponge_damp_settings
wp3_sponge_damp_settings
rtm_sponge_damp_profile
up2_vp2_sponge_damp_settings
uv_sponge_damp_profile
thlm_sponge_damp_settings
wp2_sponge_damp_settings
wp3_sponge_damp_profile
up2_vp2_sponge_damp_profile
wp2_sponge_damp_profile
File /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90 failed:
Missing threadprivate for:
ret_code
[91mcheck_for_missing_threadprivate.py did not pass.[0m
[91mfailedFiles = ['/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/corr_varnce_module.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/error_code.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/sponge_layer_damping.F90', '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90'][0m
[91m===============================================================
[91mTHREADPRIVATE CHECK FAILED
[91m THIS IS PRINTED IN ALL RED, CAPITAL LETTERS, AND USES
[91m AN EXCLAMATION MARK TO ENSURE THE DEVELOPERS FEEL SHAME!
[91m IF YOU ARE ONE OF THESE "DEVELOPERS" CHECK THE
[91m LOG FILE FOR DETAILS: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gfortran_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
CLUBBStandardsCheck.py has begun.
--------------------- 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) )
----------------------------------------------------------------------------------------------------
--------------------- 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 914 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1084 : 101 chars : rtm_ref_flip, thlm_ref_flip, um_ref_flip, vm_ref_flip, ug_flip, vg_flip, & ! Intent(in)
line 1098 : 111 chars : um_flip, vm_flip, upwp_flip, vpwp_flip, up2_flip, vp2_flip, up3_flip, vp3_flip, & ! Intent(inout)
line 1100 : 111 chars : wp2_flip, wp3_flip, rtp2_flip, rtp3_flip, thlp2_flip, thlp3_flip, rtpthlp_flip, & ! Intent(inout)
line 1106 : 109 chars : wp2rtp_flip, wp2thlp_flip, uprcp_flip, vprcp_flip, rc_coef_zm_flip, wp4_flip, & ! intent(inout)
line 1107 : 105 chars : wpup2_flip, wpvp2_flip, wp2up2_flip, wp2vp2_flip, ice_supersat_frac_flip, & ! intent(inout)
line 1594 : 115 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 1621 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2up_implicit, nzt, ngrdcol, &
line 1625 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2up_explicit, nzt, ngrdcol, &
line 1629 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wp2vp_implicit, nzt, ngrdcol, &
line 1633 : 102 chars : pdf_implicit_coefs_terms_flip%term_wp2vp_explicit, nzt, ngrdcol, &
line 1667 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpup2_implicit, nzt, ngrdcol, &
line 1671 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpup2_explicit, nzt, ngrdcol, &
line 1675 : 102 chars : pdf_implicit_coefs_terms_flip%coef_wpvp2_implicit, nzt, ngrdcol, &
line 1679 : 102 chars : pdf_implicit_coefs_terms_flip%term_wpvp2_explicit, nzt, ngrdcol, &
line 2504 : 105 chars : tolerance = max( 1.e-8_core_rknd, epsilon(tolerance) ) ! max statement for single precision runs
----------------------------------------------------------------------------------------------------
--------------------- Implicit None Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver.F90 ---------------------
WARNING: 'use' statements without 'only' found in the following lines:
WARNING: Check that comma is on same line as 'use', as CLUBB requires.
3150 : 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 201 : 110 chars : l_modify_ic_with_cubic_int, & ! Flag for interpolating the sounding profile with Steffen's monotone cubic
line 202 : 110 chars : ! method to obtain smoother initial condition profile, which is found to be
line 203 : 111 chars : ! beneficial to achive a better numerical solution convergence. If this flag
line 204 : 118 chars : ! is turned off, the initial conditions will be generated with linear interpolation.
line 205 : 111 chars : ! This is done on a case-by-case basis, since using the monotone cubic method
line 206 : 109 chars : ! requires a special sounding.in file with many additional sounding levels.
line 207 : 109 chars : l_modify_bc_for_cnvg_test, & ! Flag to activate modifications on boundary condition for convergence test
line 744 : 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 760 : 102 chars : !$omp rtp2_mc, thlp2_mc, rtpthlp_mc, Ncm_mc, Ncm, Nccnm, em, tau_zm, tau_zt, Kh_zt, Kh_zm, rfrzm, &
line 761 : 105 chars : !$omp rrm, Nrm, Nc_in_cloud, wpNcp, radht, Frad, Frad_SW_up, Frad_LW_up, Frad_SW_down, Frad_LW_down, &
line 762 : 102 chars : !$omp wpsclrp_sfc, wpedsclrp_sfc, sclrm, wpsclrp, sclrp2, sclrp3, sclrprtp, sclrpthlp, sclrpthvp, &
line 807 : 101 chars : !----------------------------------- Optional Input Variables -----------------------------------
line 1029 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 1721 : 102 chars : open(unit=iunit_grid_adaptation, file=fname_grid_adaptation, status='replace', action='write')
line 2036 : 105 chars : allocate( um_ref(ngrdcol, gr%nzt) ) ! Reference u wind for nudging; Michael Falk, 17 Oct 2007
line 2037 : 105 chars : allocate( vm_ref(ngrdcol, gr%nzt) ) ! Reference v wind for nudging; Michael Falk, 17 Oct 2007
line 2529 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 2560 : 106 chars : ! Similar call to above, checking parameter values, but because we've placed this in "init_clubb_case"
line 2568 : 116 chars : write(fstderr, *) "Fatal error calling check_clubb_settings_api and/or check_parameters_api in clubb_driver"
line 2611 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 2623 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 2624 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 2627 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 2629 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 2635 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 2636 : 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 2637 : 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 2638 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 2639 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 2640 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 2641 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 2642 : 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 2643 : 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 2644 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 2759 : 103 chars : invrs_rho_ds_zt(i,k) = invrs_rho_ds_zt_init(k) ! inv. dry, static density: t-levs
line 2765 : 103 chars : um_pert(i,k) = zero ! Variables used to track perturbed version of winds.
line 2767 : 105 chars : wm_zt(i,k) = wm_zt_init(k) ! Imposed large scale w - Thermodynamic levels
line 2961 : 113 chars : if ( clubb_config_flags%iiPDF_type == iiPDF_new .or. clubb_config_flags%iiPDF_type == iiPDF_new_hybrid ) then
line 3171 : 101 chars : !----------------------------------- Optional Input Variables -----------------------------------
line 3328 : 103 chars : um(i,:), vm(i,:), rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), &
line 3666 : 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 3668 : 104 chars : !$acc lh_rt_clipped, lh_thl_clipped, lh_rc_clipped, lh_rv_clipped, lh_Nc_clipped, &
line 3672 : 116 chars : !$acc pdf_params%cloud_frac_1, pdf_params%cloud_frac_2, pdf_params%cthl_1, pdf_params%cthl_2, &
line 3673 : 112 chars : !$acc pdf_params%crt_1, pdf_params%crt_2, pdf_params%stdev_chi_1, pdf_params%stdev_chi_2, &
line 3689 : 102 chars : wp2_zt = zm2zt_api( gr%nzm, gr%nzt, ngrdcol, gr, wp2, w_tol_sqd ) ! Positive definite quantity
line 3695 : 109 chars : thlm(i,:), p_in_Pa(i,:), exner(i,:), rho(i,:), rho_zm(i,:), rtm(i,:), &! In
line 3696 : 109 chars : rcm(i,:), cloud_frac(i,:), wm_zt(i,:), wm_zm(i,:), wp2_zt(i,:), & ! In
line 3697 : 106 chars : hydromet(i,:,:), Nc_in_cloud(i,:), & ! In
line 3713 : 103 chars : stats_zt(i), stats_zm(i), & ! Inout
line 3714 : 103 chars : stats_sfc(i), stats_lh_zt(i), & ! Inout
line 3715 : 102 chars : Nccnm(i,:), & ! Inout
line 3716 : 117 chars : hydromet_mc(i,:,:), Ncm_mc(i,:), rcm_mc(i,:), rvm_mc(i,:), & ! Out
line 3717 : 107 chars : thlm_mc(i,:), hydromet_vel_zt(i,:,:), & ! Out
line 3718 : 102 chars : hydromet_vel_covar_zt_impc(i,:,:), & ! Out
line 3719 : 102 chars : hydromet_vel_covar_zt_expc(i,:,:), & ! Out
line 3720 : 110 chars : wprtp_mc(i,:), wpthlp_mc(i,:), rtp2_mc(i,:), & ! Out
line 3721 : 105 chars : thlp2_mc(i,:), rtpthlp_mc(i,:) ) ! Out
line 3746 : 103 chars : stats_zt(i), stats_zm(i), stats_sfc(i), & ! Inout
line 3747 : 103 chars : hydromet(i,:,:), hydromet_vel_zt(i,:,:), hydrometp2(i,:,:), & ! Inout
line 3748 : 103 chars : K_hm(i,:,:), Ncm(i,:), Nc_in_cloud(i,:), rvm_mc(i,:), & ! Inout
line 3749 : 103 chars : thlm_mc(i,:), err_info, & ! Inout
line 3750 : 101 chars : wphydrometp(i,:,:), wpNcp(i,:) ) ! Out
line 3759 : 106 chars : !$acc update device( rcm_mc, rvm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3817 : 104 chars : !$acc update host( rvm_mc, rcm_mc, thlm_mc, wprtp_mc, wpthlp_mc, rtp2_mc, thlp2_mc, rtpthlp_mc )
line 3823 : 105 chars : call stat_update_var( stats_metadata%irtm_mc, rvm_mc(i,:)+rcm_mc(i,:), stats_zt(i) ) ! kg/kg/s
line 3825 : 106 chars : call stat_update_var( stats_metadata%iwprtp_mc, wprtp_mc(i,:), stats_zm(i) ) ! m*(kg/kg)/s^2
line 3827 : 104 chars : call stat_update_var( stats_metadata%irtp2_mc, rtp2_mc(i,:), stats_zm(i) ) ! (kg/kg)^2/s
line 3829 : 104 chars : call stat_update_var( stats_metadata%irtpthlp_mc, rtpthlp_mc(i,:), stats_zm(i) ) ! K*(kg/kg)/s
line 3846 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, X_nl_all_levs, &
line 3872 : 105 chars : !$acc update host( rho, rho_zm, p_in_Pa, exner, cloud_frac, ice_supersat_frac, thlm, rtm, rcm )
line 3918 : 106 chars : call update_radiation_variables( gr%nzm, gr%nzt, radht(i,:), Frad_SW_up(i,:), Frad_LW_up(i,:), &
line 4017 : 104 chars : clubb_config_flags%grid_remap_method, & ! intent(in)
line 4086 : 101 chars : .and. (stats_metadata%l_stats_last) ) then ! only adapt grid when the average of the last
line 4191 : 103 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_scheme = ', time_microphys_scheme
line 4192 : 104 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_microphys_advance = ', time_microphys_advance
line 4194 : 103 chars : write(unit=fstdout, fmt='(a,f10.4)') 'CLUBB-TIMER time_output_multi_col = ', time_output_multi_col
line 4321 : 111 chars : !$acc rcm_in_layer, cloud_cover, wprcp, w_up_in_cloud, w_down_in_cloud, cloudy_updraft_frac, &
line 4333 : 101 chars : !$acc pdf_params%corr_w_thl_2, pdf_params%corr_rt_thl_1, pdf_params%corr_rt_thl_2, &
line 4334 : 103 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 4337 : 104 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 4339 : 103 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
line 4345 : 101 chars : !$acc pdf_params_zm%thl_2, pdf_params_zm%varnce_thl_1, pdf_params_zm%varnce_thl_2, &
line 4346 : 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 4347 : 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 4348 : 115 chars : !$acc pdf_params_zm%alpha_thl, pdf_params_zm%alpha_rt, pdf_params_zm%crt_1, pdf_params_zm%crt_2, &
line 4349 : 110 chars : !$acc pdf_params_zm%cthl_1, pdf_params_zm%cthl_2, pdf_params_zm%chi_1, pdf_params_zm%chi_2, &
line 4350 : 105 chars : !$acc pdf_params_zm%stdev_chi_1, pdf_params_zm%stdev_chi_2, pdf_params_zm%stdev_eta_1, &
line 4351 : 113 chars : !$acc pdf_params_zm%stdev_eta_2, pdf_params_zm%covar_chi_eta_1, pdf_params_zm%covar_chi_eta_2, &
line 4352 : 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 4353 : 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 4354 : 110 chars : !$acc pdf_params_zm%rsatl_1, pdf_params_zm%rsatl_2, pdf_params_zm%rc_1, pdf_params_zm%rc_2, &
line 4383 : 119 chars : if ( thlm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( thlm_sponge_damp_profile )
line 4384 : 118 chars : if ( rtm_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( rtm_sponge_damp_profile )
line 4385 : 117 chars : if ( uv_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( uv_sponge_damp_profile )
line 4386 : 118 chars : if ( wp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp2_sponge_damp_profile )
line 4387 : 118 chars : if ( wp3_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( wp3_sponge_damp_profile )
line 4388 : 122 chars : if ( up2_vp2_sponge_damp_settings%l_sponge_damping ) call finalize_tau_sponge_damp_api( up2_vp2_sponge_damp_profile )
line 5071 : 104 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 5079 : 106 chars : clubb_config_flags%l_add_dycore_grid, clubb_config_flags%grid_adapt_in_time_method )
line 5093 : 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/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_tuner.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 654 : 102 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/clubb_driver_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 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"
----------------------------------------------------------------------------------------------------
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/advance_xm_wpxp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 793 : 101 chars : call wpxp_term_pr1_lhs( nzm, ngrdcol, gr, C6rt_Skw_fnc, C6thl_Skw_fnc, C7_Skw_fnc, & ! Intent(in)
line 794 : 101 chars : invrs_tau_C6_zm, l_scalar_calc, & ! Intent(in)
line 795 : 102 chars : lhs_pr1_wprtp, lhs_pr1_wpthlp, lhs_pr1_wpsclrp ) ! Intent(out)
line 874 : 101 chars : stats_zt, stats_zm, stats_sfc, & ! InOut
line 875 : 101 chars : rtm, wprtp, thlm, wpthlp, sclrm, wpsclrp, & ! InOut
line 876 : 101 chars : err_info ) ! InOut
----------------------------------------------------------------------------------------------------
--------------------- 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/grid_adaptation_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 295 : 128 chars : write(fstderr, *) "The fixed minimum grid density momentum heights need to cover at least the domain of the physics grid."
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/pdf_parameter_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 295 : 101 chars : !$acc pdf_params%alpha_thl, pdf_params%alpha_rt, pdf_params%crt_1, pdf_params%crt_2, &
line 298 : 102 chars : !$acc pdf_params%stdev_eta_2, pdf_params%covar_chi_eta_1, pdf_params%covar_chi_eta_2, &
line 300 : 101 chars : !$acc pdf_params%corr_w_eta_2, pdf_params%corr_chi_eta_1, pdf_params%corr_chi_eta_2, &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/clubb_api_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 2540 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
line 2607 : 107 chars : ! calc_derrived_params_api - Adjusts values of background eddy diffusivity based on vertical grid spacing
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/stats_clubb_utilities.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 3031 : 103 chars : call stat_update_var( stats_metadata%icovar_chi_eta_1, pdf_params%covar_chi_eta_1(icol,:), & ! In
line 3033 : 103 chars : call stat_update_var( stats_metadata%icovar_chi_eta_2, pdf_params%covar_chi_eta_2(icol,:), & ! In
line 3043 : 101 chars : call stat_update_var( stats_metadata%icorr_chi_eta_1, pdf_params%corr_chi_eta_1(icol,:), & ! In
line 3045 : 101 chars : call stat_update_var( stats_metadata%icorr_chi_eta_2, pdf_params%corr_chi_eta_2(icol,:), & ! In
line 3106 : 101 chars : chip2 = compute_variance_binormal( chi, pdf_params%chi_1(icol,:), pdf_params%chi_2(icol,:), &
line 3107 : 106 chars : pdf_params%stdev_chi_1(icol,:), pdf_params%stdev_chi_2(icol,:), &
line 3229 : 101 chars : call stat_update_var( stats_metadata%ivarnce_w_1_zm, pdf_params_zm%varnce_w_1(icol,:), & ! In
line 3231 : 101 chars : call stat_update_var( stats_metadata%ivarnce_w_2_zm, pdf_params_zm%varnce_w_2(icol,:), & ! In
line 4130 : 108 chars : 'floor stuff= ', floor( stats_metadata%stats_tout / stats_metadata%stats_tsamp )
----------------------------------------------------------------------------------------------------
--------------------- 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/parameters_tunable.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 839 : 106 chars : / sqrt( 4.0_core_rknd * ( 1.0_core_rknd - 0.4_core_rknd )**3 &
line 840 : 118 chars : + clubb_params(ngrdcol,iSkw_max_mag)**2 ) ) ) ! Known magic number
line 902 : 117 chars : mult_factor_zt(i) = 1.0_core_rknd + clubb_params(i,imult_coef) * log( avg_deltaz(i) / grid_spacing_thresh )
line 917 : 104 chars : nu_vert_res_dep%nu10(i) = clubb_params(i,inu10) * mult_factor_zt(i) !We're unsure of the grid
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/mono_flux_limiter.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 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/model_flags.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 315 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 541 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
line 845 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/err_info_type_module.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 4: module err_info_type_module
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/err_info_type_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 265 : 102 chars : write(err_info%err_header_global,'(A38, I3, A15, I3, A, A28, A33, F6.2, A1, F6.2, A, A39, A40)') &
line 276 : 105 chars : write(err_info%err_header(i),'(A38, I3, A15, I3, A, A22, I2, A11, F6.2, A13, F6.2, A, A39, A40)') &
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/penta_bicgstab_solver.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 419 : 110 chars : ! It's interesting to see how close E_inv is to the lhs diagonal, it's usually (visually) very close, but
line 420 : 104 chars : ! using it instead of E_inv requires so many more iterations (use it would be jacobi preconditioner)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/CLUBB_core/fill_holes.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 234 : 113 chars : write(*,'(A6,I5,1X,I5,A4,E30.20,A3,E30.20)') "field(", i, k, ") = ", field(i,k), " < ", threshold
line 240 : 105 chars : write(fstderr, *) "column", i, " field average = ", sum(field(i,:) * rho_ds(i,:) * dz(i,:)) &
line 584 : 102 chars : ! rho_ds_dz(i,k-num_hf_draw_points:k+num_hf_draw_points))
line 823 : 110 chars : * min( threshold, field(i,k_start:k_end:grid_dir_indx) ) ) * invrs_denom
line 828 : 106 chars : + mass_fraction * ( min( threshold, field(i,k_start:k_end:grid_dir_indx) ) - threshold )
line 833 : 110 chars : * max( threshold, field(i,k_start:k_end:grid_dir_indx) ) ) * invrs_denom
line 840 : 108 chars : + mass_fraction * ( max( threshold, field(i,k_start:k_end:grid_dir_indx) ) - threshold )
----------------------------------------------------------------------------------------------------
--------------------- 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 2485 : 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/advance_xp2_xpyp_module.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 1394 : 101 chars : sclrprtp(:,:,sclr) ) ! Intent(inout)
line 1429 : 102 chars : sclrpthlp(:,:,sclr) ) ! Intent(inout)
----------------------------------------------------------------------------------------------------
--------------------- 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 1748 : 105 chars : status = nf90_def_var( ncid_zm, "columns", NF90_INT, (/column_id_zm/), column_var_id_zm )
line 1749 : 104 chars : status = nf90_def_var( ncid_zm, "altitude", NF90_PREC, (/vertical_id_zm/), vertical_var_id_zm )
line 1750 : 104 chars : status = nf90_def_var( ncid_zm, "time", NF90_PREC, (/time_id_zm/), time_var_id_zm )
line 1802 : 105 chars : status = nf90_def_var( ncid_zt, "columns", NF90_INT, (/column_id_zt/), column_var_id_zt )
line 1803 : 104 chars : status = nf90_def_var( ncid_zt, "altitude", NF90_PREC, (/vertical_id_zt/), vertical_var_id_zt )
line 1804 : 104 chars : status = nf90_def_var( ncid_zt, "time", NF90_PREC, (/time_id_zt/), time_var_id_zt )
line 1850 : 101 chars : status = nf90_put_var( ncid_zm, column_var_id_zm, column_list, start=(/1/), count=(/ngrdcol/) )
line 1851 : 101 chars : status = nf90_put_var( ncid_zt, column_var_id_zt, column_list, start=(/1/), count=(/ngrdcol/) )
line 1853 : 101 chars : !=================================== Allocate Output Arrays ===================================
line 1980 : 116 chars : !$acc upwp_out, vpwp_out, up2_out, vp2_out, wp3_out, rcm_out, cloud_frac_out, rtm_out, thlm_out )
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 17 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/mpace_a.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 272 : 103 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 273 : 102 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/prescribe_forcings.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 1: module prescribe_forcings_module
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/prescribe_forcings.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 303 : 104 chars : !$acc enter data create( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 477 : 118 chars : !$acc update device( wm_zt, wm_zm, thlm_forcing, rtm_forcing, um_ref, vm_ref, sclrm_forcing, edsclrm_forcing )
line 908 : 105 chars : call stat_update_var_pt( stats_metadata%ish, 1, wpthlp_sfc(i)*rho_zm(i,1)*Cp,& ! intent(in)
line 911 : 105 chars : call stat_update_var_pt( stats_metadata%ilh, 1, wprtp_sfc(i)*rho_zm(i,1)*Lv, & ! intent(in)
line 914 : 103 chars : call stat_update_var_pt( stats_metadata%iwpthlp_sfc, 1, wpthlp_sfc(i), & ! intent(in)
line 917 : 103 chars : call stat_update_var_pt( stats_metadata%iwprtp_sfc, 1, wprtp_sfc(i), & ! intent(in)
line 920 : 103 chars : call stat_update_var_pt( stats_metadata%iupwp_sfc, 1, upwp_sfc(i), & ! intent(in)
line 923 : 103 chars : call stat_update_var_pt( stats_metadata%ivpwp_sfc, 1, vpwp_sfc(i), & ! intent(in)
line 926 : 102 chars : call stat_update_var_pt( stats_metadata%iustar, 1, ustar(i), & ! intent(in)
line 929 : 103 chars : call stat_update_var_pt( stats_metadata%iT_sfc, 1, T_sfc(i), & ! intent(in)
line 935 : 103 chars : !$acc exit data delete( um_bot, vm_bot, rtm_bot, thlm_bot, rho_bot, exner_bot, z_bot, ustar, ubar )
line 1138 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1142 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1146 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1147 : 112 chars : exner_zm(i,km1), exner_zm(i,k00), exner_zm(i,kp1), exner_zm(i,kp2) )
line 1150 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1151 : 108 chars : thlm_zm(i,km1), thlm_zm(i,k00), thlm_zm(i,kp1), thlm_zm(i,kp2) )
line 1154 : 101 chars : gr%zm(i,km1), gr%zm(i,k00), gr%zm(i,kp1), gr%zm(i,kp2), &
line 1155 : 104 chars : rtm_zm(i,km1), rtm_zm(i,k00), rtm_zm(i,kp1), rtm_zm(i,kp2) )
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/rico.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 159 : 103 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 160 : 102 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/dycoms2_rf02.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 121 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 122 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/bomex.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 127 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 128 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/dycoms2_rf01.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 91 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 92 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/gabls2.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 152 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 153 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/mpace_b.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 157 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 158 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/Benchmark_cases/lba.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 96 : 105 chars : if ( sclr_idx%iiedsclr_thl > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_thl) = thlm_forcing(i,k)
line 97 : 104 chars : if ( sclr_idx%iiedsclr_rt > 0 ) edsclrm_forcing(i,k,sclr_idx%iiedsclr_rt) = rtm_forcing(i,k)
line 176 : 101 chars : wpthlp_sfc(i) = convert_sens_ht_to_km_s( ( 270._core_rknd * ft**1.5_core_rknd ), rho_sfc(i) )
line 177 : 101 chars : wprtp_sfc(i) = convert_latent_ht_to_m_s( ( 554._core_rknd * ft**1.3_core_rknd ), rho_sfc(i) )
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 11 check(s) failed.
No matches for /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/KK_microphys/
CLUBBStandardsCheck.py has begun.
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/pdf_parameter_tests.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 542 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/rev_direction_grid_test.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 1: module rev_direction_grid_test
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/rev_direction_grid_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 233 : 113 chars : !seed used = -1102123834 1480475276 -248259584 525496705 1614213535 -1267052688 6290033 -78466652
line 234 : 127 chars : !seed_output = (/ -1102123834, 1480475276, -248259584, 525496705, 1614213535, -1267052688, 6290033, -78466652 /)
line 750 : 101 chars : call xpyp_term_ta_pdf_lhs_godunov( gr_ascending%nzm, gr_ascending%nzt, 1, gr_ascending, & ! In
line 754 : 104 chars : call xpyp_term_ta_pdf_lhs_godunov( gr_descending%nzm, gr_descending%nzt, 1, gr_descending, & ! In
line 795 : 101 chars : call xpyp_term_ta_pdf_rhs_godunov( gr_ascending%nzm, gr_ascending%nzt, 1, gr_ascending, & ! In
line 801 : 104 chars : call xpyp_term_ta_pdf_rhs_godunov( gr_descending%nzm, gr_descending%nzt, 1, gr_descending, & ! In
----------------------------------------------------------------------------------------------------
--------------------- Default Private Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90 ---------------------
WARNING: Missing default 'private' in the following modules:
line 19: module fill_holes_tests
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/fill_holes_tests.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 162 : 102 chars : field_initial(i,k) = 10_core_rknd * ( sin( 8.0 * pi * ( k + i ) / nz ) + one ) &
line 180 : 110 chars : if ( l_print_name ) print *, "filling holes of: easy_fill_test_"//trim(fill_type_names(fill_type))
line 195 : 121 chars : write(*,'(A6,I5,1X,I5,A4,E30.20,A3,E30.20)') "field(", i, k, ") = ", field(i,k), " < ", threshold
line 241 : 118 chars : if ( l_print_name ) print *, "filling holes of: reverse_easy_fill_test_"//trim(fill_type_names(fill_type))
line 344 : 102 chars : field_initial(i,k) = 10_core_rknd * ( sin( 8.0 * pi * ( k + i ) / nz ) + one ) &
line 363 : 113 chars : if ( l_print_name ) print *, "filling holes of: below_thresh_test_"//trim(fill_type_names(fill_type))
line 385 : 103 chars : print *, "-- error", two * abs( initial_mass - new_mass ) / ( initial_mass + new_mass )
line 419 : 121 chars : if ( l_print_name ) print *, "filling holes of: reverse_below_thresh_test_"//trim(fill_type_names(fill_type))
line 430 : 116 chars : print *, "ERROR in below_thresh_test REVERSE MODE: method was not conservative in grid reverse mode"
----------------------------------------------------------------------------------------------------
--------------------- Line Length Check: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/src/G_unit_test_types/spurious_source_test.F90 ---------------------
WARNING: Lines exceed 100 characters (limit 100):
line 508 : 104 chars : l_prescribed_avg_deltaz, & ! used in calc_derrived_params_api. If .true., avg_deltaz = deltaz
----------------------------------------------------------------------------------------------------
CLUBBStandardsCheck.py has finished.
FAIL: 6 check(s) failed.
[91m===============================================================
[91mCLUBBStandardsCheck FAILED
[91m THIS IS PRINTED IN ALL RED, CAPITAL LETTERS, AND USES
[91m AN EXCLAMATION MARK TO ENSURE THE DEVELOPERS FEEL SHAME!
[91m IF YOU ARE ONE OF THESE "DEVELOPERS" CHECK THE
[91m LOG FILE FOR DETAILS: /home/jenkins/workspace/clubb_G_unit_gfortran_tests/build/gfortran_DEBUG_PRECdouble/cmake_build_output.txt
[91m===============================================================[0m
[93mBuild completed successfully, but some source code checks have failed.[0m
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Run CLUBB's G Unit Tests)
[Pipeline] sh
+ source /etc/profile.d/larson-group.sh
++ export GIT_EDITOR=vi
++ GIT_EDITOR=vi
++ export SVN_EDITOR=vi
++ SVN_EDITOR=vi
++ export OMP_STACKSIZE=1048579
++ OMP_STACKSIZE=1048579
++ export LMOD_ROOT=/opt/lmod/
++ LMOD_ROOT=/opt/lmod/
++ source /opt/lmod//lmod/lmod/init/bash
+++ '[' -z '' ']'
+++ case "$-" in
+++ __lmod_vx=x
+++ '[' -n x ']'
+++ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for this output (/usr/local/spack/opt/spack/linux-pop22-cascadelake/gcc-12.2.0/lmod-8.7.37-bi3kyxcdrfgw3y7vv2k7c5rjxg75qzju/lmod/lmod/init/bash)
Shell debugging restarted
+++ unset __lmod_vx
+++ find /usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core -print -quit
++ export MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
++ MODULEPATH=/usr/local/spack/share/spack/lmod/linux-pop22-x86_64/Core
+ module load gcc netcdf-fortran
+ '[' -z '' ']'
+ case "$-" in
+ __lmod_sh_dbg=x
+ '[' -n x ']'
+ set +x
Shell debugging temporarily silenced: export LMOD_SH_DBG_ON=1 for Lmod's output
Lmod Warning:
-------------------------------------------------------------------------------
The following dependent module(s) are not currently loaded:
gcc-runtime/13.1.0-gzedt3f (required by: mpich/4.2.1, hdf5/1.14.3-42x3zrv,
netcdf-c/4.9.2, netcdf-fortran/4.5.3)
-------------------------------------------------------------------------------
Shell debugging restarted
+ unset __lmod_sh_dbg
+ return 0
+ python run_scripts/run_G_unit_tests.py
Running G_unit_tests
Performing KK integrals tests
Percent Difference Tolerance (maximum acceptable percent difference between the CLUBB results and the MATLAB results using the same equations): 1.50000E-12
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 1
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: 3.940566852245362E-04
MATLAB result: 3.940566852245361E-04
Percent difference: 1.37569E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: 3.927115990964176E-04
MATLAB result: 3.927115990964175E-04
Percent difference: 2.76081E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: 3.959157231937519E-04
MATLAB result: 3.959157231937519E-04
Percent difference: 1.36923E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: 3.961224408774508E-04
MATLAB result: 3.961224408774508E-04
Percent difference: 1.36852E-14%
Agreement with regards to integral: Trivariate NNL
Obtained result: 5.957899876478518E-10
MATLAB result: 5.957899876478559E-10
Percent difference: 6.94188E-13%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 5.299646183560617E-10
MATLAB result: 5.299646183560650E-10
Percent difference: 6.24329E-13%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.543555152546510E-10
MATLAB result: 5.543555152546526E-10
Percent difference: 2.98430E-13%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.286452253001267E-10
MATLAB result: 5.286452253001284E-10
Percent difference: 3.12944E-13%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -9.854231981928352E-04
MATLAB result: -9.854231981928350E-04
Percent difference: 2.20048E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549015E-04
MATLAB result: -9.922448817549015E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 8.940070763287877E-09
MATLAB result: 8.940070763287870E-09
Percent difference: 7.40201E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 8.942709549399747E-09
MATLAB result: 8.942709549399743E-09
Percent difference: 3.69991E-14%
Agreement with regards to integral: Bivariate LL
Obtained result: 1.007487885941809E-02
MATLAB result: 1.007487885941809E-02
Percent difference: 1.72183E-14%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 2
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -4.441972023148827E-05
MATLAB result: -4.441972023148829E-05
Percent difference: 4.57652E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -4.559373570842309E-05
MATLAB result: -4.559373570842311E-05
Percent difference: 4.45868E-14%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -5.040842768062482E-05
MATLAB result: -5.040842768062482E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -5.038775591225493E-05
MATLAB result: -5.038775591225493E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.629614112194526E-09
MATLAB result: 4.629614112194527E-09
Percent difference: 1.78672E-14%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.604609365066508E-09
MATLAB result: 4.604609365066509E-09
Percent difference: 1.79642E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 4.685272472815394E-09
MATLAB result: 4.685272472815394E-09
Percent difference: 1.76549E-14%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 4.683452323676534E-09
MATLAB result: 4.683452323676534E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -1.088125285831538E-04
MATLAB result: -1.088125285831538E-04
Percent difference: 2.49099E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: -9.922448817549016E-05
MATLAB result: -9.922448817549016E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 7.907812698669825E-10
MATLAB result: 7.907812698669821E-10
Percent difference: 5.23015E-14%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 6.330953526469324E-10
MATLAB result: 6.330953526469324E-10
Percent difference: 0.00000E+00%
==================================================
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals): 3
Agreement with regards to integral: Quadrivariate NNLL
Obtained result: -7.622010070476262E-05
MATLAB result: -7.622010070476262E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1)
Obtained result: -7.695119321547280E-05
MATLAB result: -7.695119321547280E-05
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Quadrivariate NNLL (const. x1 and x2)
Obtained result: -1.000000000000000E-04
MATLAB result: -1.000000000000000E-04
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL
Obtained result: 4.851589269710427E-09
MATLAB result: 4.851589269710427E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1)
Obtained result: 4.837259568295437E-09
MATLAB result: 4.837259568295437E-09
Percent difference: 1.71002E-14%
Agreement with regards to integral: Trivariate NNL (const. x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
Agreement with regards to integral: Trivariate NNL (const. x1 and x2)
Obtained result: 5.000000000000000E-09
MATLAB result: 5.000000000000000E-09
Percent difference: 0.00000E+00%
--------------------------------------------------
Agreement with regards to integral: Trivariate NLL
Obtained result: -4.609761356905442E-05
MATLAB result: -4.609761356905441E-05
Percent difference: 1.46998E-14%
Agreement with regards to integral: Trivariate NLL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL
Obtained result: 3.254808634091251E-10
MATLAB result: 3.254808634091251E-10
Percent difference: 0.00000E+00%
Agreement with regards to integral: Bivariate NL (const. x1)
Obtained result: 0.000000000000000E+00
MATLAB result: 0.000000000000000E+00
Percent difference: 0.00000E+00%
==================================================
Success!
==================================================
Performing correlation cholesky matrix setup tests
correlation matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
--
correlation matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 1.0000000000000000 0.86602540378443860
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
--
==================================================
Performing correlation matrix approximation tests
correlation cholesky matrix :
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
Agreement
Result:
1.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 1.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 1.0000000000000000 1.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
Agreement
Result:
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
--
correlation cholesky matrix :
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.0000000000000000 0.50000000000000011 0.43301270189221941
0.0000000000000000 0.0000000000000000 0.25000000000000011
Agreement
Result:
1.0000000000000000 0.86602540378443860 0.86602540378443860
0.86602540378443860 1.0000000000000000 0.96650635094610959
0.86602540378443860 0.96650635094610959 1.0000000000000000
--
==================================================
Success!
==================================================
Performing hole_filling_hm_one_lev_tests
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
==================================================
Performing fill_holes_hydromet_tests
testset_in(:, 1 ) = 1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 6.0000000000000000
testset_in(:, 2 ) = 6.0000000000000000 5.0000000000000000 4.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000
testset_in(:, 3 ) = 3.0000000000000000 3.0000000000000000 3.0000000000000000 4.0000000000000000 4.0000000000000000 4.0000000000000000
testset_in(:, 4 ) = 1.0000000000000000 3.0000000000000000 2.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 5 ) = 1.0000000000000000 -2.0000000000000000 -4.0000000000000000 -3.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 6 ) = 2.0000000000000000 -1.0000000000000000 -1.0000000000000000 1.0000000000000000 -1.0000000000000000 0.0000000000000000
testset_in(:, 7 ) = 1.0000000000000000 4.0000000000000000 2.0000000000000000 3.0000000000000000 -1.0000000000000000 0.0000000000000000
Hole-filling was conservative
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
Hole-filling was conservative
Hole not entirely filled, as expected
Hole-filling was conservative
==================================================
Success!
Nc-Ncn equations unit test
==========================
Iteration = 1
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.15865525393145707
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 4.7896918412674905E-002
Nc_in_cloud (from Ncn) = 66978771.649335034
Ncm (from Ncn) = 3208076.7610693830
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 1 is a success!
======================================================
Iteration = 2
mu_chi_1 = -1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.0000000000000000
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 0.0000000000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 2 is a success!
======================================================
Iteration = 3
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.40000000000000002
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 0.99957093966680322
cloud fraction = 0.95210308158732504
Nc_in_cloud (from Ncn) = 51650441.529586732
Ncm (from Ncn) = 49176544.545665480
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 3 is a success!
======================================================
Iteration = 4
mu_chi_1 = 1.0000000000000000E-008
mu_chi_2 = -1.0000000000000000E-008
sigma_chi_1 = 5.0000000000000002E-005
sigma_chi_2 = 5.0000000000000002E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.10000000000000001
mixture fraction = 0.50000000000000000
cloud fraction (1st PDF component) = 0.50007978845554835
cloud fraction (2nd PDF component) = 0.49992021154445165
cloud fraction = 0.50000000000000000
Nc_in_cloud (from Ncn) = 53982783.648312390
Ncm (from Ncn) = 26991391.824156195
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 4 is a success!
======================================================
Iteration = 5
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = 5.0000000000000002E-005
sigma_chi_1 = 9.9999999999999995E-007
sigma_chi_2 = 1.5000000000000000E-006
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 25000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.50000000000000000
mixture fraction = 0.29999999999999999
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 1.0000000000000000
cloud fraction = 1.0000000000000000
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 50000000.000000000
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 5 is a success!
======================================================
Iteration = 6
mu_chi_1 = -9.9999999999999995E-008
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 50000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.0000000000000000
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.49601064368536840
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 9.9545377003631155E-002
Nc_in_cloud (from Ncn) = 49999999.999999993
Ncm (from Ncn) = 4977268.8501815572
Ncnm (from Nc_in_cloud) = 49999999.999999993
Ncnm (from Ncm) = 49999999.999999993
Test iteration 6 is a success!
======================================================
Iteration = 7
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 0.0000000000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.48670438618290790
cloud fraction = 0.49732907640173213
Nc_in_cloud (from Ncn) = 50000000.000000000
Ncm (from Ncn) = 24866453.820086606
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 7 is a success!
======================================================
Iteration = 8
mu_chi_1 = 9.9999999999999995E-007
mu_chi_2 = -4.9999999999999998E-007
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-009
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.53982783727702899
cloud fraction (2nd PDF component) = 0.0000000000000000
cloud fraction = 0.10796556745540581
Nc_in_cloud (from Ncn) = 81037694.237492040
Ncm (from Ncn) = 8749280.6436284967
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 8 is a success!
======================================================
Iteration = 9
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-009
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 1.0000000000000000
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.20034324826655747
Nc_in_cloud (from Ncn) = 52151569.416871324
Ncm (from Ncn) = 10448214.819174858
Ncnm (from Nc_in_cloud) = 50000000.000000000
Ncnm (from Ncm) = 50000000.000000000
Test iteration 9 is a success!
======================================================
Iteration = 10
mu_chi_1 = 1.0000000000000001E-005
mu_chi_2 = -5.0000000000000002E-005
sigma_chi_1 = 1.0000000000000001E-005
sigma_chi_2 = 1.5000000000000000E-005
mu_Ncn_1 = mu_Ncn_2 = 50000000.000000000
sigma_Ncn_1 = sigma_Ncn_2 = 75000000.000000000
corr_chi_Ncn_1 = corr_chi_Ncn_2 = 0.69999999999999996
mixture fraction = 0.20000000000000001
cloud fraction (1st PDF component) = 0.84134474606854293
cloud fraction (2nd PDF component) = 4.2906033319683735E-004
cloud fraction = 0.16861219748026607
Nc_in_cloud (from Ncn) = 60768988.354576342
Ncm (from Ncn) = 10246392.665117815
Ncnm (from Nc_in_cloud) = 50000000.000000007
Ncnm (from Ncm) = 50000000.000000007
Test iteration 10 is a success!
======================================================
Success!
Performing corr_mtx_test
==================================================
Expected:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Actual:
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 1.00 1.10 1.20
0.00 1.00 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 0.10 0.11
0.00 0.00 1.00 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.01
0.00 0.00 0.00 1.00 -1.00 -2.00 -3.00 -4.00 -5.00 -6.00 -7.00 -8.00
0.00 0.00 0.00 0.00 1.00 -0.10 -0.20 -0.30 -0.40 -0.50 -0.60 -0.70
0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 1.00 -0.00 -0.00 -0.00 -0.00 -0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 2.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 3.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 4.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
Success!
==================================================
==================================================
Performing SILHS category test
==================================================
Success!
Performing hydrometeor mu and sigma values unit test
====================================================
PDF parameter set 1:
hmm = 2.0000000000000002E-005
mixt_frac = 2.7000000000000000E-002
precip_frac = 8.0000000000000002E-002
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 5.4470709146968145E-002
hmp2_ip_on_hmm2_ip = 3.6000000000000001
hmm_ip = 2.5000000000000001E-004
hmp2 = 2.2599999999999997E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 2:
hmm = 5.0000000000000004E-006
mixt_frac = 0.50000000000000000
precip_frac = 0.50000000000000000
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 0.50000000000000000
hmp2_ip_on_hmm2_ip = 1.0000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 7.5000000000000012E-011
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 3:
hmm = 1.0000000000000000E-010
mixt_frac = 0.25000000000000000
precip_frac = 0.40000000000000002
precip_frac_1 = 0.75000000000000000
precip_frac_2 = 0.28333333333333338
hmp2_ip_on_hmm2_ip = 0.50000000000000000
hmm_ip = 2.5000000000000002E-010
hmp2 = 2.7500000000000010E-020
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 4:
hmm = 4.0000000000000003E-005
mixt_frac = 1.0000000000000000E-002
precip_frac = 0.10000000000000001
precip_frac_1 = 1.0000000000000000E-002
precip_frac_2 = 0.10090909090909091
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 5:
hmm = 4.0000000000000003E-005
mixt_frac = 0.98999999999999999
precip_frac = 0.10000000000000001
precip_frac_1 = 0.10090909090909091
precip_frac_2 = 1.0000000000000278E-002
hmp2_ip_on_hmm2_ip = 2.5000000000000000
hmm_ip = 4.0000000000000002E-004
hmp2 = 5.4400000000000004E-008
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 6:
hmm = 1.0000000000000001E-005
mixt_frac = 0.12500000000000000
precip_frac = 1.0000000000000000
precip_frac_1 = 1.0000000000000000
precip_frac_2 = 1.0000000000000000
hmp2_ip_on_hmm2_ip = 1.5000000000000000
hmm_ip = 1.0000000000000001E-005
hmp2 = 1.5000000000000005E-010
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 7:
hmm = 6.0000000000000002E-005
mixt_frac = 5.0000000000000003E-002
precip_frac = 0.11000000000000000
precip_frac_1 = 0.90000000000000002
precip_frac_2 = 6.8421052631578952E-002
hmp2_ip_on_hmm2_ip = 5.0000000000000000
hmm_ip = 5.4545454545454548E-004
hmp2 = 1.9276363636363637E-007
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 8:
hmm = 7.9999999999999996E-007
mixt_frac = 0.80000000000000004
precip_frac = 0.25000000000000000
precip_frac_1 = 0.10000000000000001
precip_frac_2 = 0.85000000000000009
hmp2_ip_on_hmm2_ip = 0.75000000000000000
hmm_ip = 3.1999999999999999E-006
hmp2 = 3.8399999999999998E-012
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 9:
hmm = 1.3499999999999999E-005
mixt_frac = 0.10000000000000001
precip_frac = 0.10000000000000001
precip_frac_1 = 0.50000000000000000
precip_frac_2 = 5.5555555555555559E-002
hmp2_ip_on_hmm2_ip = 2.0000000000000000
hmm_ip = 1.3499999999999997E-004
hmp2 = 5.2852499999999985E-009
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
PDF parameter set 10 (randomly generated):
Random seed values = 2026 1 9 -360 3 11 49 778
hmm = 2.4553947299393127E-005
mixt_frac = 0.12496832104134623
precip_frac = 0.74387234359252541
precip_frac_1 = 0.38185109534289252
precip_frac_2 = 0.79557468610530058
hmp2_ip_on_hmm2_ip = 4.8191828596870154
hmm_ip = 3.3008280938111018E-005
hmp2 = 4.1134558904837485E-009
hm_tol = 1.0000000000000000E-010
precip_frac_tol = 1.0000000000000000E-002
Running tests for the above parameter set for all combinations of omicron and zeta. Omicron values are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Checking values with precipitation in 1st PDF component only.
Checking values with precipitation in 2nd PDF component only.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for ADG1
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 1 9 -360 3 11 49 779
wm = -2.8921054012137515E-002
rtm = 3.3462922661499296E-003
thlm = 304.14586551398270
wp2 = 0.37560716697209584
rtp2 = 2.6632079175432636E-006
thlp2 = 0.37691297220197240
Skw = 2.0571237211799849
wp3 = 0.47354495905731703
Skrt = -2.9766950162043662
rtp3 = -1.2937249395249737E-008
Skthl = -3.7229346514925865
thlp3 = -0.86148365179775666
wprtp = 9.1421934468322139E-004
wpthlp = -2.6241915768951082E-003
rtpthlp = -9.0170730592159238E-004
Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w). Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for LY93
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
PDF parameter set 10 (randomly generated):
Random seed values = 2026 1 9 -360 3 11 49 780
wm = -2.8921054012137515E-002
rtm = 3.3462922661499118E-003
thlm = 304.14586551398270
wp2 = 0.37560716676745731
rtp2 = 2.5094312957645965E-006
thlp2 = 0.37689681694308474
Skw = 1.9581984893995479
wp3 = 0.45077260700639332
Skrt = -2.2738114139106704
rtp3 = -9.0389379904481823E-009
Skthl = 1.1306566978871890
thlp3 = 0.26161609824349213
wprtp = -4.7698084118967627E-005
wpthlp = -0.24592837307315576
rtpthlp = 8.7526870208293184E-004
Running tests for the above parameter set for the setting variable (w is used here).
Running tests for the above parameter set for the full PDF.
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for TSDADG
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
PDF parameter set 10 (randomly generated):
Random seed values = 2026 1 9 -360 3 11 49 780
wm = -2.8921054012137515E-002
rtm = 3.3462922661499118E-003
thlm = 304.14586551398270
wp2 = 0.37560716676745731
rtp2 = 2.5094312957645965E-006
thlp2 = 0.37689681694308474
Skw = 1.9581984893995479
wp3 = 0.45077260700639332
Skrt = -2.2738114139106704
rtp3 = -9.0389379904481823E-009
Skthl = 1.1306566978871890
thlp3 = 0.26161609824349213
wprtp = -4.7698084118967627E-005
wpthlp = -0.24592837307315576
rtpthlp = 8.7526870208293184E-004
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 1 9 -360 3 11 49 784
wm = -2.8921054012137515E-002
rtm = 3.3462922661499851E-003
thlm = 304.14586551398253
wp2 = 0.37560716660374205
rtp2 = 5.7116239821769308E-007
thlp2 = 0.37683118909689095
Skw = 3.1437563213304536
wp3 = 0.72368518283709826
Skrt = 0.53772299821087532
rtp3 = 2.3211223006880291E-010
Skthl = 0.29559393058742778
thlp3 = 6.8377902338408611E-002
wprtp = -4.2817456058446602E-004
wpthlp = 0.27755746309631535
rtpthlp = -4.1753773794367643E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing PDF parameter values unit test
=========================================
Performing PDF parameter unit tests for the new hybrid PDF
PDF parameter set 1:
wm = -5.0000000000000001E-003
rtm = 4.2868400000000001E-003
thlm = 309.13099999999997
wp2 = 3.3741899999999998E-002
rtp2 = 1.1786499999999999E-006
thlp2 = 5.9920000000000001E-002
Skw = 8.6882480819056749
wp3 = 5.3850099999999998E-002
Skrt = 1.4880899999999999
rtp3 = 1.9041733779221134E-009
Skthl = -1.4580900000000001
thlp3 = -2.1386614365822010E-002
wprtp = 2.3725400000000002E-005
wpthlp = -5.5512699999999996E-003
rtpthlp = -2.3917778634313013E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 2:
wm = 1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.10000000000000001
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = -4.5000000000000000
wp3 = -0.14230249470757708
Skrt = 1.0000000000000000
rtp3 = 1.1180339887498949E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = -5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 3:
wm = -1.0000000000000000E-002
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.50000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.10000000000000001
Skw = 4.5000000000000000
wp3 = 1.5909902576697321
Skrt = 2.5000000000000000
rtp3 = 2.7950849718747374E-008
Skthl = -2.5000000000000000
thlp3 = -7.9056941504209499E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -6.3639610306789288E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 4:
wm = 0.10000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 5.0000000000000000
rtp2 = 6.9999999999999999E-006
thlp2 = 1.0000000000000000
Skw = 2.0000000000000000
wp3 = 22.360679774997898
Skrt = 1.5000000000000000
rtp3 = 2.7780388766178202E-008
Skthl = -1.5000000000000000
thlp3 = -1.5000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -1.0000000000000000E-002
rtpthlp = -2.3811761799581318E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 5:
wm = -1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.29999999999999999
rtp2 = 1.9999999999999999E-006
thlp2 = 0.20000000000000001
Skw = -2.0000000000000000
wp3 = -0.32863353450309962
Skrt = -0.25000000000000000
rtp3 = -7.0710678118654745E-010
Skthl = -3.5000000000000000
thlp3 = -0.31304951684997062
wprtp = 5.0000000000000002E-005
wpthlp = 5.0000000000000001E-003
rtpthlp = 5.6920997883030825E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 6:
wm = -1.0000000000000000E-002
rtm = 5.0000000000000001E-003
thlm = 305.00000000000000
wp2 = 1.0000000000000000
rtp2 = 4.9999999999999998E-007
thlp2 = 0.10000000000000001
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -2.0124611797498108E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 7:
wm = 0.20000000000000001
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.25000000000000000
rtp2 = 5.0000000000000004E-006
thlp2 = 0.14999999999999999
Skw = 0.75000000000000000
wp3 = 9.3750000000000000E-002
Skrt = 3.0000000000000000
rtp3 = 3.3541019662496846E-008
Skthl = -0.50000000000000000
thlp3 = -2.9047375096555625E-002
wprtp = 5.0000000000000002E-005
wpthlp = -5.0000000000000001E-003
rtpthlp = -7.7942286340599491E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 8:
wm = -2.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 0.0000000000000000
thlp2 = 0.0000000000000000
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = 0.0000000000000000
rtp3 = 0.0000000000000000
Skthl = 0.0000000000000000
thlp3 = 0.0000000000000000
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 0.0000000000000000
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 9:
wm = 1.0000000000000000E-003
rtm = 1.0000000000000000E-002
thlm = 305.00000000000000
wp2 = 0.0000000000000000
rtp2 = 1.9999999999999999E-006
thlp2 = 0.69999999999999996
Skw = 0.0000000000000000
wp3 = 0.0000000000000000
Skrt = -1.0000000000000000
rtp3 = -2.8284271247461898E-009
Skthl = -2.0000000000000000
thlp3 = -1.1713240371477056
wprtp = 0.0000000000000000
wpthlp = 0.0000000000000000
rtpthlp = 1.0648943609579309E-003
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
PDF parameter set 10 (randomly generated):
Random seed values = 2026 1 9 -360 3 11 49 786
wm = -2.8921054012137515E-002
rtm = 3.3462922661499743E-003
thlm = 304.14586551398253
wp2 = 0.37560716668559724
rtp2 = 6.3267304702080520E-007
thlp2 = 0.37679887661460276
Skw = 2.9459058562918887
wp3 = 0.67814047942053679
Skrt = -1.4309304415278978
rtp3 = -7.2009123162289739E-010
Skthl = 3.7541386822054772
thlp3 = 0.86830981654422656
wprtp = 1.2287517829601100E-004
wpthlp = -0.16402460445366929
rtpthlp = -4.3942724043832273E-004
Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable. F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1. Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
Success!
Performing spurious source unit test
====================================
Test 1
Vertical integral of <rt> before advance_xm_wpxp = 28.316670259590918
Vertical integral of <rt> after advance_xm_wpxp = 28.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 = 16.333839668940460
Vertical integral of <rt> after advance_xm_wpxp = 16.425848862299610
Vertical integral of <rt> forcing = -1.9277241289599983E-005
Flux of <rt> at the surface = 3.2597455248679940E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -3.2658202314336804E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 788368.99746638816
Vertical integral of <thl> after advance_xm_wpxp = 788363.77899833233
Vertical integral of <thl> forcing = -2.3914384568770465E-002
Flux of <thl> at the surface = 6.5194910497359884E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -3.9066666568388086E-013
The spurious source of <thl> is within acceptable limits.
Test 3
Vertical integral of <rt> before advance_xm_wpxp = 40.585644382337009
Vertical integral of <rt> after advance_xm_wpxp = 40.819483855011917
Vertical integral of <rt> forcing = -4.3832114709578597E-005
Flux of <rt> at the surface = 8.2329702362598698E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = -4.6173460020726420E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 796420.63966757001
Vertical integral of <thl> after advance_xm_wpxp = 796407.64628708782
Vertical integral of <thl> forcing = -5.9777208744406000E-002
Flux of <thl> at the surface = 1.6465940472519740E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -2.0875731698843936E-012
The spurious source of <thl> is within acceptable limits.
Test 4
Vertical integral of <rt> before advance_xm_wpxp = 39.607074590293301
Vertical integral of <rt> after advance_xm_wpxp = 39.834583719883319
Vertical integral of <rt> forcing = -5.9083873477513284E-005
Flux of <rt> at the surface = 8.1744763877753673E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 3.7357541105703662E-017
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 797780.91520164674
Vertical integral of <thl> after advance_xm_wpxp = 797767.93203389691
Vertical integral of <thl> forcing = -5.9626178607822038E-002
Flux of <thl> at the surface = 1.6348952775550735E-002
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -4.9400067370086731E-013
The spurious source of <thl> is within acceptable limits.
Test 5
Vertical integral of <rt> before advance_xm_wpxp = 11.286148392439264
Vertical integral of <rt> after advance_xm_wpxp = 11.326605380918174
Vertical integral of <rt> forcing = -1.3924236109356475E-005
Flux of <rt> at the surface = 1.4878086437238036E-004
Flux of <rt> at the top of the domain = 0.0000000000000000
Spurious source of <rt> = 9.7900068043652033E-018
The spurious source of <rt> is within acceptable limits.
Vertical integral of <thl> before advance_xm_wpxp = 787680.05804352730
Vertical integral of <thl> after advance_xm_wpxp = 787675.41300796927
Vertical integral of <thl> forcing = -1.8459069146791868E-002
Flux of <thl> at the surface = 2.9756172874476072E-003
Flux of <thl> at the top of the domain = 0.0000000000000000
Spurious source of <thl> = -7.5887560124776599E-013
The spurious source of <thl> is within acceptable limits.
Success!
Running tuner tests
-------------------
goldstein_price pass
rastrigin pass
himmelblau pass
eggholder pass
schaffer pass
==================================================
Performing 'mean vertical velocity in clouds' tests
Success!
==================================================
==================================================
Performing 'Smooth Heaviside Peskin' tests
Case 1: heaviside_smth_range = 0, smooth Heaviside collapses to Heaviside
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000
Case 2: heaviside_smth_range = 2.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Expected outcome: 0.0000000000000000 9.0845056908100005E-002 0.90915494309000000 1.0000000000000000
True outcome: 0.0000000000000000 9.0845056909999999E-002 0.90915494309000000 1.0000000000000000
Success!
==================================================
==================================================
Performing 'Smooth Min / Max' tests
Input: -1.0000000000000000 -1.0000000000000000E-010 1.0000000000000000E-010 1.0000000000000000
(Simple) min:
Expected outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
True outcome: -1.0000000000000000 -1.0000000000000000E-010 0.0000000000000000 0.0000000000000000
(Simple) max:
Expected outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
True outcome: 0.0000000000000000 0.0000000000000000 1.0000000000000000E-010 1.0000000000000000
Input: -2.0000000001000000 -1.0000000000000000 1.0000000000000000 2.0000000001000000
Smooth min:
Expected outcome: -2.2247448714824141 -1.3660254037844388 -0.36602540378443871 -0.22474487138241406
True outcome: -2.2247448714824136 -1.3660254037844388 -0.36602540378443871 -0.22474487138241384
Smooth max:
Expected outcome: 0.22474487138241406 0.36602540378443871 1.3660254037844388 2.2247448714824141
True outcome: 0.22474487138241384 0.36602540378443871 1.3660254037844388 2.2247448714824136
Testing smooth min and max with min_max_smth_mag=1e-7, input_var1 = equidistant grid on [-1, 1), input_var2 = 0
Input: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
True outcome: -1.1180339887498949 -1.0226805085936306 -1.0024751868658441 -1.0002497501874359
Smooth max:
Expected outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
True outcome: 1.1180339887498949 1.0226805085936306 1.0024751868658441 1.0002497501874359
Input_var1: -1000000.0000000000 -10000000.000000000 -100000000.00000000 -1000000000.0000000
Input_var2: 1.0000000000000000 10.000000000000000 100.00000000000000 1000.0000000000000
Smooth min:
Expected outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
True outcome: -1000000.0000002501 -10000000.000000026 -100000000.00000000 -1000000000.0000000
Smooth max:
Expected outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
True outcome: 1.0000002500019036 10.000000025145710 100.00000000000000 1000.0000000000000
Input_var1: -1000000.0000000000
Input_var2: 1.0000000000000000
Smooth min:
Expected outcome: -1000000.0000002501
True outcome: -1000000.0000002501
Smooth max:
Expected outcome: 1.0000002500019036
True outcome: 1.0000002500019036
Success!
==================================================
====================== Running fill holes tests ======================
-------------- easy_fill_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: easy_fill_test
-----------------------------------------------
-------------- below_thresh_test --------------
--- Testing fill type: global ---
--- Testing fill type: sliding_window ---
--- Testing fill type: widening_windows ---
--- Testing fill type: smart_window ---
--- Testing fill type: smart_window_smooth ---
--- Testing fill type: parallel ---
PASS: below_thresh_test
-----------------------------------------------
====================== End fill holes tests ======================
Running /home/jenkins/workspace/clubb_G_unit_gfortran_tests/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
