diff --git a/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f b/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f index 104eee42b..7f7153cca 100644 --- a/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f +++ b/sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f @@ -987,12 +987,12 @@ SUBROUTINE INITPOST_GFS_NEMS_MPIIO(iostatusAER) do j=jsta,jend do k=1,lm do i=1,im - p2d(i,k) = pmid(i,j,k)*0.01 - t2d(i,k) = t(i,j,k) - q2d(i,k) = q(i,j,k) - cw2d(i,k) = cwm(i,j,k) - es = min(fpvsnew(t(i,j,k)),pmid(i,j,k)) - qs2d(i,k) = eps*es/(pmid(i,j,k)+epsm1*es)!saturation q for GFS + p2d(i,k) = pmid(i,j,ll)*0.01 + t2d(i,k) = t(i,j,ll) + q2d(i,k) = q(i,j,ll) + cw2d(i,k) = cwm(i,j,ll) + es = min(fpvsnew(t(i,j,ll)),pmid(i,j,ll)) + qs2d(i,k) = eps*es/(pmid(i,j,ll)+epsm1*es)!saturation q for GFS enddo enddo call progcld1 & diff --git a/sorc/ncep_post.fd/MDL2P.f b/sorc/ncep_post.fd/MDL2P.f index 76e929114..3f4ce004a 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -2092,7 +2092,7 @@ SUBROUTINE MDL2P(iostatusD3D) if(grib == 'grib1')then ID(1:25)=0 ID(02)=141 ! Parameter Table 141 - !ID(36)=2 + CALL GRIBIT(IGET(738),LP,GRID1,IM,JM) elseif(grib == 'grib2') then cfld = cfld + 1 diff --git a/sorc/ncep_post.fd/grib2_module.f b/sorc/ncep_post.fd/grib2_module.f index 89c15be35..51d6d7916 100644 --- a/sorc/ncep_post.fd/grib2_module.f +++ b/sorc/ncep_post.fd/grib2_module.f @@ -1212,7 +1212,7 @@ subroutine getgds(ldfgrd,len3,ifield3len,igds,ifield3) ! !***** set up gds kpds to call Boi's code ! - use CTLBLK_mod, only : im,jm,gdsdegr + use CTLBLK_mod, only : im,jm,gdsdegr,modelname use gridspec_mod, only: DXVAL,DYVAL,CENLAT,CENLON,LATSTART,LONSTART,LATLAST, & & LONLAST,MAPTYPE,STANDLON,latstartv,cenlatv,lonstartv, & cenlonv,TRUELAT1,TRUELAT2,LATSTART_R,LONSTART_R, & @@ -1244,9 +1244,11 @@ subroutine getgds(ldfgrd,len3,ifield3len,igds,ifield3) ifield3(9) = jm !number of points along the y-axis ifield3(10) = latstart !latitude of first grid point ifield3(11) = lonstart !longitude of first grid point -! ifield3(12) = 8 !Resolution and component flags + ifield3(12) = 8 !Resolution and component flags ! Jili Dong change grid to earth relative - ifield3(12) = 0 !Resolution and component flags + if (modelname == 'FV3R') then + ifield3(12) = 0 !Resolution and component flags + endif ifield3(13) = TRUELAT1 ifield3(14) = STANDLON !longitude of meridian parallel to y-axis along which latitude increases @@ -1367,9 +1369,11 @@ subroutine getgds(ldfgrd,len3,ifield3len,igds,ifield3) endif ifield3(12) = latstart_r !latitude of first grid point ifield3(13) = lonstart_r !longitude of first grid point -! ifield3(14) = 56 !Resolution and component flags + ifield3(14) = 56 !Resolution and component flags ! Jili Dong change grid to earth relative (Matt Pyle) - ifield3(14) = 48 !Resolution and component flags + if(modelname=='FV3R') then + ifield3(14) = 48 !Resolution and component flags + endif ifield3(15) = latlast_r !latitude of last grid point ifield3(16) = lonlast_r !longitude of last grid point