From 5812151340d1263f1a7e0372dc5fb3b1c816e35c Mon Sep 17 00:00:00 2001
From: Dustin Swales <dustin.swales@noaa.gov>
Date: Tue, 7 Jan 2020 11:21:44 -0700
Subject: [PATCH] Added some more diagnostics.

---
 physics/radlw_main.f               |  4 ++++
 physics/rrtmgp_lw_cloud_optics.F90 |  2 +-
 physics/rrtmgp_lw_rte.F90          | 10 +++++++++-
 3 files changed, 14 insertions(+), 2 deletions(-)

diff --git a/physics/radlw_main.f b/physics/radlw_main.f
index d5c7eb460..963372526 100644
--- a/physics/radlw_main.f
+++ b/physics/radlw_main.f
@@ -1183,6 +1183,10 @@ subroutine rrtmg_lw_run                                           &
            write(47,"(23f8.2)")  plyr(iplon,k),clwp(k),relw(k),ciwp(k), &
      &          reiw(k),taucld(:,k),lon(iplon),lat(iplon)
         enddo
+        do k=1,nlay+1
+           write(47,"(5f8.2)") plvl(iplon,k),totuflux(k-1),totdflux(k-1)&
+     &         ,totuclfl(k-1),totdclfl(k-1)
+        enddo
         
 
         topflx(iplon)%upfxc = totuflux(nlay)
diff --git a/physics/rrtmgp_lw_cloud_optics.F90 b/physics/rrtmgp_lw_cloud_optics.F90
index e6cff091b..58f3eb9dd 100644
--- a/physics/rrtmgp_lw_cloud_optics.F90
+++ b/physics/rrtmgp_lw_cloud_optics.F90
@@ -447,8 +447,8 @@ subroutine rrtmgp_lw_cloud_optics_run(doLWrad, nCol, nLev, cld_optics_scheme, nr
                cld_reliq, cld_iwp, cld_reice, cld_rwp, cld_rerain, cld_swp, cld_resnow, &
                cld_frac, tau_cld)
        endif
-       lw_optical_props_cloudsByBand%tau = tau_cld
     endif
+    lw_optical_props_cloudsByBand%tau = tau_cld
     
     write(47,*) "In rrtmgp_lw_cloud_optics: "
     write(47,*),"nCol: ",nCol
diff --git a/physics/rrtmgp_lw_rte.F90 b/physics/rrtmgp_lw_rte.F90
index a4989237e..568b48ad7 100644
--- a/physics/rrtmgp_lw_rte.F90
+++ b/physics/rrtmgp_lw_rte.F90
@@ -76,7 +76,7 @@ subroutine rrtmgp_lw_rte_run(doLWrad, nCol, nLev, p_lay, t_lay, p_lev, skt, lw_g
 
     ! Local variables
     integer :: &
-         iCol, iBand
+         iCol, iBand, iLay
     type(ty_fluxes_byband) :: &
          flux_allsky, flux_clrsky
     real(kind_phys), dimension(ncol,nLev+1,lw_gas_props%get_nband()),target :: &
@@ -157,6 +157,14 @@ subroutine rrtmgp_lw_rte_run(doLWrad, nCol, nLev, p_lay, t_lay, p_lev, skt, lw_g
     !if (l_AllSky_HR_byband) then
     !endif
 
+    write(47,*) "In rrtmgp_lw_rte: "
+    do iCol=1,nCol
+       do iLay=1,nLev+1
+          write(47,"(5f8.2)") p_lev(iCol,iLay)/100.,fluxlwUP_allsky(iCol,iLay),fluxlwDOWN_allsky(iCol,iLay),&
+               fluxlwUP_clrsky(iCol,iLay),fluxlwDOWN_clrsky(iCol,iLay)
+       enddo
+    enddo
+
   end subroutine rrtmgp_lw_rte_run
   
   ! #########################################################################################