getgbm.f File Reference

Find and unpack a grib message. More...

Go to the source code of this file.

Functions/Subroutines
subroutine	getgbm (lugb, lugi, jf, j, jpds, jgds, mbuf, cbuf, nlen, nnum, mnum, kf, k, kpds, kgds, lb, f, iret)
	Read a grib index file (or optionally the grib file itself) to get the index buffer (i.e.

Detailed Description

Find and unpack a grib message.

Author

Mark Iredell

Date

1994-04-01

Definition in file getgbm.f.

Function/Subroutine Documentation

getgbm()

subroutine getgbm	(		lugb,
			lugi,
			jf,
			j,
		integer, dimension(200)	jpds,
		integer, dimension(200)	jgds,
			mbuf,
		character, dimension(mbuf)	cbuf,
			nlen,
			nnum,
			mnum,
			kf,
			k,
		integer, dimension(200)	kpds,
		integer, dimension(200)	kgds,
		logical*1, dimension(jf)	lb,
		real, dimension(jf)	f,
			iret
	)

Read a grib index file (or optionally the grib file itself) to get the index buffer (i.e.

table of contents) for the grib file. Find in the index buffer a reference to the grib message requested. The grib message request specifies the number of messages to skip and the unpacked pds and gds parameters. A requested parameter of -1 means to allow any value of this parameter to be found.) If the requested grib message is found, then it is read from the grib file and unpacked. Its message number is returned along with the unpacked pds and gds parameters, the unpacked bitmap (if any), and the unpacked data. If the grib message is not found, then the return code will be nonzero.

Program history log:

• Mark Iredell 1994-04-01
• Mark Iredell 1995-10-31 Modularized portions of code into subprograms and allowed for unspecified index file.
• Chuang 2004-07-22 Add nbitss to the argument list of getgb1r that is called in this subroutine.
• Wang 2010-03-02 wang Increase msk1 to 256000000 for nemsio files.

Parameters

[in]	lugb	integer unit of the unblocked grib data file.
[in]	lugi	integer unit of the unblocked grib index file (=0 to get index buffer from the grib file).
[in]	jf	integer maximum number of data points to unpack.
[in]	j	integer number of messages to skip (=0 to search from beginning) (<0 to read index buffer and skip -1-j messages).
[in]	jpds	integer (200) pds parameters for which to search (=-1 for wildcard). 1: id of center. 2: generating process id number. 3: grid definition. 4: gds/bms flag (right adj copy of octet 8). 5: indicator of parameter. 6: type of level. 7: height/pressure , etc of level. 8: year including (century-1). 9: month of year. 10: day of month. 11: hour of day. 12: minute of hour. 13: indicator of forecast time unit. 14: time range 1. 15: time range 2. 16: time range flag. 17: number included in average. 18: version nr of grib specification. 19: version nr of parameter table. 20: nr missing from average/accumulation. 21: century of reference time of data. 22: units decimal scale factor. 23: subcenter number. 24: pds byte 29, for nmc ensemble products. 128 if forecast field error. 64 if bias corrected fcst field. 32 if smoothed field. warning: can be combination of more than 1. 25: pds byte 30, not used.
[in]	jgds	integer (200) gds parameters for which to search (only searched if jpds(3)=255) (=-1 for wildcard). 1): data representation type. 19: number of vertical coordinate parameters. 20: octet number of the list of vertical coordinate parameters or. octet number of the list of numbers of points in each row or. 255 if neither are present. 21: for grids with pl, number of points in grid. 22: number of words in each row. Latitude/longitude grids. 2): n(i) nr points on latitude circle. 3): n(j) nr points on longitude meridian. 4): la(1) latitude of origin. 5): lo(1) longitude of origin. 6): resolution flag (right adj copy of octet 17). 7): la(2) latitude of extreme point. 8): lo(2) longitude of extreme point. 9): di longitudinal direction of increment. 10: dj latitudinal direction increment. 11: scanning mode flag (right adj copy of octet 28). Gaussian grids. 2): n(i) nr points on latitude circle. 3): n(j) nr points on longitude meridian. 4): la(1) latitude of origin. 5): lo(1) longitude of origin. 6): resolution flag (right adj copy of octet 17). 7): la(2) latitude of extreme point. 8): lo(2) longitude of extreme point. 9): di longitudinal direction of increment. 10: n - nr of circles pole to equator. 11: scanning mode flag (right adj copy of octet 28). 12: nv - nr of vert coord parameters. 13: pv - octet nr of list of vert coord parameters or pl - location of the list of numbers of points in each row (if no vert coord parameters are present) or 255 if neither are present Polar stereographic grids. 2): n(i) nr points along lat circle. 3): n(j) nr points along lon circle. 4): la(1) latitude of origin. 5): lo(1) longitude of origin. 6): resolution flag (right adj copy of octet 17). 7): lov grid orientation. 8): dx - x direction increment. 9): dy - y direction increment. 10: projection center flag. 11: scanning mode (right adj copy of octet 28). Spherical harmonic coefficients. 2): j pentagonal resolution parameter. 3): k pentagonal resolution parameter. 4): m pentagonal resolution parameter. 5): representation type. 6): coefficient storage mode. Mercator grids. 2): n(i) nr points on latitude circle. 3): n(j) nr points on longitude meridian. 4): la(1) latitude of origin. 5): lo(1) longitude of origin. 6): resolution flag (right adj copy of octet 17). 7): la(2) latitude of last grid point. 8): lo(2) longitude of last grid point. 9): latit - latitude of projection intersection. 10: reserved. 11: scanning mode flag (right adj copy of octet 28). 12: longitudinal dir grid length. 13: latitudinal dir grid length. Lambert conformal grids. 2): nx nr points along x-axis. 3): ny nr points along y-axis. 4): la1 lat of origin (lower left). 5): lo1 lon of origin (lower left). 6): resolution (right adj copy of octet 17). 7): lov - orientation of grid. 8): dx - x-dir increment. 9): dy - y-dir increment. 10: projection center flag. 11: scanning mode flag (right adj copy of octet 28). 12: latin 1 - first lat from pole of secant cone inter. 13: latin 2 - second lat from pole of secant cone inter.
[in]	mbuf	integer length of index buffer in bytes.
[in,out]	cbuf	character*1 (mbuf) index buffer (initialize by setting j=-1).
[in,out]	nlen	integer length of each index record in bytes (initialize by setting j=-1).
[in,out]	nnum	integer number of index records (initialize by setting j=-1).
[in,out]	mnum	integer number of index records skipped (initialize by setting j=-1).
[out]	kf	integer number of data points unpacked.
[out]	k	integer message number unpacked (can be same as j in calling program in order to facilitate multiple searches).
[out]	kpds	integer (200) unpacked pds parameters.
[out]	kgds	integer (200) unpacked gds parameters.
[out]	lb	logical*1 (kf) unpacked bitmap if present.
[out]	f	real (kf) unpacked data.
[out]	iret	integer return code. 0: all ok. 96: error reading index file. 97: error reading grib file. 98: number of data points greater than jf. 99: request not found. other w3fi63 grib unpacker return code.

Note

Specify an index file if feasible to increase speed. Subprogram can be called from a multiprocessing environment. Do not engage the same logical unit from more than one processor.

Author

Mark Iredell

Date

1994-04-01

Definition at line 173 of file getgbm.f.