Skip to content
Success

Console Output

Skipping 211 KB.. Full Log
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -Wno-conversion -c	../src/error.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -Wno-conversion -c -I../src	../src/enhanced_simann.F90
ar rv /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib/libclubb_other.a arm_97.o cloud_feedback.o parameters_microphys.o silhs_category_variance_module.o gabls2.o error.o clex9_oct14.o ice_dfsn_module.o wangara.o stat_file_utils.o clubb_model_settings.o lba.o jun25.o spec_hum_to_mixing_ratio.o arm.o mixed_moment_PDF_integrals.o microphys_driver.o clex9_nov02.o arm_0003.o enhanced_simann.o extended_atmosphere_module.o clubb_driver.o input_grads.o gfdl_activation.o extrapolation.o fire.o input_netcdf.o mpace_a.o bomex.o coamps_microphys_driver_module.o rico.o rad_lwsw_module.o bugsrad_driver.o cloud_sed_module.o sounding.o time_dependent_input.o morrison_microphys_module.o gabls3.o simple_rad_module.o estimate_scm_microphys_module.o quicksort.o gabls3_night.o astex_a209.o dycoms2_rf02.o twp_ice.o advance_microphys_module.o atex.o lh_microphys_driver_module.o hydrostatic_module.o mpace_b.o KK_microphys_module.o soil_vegetation.o arm_3year.o parameters_radiation.o microphys_init_cleanup.o diag_ustar_module.o input_interpret.o sfc_flux.o cobra.o input_fields.o dycoms2_rf01.o nov11.o cos_solar_zen_module.o
ar: creating /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib/libclubb_other.a
a - arm_97.o
a - cloud_feedback.o
a - parameters_microphys.o
a - silhs_category_variance_module.o
a - gabls2.o
a - error.o
a - clex9_oct14.o
a - ice_dfsn_module.o
a - wangara.o
a - stat_file_utils.o
a - clubb_model_settings.o
a - lba.o
a - jun25.o
a - spec_hum_to_mixing_ratio.o
a - arm.o
a - mixed_moment_PDF_integrals.o
a - microphys_driver.o
a - clex9_nov02.o
a - arm_0003.o
a - enhanced_simann.o
a - extended_atmosphere_module.o
a - clubb_driver.o
a - input_grads.o
a - gfdl_activation.o
a - extrapolation.o
a - fire.o
a - input_netcdf.o
a - mpace_a.o
a - bomex.o
a - coamps_microphys_driver_module.o
a - rico.o
a - rad_lwsw_module.o
a - bugsrad_driver.o
a - cloud_sed_module.o
a - sounding.o
a - time_dependent_input.o
a - morrison_microphys_module.o
a - gabls3.o
a - simple_rad_module.o
a - estimate_scm_microphys_module.o
a - quicksort.o
a - gabls3_night.o
a - astex_a209.o
a - dycoms2_rf02.o
a - twp_ice.o
a - advance_microphys_module.o
a - atex.o
a - lh_microphys_driver_module.o
a - hydrostatic_module.o
a - mpace_b.o
a - KK_microphys_module.o
a - soil_vegetation.o
a - arm_3year.o
a - parameters_radiation.o
a - microphys_init_cleanup.o
a - diag_ustar_module.o
a - input_interpret.o
a - sfc_flux.o
a - cobra.o
a - input_fields.o
a - dycoms2_rf01.o
a - nov11.o
a - cos_solar_zen_module.o
make[1]: Leaving directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
cd /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../obj; make -f Make.clubb_standalone
cd /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../obj; make -f Make.clubb_tuner
cd /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../obj; make -f Make.jacobian
cd /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../obj; make -f Make.G_unit_tests
cd /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../obj; make -f Make.int2txt
make[1]: Entering directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/clubb_standalone.F90
make[1]: Entering directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/mu_sigma_hm_tests.F90
make[1]: Entering directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran  -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c -I../src/CLUBB_core	../src/int2txt.F90
make[1]: Entering directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/jacobian.F90
make[1]: Entering directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/nrtype.f90
gfortran int2txt.o -o /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../bin/int2txt  -L/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib -lclubb_other -lclubb_coamps -lclubb_gfdlact -lsilhs -lclubb_bugsrad -lclubb_KK_microphys -lclubb_parabolic -lclubb_morrison -lmicrophys_utils -lclubb_param -lclubb_lapack -march=native -msse3 -mfpmath=sse -fopenmp -L/usr/local/NETCDF-Fortran-4.4.5-gfortran/lib -lnetcdff 
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/nrutil.f90
gfortran clubb_standalone.o -o /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../bin/clubb_standalone  -L/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib -lclubb_other -lclubb_coamps -lclubb_gfdlact -lsilhs -lclubb_bugsrad -lclubb_KK_microphys -lclubb_parabolic -lclubb_morrison -lmicrophys_utils -lclubb_param -lclubb_lapack -march=native -msse3 -mfpmath=sse -fopenmp -L/usr/local/NETCDF-Fortran-4.4.5-gfortran/lib -lnetcdff 
make[1]: Leaving directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/nr.f90
cd /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../obj; make -f Make.clubb_thread_test
make[1]: Entering directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/clubb_thread_test.F90
gfortran clubb_thread_test.o -o /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../bin/clubb_thread_test  -L/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib -lclubb_other -lclubb_coamps -lclubb_gfdlact -lsilhs -lclubb_bugsrad -lclubb_KK_microphys -lclubb_parabolic -lclubb_morrison -lmicrophys_utils -lclubb_param -lclubb_lapack -march=native -msse3 -mfpmath=sse -fopenmp -L/usr/local/NETCDF-Fortran-4.4.5-gfortran/lib -lnetcdff 
make[1]: Leaving directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	../src/clubb_tuner.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/w_up_in_cloud_tests.F90
../src/G_unit_test_types/w_up_in_cloud_tests.F90:24:8:

   24 |     use clubb_precision, only: &
      |        1
Warning: Unused parameter ‘core_rknd’ which has been explicitly imported at (1) [-Wunused-parameter]
make[1]: Leaving directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/tuner_tests.F90
gfortran jacobian.o -o /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../bin/jacobian  -L/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib -lclubb_other -lclubb_coamps -lclubb_gfdlact -lsilhs -lclubb_bugsrad -lclubb_KK_microphys -lclubb_parabolic -lclubb_morrison -lmicrophys_utils -lclubb_param -lclubb_lapack -march=native -msse3 -mfpmath=sse -fopenmp -L/usr/local/NETCDF-Fortran-4.4.5-gfortran/lib -lnetcdff 
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/hole_filling_tests.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/KK_integrals_tests.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/smooth_heaviside_tests.F90
../src/G_unit_test_types/smooth_heaviside_tests.F90:98:7:

   98 |       i
      |       1
Warning: Unused variable ‘i’ declared at (1) [-Wunused-variable]
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/spurious_source_test.F90
../src/G_unit_test_types/smooth_heaviside_tests.F90:120:80:

  120 |     total_mismatches = total_mismatches + COUNT(abs(result - result_cmp) >= eps)
      |                                                                                ^
Warning: ‘total_mismatches’ is used uninitialized in this function [-Wuninitialized]
../src/G_unit_test_types/smooth_heaviside_tests.F90:32:31:

   32 |     integer :: total_mismatches
      |                               ^
note: ‘total_mismatches’ was declared here
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/silhs_category_test.F90
make[1]: Leaving directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/corr_cholesky_mtx_tests.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/Nc_Ncn_test.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/read_corr_mtx_test.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/pdf_parameter_tests.F90
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_test_types/smooth_min_max_tests.F90
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/ran_state.f90
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/amoeba.f90
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/amebsa.f90
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/ran1_v.f90
gfortran -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2  -c	/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../src/Numerical_recipes/ran1_s.f90
gfortran ran_state.o amoeba.o ran1_v.o nrtype.o amebsa.o nr.o nrutil.o clubb_tuner.o ran1_s.o -o /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../bin/clubb_tuner  -L/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib -lclubb_other -lclubb_coamps -lclubb_gfdlact -lsilhs -lclubb_bugsrad -lclubb_KK_microphys -lclubb_parabolic -lclubb_morrison -lmicrophys_utils -lclubb_param -lclubb_lapack -march=native -msse3 -mfpmath=sse -fopenmp -L/usr/local/NETCDF-Fortran-4.4.5-gfortran/lib -lnetcdff 
make[1]: Leaving directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
gfortran -DNETCDF -D__GFORTRAN__ -DCLUBB_REAL_TYPE=8 -DCOAMPS_MICRO -DTUNER -DSILHS -Dradoffline -Dnooverlap -DCLUBB -I/usr/local/NETCDF-Fortran-4.4.5-gfortran/include/ -march=native -msse3 -mfpmath=sse -fopenmp -g -fbounds-check -mieee-fp -ffpe-trap=invalid,zero,overflow -finit-real=nan -finit-integer=-99999 -finit-logical=false -fall-intrinsics -fbacktrace -O2 -Wall -Wextra -Wconversion -Wunderflow -Wcharacter-truncation -std=f2008 -pedantic -c	../src/G_unit_tests.F90
../src/G_unit_tests.F90:176:13:

  176 |     stats_sfc
      |             1
Warning: Unused variable ‘stats_sfc’ declared at (1) [-Wunused-variable]
../src/G_unit_tests.F90:175:12:

  175 |     stats_zm, &
      |            1
Warning: Unused variable ‘stats_zm’ declared at (1) [-Wunused-variable]
../src/G_unit_tests.F90:174:12:

  174 |     stats_zt, &
      |            1
Warning: Unused variable ‘stats_zt’ declared at (1) [-Wunused-variable]
gfortran mu_sigma_hm_tests.o w_up_in_cloud_tests.o tuner_tests.o hole_filling_tests.o KK_integrals_tests.o smooth_heaviside_tests.o spurious_source_test.o silhs_category_test.o corr_cholesky_mtx_tests.o G_unit_tests.o pdf_parameter_tests.o Nc_Ncn_test.o smooth_min_max_tests.o read_corr_mtx_test.o -o /home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../bin/G_unit_tests  -L/home/jenkins/workspace/clubb_G_unit_gfortran_tests/compile/../lib -lclubb_other -lclubb_coamps -lclubb_gfdlact -lsilhs -lclubb_bugsrad -lclubb_KK_microphys -lclubb_parabolic -lclubb_morrison -lmicrophys_utils -lclubb_param -lclubb_lapack -march=native -msse3 -mfpmath=sse -fopenmp -L/usr/local/NETCDF-Fortran-4.4.5-gfortran/lib -lnetcdff 
make[1]: Leaving directory '/home/jenkins/workspace/clubb_G_unit_gfortran_tests/obj'
python ../utilities/check_for_missing_threadprivate.py ../src/CLUBB_core/ ../src/SILHS/
check_for_missing_threadprivate.py has begun.
check_for_missing_threadprivate.py passed.
perl ../utilities/CLUBBStandardsCheck.pl ../src/*.F90
CLUBBStandardsCheck.pl has begun.
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2324 :          call setup_pdf_parameters_api( gr(1), gr(1)%nz, pdf_dim, dt_main,                & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2349 :          call hydrometeor_mixed_moments( gr(1), gr(1)%nz, pdf_dim, hydromet,               & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/clubb_driver.F90
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl has finished.
perl ../utilities/CLUBBStandardsCheck.pl ../src/CLUBB_core/*.F90
CLUBBStandardsCheck.pl has begun.
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
678 :         Richardson_num(i,:) = Lscale_width_vert_avg( gr(i), Richardson_num(i,:), Lscale_zm(i,:), rho_ds_zm(i,:), &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
693 :         Cx_fnc_Richardson(i,:) = Lscale_width_vert_avg( gr(i), Cx_fnc_Richardson(i,:), Lscale_zm(i,:), rho_ds_zm(i,:), &
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/advance_helper_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1634 :                lhs(kp1_tdiag,i,:), lhs(k_tdiag,i,:), lhs(km1_tdiag,i,:), rhs(i,:,:), & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1647 :                            lhs(kp1_tdiag,i,:),  lhs(k_tdiag,i,:), lhs(km1_tdiag,i,:), rhs(i,:,:), & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1648 :                            solution(i,:,:) )                                                  ! intent(out)
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/advance_windm_edsclrm_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
574 :         a1_zt(:,:) = max( zm2zt( nz, ngrdcol, gr, a1(:,:) ), zero_threshold )  ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
807 :                            lhs_ma_zt(:,i,:) )                                           ! intent(out)
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/advance_wp2_wp3_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
860 :     write(iunit,*) "l_fix_w_chi_eta_correlations = ", clubb_config_flags%l_fix_w_chi_eta_correlations
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
877 :     write(iunit,*) "l_use_tke_in_wp3_pr_turb_term = ", clubb_config_flags%l_use_tke_in_wp3_pr_turb_term
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/model_flags.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1299 :                                               ! https://github.com/larson-group/clubb/issues/965#issuecomment-1119816722
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1300 :                                               ! for a plot on how output behaves with varying min_max_smth_mag
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1303 :     logical, parameter :: l_smooth_min_max = .false.  ! whether to apply smooth min and max functions
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1338 :         ustar(i) = smooth_max( ( upwp_sfc(i)**2 + vpwp_sfc(i)**2 )**(one_fourth), ufmin, one * min_max_smth_mag)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1348 :                                                             sqrt( (ddzt( nz, ngrdcol, gr, um ))**2 & 
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1409 :                                          one - zt2zm(nz, ngrdcol, gr, ice_supersat_frac) / 0.007_core_rknd, &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1447 :                             * sqrt(em(i,k)) / ( gr(i)%zm(k) - sfc_elevation(i) + z_displace )  ! small
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1456 :                           / smooth_max( nz, ngrdcol, 1.0e-7_core_rknd, brunt_vaisala_freq_sqd_smth, &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1458 :                       0.3_core_rknd, 0.3_core_rknd * min_max_smth_mag ), one, one * min_max_smth_mag ) &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1550 :       tau_zm           = smooth_min( nz, ngrdcol, one / invrs_tau_zm, tau_max_zm, 1.0e3_core_rknd * min_max_smth_mag )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1551 :       tau_zt           = smooth_min( nz, ngrdcol, zm2zt( nz, ngrdcol, gr, tau_zm ), tau_max_zt, 1.0e3_core_rknd * min_max_smth_mag )
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/mixing_length.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
86 :     real( kind = core_rknd ), dimension(ngrdcol,nz), intent(out) :: sigma_sqd_w ! PDF width parameter      [-]
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/sigma_sqd_w_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1198 :       um_col, vm_col, upwp_col, vpwp_col, up2_col, vp2_col, up3_col, vp3_col, &               ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1200 :       wp2_col, wp3_col, rtp2_col, rtp3_col, thlp2_col, thlp3_col, rtpthlp_col, &          ! intent(inout)
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/clubb_api_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
908 :              ( gr(i), thlm_forcing(i,:), rtm_forcing(i,:), um_forcing(i,:),                         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
909 :                vm_forcing(i,:), wm_zm(i,:), wm_zt(i,:), p_in_Pa(i,:),                                 & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
910 :                rho_zm(i,:), rho(i,:), exner(i,:), rho_ds_zm(i,:),                                     & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
911 :                rho_ds_zt(i,:), invrs_rho_ds_zm(i,:), invrs_rho_ds_zt(i,:),                       & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
912 :                thv_ds_zm(i,:), thv_ds_zt(i,:), wpthlp_sfc(i), wprtp_sfc(i), upwp_sfc(i),             & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
913 :                vpwp_sfc(i), um(i,:), upwp(i,:), vm(i,:), vpwp(i,:), up2(i,:), vp2(i,:),                            & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
914 :                rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), wp2(i,:), wp3(i,:),                                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
915 :                rtp2(i,:), thlp2(i,:), rtpthlp(i,:),                                              & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
918 :                wpsclrp_sfc(i,:), wpedsclrp_sfc(i,:), sclrm(i,:,:), wpsclrp(i,:,:), sclrp2(i,:,:),                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
919 :                sclrprtp(i,:,:), sclrpthlp(i,:,:), sclrm_forcing(i,:,:), edsclrm(i,:,:), edsclrm_forcing(i,:,:) )       ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1221 :     thlp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, thlp2(:,:) ), thl_tol**2 ) ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1222 :     rtp2_zt(:,:)    = max( zm2zt( nz, ngrdcol, gr, rtp2(:,:) ), rt_tol**2 )   ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1448 :         call calc_sfc_varnce( upwp_sfc(i), vpwp_sfc(i), wpthlp_sfc(i), wprtp_sfc(i),     & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1449 :                              um(i,2), vm(i,2), Lscale_up(i,2), wpsclrp_sfc(i,:),        & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1665 :                                       stats_zm,                                      & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1688 :       call advance_xm_wpxp( nz, ngrdcol, gr, dt_advance, sigma_sqd_w, wm_zm, wm_zt, wp2, & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1746 :                                             sclrm(i,:,:), sclrm_trsport_only(i,:,:), & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1855 :         call clip_covars_denom( gr(i), dt, rtp2(i,:), thlp2(i,:), up2(i,:), vp2(i,:), wp2(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1862 :                                 wprtp(i,:), wpthlp(i,:), upwp(i,:), vpwp(i,:), wpsclrp(i,:,:), & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1961 :         call clip_covars_denom( gr(i), dt, rtp2(i,:), thlp2(i,:), up2(i,:), vp2(i,:), wp2(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1968 :                                 wprtp(i,:), wpthlp(i,:), upwp(i,:), vpwp(i,:), wpsclrp(i,:,:), & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2060 :         call advance_xp3( gr(i), dt, rtm(i,:), thlm(i,:), rtp2(i,:), thlp2(i,:), wprtp(i,:),          & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2061 :                           wpthlp(i,:), wprtp2(i,:), wpthlp2(i,:), rho_ds_zm(i,:),         & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2063 :                           sclrm(i,:,:), sclrp2(i,:,:), wpsclrp(i,:,:), wpsclrp2(i,:,:),           & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2066 :                           rtp3(i,:), thlp3(i,:), sclrp3(i,:,:) )                         ! Intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2120 :       thlp2_zt(:,:)  = max( zm2zt( nz, ngrdcol, gr, thlp2(:,:) ), thl_tol**2 ) ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2121 :       rtp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, rtp2(:,:) ), rt_tol**2 )   ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2410 :                gr(i), um(i,:), vm(i,:), upwp(i,:), vpwp(i,:), up2(i,:), vp2(i,:),                      & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2411 :                thlm(i,:), rtm(i,:), wprtp(i,:), wpthlp(i,:),                              & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2412 :                wp2(i,:), wp3(i,:), rtp2(i,:), rtp3(i,:), thlp2(i,:), thlp3(i,:), rtpthlp(i,:),           & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2413 :                wpthvp(i,:), wp2thvp(i,:), rtpthvp(i,:), thlpthvp(i,:),                    & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2414 :                p_in_Pa(i,:), exner(i,:), rho(i,:), rho_zm(i,:),                           & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2415 :                rho_ds_zm(i,:), rho_ds_zt(i,:), thv_ds_zm(i,:), thv_ds_zt(i,:),            & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2416 :                wm_zt(i,:), wm_zm(i,:), rcm(i,:), wprcp(i,:), rc_coef(i,:), rc_coef_zm(i,:),         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2417 :                rcm_zm(i,:), rtm_zm(i,:), thlm_zm(i,:), cloud_frac(i,:), ice_supersat_frac(i,:),& ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2420 :                thvm(i,:), ug(i,:), vg(i,:), Lscale(i,:), wpthlp2(i,:), wp2thlp(i,:), wprtp2(i,:), wp2rtp(i,:),& ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2421 :                Lscale_up(i,:), Lscale_down(i,:), tau_zt(i,:), Kh_zt(i,:), wp2rcp(i,:),         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2422 :                wprtpthlp(i,:), sigma_sqd_w_zt(i,:), rsat(i,:), wp2_zt(i,:), thlp2_zt(i,:),     & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2423 :                wpthlp_zt(i,:), wprtp_zt(i,:), rtp2_zt(i,:), rtpthlp_zt(i,:), up2_zt(i,:),      & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2424 :                vp2_zt(i,:), upwp_zt(i,:), vpwp_zt(i,:), wpup2(i,:), wpvp2(i,:),                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2427 :                rtprcp(i,:), rcp2(i,:), em(i,:), a3_coef(i,:), a3_coef_zt(i,:),                 & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2428 :                wp3_zm(i,:), wp3_on_wp2(i,:), wp3_on_wp2_zt(i,:), Skw_velocity(i,:),       & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2430 :                pdf_params_single_col(i), pdf_params_zm_single_col(i), sclrm(i,:,:), sclrp2(i,:,:),              & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2431 :                sclrprtp(i,:,:), sclrpthlp(i,:,:), sclrm_forcing(i,:,:), sclrpthvp(i,:,:),         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2432 :                wpsclrp(i,:,:), sclrprcp(i,:,:), wp2sclrp(i,:,:), wpsclrp2(i,:,:), wpsclrprtp(i,:,:),     & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2433 :                wpsclrpthlp(i,:,:), wpedsclrp(i,:,:), edsclrm(i,:,:), edsclrm_forcing(i,:,:),      & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2441 :              gr(i), thlm_forcing(i,:), rtm_forcing(i,:), um_forcing(i,:),                         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2442 :              vm_forcing(i,:), wm_zm(i,:), wm_zt(i,:), p_in_Pa(i,:),                                 & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2443 :              rho_zm(i,:), rho(i,:), exner(i,:), rho_ds_zm(i,:),                                     & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2444 :              rho_ds_zt(i,:), invrs_rho_ds_zm(i,:), invrs_rho_ds_zt(i,:),                       & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2445 :              thv_ds_zm(i,:), thv_ds_zt(i,:), wpthlp_sfc(i), wprtp_sfc(i), upwp_sfc(i),             & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2446 :              vpwp_sfc(i), um(i,:), upwp(i,:), vm(i,:), vpwp(i,:), up2(i,:), vp2(i,:),                            & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2447 :              rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), wp2(i,:), wp3(i,:),                                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2448 :              rtp2(i,:), thlp2(i,:), rtpthlp(i,:),                                              & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2451 :              wpsclrp_sfc(i,:), wpedsclrp_sfc(i,:), sclrm(i,:,:), wpsclrp(i,:,:), sclrp2(i,:,:),                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2452 :              sclrprtp(i,:,:), sclrpthlp(i,:,:), sclrm_forcing(i,:,:), edsclrm(i,:,:), edsclrm_forcing(i,:,:)       ) ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3022 :     wp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, wp2(:,:) ), w_tol_sqd ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3024 :     thlp2_zt(:,:) = max( zm2zt( nz, ngrdcol, gr, thlp2(:,:) ), thl_tol**2 ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3026 :     rtp2_zt(:,:)  = max( zm2zt( nz, ngrdcol, gr, rtp2(:,:) ), rt_tol**2 ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3028 :     up2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, up2(:,:) ), w_tol_sqd ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3030 :     vp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, vp2(:,:) ), w_tol_sqd ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3034 :       sclrp2_zt(:,:,j) = max( zm2zt( nz, ngrdcol, gr, sclrp2(:,:,j) ), sclr_tol(j)**2 ) ! Pos. def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3164 :     sigma_sqd_w_zt(:,:) = max( zm2zt( nz, ngrdcol, gr, sigma_sqd_w(:,:) ), zero_threshold )  ! Pos. def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3171 :     rtp2_zt(:,:)    = max( zm2zt( nz, ngrdcol, gr, rtp2(:,:) ), rt_tol**2 )   ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3172 :     thlp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, thlp2(:,:) ), thl_tol**2 ) ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3173 :     up2_zt(:,:)     = max( zm2zt( nz, ngrdcol, gr, up2(:,:) ), w_tol_sqd )    ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3174 :     vp2_zt(:,:)     = max( zm2zt( nz, ngrdcol, gr, vp2(:,:) ), w_tol_sqd )    ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3198 :       sclrp2_zt(:,:,j)    = max( zm2zt( nz, ngrdcol, gr, sclrp2(:,:,j) ), sclr_tol(j)**2 ) ! Pos. def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3407 :         rcp2(:,:) = max( zt2zm( nz, ngrdcol, gr, rcp2_zt(:,:) ), zero_threshold )  ! Pos. def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3479 :                rcm(i,:), wp2thvp(i,:), wpsclrprtp(i,:,:), wpsclrp2(i,:,:),          & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4780 :                        ( ( 0.5_core_rknd / gr(i)%invrs_dzm(k-1) ) / ( gr(i)%zt(k) - gr(i)%zm(k-1) ) ) &
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/advance_clubb_core_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
767 :           call error_prints_xm_wpxp( gr(i), dt, sigma_sqd_w(i,:), wm_zm(i,:), wm_zt(i,:), wp2(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
768 :                                      Lscale(i,:), wp3_on_wp2(i,:), wp3_on_wp2_zt(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
769 :                                      Kh_zt(i,:), Kh_zm(i,:), invrs_tau_C6_zm(i,:), Skw_zm(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
773 :                                      wpthlp_forcing(i,:), thlm_ref(i,:), rho_ds_zm(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
779 :                                      wp2sclrp(i,:,:), sclrpthvp(i,:,:), sclrm_forcing(i,:,:), & ! intent(in) 
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
780 :                                      sclrp2(i,:,:), exner(i,:), rcm(i,:), p_in_Pa(i,:), thvm(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
783 :                                      um_forcing(i,:), vm_forcing(i,:), ug(i,:), vg(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
784 :                                      wpthvp(i,:), fcor(i), um_ref(i,:), vm_ref(i,:), up2(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
785 :                                      vp2(i,:), uprcp(i,:), vprcp(i,:), rc_coef(i,:), rtm(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
786 :                                      wprtp(i,:), thlm(i,:), wpthlp(i,:), sclrm(i,:,:), wpsclrp(i,:,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
787 :                                      um(i,:), upwp(i,:), vm(i,:), vpwp(i,:), rtm_old(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
789 :                                      sclrm_old(i,:,:), wpsclrp_old(i,:,:), um_old(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1363 :       call wpxp_terms_ac_pr2_lhs( gr(i), C7_Skw_fnc(i,:), wm_zt(i,:), gr(i)%invrs_dzm(:), & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1416 :         call term_ma_zt_lhs( gr(i), wm_zt(i,:), gr(i)%invrs_dzt(:), gr(i)%invrs_dzm(:), & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1710 :       call wpxp_terms_bp_pr3_rhs( nz, ngrdcol, gr, (one+C7_Skw_fnc), thv_ds_zm, xpthvp, & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2639 :       call xm_wpxp_rhs( nz, ngrdcol, gr, xm_wpxp_scalar, l_iter, dt, sclrm(:,:,j), wpsclrp(:,:,j), & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2710 :           call stat_update_var( iupwp_pr4, C_uu_shr * wp2(i,:) * ddzt( gr(i), um(i,:) ), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2712 :           call stat_update_var( ivpwp_pr4, C_uu_shr * wp2(i,:) * ddzt( gr(i), vm(i,:) ), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3507 :       call xm_wpxp_rhs( nz, ngrdcol, gr, xm_wpxp_scalar, l_iter, dt, sclrm(:,:,j), wpsclrp(:,:,j),            & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4105 :         call wpxp_terms_ac_pr2_lhs( gr(i), zero_vector(i,:), wm_zt(i,:), gr(i)%invrs_dzm, & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4111 :         call wpxp_terms_ac_pr2_lhs( gr(i), (one+C7_Skw_fnc(i,:)), wm_zt(i,:), gr(i)%invrs_dzm, & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4357 :           call xm_correction_wpxp_cl( gr(i), solve_type, dt, wpxp_chnge(i,:), gr(i)%invrs_dzt, & ! intent(in)
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/advance_xm_wpxp_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
505 :     call calc_xp2_xpyp_ta_terms( nz, ngrdcol, gr, wprtp, wprtp2, wpthlp, wpthlp2, wprtpthlp,      & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
560 :       call solve_xp2_xpyp_with_multiple_lhs( nz, ngrdcol, gr, C2rt_1d, C2thl_1d, C2rtthl_1d, C2sclr_1d, & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
601 :         lhs_dp1(i,k) = lhs_dp1(i,k) + term_dp1_lhs( C14_1d(i,k), invrs_tau_C14_zm(i,k) ) * gamma_over_implicit_ts
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1225 :   subroutine solve_xp2_xpyp_with_single_lhs( nz, ngrdcol, gr, C2x, invrs_tau_xp2_zm, rtm, thlm, wprtp, &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1477 :         sclrprtp(:,:,:)   = one_half * ( sclrprtp(:,:,:)   + solution(:,:,3+sclr_dim+1:3+2*sclr_dim) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1478 :         sclrpthlp(:,:,:)  = one_half * ( sclrpthlp(:,:,:)  + solution(:,:,3+2*sclr_dim+1:3+3*sclr_dim) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1520 :   subroutine solve_xp2_xpyp_with_multiple_lhs( nz, ngrdcol, gr, C2rt_1d, C2thl_1d, C2rtthl_1d, C2sclr_1d, &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1773 :         lhs_dp1(i,k) = term_dp1_lhs( C2rtthl_1d(i,k), invrs_tau_xp2_zm(i,k) ) * gamma_over_implicit_ts
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1816 :           lhs_dp1(i,k) = term_dp1_lhs( C2sclr_1d(i,k), invrs_tau_xp2_zm(i,k) ) * gamma_over_implicit_ts
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2009 :           sclrprtp(:,:,:)  = one_half * ( sclrprtp(:,:,:)  + sclr_solution(:,:,sclr_dim+1:2*sclr_dim) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2010 :           sclrpthlp(:,:,:) = one_half * ( sclrpthlp(:,:,:) + sclr_solution(:,:,2*sclr_dim+1:3*sclr_dim) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2234 :     real( kind = core_rknd ), dimension(ngrdcol) :: rcond  ! Est. of the reciprocal of the condition # on the matrix
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2283 :                lhs(kp1_mdiag,i,:), lhs(k_mdiag,i,:), lhs(km1_mdiag,i,:), rhs(i,:,1:nrhs),  & ! Intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2284 :                xapxbp(i,:,1:nrhs), rcond(i) )                                             ! Intent(out)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2703 :         rhs(i,k) = rhs(i,k) + ( one - C_uu_shr ) * term_tp( xam(i,k+1), xam(i,k), xam(i,k+1), xam(i,k), & 
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3507 :           sgn_t_vel_rtp2(i,:) = sgn_turbulent_velocity( gr(i), term_wprtp2_explicit_zm(i,:), rtp2(i,:) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3529 :           sgn_t_vel_thlp2(i,:) = sgn_turbulent_velocity( gr(i), term_wpthlp2_explicit_zm(i,:), thlp2(i,:) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3551 :           sgn_t_vel_rtpthlp(i,:) = sgn_turbulent_velocity( gr(i), term_wprtpthlp_explicit_zm(i,:), rtpthlp(i,:) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3608 :               sgn_t_vel_sclrprtp(i,:) = sgn_turbulent_velocity( gr(i), term_wprtpsclrp_explicit_zm(i,:), &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3637 :               sgn_t_vel_sclrpthlp(i,:) = sgn_turbulent_velocity( gr(i), term_wpthlpsclrp_explicit_zm(i,:), &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3746 :           = ( one - one_third * beta ) * a1_zt(:,:)**2 * wprtp_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3749 :           = ( one - one_third * beta ) * a1_zt(:,:)**2 * wpthlp_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3752 :           = ( one - one_third * beta ) * a1_zt(:,:)**2 * wprtp_zt(:,:) * wpthlp_zt(:,:) * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3791 :           sgn_t_vel_rtp2(:,:) = ( one - one_third * beta ) * a1_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3836 :           sgn_t_vel_thlp2(:,:) = ( one - one_third * beta ) * a1_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3854 :           = ( one - one_third * beta ) * a1(:,:)**2 * wprtp(:,:) * wpthlp(:,:) * wp3_on_wp2(:,:) / wp2(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3882 :           sgn_t_vel_rtpthlp(:,:) = ( one - one_third * beta ) * a1_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3885 :             call xpyp_term_ta_pdf_rhs_godunov( gr(i), term_wprtpthlp_explicit_zm(i,:), & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3912 :               = ( one - one_third * beta ) * a1(:,:)**2 * wpsclrp(:,:,sclr)**2 * wp3_on_wp2(:,:) / wp2(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3919 :               = ( one - one_third * beta ) * a1_zt(:,:)**2 * wpsclrp_zt(:,:,sclr)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3945 :               sgn_t_vel_sclrp2(:,:) = ( one - one_third * beta ) * a1_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3948 :                 call xpyp_term_ta_pdf_rhs_godunov( gr(i), term_wpsclrp2_explicit_zm(i,:), & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3953 :                                                    rhs_ta_wpsclrp2(i,:,sclr) )          ! Intent(out)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3966 :               = ( one - one_third * beta ) * a1(:,:)**2 * wpsclrp(:,:,sclr) * wprtp(:,:) * wp3_on_wp2(:,:) / wp2(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3998 :               sgn_t_vel_sclrprtp(:,:) = ( one - one_third * beta ) * a1_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4001 :                 call xpyp_term_ta_pdf_rhs_godunov( gr(i), term_wprtpsclrp_explicit_zm(i,:), & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4006 :                                                    rhs_ta_wprtpsclrp(i,:,sclr)          ) ! Intent(out)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4019 :               = ( one - one_third * beta ) * a1(:,:)**2 * wpsclrp(:,:,sclr) * wpthlp(:,:) * wp3_on_wp2(:,:) / wp2(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4052 :               sgn_t_vel_sclrpthlp(:,:) = ( one - one_third * beta ) * a1_zt(:,:)**2 * wp3_on_wp2_zt(:,:) / wp2_zt(:,:)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4055 :                 call xpyp_term_ta_pdf_rhs_godunov( gr(i), term_wpthlpsclrp_explicit_zm(i,:), & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4057 :                                                    gr(i)%invrs_dzm(:),                 & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4060 :                                                    rhs_ta_wpthlpsclrp(i,:,sclr)          ) ! Intent(out)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4086 :             coef_wprtp2_implicit_zm(i,:) = zt2zm( gr(i), pdf_implicit_coefs_terms(i)%coef_wprtp2_implicit )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4087 :             sgn_t_vel_rtp2(i,:) = sgn_turbulent_velocity( gr(i), coef_wprtp2_implicit_zm(i,:) * rtp2(i,:), rtp2(i,:) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4109 :             coef_wpthlp2_implicit_zm(i,:) = zt2zm( gr(i), pdf_implicit_coefs_terms(i)%coef_wpthlp2_implicit )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4110 :             sgn_t_vel_thlp2(i,:) = sgn_turbulent_velocity( gr(i), coef_wpthlp2_implicit_zm(i,:) * thlp2(i,:), thlp2(i,:) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4132 :             coef_wprtpthlp_implicit_zm(i,:) = zt2zm( gr(i), pdf_implicit_coefs_terms(i)%coef_wprtpthlp_implicit )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4133 :             term_wprtpthlp_explicit_zm(i,:) = zt2zm( gr(i), pdf_implicit_coefs_terms(i)%term_wprtpthlp_explicit )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4134 :             sgn_t_vel_rtpthlp(i,:) = sgn_turbulent_velocity( gr(i), coef_wprtpthlp_implicit_zm(i,:) * rtpthlp(i,:) &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4135 :                                                         + term_wprtpthlp_explicit_zm(i,:), rtpthlp(i,:) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4519 :                 = sgn_turbulent_velocity( gr(i), coef_wprtpsclrp_implicit_zm(i,:) * sclrprtp(i,:,sclr) &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4570 :                 = sgn_turbulent_velocity( gr(i), coef_wpthlpsclrp_implicit_zm(i,:) * sclrpthlp(i,:,sclr) &
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/advance_xp2_xpyp_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
238 :     real( kind = core_rknd ), dimension(ngrdcol, nz, min(1,sclr_dim), 2 ), intent(in) ::  & ! h1g, 2010-06-15
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
575 :                            wm(i,:), rtm(i,:), thlm(i,:), wp2(i,:), rtp2(i,:), thlp2(i,:),                             & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
576 :                            Skw(i,:), Skrt(i,:), Skthl(i,:), wprtp(i,:), wpthlp(i,:),                             & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
587 :         call new_pdf_driver( gr(i), wm(i,:), rtm(i,:), thlm(i,:), wp2(i,:), rtp2(i,:), thlp2(i,:), Skw(i,:),                   & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
588 :                             wprtp(i,:), wpthlp(i,:), rtpthlp(i,:),                                     & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
593 :                             Skrt(i,:), Skthl(i,:),                                                & ! In/Out
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
602 :                             F_w(i,:), F_rt(i,:), F_thl(i,:), min_F_w(i,:), max_F_w(i,:),                         & ! Out
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
603 :                             min_F_rt(i,:), max_F_rt(i,:), min_F_thl(i,:), max_F_thl(i,:) )                    ! Out
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
608 :                           wm(i,:), rtm(i,:), thlm(i,:), wp2(i,:), rtp2(i,:), thlp2(i,:),                                & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
609 :                           Skw(i,:), Skrt(i,:), Skthl(i,:), wprtp(i,:), wpthlp(i,:),                                & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
620 :         call LY93_driver( gr(i), wm(i,:), rtm(i,:), thlm(i,:), wp2(i,:), rtp2(i,:),                          & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
621 :                           thlp2(i,:), Skw(i,:), Skrt(i,:), Skthl(i,:),                           & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
632 :         call new_hybrid_pdf_driver( gr(i), wm(i,:), rtm(i,:), thlm(i,:), um(i,:), vm(i,:),          & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
633 :                                     wp2(i,:), rtp2(i,:), thlp2(i,:), up2(i,:), vp2(i,:),         & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
634 :                                     Skw(i,:), wprtp(i,:), wpthlp(i,:), upwp(i,:), vpwp(i,:),     & ! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
639 :                                     Skrt(i,:), Skthl(i,:), Sku(i,:), Skv(i,:), Sksclr(i,:,:),      & ! I/O
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1040 :             chi_at_ice_sat1 = ( sat_mixrat_ice( p_in_Pa(i,k), tl1(i,k) ) - pdf_params%rsatl_1(i,k) ) &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1045 :                                             pdf_params%ice_supersat_frac_1(i,k), rc_1_ice(i,k) )! out
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1085 :               chi_at_ice_sat2 = ( sat_mixrat_ice( p_in_Pa(i,k), tl2(i,k) ) - pdf_params%rsatl_2(i,k) ) &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1139 :         rcm(i,k) = pdf_params%mixt_frac(i,k) * pdf_params%rc_1(i,k) + ( one - pdf_params%mixt_frac(i,k) ) &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1318 :             - ( one - pdf_params%mixt_frac(i,:) ) * corr_sclr_thl_2(i,:,j) * pdf_params%cthl_2(i,:) &
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/pdf_closure_module.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2285 :     integer :: grid_level = 1  ! grid level for stats where there is only one sensible level (eg timeseries)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2835 :           write(fstderr, *) "Warning: tot_vartn_normlzd_thlm tried to divide by zero denominator ", &
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2848 :           write(fstderr, *) "Warning: tot_vartn_normlzd_wprtp tried to divide by zero denominator ", &
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/stats_clubb_utilities.F90
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl has finished.
perl ../utilities/CLUBBStandardsCheck.pl ../src/Benchmark_cases/*.F90
CLUBBStandardsCheck.pl has begun.
CLUBBStandardsCheck.pl has finished.
perl ../utilities/CLUBBStandardsCheck.pl ../src/KK_microphys/*.F90
CLUBBStandardsCheck.pl has begun.
CLUBBStandardsCheck.pl has finished.
perl ../utilities/CLUBBStandardsCheck.pl ../src/SILHS/*.F90
CLUBBStandardsCheck.pl has begun.
CLUBBStandardsCheck.pl has finished.
perl ../utilities/CLUBBStandardsCheck.pl ../src/G_unit_test_types/*.F90
CLUBBStandardsCheck.pl has begun.
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
121 :     result = smooth_min(4, 1, input, zero, min_max_smth_mag) ! Order of nz and ngrdcol is opposite to order of indexing! This is very unintuitive!
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
174 :     ! Case 3: Make sure that on a large number of arbitrarily chosen points with small min_max_smth_mag,
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/G_unit_test_types/smooth_min_max_tests.F90
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1560 :                                       wp2(1,:), rtp2(1,:), thlp2(1,:), up2(1,:), vp2(1,:),         &! In
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1610 :                                 wp2(1,:), rtp2(1,:), thlp2(1,:), up2(1,:), vp2(1,:),              & ! In 
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/G_unit_test_types/pdf_parameter_tests.F90
--------------------------------------------------------------------------------
CLUBBStandardsCheck.pl has finished.
[Pipeline] }
[Pipeline] // stage
[Pipeline] stage
[Pipeline] { (Run CLUBB's G Unit Tests)
[Pipeline] sh
+ run_scripts/run_G_unit_tests.bash
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.940566852245361E-04
MATLAB result:    3.940566852245361E-04
Percent difference:   1.37569E-14%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x1)
Obtained result:    3.927115990964175E-04
MATLAB result:    3.927115990964175E-04
Percent difference:   0.00000E+00%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x2)
Obtained result:    3.959157231937518E-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.961224408774507E-04
MATLAB result:    3.961224408774508E-04
Percent difference:   2.73704E-14%
 
Agreement with regards to integral:  Trivariate NNL
Obtained result:    5.957899876478482E-10
MATLAB result:    5.957899876478559E-10
Percent difference:   1.30160E-12%
 
Agreement with regards to integral:  Trivariate NNL (const. x1)
Obtained result:    5.299646183560589E-10
MATLAB result:    5.299646183560650E-10
Percent difference:   1.15111E-12%
 
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.286452253001269E-10
MATLAB result:    5.286452253001284E-10
Percent difference:   2.73826E-13%
 
--------------------------------------------------
 
Agreement with regards to integral:  Trivariate NLL
Obtained result:   -9.854231981928350E-04
MATLAB result:   -9.854231981928350E-04
Percent difference:   0.00000E+00%
 
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.940070763287882E-09
MATLAB result:    8.940070763287870E-09
Percent difference:   1.29535E-13%
 
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.007487885941808E-02
MATLAB result:    1.007487885941809E-02
Percent difference:   1.03310E-13%
 
==================================================
 
Saturation option (x2_opt for covariance integrals and x1_opt for mean integrals):  2
 
Agreement with regards to integral:  Quadrivariate NNLL
Obtained result:   -4.441972023148828E-05
MATLAB result:   -4.441972023148829E-05
Percent difference:   3.05102E-14%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x1)
Obtained result:   -4.559373570842311E-05
MATLAB result:   -4.559373570842311E-05
Percent difference:   1.48623E-14%
 
Agreement with regards to integral:  Quadrivariate NNLL (const. x2)
Obtained result:   -5.040842768062482E-05
MATLAB result:   -5.040842768062482E-05
Percent difference:   1.34427E-14%
 
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:   0.00000E+00%
 
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.622010070476263E-05
MATLAB result:   -7.622010070476262E-05
Percent difference:   1.77808E-14%
 
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.609761356905441E-05
MATLAB result:   -4.609761356905441E-05
Percent difference:   0.00000E+00%
 
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.95210308158732515     
 Nc_in_cloud (from Ncn) =    51650441.529586725     
 Ncm (from Ncn) =    49176544.545665473     
 Ncnm (from Nc_in_cloud) =    50000000.000000000     
 Ncnm (from Ncm) =    49999999.999999993     
 
 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.354576334     
 Ncm (from Ncn) =    10246392.665117815     
 Ncnm (from Nc_in_cloud) =    50000000.000000000     
 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 =         2022           5          14        -300           3          29          48         607        2629         612         621         307         610         636         655        1214        2629         612         621         307         610         636         655        1214        2629         612         621         307         610         636         655        1214        2629
 hmm =    3.4504317404462319E-005
 mixt_frac =   0.27509808453504964     
 precip_frac =   0.41425751605348460     
 precip_frac_1 =   0.39572105221510401     
 precip_frac_2 =   0.42129204801871056     
 hmp2_ip_on_hmm2_ip =    1.4347299202098871     
 hmm_ip =    8.3291952631723599E-005
 hmp2 =    5.8067002385099549E-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 =         2022           5          14        -300           3          29          48         610        2632         615         624         310         613         639         658        1220        2632         615         624         310         613         639         658        1220        2632         615         624         310         613         639         658        1220        2632
 wm =  -0.43997880156777514     
 rtm =    2.0745769284101084E-002
 thlm =    293.49449296851202     
 wp2 =   0.68340536624882053     
 rtp2 =    5.3686384342290111E-005
 thlp2 =    1.3612421712729121     
 Skw =    3.6161582625188102     
 wp3 =    2.0429841922464753     
 Skrt =   -4.0789344831486005     
 rtp3 =   -1.6045119793729658E-006
 Skthl =    4.1168609083636722     
 thlp3 =    6.5383668926801564     
 wprtp =   -1.2843715689053712E-003
 wpthlp =  -0.95640945448344195     
 rtpthlp =    7.6938246676026981E-003
 
 Running tests for the above parameter set for various values of sigma_sqd_w for the setting variable (w).  Values of sigma_sqd_w are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1 (or 0.99).
 
 Running tests for the above parameter set for the full PDF (the setting of sigma_sqd_w is handled the same way as in the model code).
 
 Success!
 
 
 Performing PDF parameter values unit test
 =========================================
 
 Performing PDF parameter unit tests for LY93
 
 PDF parameter set 1:
 wm =   -5.0000000000000001E-003
 rtm =    4.2868400000000001E-003
 thlm =    309.13099999999997     
 wp2 =    3.3741899999999998E-002
 rtp2 =    1.1786499999999999E-006
 thlp2 =    5.9920000000000001E-002
 Skw =    8.6882480819056749     
 wp3 =    5.3850099999999998E-002
 Skrt =    1.4880899999999999     
 rtp3 =    1.9041733779221134E-009
 Skthl =   -1.4580900000000001     
 thlp3 =   -2.1386614365822010E-002
 wprtp =    2.3725400000000002E-005
 wpthlp =   -5.5512699999999996E-003
 rtpthlp =   -2.3917778634313013E-004
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 2:
 wm =    1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.10000000000000001     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =   -4.5000000000000000     
 wp3 =  -0.14230249470757708     
 Skrt =    1.0000000000000000     
 rtp3 =    1.1180339887498949E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =   -5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 3:
 wm =   -1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.50000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =    4.5000000000000000     
 wp3 =    1.5909902576697321     
 Skrt =    2.5000000000000000     
 rtp3 =    2.7950849718747374E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 4:
 wm =   0.10000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    5.0000000000000000     
 rtp2 =    6.9999999999999999E-006
 thlp2 =    1.0000000000000000     
 Skw =    2.0000000000000000     
 wp3 =    22.360679774997898     
 Skrt =    1.5000000000000000     
 rtp3 =    2.7780388766178202E-008
 Skthl =   -1.5000000000000000     
 thlp3 =   -1.5000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -1.0000000000000000E-002
 rtpthlp =   -2.3811761799581318E-003
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 5:
 wm =   -1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.29999999999999999     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.20000000000000001     
 Skw =   -2.0000000000000000     
 wp3 =  -0.32863353450309962     
 Skrt =  -0.25000000000000000     
 rtp3 =   -7.0710678118654745E-010
 Skthl =   -3.5000000000000000     
 thlp3 =  -0.31304951684997062     
 wprtp =    5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =    5.6920997883030825E-004
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 6:
 wm =   -1.0000000000000000E-002
 rtm =    5.0000000000000001E-003
 thlm =    305.00000000000000     
 wp2 =    1.0000000000000000     
 rtp2 =    4.9999999999999998E-007
 thlp2 =   0.10000000000000001     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -2.0124611797498108E-004
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 7:
 wm =   0.20000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.25000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.14999999999999999     
 Skw =   0.75000000000000000     
 wp3 =    9.3750000000000000E-002
 Skrt =    3.0000000000000000     
 rtp3 =    3.3541019662496846E-008
 Skthl =  -0.50000000000000000     
 thlp3 =   -2.9047375096555625E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -7.7942286340599491E-004
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 8:
 wm =   -2.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    0.0000000000000000     
 thlp2 =    0.0000000000000000     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    0.0000000000000000     
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 9:
 wm =    1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.69999999999999996     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =   -1.0000000000000000     
 rtp3 =   -2.8284271247461898E-009
 Skthl =   -2.0000000000000000     
 thlp3 =   -1.1713240371477056     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    1.0648943609579309E-003
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 PDF parameter set 10 (randomly generated):
 Random seed values =         2022           5          14        -300           3          29          48         610        2632         615         624         310         613         639         658        1220        2632         615         624         310         613         639         658        1220        2632         615         624         310         613         639         658        1220        2632
 wm =  -0.43997880156777514     
 rtm =    2.0745769284101084E-002
 thlm =    293.49449296851202     
 wp2 =   0.68340536624882053     
 rtp2 =    5.3686384342290111E-005
 thlp2 =    1.3612421712729121     
 Skw =    3.6161582625188102     
 wp3 =    2.0429841922464753     
 Skrt =   -4.0789344831486005     
 rtp3 =   -1.6045119793729658E-006
 Skthl =    4.1168609083636722     
 thlp3 =    6.5383668926801564     
 wprtp =   -1.2843715689053712E-003
 wpthlp =  -0.95640945448344195     
 rtpthlp =    7.6938246676026981E-003
 
 Running tests for the above parameter set for the setting variable (w is used here).
 
 Running tests for the above parameter set for the full PDF.
 
 Success!
 
 
 Performing PDF parameter values unit test
 =========================================
 
 Performing PDF parameter unit tests for TSDADG
 
 PDF parameter set 1:
 wm =   -5.0000000000000001E-003
 rtm =    4.2868400000000001E-003
 thlm =    309.13099999999997     
 wp2 =    3.3741899999999998E-002
 rtp2 =    1.1786499999999999E-006
 thlp2 =    5.9920000000000001E-002
 Skw =    8.6882480819056749     
 wp3 =    5.3850099999999998E-002
 Skrt =    1.4880899999999999     
 rtp3 =    1.9041733779221134E-009
 Skthl =   -1.4580900000000001     
 thlp3 =   -2.1386614365822010E-002
 wprtp =    2.3725400000000002E-005
 wpthlp =   -5.5512699999999996E-003
 rtpthlp =   -2.3917778634313013E-004
 
 PDF parameter set 2:
 wm =    1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.10000000000000001     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =   -4.5000000000000000     
 wp3 =  -0.14230249470757708     
 Skrt =    1.0000000000000000     
 rtp3 =    1.1180339887498949E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =   -5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 PDF parameter set 3:
 wm =   -1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.50000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =    4.5000000000000000     
 wp3 =    1.5909902576697321     
 Skrt =    2.5000000000000000     
 rtp3 =    2.7950849718747374E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 PDF parameter set 4:
 wm =   0.10000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    5.0000000000000000     
 rtp2 =    6.9999999999999999E-006
 thlp2 =    1.0000000000000000     
 Skw =    2.0000000000000000     
 wp3 =    22.360679774997898     
 Skrt =    1.5000000000000000     
 rtp3 =    2.7780388766178202E-008
 Skthl =   -1.5000000000000000     
 thlp3 =   -1.5000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -1.0000000000000000E-002
 rtpthlp =   -2.3811761799581318E-003
 
 PDF parameter set 5:
 wm =   -1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.29999999999999999     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.20000000000000001     
 Skw =   -2.0000000000000000     
 wp3 =  -0.32863353450309962     
 Skrt =  -0.25000000000000000     
 rtp3 =   -7.0710678118654745E-010
 Skthl =   -3.5000000000000000     
 thlp3 =  -0.31304951684997062     
 wprtp =    5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =    5.6920997883030825E-004
 
 PDF parameter set 6:
 wm =   -1.0000000000000000E-002
 rtm =    5.0000000000000001E-003
 thlm =    305.00000000000000     
 wp2 =    1.0000000000000000     
 rtp2 =    4.9999999999999998E-007
 thlp2 =   0.10000000000000001     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -2.0124611797498108E-004
 
 PDF parameter set 7:
 wm =   0.20000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.25000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.14999999999999999     
 Skw =   0.75000000000000000     
 wp3 =    9.3750000000000000E-002
 Skrt =    3.0000000000000000     
 rtp3 =    3.3541019662496846E-008
 Skthl =  -0.50000000000000000     
 thlp3 =   -2.9047375096555625E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -7.7942286340599491E-004
 
 PDF parameter set 8:
 wm =   -2.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    0.0000000000000000     
 thlp2 =    0.0000000000000000     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    0.0000000000000000     
 
 PDF parameter set 9:
 wm =    1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.69999999999999996     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =   -1.0000000000000000     
 rtp3 =   -2.8284271247461898E-009
 Skthl =   -2.0000000000000000     
 thlp3 =   -1.1713240371477056     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    1.0648943609579309E-003
 
 PDF parameter set 10 (randomly generated):
 Random seed values =         2022           5          14        -300           3          29          48         615        2637         620         629         315         618         644         663        1230        2637         620         629         315         618         644         663        1230        2637         620         629         315         618         644         663        1230        2637
 wm =   -8.1656470066179709E-002
 rtm =    1.3223577929842889E-002
 thlm =    306.68128248592961     
 wp2 =   0.58788318054315858     
 rtp2 =    1.3964194782885902E-005
 thlp2 =    2.4141300816848918     
 Skw =   -3.6953600333424723     
 wp3 =   -1.6656866695952401     
 Skrt =   -4.2910631464440607     
 rtp3 =   -2.2391786914949338E-007
 Skthl =   -3.0068017680305861     
 thlp3 =   -11.278355988124812     
 wprtp =   -1.7430600997408699E-003
 wpthlp =    2.1682368338452503E-002
 rtpthlp =   -5.2255353773959260E-003
 
 Success!
 
 
 Performing PDF parameter values unit test
 =========================================
 
 Performing PDF parameter unit tests for the new PDF
 
 PDF parameter set 1:
 wm =   -5.0000000000000001E-003
 rtm =    4.2868400000000001E-003
 thlm =    309.13099999999997     
 wp2 =    3.3741899999999998E-002
 rtp2 =    1.1786499999999999E-006
 thlp2 =    5.9920000000000001E-002
 Skw =    8.6882480819056749     
 wp3 =    5.3850099999999998E-002
 Skrt =    1.4880899999999999     
 rtp3 =    1.9041733779221134E-009
 Skthl =   -1.4580900000000001     
 thlp3 =   -2.1386614365822010E-002
 wprtp =    2.3725400000000002E-005
 wpthlp =   -5.5512699999999996E-003
 rtpthlp =   -2.3917778634313013E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 2:
 wm =    1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.10000000000000001     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =   -4.5000000000000000     
 wp3 =  -0.14230249470757708     
 Skrt =    1.0000000000000000     
 rtp3 =    1.1180339887498949E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =   -5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 3:
 wm =   -1.0000000000000000E-002
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.50000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.10000000000000001     
 Skw =    4.5000000000000000     
 wp3 =    1.5909902576697321     
 Skrt =    2.5000000000000000     
 rtp3 =    2.7950849718747374E-008
 Skthl =   -2.5000000000000000     
 thlp3 =   -7.9056941504209499E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -6.3639610306789288E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 4:
 wm =   0.10000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    5.0000000000000000     
 rtp2 =    6.9999999999999999E-006
 thlp2 =    1.0000000000000000     
 Skw =    2.0000000000000000     
 wp3 =    22.360679774997898     
 Skrt =    1.5000000000000000     
 rtp3 =    2.7780388766178202E-008
 Skthl =   -1.5000000000000000     
 thlp3 =   -1.5000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -1.0000000000000000E-002
 rtpthlp =   -2.3811761799581318E-003
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 5:
 wm =   -1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.29999999999999999     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.20000000000000001     
 Skw =   -2.0000000000000000     
 wp3 =  -0.32863353450309962     
 Skrt =  -0.25000000000000000     
 rtp3 =   -7.0710678118654745E-010
 Skthl =   -3.5000000000000000     
 thlp3 =  -0.31304951684997062     
 wprtp =    5.0000000000000002E-005
 wpthlp =    5.0000000000000001E-003
 rtpthlp =    5.6920997883030825E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 6:
 wm =   -1.0000000000000000E-002
 rtm =    5.0000000000000001E-003
 thlm =    305.00000000000000     
 wp2 =    1.0000000000000000     
 rtp2 =    4.9999999999999998E-007
 thlp2 =   0.10000000000000001     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -2.0124611797498108E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 7:
 wm =   0.20000000000000001     
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =   0.25000000000000000     
 rtp2 =    5.0000000000000004E-006
 thlp2 =   0.14999999999999999     
 Skw =   0.75000000000000000     
 wp3 =    9.3750000000000000E-002
 Skrt =    3.0000000000000000     
 rtp3 =    3.3541019662496846E-008
 Skthl =  -0.50000000000000000     
 thlp3 =   -2.9047375096555625E-002
 wprtp =    5.0000000000000002E-005
 wpthlp =   -5.0000000000000001E-003
 rtpthlp =   -7.7942286340599491E-004
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 8:
 wm =   -2.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    0.0000000000000000     
 thlp2 =    0.0000000000000000     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =    0.0000000000000000     
 rtp3 =    0.0000000000000000     
 Skthl =    0.0000000000000000     
 thlp3 =    0.0000000000000000     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    0.0000000000000000     
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 9:
 wm =    1.0000000000000000E-003
 rtm =    1.0000000000000000E-002
 thlm =    305.00000000000000     
 wp2 =    0.0000000000000000     
 rtp2 =    1.9999999999999999E-006
 thlp2 =   0.69999999999999996     
 Skw =    0.0000000000000000     
 wp3 =    0.0000000000000000     
 Skrt =   -1.0000000000000000     
 rtp3 =   -2.8284271247461898E-009
 Skthl =   -2.0000000000000000     
 thlp3 =   -1.1713240371477056     
 wprtp =    0.0000000000000000     
 wpthlp =    0.0000000000000000     
 rtpthlp =    1.0648943609579309E-003
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 PDF parameter set 10 (randomly generated):
 Random seed values =         2022           5          14        -300           3          29          48         634        2656         639         648         334         637         663         682        1268        2656         639         648         334         637         663         682        1268        2656         639         648         334         637         663         682        1268        2656
 wm =   0.15533664363192701     
 rtm =    1.5234077663377412E-002
 thlm =    302.82370489768192     
 wp2 =   0.68346999441464684     
 rtp2 =    5.2756890179689701E-005
 thlp2 =    2.3469540417162582     
 Skw =    1.2622467298625291     
 wp3 =   0.71321986966861384     
 Skrt =   0.20166567728766915     
 rtp3 =    7.7277091967535028E-008
 Skthl =    1.1260208419835247     
 thlp3 =    4.0485854226347850     
 wprtp =   -1.4282168522961945E-004
 wpthlp =   0.30524507853248506     
 rtpthlp =   -1.0014618179183810E-002
 
 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 =         2022           5          14        -300           3          29          48         654        2676         659         668         354         657         683         702        1308        2676         659         668         354         657         683         702        1308        2676         659         668         354         657         683         702        1308        2676
 wm =  -0.33496421707113988     
 rtm =    1.8983969873240230E-002
 thlm =    300.64360118752199     
 wp2 =   0.45419477467874492     
 rtp2 =    4.8452581345960511E-006
 thlp2 =   0.55324963092594270     
 Skw =    3.2733464444615477     
 wp3 =    1.0019710483531068     
 Skrt =    2.4309565840031828     
 rtp3 =    2.5927018740641533E-008
 Skthl =   0.55092936059864728     
 thlp3 =   0.22671362292815567     
 wprtp =    5.7697599724967756E-004
 wpthlp =  -0.32849366071170605     
 rtpthlp =   -1.4735386633923152E-003
 
 Running tests for the above parameter set for all combinations of F_w and zeta_w for the setting variable.  F_w values are 0 (or 1.0 x 10^-5), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.  Zeta_w values are 0, -1/2, 1, -1/3, 1/2, -1/5, 1/4, -9/10, 9, -3/4, and 3.
 
 Running tests for the above parameter set for the full PDF (the setting of F and zeta values is handled internally).
 
 Success!
 

 Performing spurious source unit test
 ====================================

 Test            1
 
 Vertical integral of <rt> before advance_xm_wpxp =    28.316670259590918     
 Vertical integral of <rt> after advance_xm_wpxp =    28.454761250379878     
 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> =    2.4909544912854464E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    793918.42642888101     
 Vertical integral of <thl> after advance_xm_wpxp =    793908.03522151255     
 Vertical integral of <thl> forcing =   -4.4637357896125940E-002
 Flux of <thl> at the surface =    1.0000000000000000E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    1.2514503322513804E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            2
 
 Vertical integral of <rt> before advance_xm_wpxp =    33.985374857862212     
 Vertical integral of <rt> after advance_xm_wpxp =    34.100497221330777     
 Vertical integral of <rt> forcing =   -5.3171450884365342E-005
 Flux of <rt> at the surface =    4.3691266244624062E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    7.9621097041904232E-018
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    795868.31370180077     
 Vertical integral of <thl> after advance_xm_wpxp =    795852.87636945362     
 Vertical integral of <thl> forcing =   -6.0196027737386831E-002
 Flux of <thl> at the surface =    8.7382532489248127E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -2.0440593662129913E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            3
 
 Vertical integral of <rt> before advance_xm_wpxp =    47.299184187250987     
 Vertical integral of <rt> after advance_xm_wpxp =    47.549037354259163     
 Vertical integral of <rt> forcing =   -6.8072410655938367E-005
 Flux of <rt> at the surface =    9.0091630068319382E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -1.9515639104739080E-018
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    803002.87671383703     
 Vertical integral of <thl> after advance_xm_wpxp =    802986.52396187163     
 Vertical integral of <thl> forcing =   -7.2527499230006037E-002
 Flux of <thl> at the surface =    1.8018326013663877E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -1.6634887911592955E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            4
 
 Vertical integral of <rt> before advance_xm_wpxp =    32.429623876547488     
 Vertical integral of <rt> after advance_xm_wpxp =    32.601659691674527     
 Vertical integral of <rt> forcing =   -5.2365219016074127E-005
 Flux of <rt> at the surface =    6.2581793610621498E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -1.1160506113022661E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    795191.62663351547     
 Vertical integral of <thl> after advance_xm_wpxp =    795180.17861771386     
 Vertical integral of <thl> forcing =   -5.0676411394001600E-002
 Flux of <thl> at the surface =    1.2516358722124300E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -1.5124013152956195E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            5
 
 Vertical integral of <rt> before advance_xm_wpxp =    30.483558549529047     
 Vertical integral of <rt> after advance_xm_wpxp =    30.757931418701155     
 Vertical integral of <rt> forcing =   -3.7018567165692812E-005
 Flux of <rt> at the surface =    9.5159479773940302E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -1.6432439176733427E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    795839.32450827875     
 Vertical integral of <thl> after advance_xm_wpxp =    795835.15876627329     
 Vertical integral of <thl> forcing =   -3.2917702637392646E-002
 Flux of <thl> at the surface =    1.9031895954788061E-002
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -2.2433374602393030E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Success!
 
 Running tuner tests
 -------------------
 goldstein_price pass
 rastrigin pass
 himmelblau pass
 eggholder pass
 schaffer pass
 ==================================================
  
 Performing 'mean vertical velocity in clouds' tests
  
 Success!
  
 ==================================================
  
 ==================================================
  
 Performing 'Smooth Heaviside Peskin' tests
  
 Case 1: heaviside_smth_range = 0, smooth Heaviside collapses to Heaviside
 Input:   -1.0000000000000000       -1.0000000000000000E-010   1.0000000000000000E-010   1.0000000000000000     
 Expected outcome:    0.0000000000000000        0.0000000000000000        1.0000000000000000        1.0000000000000000     
 True outcome:        0.0000000000000000        0.0000000000000000        1.0000000000000000        1.0000000000000000     
 
Case 2: heaviside_smth_range =    2.0000000000000000     
 Input:   -2.0000000001000000       -1.0000000000000000        1.0000000000000000        2.0000000001000000     
 Expected outcome:    0.0000000000000000        9.0845056908100005E-002  0.90915494309000000        1.0000000000000000     
 True outcome:        0.0000000000000000        9.0845056909999999E-002  0.90915494309000000        1.0000000000000000     
 Success!
  
 ==================================================
  
 ==================================================
  
 Performing 'Smooth Min / Max' tests
  
 Input:   -1.0000000000000000       -1.0000000000000000E-010   1.0000000000000000E-010   1.0000000000000000     
 (Simple) min:
 Expected outcome:   -1.0000000000000000       -1.0000000000000000E-010   0.0000000000000000        0.0000000000000000     
 True outcome:       -1.0000000000000000       -1.0000000000000000E-010   0.0000000000000000        0.0000000000000000      

 (Simple) max:
 Expected outcome:    0.0000000000000000        0.0000000000000000        1.0000000000000000E-010   1.0000000000000000     
 True outcome:        0.0000000000000000        0.0000000000000000        1.0000000000000000E-010   1.0000000000000000      


 Input:   -2.0000000001000000       -1.0000000000000000        1.0000000000000000        2.0000000001000000     
 Smooth min:
 Expected outcome:   -2.2247448714824141       -1.3660254037844388      -0.36602540378443871      -0.22474487138241406     
 True outcome:       -2.2247448714824141       -1.3660254037844388      -0.36602540378443871      -0.22474487138241406     
 Smooth max:
 Expected outcome:   0.22474487138241406       0.36602540378443871        1.3660254037844388        2.2247448714824141     
 True outcome:       0.22474487138241406       0.36602540378443871        1.3660254037844388        2.2247448714824141      


 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!
  
 ==================================================
  
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