diff --git a/physics/GFS_rrtmgp_lw_post.F90 b/physics/GFS_rrtmgp_lw_post.F90
index 0d5a9c974..b423e4d22 100644
--- a/physics/GFS_rrtmgp_lw_post.F90
+++ b/physics/GFS_rrtmgp_lw_post.F90
@@ -228,13 +228,7 @@ subroutine GFS_rrtmgp_lw_post_run (Model, Grid, Radtend, Coupling, Diag,  Statei
                 Diag%fluxr(i,46-j) = Diag%fluxr(i,46-j) + tem0d * (1.0-exp(-tem2))
              enddo
           enddo
-       endif
-       
-       if (Model%lgocart .or. Model%ldiag3d) then
-          do k = 1, Model%levs
-             Coupling%cldcovi(1:im,k) = cld_frac(1:im,k)
-          enddo
-       endif
+       endif       
     endif
 
   end subroutine GFS_rrtmgp_lw_post_run
diff --git a/physics/GFS_rrtmgp_sw_post.F90 b/physics/GFS_rrtmgp_sw_post.F90
index 1dbfe6843..ed272f530 100644
--- a/physics/GFS_rrtmgp_sw_post.F90
+++ b/physics/GFS_rrtmgp_sw_post.F90
@@ -313,12 +313,6 @@ subroutine GFS_rrtmgp_sw_post_run (Model, Grid, Diag, Radtend, Coupling, Statein
              enddo
           enddo
        endif
-
-       if (Model%lgocart .or. Model%ldiag3d) then
-          do k = 1, Model%levs
-             Coupling%cldcovi(1:im,k) = cld_frac(1:im,k)
-          enddo
-       endif
     endif