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

Make array used by grib packer for gds. More...

Go to the source code of this file.

Functions/Subroutines

subroutine w3fi71 (IGRID, IGDS, IERR)
 W3FI71 Makes a 18, 37, 55, 64, or 91 word integer array used by w3fi72() grib packer to make the grid description section (gds) - section 2. More...
 

Detailed Description

Make array used by grib packer for gds.

Author
Ralph Jones
Date
1992-02-21

Definition in file w3fi71.f.

Function/Subroutine Documentation

◆ w3fi71()

subroutine w3fi71 ( integer  IGRID,
integer, dimension (*)  IGDS,
  IERR 
)

W3FI71 Makes a 18, 37, 55, 64, or 91 word integer array used by w3fi72() grib packer to make the grid description section (gds) - section 2.

Program history log:

  • Ralph Jones 1992-02-21
  • M. Farley 1992-07-01 Added remarks for 'igds' array elements. added lambert conformal grids and enlarged idgs array from 14 to 18 words.
  • Ralph Jones 1992-10-03 Added corrections to awips grib tables
  • Ralph Jones 1992-10-16 Add gaussian grid 126 to tables
  • Ralph Jones 1992-10-18 Corrections to lambert conformal tables and other tables
  • Ralph Jones 1992-10-19 Add gaussian grid 98 to tables
  • Ralph Jones 1993-01-25 Add on84 grids 87, 106, 107 to tables
  • Ralph Jones 1993-03-10 Add on84 grids 1, 55, 56 to tables
  • Ralph Jones 1993-03-26 Add grib grids 2, 3 to tables
  • Ralph Jones 1993-03-29 Add save statement
  • Ralph Jones 1993-06-15 Add grib grids 37 to 44 to tables
  • Ralph Jones 1993-09-29 Gaussian grid document not correct, w3fi74 will be changed to agree with it. gaussian grid 98 table has wrong value.
  • Ralph Jones 1993-10-12 Changes for on388 rev. oct 8,1993 for grid 204, 208.
  • Ralph Jones 1993-10-13 Correction for grids 37-44, bytes 7-8, 24-25 set to all bits 1 for missing.
  • Ralph Jones 1993-11-23 Add grids 90-93 for eta model add grid 4 for 720*361 .5 deg. grid
  • Ralph Jones 1994-04-12 Correction for grid 28
  • Ralph Jones 1994-06-01 Add grid 45, 288*145 1.25 deg. grid
  • Ralph Jones 1994-06-22 Add grids 94, 95 for eta model
  • Ralph Jones 1995-04-11 Add grids 96, 97 for eta model
  • Ralph Jones 1995-05-19 Add from 20 km eta model awips grid 215
  • Ralph Jones 1995-10-19 Add from 20 km eta model alaska grid 216
  • Mark Iredell 1995-10-31 Removed saves and prints
  • Mark Iredell 1996-05-08 Correct first latitude for grids 27 and 28
  • Ralph Jones 1996-07-02 Add from 10 km eta model olympic grid 218
  • Ralph Jones 1996-07-02 Add 196 for eta model
  • Ralph Jones 1996-08-15 Add o.n. 84 grid 8 and 53 as grib grid 8 and 53
  • Ralph Jones 1996-11-29 Correction to tables for grid 21-26, 61-64
  • Mark Iredell 1997-01-31 Correct first latitude for grid 30
  • Mark Iredell 1997-10-20 Correct last longitude for grid 98
  • Stephen Gilbert 1998-07-07 Add grids 217 and 219 through 235
  • Baldwin 1998-09-21 Add grids 190, 192 for eta model
  • Bladwin 1999-01-20 Add grids 236, 237
  • Mark Iredell 1999-08-18 Add grid 170
  • 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
  • Boi Vuong 2001-04-02 Correct last longitude for grid 225
  • Eric Rogers 2001-05-03 Added grid 249
  • Eric Rogers 2001-10-10 Redefined 218 for 12-km eta redefined grid 192 for new 32-km eta grid
  • Boi Vuong 2002-03-27 Added rsas grid 88 and awips grids 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 Added grids 145 and 146 for cmaq and grid 175 for awips over guam.
  • Boi Vuong 2003-07-08 Corrected latitude for grid 253 and 170, add grid 110, 127, 171 and 172
  • Boi Vuong 2004-08-05 Corrected latitude for grid 253
  • Stephen Gilbert 2004-09-01 Corrected the orientation and projection center flag for southern hemisphere grids 28, 172, 220 and 224
  • Boi Vuong 2004-09-02 Added grids 147, 148, 173 and 254
  • Matt Cooke 2005-01-04 Added grids 160, 161 and corrected longitude of orientation for grid 172
  • Boi Vuong 2005-03-03 Moved grid 170 to grid 174 and add grid 170
  • Boi Vuong 2005-03-21 Added grids 130
  • Boi Vuong 2005-09-12 Added grids 163
  • Boi Vuong 2006-10-27 Corrected x and y-direction grid length for grids 252
  • Boi Vuong 2006-11-16 Changed the longitude from negative to positive degree for grids 252
  • Boi Vuong 2006-12-12 Changed data representation type (octet 6) from 0 to 1 for grid 254 add grid 120 (curvilinear orthogonal grid)
  • Boi Vuong 2006-12-27 Corrected the lat/lon direction increment for grid 160
  • Boi Vuong 2007-03-21 Corrected the lat/lon direction increment, resoulution, scanning mode for grid 235 and grid type 204 for grid 120
  • Boi Vuong 2007-04-24 Corrected the lat/lon direction increment, resoulution, for grids (219,173,220,171,233,238,239,244,253) and added grid 176.
  • Boi Vuong 2007-06-11 Added new grids (11,12,13,14,15,16,18,122,123,124,125,138 180, 181, 182, 183) and corrected the lat/lon direction increment for grid 240.
  • Boi Vuong 2007-11-06 Corrected the scanning mode for grids (11,12,13,14,15,16,18) changed grid 198 from arakawa staggered e-grid to polar stereographic 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 2010-02-15 Modified to correct latitude for grid 151 and added
  • Boi Vuong 2010-06-01 Modified to correct latitude and longitude for grid 196
  • 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 grid 132, 200
  • Boi Vuong 2012-11-07 Corrected grid 174 for res. and comp. flag set to 128
  • Boi Vuong 2017-07-17 Correct grid 161 number of point nj from 102 to 103 and map size from 13974 to 14111
  • Boi Vuong 2020-06-15 Corrected grid 200,212,216 and 236 for res. and comp. flag set to 136 and south pole to -90.00
Parameters
[in]IGRIDGRIB grid number, or office note 84 grid number
[out]IGDS18, 37, 55, 64, or 91 word integer array with information to make a grib grid description section.
[out]IERR
  • 0 Correct exit
  • 1 Grid type in igrid is not in table
Note
  • 1) Office note grid type 26 is 6 in grib, 26 is an international exchange grid.
  • 2) Values returned in 18, 37, 55, 64, or 91 word integer array igds vary depending on grid representation type.
  • LAT/LON GRID:
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6)
    • IGDS( 4) = no. of points along a latitude
    • IGDS( 5) = no. of points along a longitude meridian
    • IGDS( 6) = latitude of origin (south - ive)
    • IGDS( 7) = longitude of origin (west -ive)
    • IGDS( 8) = resolution flag (code table 7)
    • IGDS( 9) = latitude of extreme point (south - ive)
    • IGDS(10) = longitude of extreme point (west - ive)
    • IGDS(11) = latitude increment
    • IGDS(12) = longitude increment
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = ... through ...
    • IGDS(18) = ... not used for this grid
    • IGDS(19) - igds(91) for grids 37-44, number of points
    • in each of 73 rows.
  • GAUSSIAN GRID:
    • IGDS( 1) = ... through ...
    • IGDS(10) = ... same as lat/lon grid
    • IGDS(11) = number of latitude lines between a pole
    • and the equator
    • IGDS(12) = longitude increment
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = ... through ...
    • IGDS(18) = ... not used for this grid
  • SPHERICAL HARMONICS:
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6)
    • IGDS( 4) = j - pentagonal resolution parameter
    • IGDS( 5) = k - pentagonal resolution parameter
    • IGDS( 6) = m - pentagonal resolution parameter
    • IGDS( 7) = representation type (code table 9)
    • IGDS( 8) = representation mode (code table 10)
    • IGDS( 9) = ... through ...
    • IGDS(18) = ... not used for this grid
  • POLAR STEREOGRAPHIC:
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6)
    • IGDS( 4) = no. of points along x-axis
    • IGDS( 5) = no. of points along y-axis
    • IGDS( 6) = latitude of origin (south -ive)
    • IGDS( 7) = longitute of origin (west -ive)
    • IGDS( 8) = resolution flag (code table 7)
    • IGDS( 9) = longitude of meridian parallel to y-axis
    • IGDS(10) = x-direction grid length (increment)
    • IGDS(11) = y-direction grid length (increment)
    • IGDS(12) = projection center flag (0=north pole on plane,
    • 1=south pole on plane,
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = ... through ...
    • IGDS(18) = .. not used for this grid
  • MERCATOR:
    • IGDS( 1) = ... through ...
    • IGDS(12) = ... same as lat/lon grid
    • IGDS(13) = latitude at which projection cylinder
    • intersects earth
    • IGDS(14) = scanning mode flags
    • IGDS(15) = ... through ...
    • IGDS(18) = .. not used for this grid
  • LAMBERT CONFORMAL:
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6)
    • IGDS( 4) = no. of points along x-axis
    • IGDS( 5) = no. of points along y-axis
    • IGDS( 6) = latitude of origin (south -ive)
    • IGDS( 7) = longitute of origin (west -ive)
    • IGDS( 8) = resolution flag (code table 7)
    • IGDS( 9) = longitude of meridian parallel to y-axis
    • IGDS(10) = x-direction grid length (increment)
    • IGDS(11) = y-direction grid length (increment)
    • IGDS(12) = projection center flag (0=north pole on plane,
    • 1=south pole on plane,
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = not used
    • IGDS(15) = first latitude from the pole at which the
    • secant cone cuts the sperical earth
    • IGDS(16) = second latitude ...
    • IGDS(17) = latitude of south pole (millidegrees)
    • IGDS(18) = longitude of south pole (millidegrees)
  • ARAKAWA SEMI-STAGGERED E-GRID ON ROTATED LAT/LON GRID
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6) [201]
    • IGDS( 4) = ni - total number of actual data points
    • included on grid
    • IGDS( 5) = nj - dummy second dimension; set=1
    • IGDS( 6) = la1 - latitude of first grid point
    • IGDS( 7) = lo1 - longitude of first grid point
    • IGDS( 8) = resolution and component flag (code table 7)
    • IGDS( 9) = la2 - number of mass points along
    • southernmost row of grid
    • IGDS(10) = lo2 - number of rows in each column
    • IGDS(11) = di - longitudinal direction increment
    • IGDS(12) = dj - latitudinal direction increment
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = ... through ...
    • IGDS(18) = ... not used for this grid (set to zero)
  • ARAKAWA FILLED E-GRID ON ROTATED LAT/LON GRID
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6) [202]
    • IGDS( 4) = ni - total number of actual data points
    • included on grid
    • IGDS( 5) = nj - dummy second dimention; set=1
    • IGDS( 6) = la1 - latitude latitude of first grid point
    • IGDS( 7) = lo1 - longitude of first grid point
    • IGDS( 8) = resolution and component flag (code table 7)
    • IGDS( 9) = la2 - number of (zonal) points in each row
    • IGDS(10) = lo2 - number of (meridional) points in each
    • column
    • IGDS(11) = di - longitudinal direction increment
    • IGDS(12) = dj - latitudinal direction increment
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = ... through ...
    • IGDS(18) = ... not used for this grid
  • ARAKAWA STAGGERED E-GRID ON ROTATED LAT/LON GRID
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6) [203]
    • IGDS( 4) = ni - number of data points in each row
    • IGDS( 5) = nj - number of rows
    • IGDS( 6) = la1 - latitude of first grid point
    • IGDS( 7) = lo1 - longitude of first grid point
    • IGDS( 8) = resolution and component flag (code table 7)
    • IGDS( 9) = la2 - central latitude
    • IGDS(10) = lo2 - central longtitude
    • IGDS(11) = di - longitudinal direction increment
    • IGDS(12) = dj - latitudinal direction increment
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = ... through ...
    • IGDS(18) = ... not used for this grid
  • CURVILINEAR ORTHOGONAL GRID
    • IGDS( 1) = number of vertical coordinates
    • IGDS( 2) = pv, pl or 255
    • IGDS( 3) = data representation type (code table 6) [204]
    • IGDS( 4) = ni - number of data points in each row
    • IGDS( 5) = nj - number of rows
    • IGDS( 6) = reserved (set to 0)
    • IGDS( 7) = reserved (set to 0)
    • IGDS( 8) = resolution and component flag (code table 7)
    • IGDS( 9) = reserved (set to 0)
    • IGDS(10) = reserved (set to 0)
    • IGDS(11) = reserved (set to 0)
    • IGDS(12) = reserved (set to 0)
    • IGDS(13) = scanning mode flags (code table 8)
    • IGDS(14) = ... through ...
    • IGDS(18) = ... not used for this grid

Definition at line 293 of file w3fi71.f.