Skip to content

Commit

Permalink
Develop mods 20200114 (NOAA-EMC#50)
Browse files Browse the repository at this point in the history
* Fix ID index in WRF-CHEM SMOKE block, which currently causes failed build in debug mode

* Propose fix for filling a handful of 2d variables currently causing a divide-by-zero error when run in debug mode.
Currently filled with 'k' which fills with zeros
Needs to fill with 'll' which is what's being filled within the larger loop and different index value than 'k'

Co-authored-by: WenMeng-NOAA <48260754+WenMeng-NOAA@users.noreply.github.com>
  • Loading branch information
hertneky and WenMeng-NOAA committed Jan 19, 2020
1 parent 183cc64 commit ce03b79
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 7 deletions.
12 changes: 6 additions & 6 deletions sorc/ncep_post.fd/INITPOST_GFS_NEMS_MPIIO.f
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
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

0 comments on commit ce03b79

Please sign in to comment.