diff --git a/tests/rte_lw_solver_unit_tests.F90 b/tests/rte_lw_solver_unit_tests.F90
index f2a834a2f..3f77aa606 100644
--- a/tests/rte_lw_solver_unit_tests.F90
+++ b/tests/rte_lw_solver_unit_tests.F90
@@ -305,10 +305,8 @@ function check_gray_rad_equil(sfc_T, lw_tau, top_at_1, up_flux, net_flux)
     ! Check top-of-atmosphere energy balance 
     !
     if(.not. allclose(up_flux(:,toa), &
-                      gray_rad_equil_olr(sfc_t, lw_tau), tol=4._wp)) then
+                      gray_rad_equil_olr(sfc_t, lw_tau), tol=8._wp)) then
       call report_err("OLR is not consistent with gray radiative equilibrium")
-      print *, maxval((up_flux(:,toa) - &
-                       gray_rad_equil_olr(sfc_t, lw_tau))/spacing(gray_rad_equil_olr(sfc_t, lw_tau)))
       check_gray_rad_equil = .false.
     end if 
     !
@@ -319,8 +317,6 @@ function check_gray_rad_equil(sfc_T, lw_tau, top_at_1, up_flux, net_flux)
     if(.not. allclose(net_flux(:,:), & 
                       spread(net_flux(:,1), dim=2, ncopies=size(net_flux,2)), &
                       tol = 100._wp)) then 
-      print *, maxval((net_flux(:,:) - &
-                       spread(net_flux(:,1), dim=2, ncopies=size(net_flux,2)))/spacing(net_flux(:,:)))
       call report_err("Net flux not constant with tau in gray radiative equilibrium")
       check_gray_rad_equil = .false.
     end if