W3AI08

The W3AI08 routine unpacks a grib field to the exact grid specified
           in the message, isolate the bit map and make the values 
           of the product description section (PDS) and the grid
           description section (GDS) available in return arrays.

USAGE:    CALL W3AI08(MSGA,KPDS,KGDS,KBMS,DATA,KPTR,KRET)

   Input argument list:
      MSGA     - GRIB FIELD - "GRIB" THRU "7777"   CHAR*1
 
   Output argument list:
      DATA     - ARRAY CONTAINING DATA ELEMENTS
      KPDS     - ARRAY CONTAINING PDS ELEMENTS. (VERSION 0)
           (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 INCLUDING 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)  - GRIB SPECIFICATION EDITION NUMBER
      KPDS     - ARRAY CONTAINING PDS ELEMENTS.  (VERSION 1)
           (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 INCLUDING 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)  - VERSION NR OF GRIB SPECIFICATION
           (19)  - VERSION NR OF PARAMETER TABLE
           (20)  - TOTAL LENGTH OF GRIB MESSAGE (INCLUDING SECTION 0)
      KGDS     - ARRAY CONTAINING GDS ELEMENTS.
           (1)   - DATA REPRESENTATION TYPE
        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
           (7)   - LA(2) LATITUDE OF EXTREME POINT
           (8)   - LO(2) LONGITUDE OF EXTREME POINT
           (9)   - DI LONGITUDINAL DIRECTION OF INCREMENT
           (10)  - DJ LATITUNDINAL DIRECTION OF 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      "          "         "
           (4)   - M      "          "         "
           (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)   - LONGIT DIR INCREMENT
           (10)  - LATIT DIR INCREMENT
           (11)  - SCANNING MODE FLAG
           (12)  - LATITUDE INTERSECTION
        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)   - RESERVED
           (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
      KBMS       - BITMAP DESCRIBING LOCATION OF OUTPUT ELEMENTS.
      KPTR       - ARRAY CONTAINING STORAGE FOR FOLLOWING PARAMETERS
           (1)   - UNUSED
           (2)   - UNUSED
           (3)   - LENGTH OF PDS
           (4)   - LENGTH OF GDS
           (5)   - LENGTH OF BMS
           (6)   - LENGTH OF BDS
           (7)   - VALUE OF CURRENT BYTE
           (8)   - UNUSED
           (9)   - GRIB START BYTE NR
          (10)   - GRIB/GRID ELEMENT COUNT
      KRET       - FLAG INDICATING QUALITY OF COMPLETION
 
REMARKS: 
      Values for return flag (kret)
      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 32768
           = 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

 


  
W3lib.tar
  
Library contains Fortran 90 decoder/encoder
   routines for GRIB edition 1. (Fortran90) 

   Date posted: 2/22/2007