Added some more diagnostics.

NCAR · Jan 7, 2020 · 5812151 · 5812151
1 parent 203cd4a
commit 5812151
Show file tree

Hide file tree

Showing 3 changed files with 14 additions and 2 deletions.
diff --git 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
@@ -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
@@ -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
 
   ! #########################################################################################