diff --git a/io/post_fv3.F90 b/io/post_fv3.F90
index d754dc45d..17b1572f1 100644
--- a/io/post_fv3.F90
+++ b/io/post_fv3.F90
@@ -4758,12 +4758,18 @@ subroutine set_postvars_fv3(wrt_int_state,grid_id,mype,mpicomp)
           enddo
         enddo ! do loop for l
 
-        l=lm
-        do j=jsta,jend
+        do l=1,lm
+         do j=jsta,jend
           do i=ista,iend
-
             tv = t(i,j,l) * (h1+d608*MAX(q(I,J,L),qmin))
             rhomid(i,j,l) = pmid(i,j,l) / (rd*tv)
+          enddo
+         enddo
+        enddo
+
+        l=lm
+        do j=jsta,jend
+          do i=ista,iend
 
             dustcb(i,j) = MAX(dustcb(i,j), 0.0)
             dustallcb(i,j) = MAX(dustallcb(i,j), 0.0)