GEMPAK Manual |
Programs
NAGRIB converts gridded data in GRIB files to GEMPAK gridded data.
INPUT PARAMETERS
GBFILE GRIB data file name
INDXFL GRIB index file name
GDOUTF Output grid file
PROJ Map projection/angles/margins|drop flag
GRDAREA Area covered by grid
KXKY Number of grid points in x;y
MAXGRD Maximum number of grids
CPYFIL Grid file whose navigation is to be used in new grid file | subare
GAREA Graphics area
OUTPUT Output device/filename
GBTBLS Input GRIB decoding tables
GBDIAG GRIB diagnostic elements
PDSEXT Y or N, add PDS extension if found
OVERWR Overwr flag
PROGRAM DESCRIPTION
NAGRIB will convert gridded data which is in GRIB files to
gridded data in a GEMPAK file.
The GRIB input file is specified in GBFILE, and the index
file associated with the GRIB file is in INDXFL.
The GEMPAK output file is given in GDOUTF. This file will
be opened if it exists and the projection information
matches the GRIB message projection. If the file does not
exist, it will be created using information given by the
user. If the output filename is "LIST", the GRIB file will
be scanned and the decoded GEMPAK header information will be
printed out.
The navigation information given by the user may take several
forms. The first possibility is for the user to enter
CPYFIL=#NNN, where NNN is a grid projection number listed in
GRDNAV.TBL. Second, the user may enter CPYFIL=GDS. This will
read the navigation from the Grid Definition Section of the
first GRIB message. Third, the user could enter CPYFIL=OLDFILE,
where OLDFILE is the name of an existing GEMPAK file from
which to copy the navigation information. Last, the user
may choose to use the PROJ, GRDAREA, and KXKY parameters to
specify the navigation.
GAREA is used to subset the grid. It may be specified as
any valid input for GAREA which are equivalent to lat/lon bounds
or as a single @ followed by the lower-left and upper-right grid
point coordinates separated by semicolons. In the former case,
the locations are rounded to the nearest grid points. The
projection type and angles entered as described above along with
the lat/lon coordinates of the corners of the subset define a
new grid navigation different from that of the full grid. The
navigation entered as described in the preceding paragraph is
always that of the full grid. To get the full grid, set GAREA to
blank, GRID or DSET.
OUTPUT defines the direction and destination of printed output.
GBTBLS allows the user to override the default GRIB decoding tables
by listing the filenames of the WMO and NCEP parameter tables,
the vertical coordinate table and the originating center table.
GBDIAG allows for detailed GRIB message section information,
byte-by-byte, to be printed out for selected GRIB messages.
Simply list those GRIB sections to be examined (IDS, PDS, GDS,
BDS, BMS, END or ALL for all sections), e.g., pds;gds.
Selected GRIB messages for examination may also be identified by
number in list and/or range format, e.g., 2;4;5-9.
PDSEXT is a logical flag which only becomes applicable when a PDS
extension exists in the GRIB message. If PDSEXT is YES, then a
sequence of characters specifying the extension information will be
appended onto the standard GEMPAK parameter name. If PDSEXT is NO,
then the append will not be performed. For instance, if the parameter
is HGHT and the PDS extension is an ensemble extension indicating that
this particular height field is an ensemble mean average, then the PDS
extension suffix will be "ENMA". The final GEMPAK parameter name
will be "HGHTENMA" and must be referenced as such in any GEMPAK program.
OVERWR is a logical flag that allows the user to overwrite the existing
GEMPAK grid file. If OVERWR is YES, then the grid file is being
overwritten. Otherwise, it remains the same.
EXAMPLES
1. Convert the data in the NGM GRIB file for 2 Sep 1993.
Create the GEMPAK file using projection #105 from the
grid navigation table. Define a subset using lat/lon.
GBFILE = ngm_930902_00_000.grib
INDXFL =
GDOUTF = ngm_930902_00.grd
PROJ =
GRDAREA =
KXKY =
MAXGRD = 3000
CPYFIL = #105
GAREA = ks-
OUTPUT = t
2. Convert the data in the MRF GRIB and INDEX files for
2 Sep 1993. Create the GEMPAK file using the navigation
information given the Grid Definition Section of the
GRIB message. Define a subset using grid coordinates.
GBFILE = mrf_930902_00_000.grib
INDXFL = mrf_930902_00_000.indx
GDOUTF = mrf_930902_00.grd
PROJ =
GRDAREA =
KXKY =
MAXGRD = 3000
CPYFIL = gds
GAREA = @230;110;290;140
OUTPUT = t
3. Convert the data in the GRIB file for North America. Create
the GEMPAK file using the user input for PROJ, GRDAREA and
KXKY. Do not do a subset, use the full grid.
GBFILE = north_amer.grib
INDXFL =
GDOUTF = north_amer.grd
PROJ = ced
GRDAREA = 0;-180;90;0
KXKY = 181;91
MAXGRD = 3000
CPYFIL =
GAREA = grid
OUTPUT = t
4. Convert the data in the GRIB file for the MESO model.
Override the default decoding tables.
GBFILE = /tmp/model/meso.mdl
INDXFL =
GDOUTF = /tmp/model/meso.mdl.gem
PROJ =
GRDAREA =
KXKY =
MAXGRD = 3000
CPYFIL = gds
GAREA = grid
OUTPUT = t
GBTBLS = /tmp/wmogrib2.tbl;/tmp/nmcgrib2.tbl;/tmp/vcrdgrib1.tbl
GBDIAG =
5. Scan a GRIB file and print out the PDS and GDS sections
for GRIB messages 1, 2, 6 through 12, and 88.
GBFILE = ensemble.grib
INDXFL =
GDOUTF = LIST
PROJ = ced
GRDAREA = 0;-180;90;0
KXKY = 181;91
MAXGRD = 3999
CPYFIL =
GAREA = grid
OUTPUT = t
GBTBLS =
GBDIAG = pds;gds|1;2;6-12;88
6. Convert the data in the GRIB file for the MESO model.
Overwrite the existing GEMPAK grid file.
GBFILE = /tmp/model/meso.mdl
INDXFL =
GDOUTF = /tmp/model/meso.mdl.gem
PROJ =
GRDAREA =
KXKY =
MAXGRD = 3000
CPYFIL = gds
GAREA = grid
OUTPUT = t
GBTBLS =
GBDIAG =
OVERWR = yes
ERROR MESSAGES
[NAGRIB +5] GRIB version 0 cannot be decoded.
[NAGRIB +4] Invalid parameter code table version.
[NAGRIB +3] No valid parameter found for ....
[NAGRIB +2] No valid vertical coordinate found for ....
[NAGRIB +1] WARNING: This grid is too large for GEMPAK programs.
[NAGRIB -1] Fatal error initializing TAE.
[NAGRIB -2] Fatal error reading TAE parameters.
[NAGRIB -3] Fatal error initializing GEMPLT.
[NAGRIB -4] Navigation information is invalid.
[NAGRIB -5] Grid area ... is invalid.
[NAGRIB -6] Grid size is invalid.
[NAGRIB -7] The grid file name may not be blank.
[NAGRIB -8] Navigation table cannot be read.
[NAGRIB -9] Grid name ... cannot be found in grid table.
[NAGRIB -10] Input for GEMPAK output file is blank.
[NAGRIB -11] Error opening GEMPAK grid file.
[NAGRIB -12] Error creating GEMPAK grid file.
[NAGRIB -13] Cannot open GRIB vertical coordinate table.
[NAGRIB -14] Cannot open GRIB parameter table.
[NAGRIB -15] Error opening GRIB file.
[NAGRIB -16] Error getting next message.
[NAGRIB -17] Error setting date/time.
[NAGRIB -18] Invalid bitmap specification.
[NAGRIB -19] Error reading GRIB file.
[NAGRIB -20] Invalid input for GAREA.
[NAGRIB -21] Grid is too large.
[NAGRIB -22] Length of BDS section less than or equal to 0