diff --git a/physics/rrtmgp_lw_gas_optics.F90 b/physics/rrtmgp_lw_gas_optics.F90
index 45902b967..12fcb4552 100644
--- a/physics/rrtmgp_lw_gas_optics.F90
+++ b/physics/rrtmgp_lw_gas_optics.F90
@@ -10,9 +10,6 @@ module rrtmgp_lw_gas_optics
   use rrtmgp_aux,            only: check_error_msg
   use netcdf
 
-  ! Parameters
-
-
 contains
 
   ! #########################################################################################
@@ -453,14 +450,6 @@ subroutine rrtmgp_lw_gas_optics_run(Model, lw_gas_props, ncol, p_lay, p_lev, t_l
     call check_error_msg('rrtmgp_lw_gas_optics_run',sources%init(lw_gas_props))
     call check_error_msg('rrtmgp_lw_gas_optics_run',sources%alloc(ncol, Model%levs))
 
-    ! Print some diagnostics to test MPI
-    write(*,*) "lw_gas_props%get_ngas():      ",   lw_gas_props%get_ngas()
-    write(*,*) "lw_gas_props%get_gases():     ",   lw_gas_props%get_gases()
-    write(*,*) "lw_gas_props%get_press_min(): ",   lw_gas_props%get_press_min()
-    write(*,*) "lw_gas_props%get_press_max(): ",   lw_gas_props%get_press_max()
-    write(*,*) "lw_gas_props%get_temp_min():  ",   lw_gas_props%get_temp_min()
-    write(*,*) "lw_gas_props%get_temp_max():  ",   lw_gas_props%get_temp_max()
-
     ! Gas-optics 
     call check_error_msg('rrtmgp_lw_gas_optics_run',lw_gas_props%gas_optics(&
          p_lay,                   & ! IN  - Pressure @ layer-centers (Pa)
@@ -472,6 +461,8 @@ subroutine rrtmgp_lw_gas_optics_run(Model, lw_gas_props, ncol, p_lay, p_lev, t_l
          sources,                 & ! OUT - RRTMGP DDT: source functions
          tlev=t_lev))               ! IN  - Temperature @ layer-interfaces (K) (optional)
 
+    write(*,*) "lw_optical_props_clrsky(gas_optics): ",lw_optical_props_clrsky%tau(:,:,1)
+
   end subroutine rrtmgp_lw_gas_optics_run
 
   ! #########################################################################################
diff --git a/physics/rrtmgp_lw_rte.F90 b/physics/rrtmgp_lw_rte.F90
index 7265da45b..9864da635 100644
--- a/physics/rrtmgp_lw_rte.F90
+++ b/physics/rrtmgp_lw_rte.F90
@@ -109,6 +109,8 @@ subroutine rrtmgp_lw_rte_run(Model, Statein, Interstitial, ncol, lw_gas_props, p
     !   flux_allsky%bnd_flux_dn => fluxLWBB_dn_allsky
     !endif
 
+    write(*,*) "lw_optical_props_clrsky(rte_lw): ",lw_optical_props_clrsky%tau(:,:,1)
+
     ! Compute clear-sky fluxes (if requested)
     ! Clear-sky fluxes are gas+aerosol
     call check_error_msg('rrtmgp_lw_rte_run',lw_optical_props_aerosol%increment(lw_optical_props_clrsky))