From 2f94bf1a1781de0bfc63278ecf5ad2a4666fa51f Mon Sep 17 00:00:00 2001
From: Wen Meng <Wen.Meng@noaa.gov>
Date: Fri, 19 Jul 2024 15:24:48 +0000
Subject: [PATCH] Add aerosol fix from Li Pan.

---
 io/post_fv3.F90 | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

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)