[gembud] GDBIINT bug or corrupt NOAA grib2 file?

• To: gembud@xxxxxxxxxxxxxxxx
• Subject: [gembud] GDBIINT bug or corrupt NOAA grib2 file?
• From: elowery@xxxxxxxxxxxx
• Date: Thu, 20 Jun 2019 19:51:24 -0400 (EDT)

#!/bin/csh

source $NAWIPS/Gemenviron

echo Script Name: "$0"
echo Total Number of Arguments Passed: "$#"

if ( $# != 2 ) then
    echo not two arguments, leaving
    exit 1
endif

## initializing variables
set res="$1"
set hhh="$2"
set grib2_in="gfs.t00z.pgrb2.$res.f$hhh"
set grid_in="gfs.t00z.pgrob2.$res.f$hhh.grd"
set grid_2p50_out="gfs.t00z.pgrb2.2p50.f$hhh.grd"
set grid_2p50_template="gfs_25deg_12hr.grd"


## make sure files exist
if ( ! -f $grib2_in ) then
    echo "file($grib2_in) does not exist, exiting"
    exit 1
endif
if ( ! -f $grid_2p50_template ) then
    echo "file($grid_2p50_template) does not exist, exiting"
    exit 1
endif


## cleanup files to be created if they exist
if ( -f $grid_in ) then
    echo "removing residual grid file($grid_in)"
    rm "$grid_in"
endif
if ( -f $grid_2p50_out ) then
    echo "removing residual grid_file($grid_2p50_out)"
    rm "$grid_2p50_out"
endif

echo "using nagrib2 to convert $grib2_in to $grid_in"

  nagrib2 << NAGRIB2_EOF
  GBFILE=$grib2_in
  GDOUTF=$grid_in
  PROJ=
  GRDAREA=
  KXKY=
  MAXGRD=1500
  CPYFIL=gds
  GAREA=grid
  OUTPUT=n
  G2TBLS=
  G2DIAG=
r

e
NAGRIB2_EOF

gpend


rm *nts
echo "using gdcfil to create $grid_2p50_out"

  gdcfil << GDCFIL_EOF
  PROJ=
  GRDAREA=
  KXKY=
  MAXGRD=1500
  ANLYSS=
  GDOUTF=$grid_2p50_out
  CPYFIL=$grid_2p50_template
r

e
GDCFIL_EOF

gpend


rm *.nts
echo "using gdbiint to convert $grid_in to $grid_2p50_out"

  gdbiint << GDBIINT_EOF
  GDFILE=$grid_in
  GDOUTF=$grid_2p50_out
  GFUNC=
  GLEVEL=
  GVCORD=
  GDATTIM=
  GDNUM=ALL
r

e
GDBIINT_EOF

gpend


#rm *.nts

• Follow-Ups:
  • Re: [gembud] GDBIINT bug or corrupt NOAA grib2 file?
    • From: Steven Decker
  • Re: [gembud] GDBIINT bug or corrupt NOAA grib2 file?
    • From: Steven Decker