Skip to content

Commit

Permalink
Merge pull request NOAA-EMC#21 from EricJames-NOAA/missing_n_parm
Browse files Browse the repository at this point in the history
GRIB2 encoding fix for 6 variables in RRFS
  • Loading branch information
EricJames-NOAA authored Feb 3, 2022
2 parents b81fa4f + bd86c3f commit ac752ec
Showing 1 changed file with 36 additions and 12 deletions.
48 changes: 36 additions & 12 deletions sorc/ncep_post.fd/SURFCE.f
Original file line number Diff line number Diff line change
Expand Up @@ -791,8 +791,12 @@ SUBROUTINE SURFCE
if(grib=='grib2') then
cfld=cfld+1
fld_info(cfld)%ifld=IAVBLFLD(IGET(501))
fld_info(cfld)%ntrange=IFHR-ID(18)
fld_info(cfld)%tinvstat=1
if(ITSRFC>0) then
fld_info(cfld)%ntrange=1
else
fld_info(cfld)%ntrange=0
endif
fld_info(cfld)%tinvstat=IFHR-ID(18)
!$omp parallel do private(i,j,jj)
do j=1,jend-jsta+1
jj = jsta+j-1
Expand Down Expand Up @@ -820,8 +824,12 @@ SUBROUTINE SURFCE
if(grib=='grib2') then
cfld=cfld+1
fld_info(cfld)%ifld=IAVBLFLD(IGET(502))
fld_info(cfld)%ntrange=IFHR-ID(18)
fld_info(cfld)%tinvstat=1
if(ITSRFC>0) then
fld_info(cfld)%ntrange=1
else
fld_info(cfld)%ntrange=0
endif
fld_info(cfld)%tinvstat=IFHR-ID(18)
!$omp parallel do private(i,j,jj)
do j=1,jend-jsta+1
jj = jsta+j-1
Expand Down Expand Up @@ -2242,8 +2250,12 @@ SUBROUTINE SURFCE
if(grib=='grib2') then
cfld=cfld+1
fld_info(cfld)%ifld=IAVBLFLD(IGET(506))
fld_info(cfld)%ntrange=IFHR-ID(18)
fld_info(cfld)%tinvstat=1
if(ITSRFC>0) then
fld_info(cfld)%ntrange=1
else
fld_info(cfld)%ntrange=0
endif
fld_info(cfld)%tinvstat=IFHR-ID(18)
!$omp parallel do private(i,j,jj)
do j=1,jend-jsta+1
jj = jsta+j-1
Expand All @@ -2253,8 +2265,12 @@ SUBROUTINE SURFCE
enddo
cfld=cfld+1
fld_info(cfld)%ifld=IAVBLFLD(IGET(507))
fld_info(cfld)%ntrange=IFHR-ID(18)
fld_info(cfld)%tinvstat=1
if(ITSRFC>0) then
fld_info(cfld)%ntrange=1
else
fld_info(cfld)%ntrange=0
endif
fld_info(cfld)%tinvstat=IFHR-ID(18)
!$omp parallel do private(i,j,jj)
do j=1,jend-jsta+1
jj = jsta+j-1
Expand Down Expand Up @@ -6289,8 +6305,12 @@ SUBROUTINE SURFCE
if(grib=='grib2') then
cfld=cfld+1
fld_info(cfld)%ifld=IAVBLFLD(IGET(503))
fld_info(cfld)%ntrange=IFHR-ID(18)
fld_info(cfld)%tinvstat=1
if(ITSRFC>0) then
fld_info(cfld)%ntrange=1
else
fld_info(cfld)%ntrange=0
endif
fld_info(cfld)%tinvstat=IFHR-ID(18)
datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
endif
ENDIF
Expand All @@ -6314,8 +6334,12 @@ SUBROUTINE SURFCE
if(grib=='grib2') then
cfld=cfld+1
fld_info(cfld)%ifld=IAVBLFLD(IGET(504))
fld_info(cfld)%ntrange=IFHR-ID(18)
fld_info(cfld)%tinvstat=1
if(ITSRFC>0) then
fld_info(cfld)%ntrange=1
else
fld_info(cfld)%ntrange=0
endif
fld_info(cfld)%tinvstat=IFHR-ID(18)
datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
endif

Expand Down

0 comments on commit ac752ec

Please sign in to comment.