diff --git a/parm/post_avblflds.xml b/parm/post_avblflds.xml
index ff6a16031..815183a6f 100755
--- a/parm/post_avblflds.xml
+++ b/parm/post_avblflds.xml
@@ -6219,16 +6219,17 @@
 
     <param>
        <post_avblfldidx>685</post_avblfldidx>
-       <shortname>DU_CR_AER_SFC_MASS_CON</shortname>
+       <shortname>DUST10_SFC_MASS_CON</shortname>
        <stats_proc>AVE</stats_proc>
        <pdstmpl>tmpl4_48</pdstmpl>
-       <pname>MASSDEN</pname>
+       <pname>PMTC</pname>
+       <fixed_sfc1_type>surface</fixed_sfc1_type>
        <aerosol_type>dust_dry</aerosol_type>
        <typ_intvl_size>smaller_than_first_limit</typ_intvl_size>
        <scale_fact_1st_size>6</scale_fact_1st_size>
        <scale_val_1st_size>10</scale_val_1st_size>
        <fixed_sfc1_type>surface</fixed_sfc1_type>
-       <scale> 9.0</scale>
+       <scale>9.0</scale>
     </param>
 
     <param>
diff --git a/parm/postcntrl_gefs_chem.xml b/parm/postcntrl_gefs_chem.xml
index c4ce522cf..5abf7128b 100755
--- a/parm/postcntrl_gefs_chem.xml
+++ b/parm/postcntrl_gefs_chem.xml
@@ -225,6 +225,12 @@
        <scale>9.0</scale>
     </param>
 
+    <param>
+       <shortname>DUST10_SFC_MASS_CON</shortname>
+       <table_info>NCEP</table_info>
+       <scale>9.0</scale>
+    </param>
+    
     <param>
        <shortname>DUST25_SFC_MASS_CON</shortname>
        <table_info>NCEP</table_info>
diff --git a/sorc/ncep_post.fd/ALLOCATE_ALL.f b/sorc/ncep_post.fd/ALLOCATE_ALL.f
index afe9c7eed..75d6bfc64 100644
--- a/sorc/ncep_post.fd/ALLOCATE_ALL.f
+++ b/sorc/ncep_post.fd/ALLOCATE_ALL.f
@@ -1260,7 +1260,9 @@ SUBROUTINE ALLOCATE_ALL()
         allocate(dustallcb(im,jsta_2l:jend_2u))
         allocate(ssallcb(im,jsta_2l:jend_2u))
         allocate(dustpm(im,jsta_2l:jend_2u))
+        allocate(dustpm10(im,jsta_2l:jend_2u))	
         allocate(sspm(im,jsta_2l:jend_2u))
+	allocate(maod(im,jsta_2l:jend_2u))
 !Initialization
 !$omp parallel do private(i,j)
        do j=jsta_2l,jend_2u
@@ -1295,7 +1297,9 @@ SUBROUTINE ALLOCATE_ALL()
            dustallcb(i,j)=spval
            ssallcb(i,j)=spval
            dustpm(i,j)=spval
+	   dustpm10(i,j)=spval
            sspm(i,j)=spval
+	   maod(i,j)=spval
          enddo
        enddo
       endif
diff --git a/sorc/ncep_post.fd/CALPW.f b/sorc/ncep_post.fd/CALPW.f
index 51aa4c457..b63e9ce1f 100644
--- a/sorc/ncep_post.fd/CALPW.f
+++ b/sorc/ncep_post.fd/CALPW.f
@@ -273,12 +273,10 @@ SUBROUTINE CALPW(PW,IDECID)
           DO I=1,IM
             if(PINT(I,J,L+1) <spval .and. Qdum(I,J) < spval) then
              DP      = PINT(I,J,L+1) - PINT(I,J,L)
-             PW(I,J) = PW(I,J) + Qdum(I,J)*DP*GI*HTM(I,J,L)
-            IF (IDECID == 17 .or. IDECID == 20 .or. IDECID == 21) THEN
-             PW(I,J) = PW(I,J) + Qdum(I,J)*MAX(DP,0.)*GI*HTM(I,J,L)
-            ENDIF
             IF (IDECID == 19) THEN
              PW(I,J) = PW(I,J) + Qdum(I,J)
+	    ELSE
+	     PW(I,J) = PW(I,J) + Qdum(I,J)*MAX(DP,0.)*GI*HTM(I,J,L)
             ENDIF
             IF (IDECID == 14) PWS(I,J) = PWS(I,J) + QS(I,J)*DP*GI*HTM(I,J,L)
             else
diff --git a/sorc/ncep_post.fd/CLDRAD.f b/sorc/ncep_post.fd/CLDRAD.f
index 1eabc8610..26c375205 100644
--- a/sorc/ncep_post.fd/CLDRAD.f
+++ b/sorc/ncep_post.fd/CLDRAD.f
@@ -90,7 +90,7 @@ SUBROUTINE CLDRAD
                          AVGCPRATE,                                           &
                          DUSTCB,SSCB,BCCB,OCCB,SULFCB,DUSTPM,SSPM,aod550,     &
                          du_aod550,ss_aod550,su_aod550,oc_aod550,bc_aod550,   &
-                         PWAT
+                         PWAT,DUSTPM10,MAOD
       use masks,    only: LMH, HTM
       use params_mod, only: TFRZ, D00, H99999, QCLDMIN, SMALL, D608, H1, ROG, &
                             GI, RD, QCONV, ABSCOEFI, ABSCOEF, STBOL, PQ0, A2, &
@@ -5160,6 +5160,26 @@ SUBROUTINE CLDRAD
           datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
          endif
       ENDIF
+      
+!! Multiply by 1.E-6 to revert these fields back
+      IF (IGET(667)>0) THEN
+         GRID1=SPVAL
+!$omp parallel do private(i,j)
+         DO J = JSTA,JEND
+            DO I = 1,IM
+               IF(BCEM(I,J,1)<SPVAL) GRID1(I,J) = BCEM(I,J,1)
+               DO K=2,NBIN_BC
+               IF(BCEM(I,J,K)<SPVAL)&
+                GRID1(I,J) = GRID1(I,J) + BCEM(I,J,K)
+               END DO
+            END DO
+         END DO
+         if(grib=='grib2') then
+          cfld=cfld+1
+          fld_info(cfld)%ifld=IAVBLFLD(IGET(667))
+          datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
+         endif
+      ENDIF
 
       IF (IGET(660)>0) THEN
          GRID1=SPVAL
@@ -5179,6 +5199,21 @@ SUBROUTINE CLDRAD
           datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
          endif
       ENDIF
+      
+      IF (IGET(699)>0) THEN
+         GRID1=SPVAL
+!$omp parallel do private(i,j)
+         DO J = JSTA,JEND
+            DO I = 1,IM
+                GRID1(I,J) = MAOD(I,J)
+            END DO
+         END DO
+         if(grib=='grib2') then
+          cfld=cfld+1
+          fld_info(cfld)%ifld=IAVBLFLD(IGET(699))
+          datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
+         endif
+      ENDIF
 !! ADD DUST DRY DEPOSITION FLUXES (kg/m2/sec)
 !
 !      IF (IGET(661)>0) THEN
@@ -5217,6 +5252,20 @@ SUBROUTINE CLDRAD
          endif
       ENDIF
 
+      IF (IGET(685)>0 ) THEN
+!$omp parallel do private(i,j)
+         DO J = JSTA,JEND
+            DO I = 1,IM
+               GRID1(I,J) = DUSTPM10(I,J)   !ug/m3 
+            END DO
+         END DO
+         if(grib=='grib2') then
+           cfld=cfld+1
+           fld_info(cfld)%ifld=IAVBLFLD(IGET(685))
+           datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
+         endif
+      ENDIF
+      
 !! ADD DUST WET DEPOSITION FLUXES (kg/m2/sec)
 !      IF (IGET(662)>0) THEN
 !         DO J = JSTA,JEND
@@ -5430,7 +5479,10 @@ SUBROUTINE CLDRAD
       IF (IGET(674)>0) call wrt_aero_diag(674,nbin_oc,ocwt)
       IF (IGET(682)>0) call wrt_aero_diag(682,nbin_oc,ocsv)
 !      print *,'aft wrt disg ocwt'
-
+!! wrt MIE AOD at 550nm
+      IF (IGET(699).GT.0) call wrt_aero_diag(699,1,maod)
+      print *,'aft wrt disg maod'
+      
 !! wrt SU diag field
 !      IF (IGET(675)>0) call wrt_aero_diag(675,nbin_su,suem)
 !      IF (IGET(676)>0) call wrt_aero_diag(676,nbin_su,susd)
diff --git a/sorc/ncep_post.fd/DEALLOCATE.f b/sorc/ncep_post.fd/DEALLOCATE.f
index 1bcfdd788..10fac93fa 100644
--- a/sorc/ncep_post.fd/DEALLOCATE.f
+++ b/sorc/ncep_post.fd/DEALLOCATE.f
@@ -532,7 +532,9 @@ SUBROUTINE DE_ALLOCATE
         deallocate(dustallcb)
         deallocate(ssallcb)
         deallocate(dustpm)
+	deallocate(dustpm10)
         deallocate(sspm)
+	deallocate(maod)
       endif
 !
 ! HWRF RRTMG output 
diff --git a/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f b/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f
index 0d9b4da88..e8fcdcc28 100644
--- a/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f
+++ b/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f
@@ -55,8 +55,8 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
               avgedir,avgecan,avgetrans,avgesnow,avgprec_cont,avgcprate_cont, &
               avisbeamswin,avisdiffswin,airbeamswin,airdiffswin, &
               alwoutc,alwtoac,aswoutc,aswtoac,alwinc,aswinc,avgpotevp,snoavg, &
-              dustcb,bccb,occb,sulfcb,sscb,dustallcb,ssallcb,dustpm,sspm,pp25cb,pp10cb, &
-              ti 
+              dustcb,bccb,occb,sulfcb,sscb,dustallcb,ssallcb,dustpm,dustpm10,sspm,pp25cb, &
+              pp10cb,maod,ti 
       use soil,  only: sldpth, sh2o, smc, stc
       use masks, only: lmv, lmh, htm, vtm, gdlat, gdlon, dx, dy, hbm2, sm, sice
 !     use kinds, only: i_llong
@@ -1413,7 +1413,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
 
            dustallcb(1:im,jsta_2l:jend_2u)=dustallcb(1:im,jsta_2l:jend_2u)+ &
            (dust(1:im,jsta_2l:jend_2u,ll,1)+dust(1:im,jsta_2l:jend_2u,ll,2)+ &
-           dust(1:im,jsta_2l:jend_2u,ll,3)+0.67*dust(1:im,jsta_2l:jend_2u,ll,4))* &
+           dust(1:im,jsta_2l:jend_2u,ll,3)+0.74*dust(1:im,jsta_2l:jend_2u,ll,4))* &
            dpres(1:im,jsta_2l:jend_2u,ll)/grav
 
 !         if(debugprint)print*,'sample l ',VarName,' = ',ll,dust(isa,jsa,ll,5)
@@ -1470,7 +1470,8 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
           ,salt(1:im,jsta_2l:jend_2u,ll,3))
 
             sscb(1:im,jsta_2l:jend_2u)=sscb(1:im,jsta_2l:jend_2u)+ &
-         (salt(1:im,jsta_2l:jend_2u,ll,2)+0.75*salt(1:im,jsta_2l:jend_2u,ll,3))* &
+	    (salt(1:im,jsta_2l:jend_2u,ll,1)+ &
+	    salt(1:im,jsta_2l:jend_2u,ll,2)+0.83*salt(1:im,jsta_2l:jend_2u,ll,3))* &
            dpres(1:im,jsta_2l:jend_2u,ll)/grav
      
 !         if(debugprint)print*,'sample l ',VarName,' = ',ll,salt(isa,jsa,ll,3)
@@ -1501,7 +1502,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
           ssallcb(1:im,jsta_2l:jend_2u)=ssallcb(1:im,jsta_2l:jend_2u)+ &
          (salt(1:im,jsta_2l:jend_2u,ll,1)+salt(1:im,jsta_2l:jend_2u,ll,2)+ &
           salt(1:im,jsta_2l:jend_2u,ll,3)+ &
-          salt(1:im,jsta_2l:jend_2u,ll,4)*0.83)* &
+          salt(1:im,jsta_2l:jend_2u,ll,4))* &
            dpres(1:im,jsta_2l:jend_2u,ll)/grav
 
 !         if(debugprint)print*,'sample l ',VarName,' = ',ll,salt(isa,jsa,ll,5)
@@ -1728,17 +1729,18 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
             bccb(i,j) = MAX(bccb(i,j), 0.0)
             occb(i,j) = MAX(occb(i,j), 0.0)
             sulfcb(i,j) = MAX(sulfcb(i,j), 0.0)
-            pp25cb(i,j) = MAX(sulfcb(i,j), 0.0)
-            pp10cb(i,j) = MAX(sulfcb(i,j), 0.0)
+            pp25cb(i,j) = MAX(pp25cb(i,j), 0.0)
+            pp10cb(i,j) = MAX(pp10cb(i,j), 0.0)
 !      PM10 concentration
        dusmass(i,j)=(dust(i,j,l,1)+dust(i,j,l,2)+dust(i,j,l,3)+ &
        0.74*dust(i,j,l,4)+salt(i,j,l,1)+salt(i,j,l,2)+salt(i,j,l,3)+ &
-       salt(i,j,l,4) + &
-       salt(i,j,l,5)+soot(i,j,l,1)+soot(i,j,l,2)+waso(i,j,l,1)+ &
+       salt(i,j,l,4) + soot(i,j,l,1)+soot(i,j,l,2)+waso(i,j,l,1)+ &
        waso(i,j,l,2) +suso(i,j,l,1)+pp25(i,j,l,1)+pp10(i,j,l,1)) &
        *RHOMID(i,j,l)  !ug/m3
 !      PM25 dust and seasalt      
        dustpm(i,j)=(dust(i,j,l,1)+0.38*dust(i,j,l,2))*RHOMID(i,j,l) !ug/m3
+       dustpm10(i,j)=(dust(i,j,l,1)+dust(i,j,l,2)+dust(i,j,l,3)+ &
+       0.74*dust(i,j,l,4))*RHOMID(i,j,l) !ug/m3
        sspm(i,j)=(salt(i,j,l,1)+salt(i,j,l,2)+ &
        0.83*salt(i,j,l,3))*RHOMID(i,j,l)  !ug/m3 
 !      PM25 concentration       
@@ -3673,43 +3675,6 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
         enddo
       enddo
 
-! done with flux file, close it for now
-      call nemsio_close(ffile,iret=status)
-      deallocate(tmp,recname,reclevtyp,reclev)
-
-
-! Retrieve aer fields if it's listed (GOCART)
-      print *, 'iostatus for aer file=', iostatusAER
-      if(iostatusAER == 0) then ! start reading aer file
-       call nemsio_open(rfile,trim(fileNameAER),'read',mpi_comm_comp &
-                       ,iret=status)
-       if ( Status /= 0 ) then
-        print*,'error opening ',fileNameAER, ' Status = ', Status
-       endif
-       call nemsio_getfilehead(rfile,iret=status,nrec=nrec)
-       print*,'nrec for aer file=',nrec
-       allocate(recname(nrec),reclevtyp(nrec),reclev(nrec))
-       call nemsio_getfilehead(rfile,iret=iret,recname=recname       &
-                              ,reclevtyp=reclevtyp,reclev=reclev)
-       if(debugprint)then
-         if (me == 0)then
-           do i=1,nrec
-             print *,'recname,reclevtyp,reclev=',trim(recname(i)),' ', &
-                      trim(reclevtyp(i)),reclev(i)
-           end do
-         end if
-       end if
-! start reading nemsio aer files using parallel read
-      fldsize=(jend-jsta+1)*im
-      allocate(tmp(fldsize*nrec))
-      print*,'allocate tmp successfully'
-      tmp=0.
-      call nemsio_denseread(rfile,1,im,jsta,jend,tmp,iret=iret)
-      if(iret/=0)then
-        print*,"fail to read aer file using mpi io read, stopping"
-        stop 
-      end if
-
 ! retrieve dust emission fluxes
       do K = 1, nbin_du
        if ( K == 1) VarName='duem001'
@@ -3717,7 +3682,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='duem003'
        if ( K == 4) VarName='duem004'
        if ( K == 5) VarName='duem005'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3728,12 +3693,12 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
 
 ! retrieve dust sedimentation fluxes
       do K = 1, nbin_du
-       if ( K == 1) VarName='dust1SD'
-       if ( K == 2) VarName='dust2SD'
-       if ( K == 3) VarName='dust3SD'
-       if ( K == 4) VarName='dust4SD'
-       if ( K == 5) VarName='dsut5SD'
-       VcoordName='atmos col'
+       if ( K == 1) VarName='dust1sd'
+       if ( K == 2) VarName='dust2sd'
+       if ( K == 3) VarName='dust3sd'
+       if ( K == 4) VarName='dust4sd'
+       if ( K == 5) VarName='dsut5sd'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3749,7 +3714,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='dust3dp'
        if ( K == 4) VarName='dust4dp'
        if ( K == 5) VarName='dust5dp'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3767,7 +3732,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='dust3wtl'
        if ( K == 4) VarName='dust4wtl'
        if ( K == 5) VarName='dust5wtl'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3781,7 +3746,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='dust3wtc'
        if ( K == 4) VarName='dust4wtc'
        if ( K == 5) VarName='dust5wtc'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3796,13 +3761,29 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='ssem003'
        if ( K == 4) VarName='ssem004'
        if ( K == 5) VarName='ssem005'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
                            ,recname,reclevtyp,reclev,VarName,VcoordName&
                            ,ssem(1,jsta_2l,K))
       enddo
+      
+! retrieve seasalt emission fluxes
+      do K = 1, nbin_ss
+       if ( K == 1) VarName='seas1sd'
+       if ( K == 2) VarName='seas2sd'
+       if ( K == 3) VarName='seas3sd'
+       if ( K == 4) VarName='seas4sd'
+       if ( K == 5) VarName='seas5sd'
+       VcoordName='sfc'
+       l=1
+       call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
+                           ,l,nrec,fldsize,spval,tmp                   &
+                           ,recname,reclevtyp,reclev,VarName,VcoordName&
+                           ,sssd(1,jsta_2l,K))
+      enddo
+      
 
 ! retrieve seasalt dry deposition fluxes
       do K = 1, nbin_ss
@@ -3811,7 +3792,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='seas3dp'
        if ( K == 4) VarName='seas4dp'
        if ( K == 5) VarName='seas5dp'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3826,7 +3807,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='seas3wtl'
        if ( K == 4) VarName='seas4wtl'
        if ( K == 5) VarName='seas5wtl'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3841,7 +3822,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
        if ( K == 3) VarName='seas1wtc'
        if ( K == 4) VarName='seas1wtc'
        if ( K == 5) VarName='seas1wtc'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3853,7 +3834,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_bc
        if ( K == 1) VarName='bceman'
        if ( K == 2) VarName='bcembb'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3865,7 +3846,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_bc
        if ( K == 1) VarName='bc1sd'
        if ( K == 2) VarName='bc2sd'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3877,7 +3858,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_bc
        if ( K == 1) VarName='bc1dp'
        if ( K == 2) VarName='bc2dp'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3889,7 +3870,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_bc
        if ( K == 1) VarName='bc1wtl'
        if ( K == 2) VarName='bc2wtl'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3901,7 +3882,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_bc
        if ( K == 1) VarName='bc1wtc'
        if ( K == 2) VarName='bc2wtc'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3913,7 +3894,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_oc
        if ( K == 1) VarName='oceman'
        if ( K == 2) VarName='ocembb'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3925,7 +3906,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_oc
        if ( K == 1) VarName='oc1sd'
        if ( K == 2) VarName='oc2sd'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3935,9 +3916,9 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
 
 ! retrieve oc dry deposition fluxes
       do K = 1, nbin_oc
-       if ( K == 1) VarName='c1dp'
-       if ( K == 2) VarName='c2dp'
-       VcoordName='atmos col'
+       if ( K == 1) VarName='oc1dp'
+       if ( K == 2) VarName='oc2dp'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3949,7 +3930,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_oc
        if ( K == 1) VarName='oc1wtl'
        if ( K == 2) VarName='oc2wtl'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3961,7 +3942,7 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
       do K = 1, nbin_oc
        if ( K == 1) VarName='oc1wtc'
        if ( K == 2) VarName='oc2wtc'
-       VcoordName='atmos col'
+       VcoordName='atmos sfc'
        l=1
        call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
                            ,l,nrec,fldsize,spval,tmp                   &
@@ -3969,8 +3950,20 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
                            ,ocsv(1,jsta_2l,K))
       enddo
 
+! retrieve MIE AOD
+       VarName='maod'
+       VcoordName='sfc'
+       l=1
+       call assignnemsiovar(im,jsta,jend,jsta_2l,jend_2u               &
+                           ,l,nrec,fldsize,spval,tmp                   &
+                           ,recname,reclevtyp,reclev,VarName,VcoordName&
+                           ,maod(1,jsta_2l))
 
 
+! done with flux file, close it for now
+      call nemsio_close(ffile,iret=status)
+      deallocate(tmp,recname,reclevtyp,reclev)
+      
 !lzhang
 !! retrieve sfc mass concentration
 !      VarName='DUSMASS'
@@ -4015,11 +4008,6 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER)
 !                          ,ducmass25)
 !     if(debugprint)print*,'sample ',VarName,' = ',ducmass25(isa,jsa)
 
-        if (me == 0) print *,'after aer files reading,mype=',me
-       call nemsio_close(rfile,iret=status)
-       deallocate(tmp,recname,reclevtyp,reclev)
-      end if ! end of aer file read
-
 ! pos east
        call collect_loc(gdlat,dummy)
        if(me == 0)then
diff --git a/sorc/ncep_post.fd/MPI_FIRST.f b/sorc/ncep_post.fd/MPI_FIRST.f
index 920f78e18..f24dd90fe 100644
--- a/sorc/ncep_post.fd/MPI_FIRST.f
+++ b/sorc/ncep_post.fd/MPI_FIRST.f
@@ -79,7 +79,7 @@ SUBROUTINE MPI_FIRST()
               ocsmass25, occmass25, bcsmass, bccmass, bcsmass25, bccmass25,            &
               sssmass, sscmass, sssmass25, sscmass25, ducmass25,                       &
               dustcb, sscb, bccb, occb, sulfcb, dustallcb, ssallcb,dustpm,sspm, pp25cb,&
-              pp10cb, ti
+              dustpm10, pp10cb, maod, ti
       use soil, only:  smc, stc, sh2o, sldpth, rtdpth, sllevel
       use masks, only: htm, vtm, hbm2, sm, sice, lmh, gdlat, gdlon, dx, dy, lmv
       use ctlblk_mod, only: me, num_procs, jm, jsta, jend, jsta_m, jsta_m2,           &
diff --git a/sorc/ncep_post.fd/VRBLS2D_mod.f b/sorc/ncep_post.fd/VRBLS2D_mod.f
index 134d014f1..569f34ea5 100644
--- a/sorc/ncep_post.fd/VRBLS2D_mod.f
+++ b/sorc/ncep_post.fd/VRBLS2D_mod.f
@@ -95,7 +95,7 @@ module vrbls2d
       ,SSSMASS(:,:),SSCMASS(:,:),SSSMASS25(:,:),SSCMASS25(:,:)               &
       ,DUSTCB(:,:),SSCB(:,:),OCCB(:,:),BCCB(:,:),SULFCB(:,:)                 &
       ,DUSTALLCB(:,:),SSALLCB(:,:),DUSTPM(:,:),SSPM(:,:),PP25CB(:,:)         &
-      ,PP10CB(:,:)!lzhang, add for FV3-Chem
+      ,DUSTPM10(:,:),PP10CB(:,:),maod(:,:)!lzhang, add for FV3-Chem
  
 !
       end module vrbls2d