NCEPLIBS-w3emc  2.9.3
Functions/Subroutines
w3fi63.f File Reference

Unpk grib field to grib grid. More...

Go to the source code of this file.

Functions/Subroutines

subroutine fi631 (MSGA, KPTR, KPDS, KRET)
 Find 'grib' chars & reset pointers. More...
 
subroutine fi632 (MSGA, KPTR, KPDS, KRET)
 Gather info from product definition sec. More...
 
subroutine fi633 (MSGA, KPTR, KGDS, KRET)
 Extract info from grib-gds. More...
 
subroutine fi634 (MSGA, KPTR, KPDS, KGDS, KBMS, KRET)
 Extract or generate bit map for output. More...
 
subroutine fi634x (NPTS, NSKP, MSGA, KBMS)
 Extract bit map. More...
 
subroutine fi635 (MSGA, KPTR, KPDS, KGDS, KBMS, DATA, KRET)
 Extract grib data elements from bds. More...
 
subroutine fi636 (DATA, MSGA, KBMS, REFNCE, KPTR, KPDS, KGDS)
 Process second order packing. More...
 
subroutine fi637 (J, KPDS, KGDS, KRET)
 Grib grid/size test. More...
 
subroutine w3fi63 (MSGA, KPDS, KGDS, KBMS, DATA, KPTR, KRET)
 Unpack a grib (edition 1) field to the exact grid specified in the grib message, isolate the bit map, and make the values of the product descripton section (pds) and the grid description section (gds) available in return arrays. More...
 

Detailed Description

Unpk grib field to grib grid.

Author
Bill Cavanaugh
Date
1991-09-13

Definition in file w3fi63.f.

Function/Subroutine Documentation

◆ fi631()

subroutine fi631 ( character*1, dimension(*)  MSGA,
integer, dimension(*)  KPTR,
integer, dimension(*)  KPDS,
integer  KRET 
)

Find 'grib' chars & reset pointers.

Author
Bill Cavanaugh
Date
1991-09-13 Find 'grib; characters and set pointers to the next byte following 'grib'. If they exist extract counts from gds and bms. Extract count from bds. Determine if sum of counts actually places terminator '7777' at the correct location.

Program history log:

  • Bill Cavanaugh 1991-09-13
  • Mark Iredell 1995-10-31 Removed saves and prints.
Parameters
[in]MSGAGrib field - "grib" thru "7777"
[in,out]KPTRArray containing storage for following parameters
  • 1 Total length of grib message
  • 2 Length of indicator (section 0)
  • 3 Length of pds (section 1)
  • 4 Length of gds (section 2)
  • 5 Length of bms (section 3)
  • 6 Length of bds (section 4)
  • 7 Value of current byte
  • 8 Bit pointer
  • 9 Grib start bit nr
  • 10 Grib/grid element count
  • 11 Nr unused bits at end of section 3
  • 12 Bit map flag
  • 13 Nr unused bits at end of section 2
  • 14 Bds flags
  • 15 Nr unused bits at end of section 4
[out]KPDSArray containing pds elements.
  • 1 Id of center
  • 2 Model identification
  • 3 Grid identification
  • 4 Gds/bms flag
  • 5 Indicator of parameter
  • 6 Type of level
  • 7 Height/pressure , etc of level
  • 8 Year of century
  • 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
[out]KRETError return
Note
ERROR RETURNS KRET:
  • 1 NO 'GRIB'
  • 2 NO '7777' OR MISLOCATED (BY COUNTS)
Author
Bill Cavanaugh
Date
1991-09-13

Definition at line 848 of file w3fi63.f.

◆ fi632()

subroutine fi632 ( character*1, dimension(*)  MSGA,
integer, dimension(*)  KPTR,
integer, dimension(*)  KPDS,
integer  KRET 
)

Gather info from product definition sec.

Author
Bill Cavanaugh
Date
1991-09-13 Extract information from the product description sec , and generate label information to permit storage in office note 84 format.

Program history log:

  • Bill Cavanaugh 1991-09-13
  • Bill Cavanaugh 1993-12-08 Corrected test for edition number instead of version number.
  • Mark Iredell 1995-10-31 Removed saves and prints.
  • M. Baldwin 1999-01-20 Modified to handle grid 237.
Parameters
[in]MSGAArray containing grib message.
[in,out]KPTRArray containing storage for following parameters.
  • 1 Total length of grib message
  • 2 Length of indicator (section 0)
  • 3 Length of pds (section 1)
  • 4 Length of gds (section 2)
  • 5 Length of bms (section 3)
  • 6 Length of bds (section 4)
  • 7 Value of current byte
  • 8 Bit pointer
  • 9 Grib start bit nr
  • 10 Grib/grid element count
  • 11 Nr unused bits at end of section 3
  • 12 Bit map flag
  • 13 Nr unused bits at end of section 2
  • 14 Bds flags
  • 15 Nr unused bits at end of section 4
[out]KPDSArray containing pds elements.
  • 1 Id of center
  • 2 Model identification
  • 3 Grid identification
  • 4 Gds/bms flag
  • 5 Indicator of parameter
  • 6 Type of level
  • 7 Height/pressure , etc of level
  • 8 Year of century
  • 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
  • 19
  • 20 Number missing from avgs/accumulations
  • 21 Century
  • 22 Units decimal scale factor
  • 23 Subcenter
[out]KRETError return.
Note
ERROR RETURN:
  • 0 - NO ERRORS
  • 8 - TEMP GDS INDICATED, BUT NO GDS
Author
Bill Cavanaugh
Date
1991-09-13

Definition at line 1005 of file w3fi63.f.

◆ fi633()

subroutine fi633 ( character*1, dimension(*)  MSGA,
integer, dimension(*)  KPTR,
integer, dimension(*)  KGDS,
integer  KRET 
)

Extract info from grib-gds.

Author
Bill Cavanaugh
Date
1991-09-13 Extract information on unlisted grid to allow conversion to office note 84 format.

Program history log:

  • Bill Cavanaugh 1991-09-13
  • M. Baldwin 1995-03-20 fi633 modification to get data rep types [kgds(1)] 201 and 202 to work.
  • Mark Iredell 1995-10-31 Removed saves and prints
  • M. Baldwin 1998-09-08 Add data rep type [kgds(1)] 203
  • Boi Vuong 2007-04-24 Add data rep type [kgds(1)] 204
  • George Gayno 2010-07-20 Add data rep type [kgds(1)] 205
Parameters
[in]MSGAArray containing grib message
[in,out]KPTRArray containing storage for following parameters
  • 1 Total length of grib message
  • 2 Length of indicator (section 0)
  • 3 Length of pds (section 1)
  • 4 Length of gds (section 2)
  • 5 Length of bms (section 3)
  • 6 Length of bds (section 4)
  • 7 Value of current byte
  • 8 Bit pointer
  • 9 Grib start bit nr
  • 10 Grib/grid element count
  • 11 Nr unused bits at end of section 3
  • 12 Bit map flag
  • 13 Nr unused bits at end of section 2
  • 14 Bds flags
  • 15 Nr unused bits at end of section 4
[out]KGDSArray containing gds elements.
  • 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
  • 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
    • 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
  • 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) Reserved
    • 7) Lov grid orientation
    • 8) Dx - x direction increment
    • 9) Dy - y direction increment
    • 10 Projection center flag
    • 11 Scanning mode
  • 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
    • 7 La(2) latitude of last grid point
    • 8 Lo(2) longitude of last grid point
    • 9 Latin - latitude of projection intersection
    • 10 Reserved
    • 11 Scanning mode flag
    • 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
    • 12 Latin 1 - first lat from pole of secant cone inter
    • 13 Latin 2 - second lat from pole of secant cone inter
  • Staggered arakawa rotated lat/lon grids (203 e stagger)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
  • Staggered arakawa rotated lat/lon grids (205 a,b,c,d staggers)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
    • 12 Latitude of last point
    • 13 Longitude of last point
      Parameters
      [out]KRETError return
      Note
  • KRET
    • 0
    • 4 - Data representation type not currently acceptable
Author
Bill Cavanaugh
Date
1991-09-13

Definition at line 1351 of file w3fi63.f.

◆ fi634()

subroutine fi634 ( character*1, dimension(*)  MSGA,
integer, dimension(*)  KPTR,
integer, dimension(*)  KPDS,
integer, dimension(*)  KGDS,
logical*1, dimension(*)  KBMS,
integer  KRET 
)

Extract or generate bit map for output.

Author
Bill Cavanaugh
Date
1991-09-13 If bit map sec is available in grib message, extract for program use, otherwise generate an appropriate bit map.

Program history log:

  • Bill Cavanaugh 1991-09-13
  • Bill Cavanaugh 1991-11-12 Modified size of ecmwf grids 5 - 8.
  • Mark Iredell 1995-10-31 removed saves and prints
  • W. Bostelman 1997-02-12 corrects ecmwf us grid 2 processing
  • Mark Iredell 1997-09-19 vectorized bitmap decoder
  • Stephen Gilbert 1998-09-02 corrected error in map size for u.s. grid 92
  • M. Baldwin 1998-09-08 add grids 190,192
  • M. Baldwin 1999-01-20 add grids 236,237
  • Eric Rogers 2001-10-02 redefined grid #218 for 12 km eta redefined grid 192 for new 32-km eta grid
  • Stephen Gilbert 2003-06-30 added grids 145 and 146 for cmaq and grid 175 for awips over guam.
  • Boi Vuong 2004-09-02 Added awips grids 147, 148, 173 and 254
  • Boi Vuong 2006-12-12 Added awips grids 120
  • Boi Vuong 2007-04-20 Added awips grids 176
  • Boi Vuong 2007-06-11 Added awips grids 11 to 18 and 122 to 125 and 180 to 183
  • Boi Vuong 2010-08-05 Added new grid 184, 199, 83 and redefined grid 90 for new rtma conus 1.27-km redefined grid 91 for new rtma alaska 2.976-km redefined grid 92 for new rtma alaska 1.488-km
  • Boi Vuong 2012-02-28 Added new grid 200
Parameters
[in]MSGABufr message
[in,out]KPTRArray containing storage for following parameters
  • 1 Total length of grib message
  • 2 Length of indicator (section 0)
  • 3 Length of pds (section 1)
  • 4 Length of gds (section 2)
  • 5 Length of bms (section 3)
  • 6 Length of bds (section 4)
  • 7 Value of current byte
  • 8 Bit pointer
  • 9 Grib start bit nr
  • 10 Grib/grid element count
  • 11 Nr unused bits at end of section 3
  • 12 Bit map flag
  • 13 Nr unused bits at end of section 2
  • 14 Bds flags
  • 15 Nr unused bits at end of section 4
[in]KPDSArray containing pds elements.
  • 1 Id of center
  • 2 Model identification
  • 3 Grid identification
  • 4 Gds/bms flag
  • 5 Indicator of parameter
  • 6 Type of level
  • 7 Height/pressure , etc of level
  • 8 Year of century
  • 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
[in]KGDSArray containing gds elements.
  • 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
  • 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
    • 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
  • 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) Reserved
    • 7) Lov grid orientation
    • 8) Dx - x direction increment
    • 9) Dy - y direction increment
    • 10 Projection center flag
    • 11 Scanning mode
  • 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
    • 7 La(2) latitude of last grid point
    • 8 Lo(2) longitude of last grid point
    • 9 Latin - latitude of projection intersection
    • 10 Reserved
    • 11 Scanning mode flag
    • 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
    • 12 Latin 1 - first lat from pole of secant cone inter
    • 13 Latin 2 - second lat from pole of secant cone inter
  • Staggered arakawa rotated lat/lon grids (203 e stagger)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
  • Staggered arakawa rotated lat/lon grids (205 a,b,c,d staggers)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
    • 12 Latitude of last point
    • 13 Longitude of last point
      Parameters
      [out]KBMSBitmap describing location of output elements.
      [out]KRETError return
      Note
  • KRET
    • 0 - No error
    • 5 - Grid not avail for center indicated
    • 10 - Incorrect center indicator
    • 12 - Bytes 5-6 are not zero in bms, predefined bit map not provided by this center
Author
Bill Cavanaugh
Date
1991-09-13

Definition at line 1897 of file w3fi63.f.

◆ fi634x()

subroutine fi634x (   NPTS,
  NSKP,
character*1, dimension(*)  MSGA,
logical*1, dimension(npts)  KBMS 
)

Extract bit map.

Author
Mark Iredell
Date
1997-09-19 Extract the packed bitmap into a logical array.

Program history log: 97-09-19 Vectorized bitmap decoder.

Parameters
[in]NPTSXInteger number of points in the bitmap field
[in]NSKPInteger number of bits to skip in grib message
[in]MSGACharacter*1 grib message
[out]KBMSLogical*1 bitmap
Note
Subprogram can be called from a multiprocessing environment.
Author
Mark Iredell
Date
1997-09-19

Definition at line 2882 of file w3fi63.f.

◆ fi635()

subroutine fi635 ( character*1, dimension(*)  MSGA,
integer, dimension(*)  KPTR,
integer, dimension(*)  KPDS,
integer, dimension(*)  KGDS,
logical*1, dimension(*)  KBMS,
real, dimension(*)  DATA,
  KRET 
)

Extract grib data elements from bds.

Author
Bill Cavanaugh
Date
1991-09-13 Extract grib data from binary data section and place into output array in proper position.

Program history log:

  • Bill Cavanaugh 1991-09-13
  • Bill Cavanaugh 1994-04-01 Modified code to include decimal scaling when calculating the value of data points specified as being equal to the reference value
  • Farley 1994-11-10 Increased mxsize from 72960 to 260000 for .5 degree sst analysis fields.
  • Mark Iredell 1995-10-31 Removed saves and prints
  • Mark Iredell 1998-08-31 Eliminated need for mxsize
Parameters
[in]MSGAArray containing grib message
[in,out]KPTRArray containing storage for following parameters
  • 1 Total length of grib message
  • 2 Length of indicator (section 0)
  • 3 Length of pds (section 1)
  • 4 Length of gds (section 2)
  • 5 Length of bms (section 3)
  • 6 Length of bds (section 4)
  • 7 Value of current byte
  • 8 Bit pointer
  • 9 Grib start bit nr
  • 10 Grib/grid element count
  • 11 Nr unused bits at end of section 3
  • 12 Bit map flag
  • 13 Nr unused bits at end of section 2
  • 14 Bds flags
  • 15 Nr unused bits at end of section 4
  • 16 Reserved
  • 17 Reserved
  • 18 Reserved
  • 19 Binary scale factor
  • 20 Num bits used to pack each datum
[in]KPDSArray containing pds elements. See initial routine
[in]KGDSArray containing gds elements.
  • 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
  • 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
    • 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
  • 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) Reserved
    • 7) Lov grid orientation
    • 8) Dx - x direction increment
    • 9) Dy - y direction increment
    • 10 Projection center flag
    • 11 Scanning mode
  • 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
    • 7 La(2) latitude of last grid point
    • 8 Lo(2) longitude of last grid point
    • 9 Latin - latitude of projection intersection
    • 10 Reserved
    • 11 Scanning mode flag
    • 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
    • 12 Latin 1 - first lat from pole of secant cone inter
    • 13 Latin 2 - second lat from pole of secant cone inter
  • Staggered arakawa rotated lat/lon grids (203 e stagger)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
  • Staggered arakawa rotated lat/lon grids (205 a,b,c,d staggers)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
    • 12 Latitude of last point
    • 13 Longitude of last point
      Parameters
      [in]KBMSBitmap describing location of output elements. -KBDS Information extracted from binary data section
    • KBDS(1) - N1
    • KBDS(2) - N2
    • KBDS(3) - P1
    • KBDS(4) - P2
    • KBDS(5) - Bit pointer to 2nd order widths
    • KBDS(6) - Bit pointer to 2nd order bit maps
    • KBDS(7) - Bit pointer to first order values
    • KBDS(8) - Bit pointer to second order values
    • KBDS(9) - Bit pointer start of bds
    • KBDS(10) - Bit pointer main bit map
    • KBDS(11) - Binary scaling
    • KBDS(12) - Decimal scaling
    • KBDS(13) - Bit width of first order values
    • KBDS(14) - Bit map flag 0 = no second order bit map 1 = second order bit map present
  • KBDS(15) - Second order bit width
  • KBDS(16) - Constant / different widths 0 = constant widths 1 = different widths
  • KBDS(17) - Single datum / matrix
    • 0 = single datum at each grid point
    • 1 = matrix of values at each grid point
    • (18-20) - Unused
      Parameters
      [out]DATAReal*4 array of gridded elements in grib message.
      [out]KRETError return
      Note
  • Error return
    • 3 = Unpacked field is larger than 65160
    • 6 = Does not match nr of entries for this grib/grid
    • 7 = Number of bits in fill too large
Author
Bill Cavanaugh
Date
1991-09-13

Definition at line 3056 of file w3fi63.f.

◆ fi636()

subroutine fi636 ( real, dimension(*)  DATA,
character*1, dimension(*)  MSGA,
logical*1, dimension(*)  KBMS,
real  REFNCE,
integer, dimension(*)  KPTR,
integer, dimension(*)  KPDS,
integer, dimension(*)  KGDS 
)

Process second order packing.

Author
Bill Cavanaugh
Date
1992-09-22 Process second order packing from the binary data section (bds) for single data items grid point data.

Program history log:

  • Bill Cavanaugh 1993-06-08
  • Bill Cavanaugh 1993-12-15 Modified second order pointers to first order values and second order values correctly.
  • Ralph Jones 1995-04-26 Fi636 corection for 2nd order complex Unpacking.
  • Mark Iredell 1995-10-31 Saves and prints.
Parameters
[in]MSGAArray containing grib message
[in]REFNCEReference value
[in]KPTRWork array
[out]DATALocation of output array
  • KBDS Working array
    • KBDS(1) N1
    • KBDS(2) N2
    • KBDS(3) P1
    • KBDS(4) P2
    • KBDS(5) Bit pointer to 2nd order widths
    • KBDS(6) Bit pointer to 2nd order bit maps
    • KBDS(7) Bit pointer to first order values
    • KBDS(8) Bit pointer to second order values
    • KBDS(9) Bit pointer start of bds
    • KBDS(10) Bit pointer main bit map
    • KBDS(11) Binary scaling
    • KBDS(12) Decimal scaling
    • KBDS(13) Bit width of first order values
    • KBDS(14) Bit map flag
      • 0 = No second order bit map
      • 1 = Second order bit map present
    • KBDS(15) Second order bit width
    • KBDS(16) Constant / different widths
      • 0 = Constant widths
      • 1 = Different widths
    • KBDS(17) Single datum / matrix
      • 0 = Single datum at each grid point
      • 1 = Matrix of values at each grid point
    • KBDS(18-20) Unused
[in]KBMS
[in]KPDS
[in]KGDSArray containing gds elements.
  • 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
  • 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
    • 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
  • 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) Reserved
    • 7) Lov grid orientation
    • 8) Dx - x direction increment
    • 9) Dy - y direction increment
    • 10 Projection center flag
    • 11 Scanning mode
  • 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
    • 7 La(2) latitude of last grid point
    • 8 Lo(2) longitude of last grid point
    • 9 Latin - latitude of projection intersection
    • 10 Reserved
    • 11 Scanning mode flag
    • 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
    • 12 Latin 1 - first lat from pole of secant cone inter
    • 13 Latin 2 - second lat from pole of secant cone inter
  • Staggered arakawa rotated lat/lon grids (203 e stagger)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
  • Staggered arakawa rotated lat/lon grids (205 a,b,c,d staggers)
    • 2 N(i) nr points on rotated latitude circle
    • 3 N(j) nr points on rotated longitude meridian
    • 4 La(1) latitude of origin
    • 5 Lo(1) longitude of origin
    • 6 Resolution flag
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag
    • 12 Latitude of last point
    • 13 Longitude of last point
Author
Bill Cavanaugh
Date
1992-09-22

Definition at line 3701 of file w3fi63.f.

◆ fi637()

subroutine fi637 ( integer  J,
integer, dimension(*)  KPDS,
integer, dimension(*)  KGDS,
  KRET 
)

Grib grid/size test.

Author
Bill Cavanaugh
Date
1991-09-13 To test when gds is available to see if size mismatch on existing grids (by center) is indicated.

Program history log:

  • Bill Cavanaugh 1991-09-13
  • Mark Iredell 1995-10-31 Removed saves and prints
  • M. Bostelman 1997-02-12 Corrects ecmwf us grid 2 processing
  • Mark Iredell 1998-06-17 Removed alternate return
  • M. Baldwin 1999-01-20 Modify to handle grid 237
  • Boi Vuong 1909-05-21 Modify to handle grid 45
Parameters
[in,out]JSize for indicated grid modified for ecmwf-us 2
[in]KPDS
[in]KGDS
[out]KRETError return (a mismatch was detected if kret is not zero)
Note
  • KRET:
    • 9 - Gds indicates size mismatch with std grid
Author
Bill Cavanaugh
Date
1991-09-13

Definition at line 3968 of file w3fi63.f.

◆ w3fi63()

subroutine w3fi63 ( character*1, dimension(*)  MSGA,
integer, dimension(*)  KPDS,
integer, dimension(*)  KGDS,
logical*1, dimension(*)  KBMS,
real, dimension(*)  DATA,
integer, dimension(*)  KPTR,
  KRET 
)

Unpack a grib (edition 1) field to the exact grid specified in the grib message, isolate the bit map, and make the values of the product descripton section (pds) and the grid description section (gds) available in return arrays.

When decoding is completed, data at each grid point has been returned in the units specified in the grib manual.

PROGRAM HISTORY LOG:

  • Bill Cavanaugh 1991-09-13
  • Bill Cavanaugh 1991-11-12 Modified size of ecmwf grids 5-8
  • Bill Cavanaugh 1991-12-22 Corrected processing of mercator projections in grid definition section (gds) in routine fi633
  • Bill Cavanaugh 1992-08-05 Corrected maximum grid size to allow for one degree by one degree global grids
  • Bill Cavanaugh 1992-08-27 Corrected typo error, added code to compare total byte size from section 0 with sum of section sizes.
  • Bill Cavanaugh 1992-10-21 Corrections were made (in fi634) to reduce processing time for international grids. removed a typographical error in fi635.
  • Bill Cavanaugh 1993-01-07 Corrections were made (in fi635) to facilitate use of these routines on a pc. a typographical error was also corrected
  • Bill Cavanaugh 1993-01-13 Corrections were made (in fi632) to properly handle condition when time range indicator = 10. added u.s.grid 87.
  • Bill Cavanaugh 1993-02-04 Added u.s.grids 85 and 86
  • Bill Cavanaugh 1993-02-26 Added grids 2, 3, 37 thru 44,and grids 55, 56, 90, 91, 92, and 93 to list of u.s. grids.
  • Bill Cavanaugh 1993-04-07 Added grids 67 thru 77 to list of u.s. grids.
  • Bill Cavanaugh 1993-04-20 Increased max size to accomodate gaussian grids.
  • Bill Cavanaugh 1993-05-26 Corrected grid range selection in fi634 for ranges 67-71 & 75-77
  • Bill Cavanaugh 1993-06-08 Corrected fi635 to accept grib messages with second order packing. added routine fi636 to process messages with second order packing.
  • Bill Cavanaugh 1993-09-22 Modified to extract sub-center number from pds byte 26
  • Bill Cavanaugh 1993-10-13 Modified fi634 to correct grid sizes for grids 204 and 208
  • Bill Cavanaugh 1993-10-14 Increased size of kgds to include entries for number of points in grid and number of words in each row
  • Bill Cavanaugh 1993-12-08 Corrected test for edition number instead of version number
  • Bill Cavanaugh 1993-12-15 Modified second order pointers to first order values and second order values correctly in routine fi636
  • Bill Cavanaugh 1994-03-02 Added call to w3fi83 within decoder. user no longer needs to make call to this routine
  • Bill Cavanaugh 1994-04-22 Modified fi635, fi636 to process row by row second order packing, added scaling correction to fi635, and corrected typographical errors in comment fields in fi634
  • Bill Cavanaugh 1994-05-17 COrrected error in fi633 to extract resolution for lambert-conformal grids. added clarifying information to docblock entries
  • Bill Cavanaugh 1994-05-25 Added code to process column by column as well as row by row ordering of second order data
  • Bill Cavanaugh 1994-06-27 Added processing for grids 45, 94 and 95. includes construction of second order bit maps for thinned grids in fi636.
  • Bill Cavanaugh 1994-07-08 Commented out print outs used for debugging
  • Bill Cavanaugh 1994-09-08 Added grids 220, 221, 223 for fnoc
  • Farley 1994-11-10 Increased mxsize from 72960 to 260000 for .5 degree sst analysis fields
  • Ralph Jones 1994-12-06 Changes in fi632 for pds greater than 28
  • Ralph Jones 1995-02-14 Correct in fi633 for navy wafs grib
  • M Baldwin 1995-03-20 Fi633 modification to get data rep types [kgds(1)] 201 and 202 to work.
  • M. Baldwin 1995-04-10 Added grids 96 and 97 for eta model in fi634.
  • Ralph Jones 1995-04-26 Fi636 corection for 2nd order complex unpacking. r
  • Ralph Jones 1995-05-19 Added grid 215, 20 km awips grid
  • Ralph Jones 1995-07-06 Added gaussian t62, t126 grid 98, 126
  • Ralph Jones 1995-10-19 Added grid 216, 45 km eta awips alaska grid
  • Mark Iredell 1995-10-31 Removed saves and prints
  • Ralph Jones 1996-03-07 Continue unpack with kret error 9 in fi631.
  • Ralph Jones 1996-08-19 Added mercator grids 8 and 53, and grid 196
  • W. Bostelman 1997-02-12 Corrects ecmwf us grid 2 processing
  • Mark Iredell 1998-06-17 Removed alternate return in fi637
  • Mark Iredell 1998-08-31 Eliminated need for mxsize
  • Stephen Gilbert 1998-09-02 Corrected error in map size for U.S. Grid 92
  • M. Baldwin 1998-09-08 Add data rep type [kgds(1)] 203
  • Eric Rogers 2001-03-08 Changed eta grids 90-97, added eta grids 194, 198. added awips grids 241,242,243, 245, 246, 247, 248, and 250
  • Boi Vuong 2001-03-19 Added awips grids 238,239,240, and 244
  • Stephen Gilbert 2001-06-06 Changed gbyte/sbyte calls to refer to Wesley Ebisuzaki's endian independent versions gbytec/sbytec. Removed equivalences.
  • Eric Rogers 2001-05-03 Added grid 249 (12km for alaska)
  • Eric Rogers 2001-10-10 Redefined grid 218 for 12 km eta redefined grid 192 for new 32-km eta grid
  • Boi Vuong 2002-03-27 VUONG Added rsas grid 88 and awips grids 219, 220, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 251, and 252
  • Eric Rogers 2002-08-06 Redefined grids 90-93,97,194,245-250 for the 8km hi-res-window model and add awips grid 253
  • Stephen Gilbert 2003-06-30 Set new values in array kptr to pass back additional packing info. kptr(19) - binary scale factor kptr(20) - num bits used to pack each datum
  • Stephen Gilbert 2003-06-30 Added grids 145 and 146 for cmaq and grid 175 for awips over guam.
  • Boi Vuong 2003-07-08 Added grids 110, 127, 171, 172 and modified grid 170
  • Boi Vuong 2004-09-02 Added awips grids 147, 148, 173 and 254 2005-01-04 COOKE Added awips grids 160 and 161
  • Boi Vuong 2005-03-03 Moved grid 170 to grid 174 and add grid 170
  • Boi Vuong 2005-03-21 Added awips grid 130
  • Boi Vuong 2005-10-11 Added awips grid 163
  • Boi Vuong 2006-12-12 Added awips grid 120
  • Boi Vuong 2007-04-12 Added awips 176 and data rep type kgds(1) 204
  • Boi Vuong 2007-06-11 Added new grids 11 to 18 and 122 to 125 and 138 and 180 to 183
  • Boi Vuong 2007-11-06 Changed grid 198 from arakawa staggered e-grid to polar stereograpgic grid added new grid 10, 99, 150, 151, 197
  • Boi Vuong 2008-01-17 Added new grid 195 and changed grid 196 (arakawa-e to mercator)
  • Boi Vuong 2009-05-21 Modified to handle grid 45
  • Boi Vuong 2010-05-11 Data rep type kgds(1) 205
  • Boi Vuong 2010-02-18 Added grid 128, 139 and 140 2010-07-20 Added rotated lat/lon "a,b,c,d" staggers -> kgds(1) 205
  • Boi Vuong 2010-08-05 Added new grid 184, 199, 83 and redefined grid 90 for new rtma conus 1.27-km redefined grid 91 for new rtma alaska 2.976-km redefined grid 92 for new rtma alaska 1.488-km
  • Eric Rogers 2010-09-08 Changed grid 94 to alaska 6km staggered b-grid changed grid 95 to puerto rico 3km staggered b-grid changed grid 96 to hawaii 3km staggered b-grid changed grid 96 to hawaii 3km staggered b-grid changed grid 97 to conus 4km staggered b-grid changed grid 99 to nam 12km staggered b-grid added grid 179 (12 km polar stereographic over north america) changed grid 194 to 3km mercator grid over puerto rico corrected latitude of sw corner point of grid 151
  • Boi Vuong 2011-10-12 Added grid 129, 187, 188, 189 and 193
  • Boi Vuong 2012-04-16 Added new grid 132, 200
  • Boi Vuong 2017-07-17 Correct grid 161 number of point nj from 102 to 103 and map size from 13974 to 14111
Parameters
[in]MSGAGrib field - "grib" thru "7777" char*1 (message can be preceded by junk chars)
[out]DATAArray containing data elements
[out]KPDSArray containing pds elements. (edition 1)
  • 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
  • 26-35 Reserved
  • 36-N Consecutive bytes extracted from program Definition section (pds) of grib message
[out]KGDSARRAY CONTAINING GDS ELEMENTS.
  • 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
  • E-STAGGERED ARAKAWA ROTATED LAT/LON GRIDS (TYPE 203)
    • 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 center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag (right adj copy of octet 28)
  • CURVILINEAR ORTHIGINAL GRID (TYPE 204)
    • 2 N(i) nr points on latitude circle
    • 3 N(j) nr points on longitude meridian
    • 4 Reserved set to 0
    • 5 Reserved set to 0
    • 6 Resolution flag (right adj copy of octet 17)
    • 7 Reserved set to 0
    • 8 Reserved set to 0
    • 9 Reserved set to 0
    • 10 Reserved set to 0
    • 11 Scanning mode flag (right adj copy of octet 28)
  • ROTATED LAT/LON A,B,C,D-STAGGERED (TYPE 205)
    • 2 N(i) nr points on latitude circle
    • 3 N(j) nr points on longitude meridian
    • 4 La(1) latitude of first point
    • 5 Lo(1) longitude of first point
    • 6 Resolution flag (right adj copy of octet 17)
    • 7 La(2) latitude of center
    • 8 Lo(2) longitude of center
    • 9 Di longitudinal direction of increment
    • 10 Dj latitudinal direction increment
    • 11 Scanning mode flag (right adj copy of octet 28)
    • 12 Latitude of last point
    • 13 Longitude of last point
      Parameters
      [out]KBMSBitmap describing location of output elements. (always constructed)
      [out]KPTRArray containing storage for following parameters
  • 1 Total length of grib message
  • 2 Length of indicator (section 0)
  • 3 Length of pds (section 1)
  • 4 Length of gds (section 2)
  • 5 Length of bms (section 3)
  • 6 Length of bds (section 4)
  • 7 Value of current byte
  • 8 Bit pointer
  • 9 Grib start bit nr
  • 10 Grib/grid element count
  • 11 Nr unused bits at end of section 3
  • 12 Bit map flag (copy of bms octets 5,6)
  • 13 Nr unused bits at end of section 2
  • 14 Bds flags (right adj copy of octet 4)
  • 15 Nr unused bits at end of section 4
  • 16 Reserved
  • 17 Reserved
  • 18 Reserved
  • 19 Binary scale factor
  • 20 Num bits used to pack each datum
    Parameters
    [out]KRETFlag indicating quality of completion
    Note
    When decoding is completed, data at each grid point has been returned in the units specified in the grib manual.
  • Values for return flag (kret)
    • 0 - Normal return, no errors
    • 1 - 'grib' not found in first 100 chars
    • 2 - '7777' not in correct location
    • 3 - Unpacked field is larger than 260000
    • 4 - Gds/ grid not one of currently accepted values
    • 5 - Grid not currently avail for center indicated
    • 8 - Temp gds indicated, but gds flag is off
    • 9 - Gds indicates size mismatch with std grid
    • 10 - Incorrect center indicator
    • 11 - Binary data section (bds) not completely processed. program is not set to process flag combinations shown in octets 4 and 14.
    • 12 - Binary data section (bds) not completely processed. program is not set to process flag combinations
Author
Bill Cavanaugh
Date
1991-09-13

Definition at line 347 of file w3fi63.f.