Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Feature/post update 20200122 #8

Merged
merged 3 commits into from
Jan 22, 2020
Merged

Feature/post update 20200122 #8

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
ce03b79
Develop mods 20200114 (#50)
hertneky Jan 19, 2020
78078f6
Lambert update (#45)
JiliDong-NOAA Jan 22, 2020
22b10da
solve merge conflict
JiliDong-NOAA Jan 22, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 6 additions & 6 deletions sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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 &
Expand Down
2 changes: 1 addition & 1 deletion sorc/ncep_post.fd/MDL2P.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
14 changes: 9 additions & 5 deletions sorc/ncep_post.fd/grib2_module.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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, &
Expand Down Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand Down