IJKGDS1

The IJKGDS1 function decodes the grib grid description section
            and returns the field position for a given grid point.
            call IJKGDS0 to set up the navigation parameter array.

USAGE:    ...IJKGDS1(I,J,IJKGDSA)

   Input argument list:
      I        - INTEGER X GRID POINT
      J        - INTEGER Y GRID POINT
      IJKGDSA  - INTEGER (20) NAVIGATION PARAMETER ARRAY
                 IJKGDSA(1) IS NUMBER OF X POINTS
                 IJKGDSA(2) IS NUMBER OF Y POINTS
                 IJKGDSA(3) IS X WRAPAROUND INCREMENT
                            (0 IF NO WRAPAROUND)
                 IJKGDSA(4) IS Y WRAPAROUND LOWER PIVOT POINT
                            (0 IF NO WRAPAROUND)
                 IJKGDSA(5) IS Y WRAPAROUND UPPER PIVOT POINT
                            (0 IF NO WRAPAROUND)
                 IJKGDSA(6) IS SCANNING MODE
                            (0 IF X FIRST THEN Y; 1 IF Y FIRST THEN X;
                             2 IF STAGGERED DIAGONAL LIKE PROJECTION 201;
                             3 IF STAGGERED DIAGONAL LIKE PROJECTION 203)
                 IJKGDSA(7) IS MASS/WIND FLAG FOR STAGGERED DIAGONAL
                            (0 IF MASS; 1 IF WIND)
                 IJKGDSA(8:20) ARE UNUSED AT THE MOMENT
 
   Output argument list:
      IJKGDS   - INTEGER POSITION IN GRIB FIELD TO LOCATE GRID POINT
                 (0 IF OUT OF BOUNDS)
 
 
iplib.tar Library contains subroutines to be used for interpolating almost any grids used at NCEP. (Fortran90)
Date posted: 2/22/2007