Update post_tag_gfs128 for chemical output

Update diagnostic output option for 2D chemical output
Update RHOMID calculation to prevent division by zero

parm/post_tag_gfs128

@@ -2,6 +2,6 @@
 MODELNAME='GFS'
 /
 &NAMPGB
-  KPO=57,PO=1000.,975.,950.,925.,900.,875.,850.,825.,800.,775.,750.,725.,700.,675.,650.,625.,600.,575.,550.,525.,500.,475.,450.,425.,400.,375.,350.,325.,300.,275.,250.,225.,200.,175.,150.,125.,100.,70.,50.,40.,30.,20.,15.,10.,7.,5.,3.,2.,1.,0.7,0.4,0.2,0.1,0.07,0.04,0.02,0.01,
+ KPO=50,PO=1000.,975.,950.,925.,900.,875.,850.,825.,800.,775.,750.,725.,700.,675.,650.,625.,600.,575.,550.,525.,500.,475.,450.,425.,400.,375.,350.,325.,300.,275.,250.,225.,200.,175.,150.,125.,100.,70.,50.,40.,30.,20.,15.,10.,7.,5.,3.,2.,1.,0.4,gccpp_on=.true., rdaod=.true.,
   /
 EOF

sorc/ncep_post.fd/CLDRAD.f

@@ -109,7 +109,7 @@ SUBROUTINE CLDRAD
                             TCLOD, ARDSW, TRDSW, ARDLW, NBIN_DU, TRDLW, IM,   &
                             NBIN_SS, NBIN_OC,NBIN_BC,NBIN_SU,NBIN_NO3,DTQ2,   &
                             JM, LM, gocart_on, gccpp_on, nasa_on, me, rdaod,  &
-                            ISTA, IEND,aqf_on
+                            ISTA, IEND,aqf_on,d2d_chem
       use rqstfld_mod, only: IGET, ID, LVLS, IAVBLFLD
       use gridspec_mod, only: dyval, gridtype
       use cmassi_mod,  only: TRAD_ice
@@ -5586,7 +5586,7 @@ SUBROUTINE CLDRAD
 
       endif  !nasa_on
 
-      if (gocart_on .or. gccpp_on ) then
+      if ((gocart_on .or. gccpp_on ) .and. d2d_chem) then
 !! ADD EMISSION FLUXES,dry depostion, wet/convective depostion (kg/m2/sec)
 !! The AER file uses 1.E6 to scale all 2d diagnosis fields
 !! Multiply by 1.E-6 to revert these fields back

sorc/ncep_post.fd/INITPOST_NETCDF.f

@@ -1413,7 +1413,7 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
           do j=jsta,jend
             do i=ista,iend
 
-              TV = T(I,J,L) * (H1+D608*MAX(Q(I,J,L),QMIN))
+              TV = MAX(T(I,J,L),QMIN) * (H1+D608*MAX(Q(I,J,L),QMIN))
               RHOMID(I,J,L) = PMID(I,J,L) / (RD*TV)
               do n = 1,  NBIN_DU
                 IF ( dust(i,j,l,n) < SPVAL) THEN