diff --git a/physics/radlw_main.f b/physics/radlw_main.f
index c34a1f77d..32dafd7e4 100644
--- a/physics/radlw_main.f
+++ b/physics/radlw_main.f
@@ -362,11 +362,7 @@ subroutine rrtmg_lw_init (mpicomm, mpirank, mpiroot)
           ! Inputs
           integer, intent(in) :: mpicomm,mpirank,mpiroot
 
-          if (mpirank .eq. mpiroot) then
-             print*,'DJS+ Opening file containing RRTMG LW cloud fields'
-             open(47,file='rrtmg_clds.txt',status='unknown')
-          endif
-
+        
         end subroutine rrtmg_lw_init
 
 !> \defgroup module_radlw_main GFS RRTMG Longwave Module 
@@ -730,10 +726,6 @@ subroutine rrtmg_lw_run                                           &
 !     endif
 
 !  --- ...  loop over horizontal npts profiles
-      write(47,*) "In radlw_main: "
-      write(47,*),"nCol: ",npts
-      write(47,*),"nLay: ",nLay
-
       lab_do_iplon : do iplon = 1, npts
 
 !> -# Read surface emissivity.
@@ -1179,16 +1171,6 @@ subroutine rrtmg_lw_run                                           &
         endif   ! end if_isubclw_block
 
 !> -# Save outputs.
-
-
-        do k=1,nlay
-           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)
@@ -1287,7 +1269,6 @@ subroutine rrtmg_lw_finalize (mpicomm, mpirank, mpiroot)
 #ifdef MPI
         call MPI_BARRIER(mpicomm, ierr)
 #endif
-        close(47)
       end subroutine rrtmg_lw_finalize
 
 
@@ -3706,12 +3687,6 @@ subroutine rtrnmc                                                 &
         totdclfl(k) = totdclfl(k) * flxfac
       enddo
 
-
-      do k=0,nlay
-         write(47,"(32f8.2)") toturad(k,:)*flxfac,                      &
-     &        totdrad(k,:)*flxfac
-      enddo
-
 !> -# Calculate net fluxes and heating rates.
       fnet(0) = totuflux(0) - totdflux(0)
 
diff --git a/physics/rrtmgp_lw_cloud_optics.F90 b/physics/rrtmgp_lw_cloud_optics.F90
index 7e546acd0..7e50f2a81 100644
--- a/physics/rrtmgp_lw_cloud_optics.F90
+++ b/physics/rrtmgp_lw_cloud_optics.F90
@@ -99,11 +99,6 @@ subroutine rrtmgp_lw_cloud_optics_init(cld_optics_scheme, nrghice, rrtmgp_root_d
     errmsg = ''
     errflg = 0
 
-    if (mpirank .eq. mpiroot) then
-       print*,'DJS+ Opening file containing RRTMGP LW cloud fields'
-       open(47,file='rrtmgp_clds.txt',status='unknown')
-    endif
-
     if (cld_optics_scheme .eq. 0) return
 
     ! Filenames are set in the physics_nml
@@ -450,19 +445,7 @@ subroutine rrtmgp_lw_cloud_optics_run(doLWrad, nCol, nLev, cld_optics_scheme, nr
        endif
        lw_optical_props_cloudsByBand%tau = tau_cld
     endif    
-    
-    write(47,*) "In rrtmgp_lw_cloud_optics: "
-    write(47,*),"nCol: ",nCol
-    write(47,*),"nLay: ",nLev
-    do iCol=1,nCol
-       do iLay=1,nLev
-          write(47,"(23f8.2)") p_lay(iCol,iLay)/100.,cld_lwp(iCol,iLay),cld_reliq(iCol,iLay),&
-               cld_iwp(iCol,iLay),cld_reice(iCol,iLay),tau_cld(iCol,iLay,:),lon(iCol), lat(iCol)
-       enddo
-    enddo
-
-    !lw_optical_props_cloudsByBand%tau(:,:,2) = lw_optical_props_cloudsByBand%tau(:,:,1)
-
+ 
     ! All-sky LW optical depth ~10microns
     cldtaulw = lw_optical_props_cloudsByBand%tau(:,:,7)
 
@@ -486,7 +469,6 @@ subroutine rrtmgp_lw_cloud_optics_finalize(mpicomm, mpirank, mpiroot)
 #ifdef MPI
     call MPI_BARRIER(mpicomm, ierr)
 #endif
-    close(47)
 
   end subroutine rrtmgp_lw_cloud_optics_finalize
 end module rrtmgp_lw_cloud_optics
diff --git a/physics/rrtmgp_lw_cloud_sampling.F90 b/physics/rrtmgp_lw_cloud_sampling.F90
index 45a1811f1..af5f77210 100644
--- a/physics/rrtmgp_lw_cloud_sampling.F90
+++ b/physics/rrtmgp_lw_cloud_sampling.F90
@@ -109,6 +109,7 @@ subroutine rrtmgp_lw_cloud_sampling_run(doLWrad, nCol, nLev, ipsdlw0, icseed_lw,
 
     ! Call McICA
     select case ( iovrlw )
+    case(0)
        ! Maximumn-random 
     case(1)
        call check_error_msg('rrtmgp_lw_cloud_sampling_run',sampled_mask_max_ran(rng3D,cld_frac_noSamp,cldfracMCICA))       
diff --git a/physics/rrtmgp_lw_rte.F90 b/physics/rrtmgp_lw_rte.F90
index 3f2651e79..cf0d0ea42 100644
--- a/physics/rrtmgp_lw_rte.F90
+++ b/physics/rrtmgp_lw_rte.F90
@@ -160,14 +160,6 @@ 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,"(35f8.2)") p_lev(iCol,iLay)/100.,flux_allsky%bnd_flux_up(iCol,iLay,:),flux_allsky%bnd_flux_dn(iCol,iLay,:),&
-               fluxlwUP_clrsky(iCol,iLay),fluxlwDOWN_clrsky(iCol,iLay)
-       enddo
-    enddo
-
   end subroutine rrtmgp_lw_rte_run
   
   ! #########################################################################################