diff --git a/physics/radlw_main.f b/physics/radlw_main.f
index ca953c1c3..34cc389cb 100644
--- a/physics/radlw_main.f
+++ b/physics/radlw_main.f
@@ -390,7 +390,7 @@ subroutine rrtmg_lw_run                                           &
      &       gasvmr_ch4, gasvmr_o2, gasvmr_co, gasvmr_cfc11,            &
      &       gasvmr_cfc12, gasvmr_cfc22, gasvmr_ccl4,                   &
      &       icseed,aeraod,aerssa,sfemis,sfgtmp,                        &
-     &       dzlyr,delpin,de_lgth,                                      &
+     &       dzlyr,delpin,de_lgth,xlon,xlat,                            &
      &       npts, nlay, nlp1, lprnt, cld_cf, lslwr,                    &
      &       hlwc,topflx,sfcflx,cldtau,                                 &   !  ---  outputs
      &       HLW0,HLWB,FLXPRF,                                          &   !  ---  optional
@@ -589,7 +589,7 @@ subroutine rrtmg_lw_run                                           &
      &       cld_od
 
       real (kind=kind_phys), dimension(npts), intent(in) :: sfemis,     &
-     &       sfgtmp, de_lgth
+     &       sfgtmp, de_lgth, xlon, xlat
 
       real (kind=kind_phys), dimension(npts,nlay,nbands),intent(in)::   &
      &       aeraod, aerssa
@@ -1168,8 +1168,8 @@ subroutine rrtmg_lw_run                                           &
 
 
         do k=1,nlay
-           write(47,"(21f8.2)")  plyr(iplon,k),clwp(k),relw(k),ciwp(k), &
-     &          reiw(k),taucld(:,k)
+           write(47,"(23f8.2)")  plyr(iplon,k),clwp(k),relw(k),ciwp(k), &
+     &          reiw(k),taucld(:,k),xlon(k),xlat(k)
         enddo
         
 
diff --git a/physics/radlw_main.meta b/physics/radlw_main.meta
index 73977e5cb..1d8a08fb3 100644
--- a/physics/radlw_main.meta
+++ b/physics/radlw_main.meta
@@ -207,6 +207,24 @@
   kind = kind_phys
   intent = in
   optional = F
+[xlon]
+  standard_name = longitude
+  long_name = longitude
+  units = radians
+  dimensions = (horizontal_dimension)
+  type = real
+  kind = kind_phys
+  intent = in
+  optional = F
+[xlat]
+  standard_name = latitude
+  long_name = latitude
+  units = radians
+  dimensions = (horizontal_dimension)
+  type = real
+  kind = kind_phys
+  intent = in
+  optional = F
 [npts]
   standard_name = horizontal_loop_extent
   long_name = horizontal dimension
diff --git a/physics/rrtmgp_lw_cloud_optics.F90 b/physics/rrtmgp_lw_cloud_optics.F90
index c7f5ee24c..1961066bc 100644
--- a/physics/rrtmgp_lw_cloud_optics.F90
+++ b/physics/rrtmgp_lw_cloud_optics.F90
@@ -359,7 +359,7 @@ end subroutine rrtmgp_lw_cloud_optics_init
 !!
   subroutine rrtmgp_lw_cloud_optics_run(doLWrad, nCol, nLev, cld_optics_scheme, nrghice,    &
        cld_frac, cld_lwp, cld_reliq, cld_iwp, cld_reice, cld_swp, cld_resnow, cld_rwp,      &
-       cld_rerain, p_lay, lw_cloud_props, lw_gas_props,                                     &
+       cld_rerain, p_lay, lw_cloud_props, lw_gas_props, lon, lat,                           &
        cldtaulw, lw_optical_props_cloudsByBand, errmsg, errflg)
     
     ! Inputs
@@ -370,6 +370,9 @@ subroutine rrtmgp_lw_cloud_optics_run(doLWrad, nCol, nLev, cld_optics_scheme, nr
          nLev,             & ! Number of vertical levels
          nrghice,          & ! Number of ice-roughness categories
          cld_optics_scheme   ! Cloud-optics scheme
+    real(kind_phys), dimension(nCol), intent(in) :: &
+         lon,              & ! Longitude
+         lat                 ! Latitude
     real(kind_phys), dimension(ncol,nLev),intent(in) :: &
          p_lay,            & ! Layer pressure (Pa)
          cld_frac,         & ! Total cloud fraction by layer
@@ -449,8 +452,8 @@ subroutine rrtmgp_lw_cloud_optics_run(doLWrad, nCol, nLev, cld_optics_scheme, nr
     write(47,*),"nLay: ",nLev
     do iCol=1,nCol
        do iLay=1,nLev
-          write(47,"(21f8.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,:)
+          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
 
diff --git a/physics/rrtmgp_lw_cloud_optics.meta b/physics/rrtmgp_lw_cloud_optics.meta
index 2582362d1..850602d64 100644
--- a/physics/rrtmgp_lw_cloud_optics.meta
+++ b/physics/rrtmgp_lw_cloud_optics.meta
@@ -226,6 +226,24 @@
   intent = in
   type = ty_cloud_optics
   optional = F
+[lon]
+  standard_name = longitude
+  long_name = longitude
+  units = radians
+  dimensions = (horizontal_dimension)
+  type = real
+  kind = kind_phys
+  intent = in
+  optional = F
+[lat]
+  standard_name = latitude
+  long_name = latitude
+  units = radians
+  dimensions = (horizontal_dimension)
+  type = real
+  kind = kind_phys
+  intent = in
+  optional = F
 [cldtaulw]
   standard_name = RRTMGP_cloud_optical_depth_layers_at_10mu_band
   long_name = approx 10mu band layer cloud optical depth