Skip to content
Success

Console Output

Skipping 206 KB.. Full Log
../src/KK_microphys_module.F90:43:35:

   43 |   subroutine KK_local_microphys( gr, dt, nz, l_latin_hypercube,             & ! In
      |                                   1
Warning: Unused dummy argument ‘gr’ at (1) [-Wunused-dummy-argument]
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/input_fields.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	../src/bugsrad_driver.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	../src/sounding.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/estimate_scm_microphys_module.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/lh_microphys_driver_module.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	../src/microphys_driver.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	../src/clubb_driver.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	../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/jacobian.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 -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
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
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 
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
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_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 -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
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 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 
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/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]
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 -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
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/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
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
../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/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.
664 :         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.
679 :         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.
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.
1335 :                                                             sqrt( (ddzt( nz, ngrdcol, gr, um ))**2 & 
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1356 :             / max( ( ddzt( nz, ngrdcol, gr, um)**2 + ddzt(nz, ngrdcol, gr, vm)**2 ), 1.0e-7_core_rknd ) )
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1391 :                             * sqrt(em(i,k)) / ( gr(i)%zm(k) - sfc_elevation(i) + z_displace )  ! small
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
1395 :       invrs_tau_xp2_zm = min( max( sqrt( ( ddzt( nz, ngrdcol, gr, um)**2 + ddzt(nz, ngrdcol, gr, vm)**2 ) &
CLUBBStandardsCheck.pl WARNING: File has lines that exceed 100 characters.
../src/CLUBB_core/mixing_length.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.
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.
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.
2000 :                                     wp2(i,:), up2(i,:), vp2(i,:), um_forcing(i,:), vm_forcing(i,:), & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2012 :                                     stats_zt(i), stats_zm(i), stats_sfc(i),         & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2013 :                                     um(i,:), vm(i,:), edsclrm(i,:,:),               & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2014 :                                     upwp(i,:), vpwp(i,:), wpedsclrp(i,:,:),         & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2015 :                                     um_pert(i,:), vm_pert(i,:), upwp_pert(i,:), vpwp_pert(i,:) ) ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2062 :         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.
2063 :                           wpthlp(i,:), wprtp2(i,:), wpthlp2(i,:), rho_ds_zm(i,:),         & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2065 :                           sclrm(i,:,:), sclrp2(i,:,:), wpsclrp(i,:,:), wpsclrp2(i,:,:),           & ! Intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2068 :                           rtp3(i,:), thlp3(i,:), sclrp3(i,:,:) )                         ! Intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2122 :       thlp2_zt(:,:)  = max( zm2zt( nz, ngrdcol, gr, thlp2(:,:) ), thl_tol**2 ) ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2123 :       rtp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, rtp2(:,:) ), rt_tol**2 )   ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2412 :                gr(i), um(i,:), vm(i,:), upwp(i,:), vpwp(i,:), up2(i,:), vp2(i,:),                      & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2413 :                thlm(i,:), rtm(i,:), wprtp(i,:), wpthlp(i,:),                              & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2414 :                wp2(i,:), wp3(i,:), rtp2(i,:), rtp3(i,:), thlp2(i,:), thlp3(i,:), rtpthlp(i,:),           & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2415 :                wpthvp(i,:), wp2thvp(i,:), rtpthvp(i,:), thlpthvp(i,:),                    & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2416 :                p_in_Pa(i,:), exner(i,:), rho(i,:), rho_zm(i,:),                           & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2417 :                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.
2418 :                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.
2419 :                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.
2422 :                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.
2423 :                Lscale_up(i,:), Lscale_down(i,:), tau_zt(i,:), Kh_zt(i,:), wp2rcp(i,:),         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2424 :                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.
2425 :                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.
2426 :                vp2_zt(i,:), upwp_zt(i,:), vpwp_zt(i,:), wpup2(i,:), wpvp2(i,:),                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2429 :                rtprcp(i,:), rcp2(i,:), em(i,:), a3_coef(i,:), a3_coef_zt(i,:),                 & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2430 :                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.
2432 :                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.
2433 :                sclrprtp(i,:,:), sclrpthlp(i,:,:), sclrm_forcing(i,:,:), sclrpthvp(i,:,:),         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2434 :                wpsclrp(i,:,:), sclrprcp(i,:,:), wp2sclrp(i,:,:), wpsclrp2(i,:,:), wpsclrprtp(i,:,:),     & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2435 :                wpsclrpthlp(i,:,:), wpedsclrp(i,:,:), edsclrm(i,:,:), edsclrm_forcing(i,:,:),      & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2443 :              gr(i), thlm_forcing(i,:), rtm_forcing(i,:), um_forcing(i,:),                         & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2444 :              vm_forcing(i,:), wm_zm(i,:), wm_zt(i,:), p_in_Pa(i,:),                                 & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2445 :              rho_zm(i,:), rho(i,:), exner(i,:), rho_ds_zm(i,:),                                     & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2446 :              rho_ds_zt(i,:), invrs_rho_ds_zm(i,:), invrs_rho_ds_zt(i,:),                       & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2447 :              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.
2448 :              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.
2449 :              rtm(i,:), wprtp(i,:), thlm(i,:), wpthlp(i,:), wp2(i,:), wp3(i,:),                                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2450 :              rtp2(i,:), thlp2(i,:), rtpthlp(i,:),                                              & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2453 :              wpsclrp_sfc(i,:), wpedsclrp_sfc(i,:), sclrm(i,:,:), wpsclrp(i,:,:), sclrp2(i,:,:),                & ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
2454 :              sclrprtp(i,:,:), sclrpthlp(i,:,:), sclrm_forcing(i,:,:), edsclrm(i,:,:), edsclrm_forcing(i,:,:)       ) ! intent(in)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3024 :     wp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, wp2(:,:) ), w_tol_sqd ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3026 :     thlp2_zt(:,:) = max( zm2zt( nz, ngrdcol, gr, thlp2(:,:) ), thl_tol**2 ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3028 :     rtp2_zt(:,:)  = max( zm2zt( nz, ngrdcol, gr, rtp2(:,:) ), rt_tol**2 ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3030 :     up2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, up2(:,:) ), w_tol_sqd ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3032 :     vp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, vp2(:,:) ), w_tol_sqd ) ! Positive definite quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3036 :       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.
3166 :     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.
3173 :     rtp2_zt(:,:)    = max( zm2zt( nz, ngrdcol, gr, rtp2(:,:) ), rt_tol**2 )   ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3174 :     thlp2_zt(:,:)   = max( zm2zt( nz, ngrdcol, gr, thlp2(:,:) ), thl_tol**2 ) ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3175 :     up2_zt(:,:)     = max( zm2zt( nz, ngrdcol, gr, up2(:,:) ), w_tol_sqd )    ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3176 :     vp2_zt(:,:)     = max( zm2zt( nz, ngrdcol, gr, vp2(:,:) ), w_tol_sqd )    ! Positive def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3200 :       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.
3409 :         rcp2(:,:) = max( zt2zm( nz, ngrdcol, gr, rcp2_zt(:,:) ), zero_threshold )  ! Pos. def. quantity
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
3481 :                rcm(i,:), wp2thvp(i,:), wpsclrprtp(i,:,:), wpsclrp2(i,:,:),          & ! intent(inout)
CLUBBStandardsCheck.pl WARNING: Line has exceeded 100 characters.
4782 :                        ( ( 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.
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.
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.
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 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.
115 :     result = smooth_min(4, 1, input, zero, smth_coef) ! Order of nz and ngrdcol is opposite to order of indexing! This is very unintuitive!
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           4        -300           3          31          29         576        2598         581         580         276         579         607         605        1152        2598         581         580         276         579         607         605        1152        2598         581         580         276         579         607         605        1152        2598
 hmm =    2.6226263990318500E-005
 mixt_frac =   0.37698867721339019     
 precip_frac =   0.88153278362204268     
 precip_frac_1 =   0.83316517708991833     
 precip_frac_2 =   0.91080037369563205     
 hmp2_ip_on_hmm2_ip =    1.8360213029926327     
 hmm_ip =    2.9750752867703914E-005
 hmp2 =    1.5249918144773755E-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           4        -300           3          31          29         579        2601         584         583         279         582         610         608        1158        2601         584         583         279         582         610         608        1158        2601         584         583         279         582         610         608        1158        2601
 wm =  -0.38446683938453019     
 rtm =    3.7524562539477460E-003
 thlm =    303.46718168330068     
 wp2 =   0.95786644163099044     
 rtp2 =    6.9729682080903355E-007
 thlp2 =    2.9616950877639718     
 Skw =    2.5085301906517712     
 wp3 =    2.3516721151523083     
 Skrt =   -1.4263720155888677     
 rtp3 =   -8.3053767425746374E-010
 Skthl =    3.7587591291675420     
 thlp3 =    19.158213707748349     
 wprtp =   -5.5822607499803533E-004
 wpthlp =  -0.65794920846699390     
 rtpthlp =    1.2933662516131598E-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           4        -300           3          31          29         579        2601         584         583         279         582         610         608        1158        2601         584         583         279         582         610         608        1158        2601         584         583         279         582         610         608        1158        2601
 wm =  -0.38446683938453019     
 rtm =    3.7524562539477460E-003
 thlm =    303.46718168330068     
 wp2 =   0.95786644163099044     
 rtp2 =    6.9729682080903355E-007
 thlp2 =    2.9616950877639718     
 Skw =    2.5085301906517712     
 wp3 =    2.3516721151523083     
 Skrt =   -1.4263720155888677     
 rtp3 =   -8.3053767425746374E-010
 Skthl =    3.7587591291675420     
 thlp3 =    19.158213707748349     
 wprtp =   -5.5822607499803533E-004
 wpthlp =  -0.65794920846699390     
 rtpthlp =    1.2933662516131598E-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           4        -300           3          31          29         584        2606         589         588         284         587         615         613        1168        2606         589         588         284         587         615         613        1168        2606         589         588         284         587         615         613        1168        2606
 wm =  -0.29831119246933535     
 rtm =    2.0892192309480643E-002
 thlm =    297.83947505837972     
 wp2 =    5.0605775304216394E-002
 rtp2 =    7.7815645837505952E-005
 thlp2 =    2.8060175317470137     
 Skw =   0.92238755722830967     
 wp3 =    1.0500586577951868E-002
 Skrt =  -0.87336062753336041     
 rtp3 =   -5.9950664157548127E-007
 Skthl =   -2.1428992375540092     
 thlp3 =   -10.072501031204647     
 wprtp =    1.2094693188474585E-003
 wpthlp =   0.15226537567260756     
 rtpthlp =    1.3299066652816249E-002
 
 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           4        -300           3          31          29         606        2628         611         610         306         609         637         635        1212        2628         611         610         306         609         637         635        1212        2628         611         610         306         609         637         635        1212        2628
 wm =   0.34585218158503550     
 rtm =    1.7370813430817626E-003
 thlm =    291.57678534027326     
 wp2 =   0.63144506646216070     
 rtp2 =    1.2510345311114163E-008
 thlp2 =   0.58183544133497300     
 Skw =    4.4916308784005894     
 wp3 =    2.2537587218737944     
 Skrt =   -4.2104323097705141     
 rtp3 =   -5.8915644859046640E-012
 Skthl =  -0.27281342889093385     
 thlp3 =  -0.12107821282718242     
 wprtp =   -6.3360662879617181E-005
 wpthlp =   0.51311807238757745     
 rtpthlp =   -7.6785151241170755E-005
 
 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           4        -300           3          31          29         626        2648         631         630         326         629         657         655        1252        2648         631         630         326         629         657         655        1252        2648         631         630         326         629         657         655        1252        2648
 wm =  -0.38862485223264809     
 rtm =    1.6272644748802335E-002
 thlm =    290.06061014835421     
 wp2 =   0.76800992557015924     
 rtp2 =    5.6736190373187382E-005
 thlp2 =   0.49739113969952520     
 Skw =   -4.0041780544386452     
 wp3 =   -2.6950301632719333     
 Skrt =    1.4898048703751730     
 rtp3 =    6.3667771133137328E-007
 Skthl =   -1.7696009440263494     
 thlp3 =  -0.62075811878431675     
 wprtp =   -6.2348696333160312E-003
 wpthlp =    3.8458373076872242E-002
 rtpthlp =   -4.7810316352710326E-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 =    37.689405075602664     
 Vertical integral of <rt> after advance_xm_wpxp =    37.766563189249773     
 Vertical integral of <rt> forcing =   -6.3764422484994575E-005
 Flux of <rt> at the surface =    3.2095813464191750E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    1.0526247842118641E-016
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    798154.94750143215     
 Vertical integral of <thl> after advance_xm_wpxp =    798135.06086753064     
 Vertical integral of <thl> forcing =   -7.2707942365705810E-002
 Flux of <thl> at the surface =    6.4191626928383492E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    1.1754763828975001E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            3
 
 Vertical integral of <rt> before advance_xm_wpxp =    35.623248502093965     
 Vertical integral of <rt> after advance_xm_wpxp =    35.704550314483619     
 Vertical integral of <rt> forcing =   -5.8907472478841943E-005
 Flux of <rt> at the surface =    3.2991351377770640E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =   -1.7611509039311413E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    798409.08604601328     
 Vertical integral of <thl> after advance_xm_wpxp =    798390.93532278435     
 Vertical integral of <thl> forcing =   -6.7100681038468668E-002
 Flux of <thl> at the surface =    6.5982702755541280E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -1.8400558854381188E-013
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            4
 
 Vertical integral of <rt> before advance_xm_wpxp =    35.602692586752227     
 Vertical integral of <rt> after advance_xm_wpxp =    35.690519625507612     
 Vertical integral of <rt> forcing =   -6.3964758636434081E-005
 Flux of <rt> at the surface =    3.5672155448770060E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    1.8282359133536819E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    799491.40785954497     
 Vertical integral of <thl> after advance_xm_wpxp =    799473.60748013237     
 Vertical integral of <thl> forcing =   -6.6469029133198548E-002
 Flux of <thl> at the surface =    7.1344310897540125E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =    1.4512835377900046E-012
 
 The spurious source of <thl> is within acceptable limits.
 
 Test            5
 
 Vertical integral of <rt> before advance_xm_wpxp =    37.996833910063842     
 Vertical integral of <rt> after advance_xm_wpxp =    38.102355089217163     
 Vertical integral of <rt> forcing =   -6.5513627892723324E-005
 Flux of <rt> at the surface =    4.1725089173708476E-004
 Flux of <rt> at the top of the domain =    0.0000000000000000     
 Spurious source of <rt> =    4.2053491765281503E-017
 
 The spurious source of <rt> is within acceptable limits.
 
 Vertical integral of <thl> before advance_xm_wpxp =    797938.11762778799     
 Vertical integral of <thl> after advance_xm_wpxp =    797919.47030101134     
 Vertical integral of <thl> forcing =   -7.0502773755517556E-002
 Flux of <thl> at the surface =    8.3450178347416956E-003
 Flux of <thl> at the top of the domain =    0.0000000000000000     
 Spurious source of <thl> =   -1.3622158956394514E-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: 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: 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 smth_coef=1e-7, input_var1 = equidistant grid on [-1, 1), input_var2 = 0

 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