IPXWAFS2

The IPXWAFS2 routine transforms between the thinned wafs grids
             as used for transmitting to the aviation community
             and their full expansion as used for general interpolation
             and graphics.  The WAFS grids are latitude-longitude grids
             thinned only in the zonal direction as indicated
             by the PL parameters in the grib grid description section.
             The PL parameters must be supplied for contraction
             (starting at KGDS1(22)).  Otherwise (if KGDS1(22)<=0)
             grid contraction is performed specifically for
             The 1.25 degree NCEP WAFS grid identifications 37-44.
             The expansion and contraction of the fields are done
             by linear interpolation, so that they are not reversible.
             This version allows a bitmap.

USAGE:    CALL IPXWAFS2(IDIR,M1,M2,KM,
                        KGDS1,IB1,L1,F1,KGDS2,IB2,L2,F2,IRET)

   Input argument list:
      IDIR     - INTEGER TRANSFORM OPTION
                 (+1 TO EXPAND THINNED FIELDS TO FULL FIELDS)
                 (-1 TO CONTRACT FULL FIELDS TO THINNED FIELDS)
      M1       - INTEGER SKIP NUMBER BETWEEN THINNED GRID FIELDS
      M2       - INTEGER SKIP NUMBER BETWEEN FULL GRID FIELDS
      KM       - INTEGER NUMBER OF FIELDS TO TRANSFORM
      KGDS1    - INTEGER (200) GDS PARMS OF THINNED GRID IF IDIR>0
                 (IF IDIR<0, THEN EITHER THE PL PARAMETERS STARTING AT
                  KGDS1(22) MUST BE SUPPLIED OR IF KGDS1(22)<=0,
                  THEN THE PL PARAMETERS DEFAULT TO THOSE FOR
                  SPECIFIC NCEP WAFS GRIDS 37-44).
      IB1      - INTEGER (KM) THINNED BITMAP FLAGS IF IDIR>0
      L1       - LOGICAL*1 (M1,KM) THINNED BITMAP FIELDS IF IDIR>0
      F1       - REAL (M1,KM) THINNED GRID FIELDS IF IDIR>0
      KGDS2    - INTEGER (200) GDS PARMS OF FULL GRID IF IDIR<0
      IB2      - INTEGER (KM) FULL BITMAP FLAGS IF IDIR<0
      L2       - LOGICAL*1 (M1,KM) FULL BITMAP FIELDS IF IDIR<0
      F2       - REAL (M2,KM) FULL GRID FIELDS IF IDIR<0
 
   Output argument list:
      KGDS1    - INTEGER (200) GDS PARMS OF THINNED GRID IF IDIR<0
      IB1      - INTEGER (KM) THINNED BITMAP FLAGS IF IDIR<0
      L1       - LOGICAL*1 (M1,KM) THINNED BITMAP FIELDS IF IDIR<0
      F1       - REAL (M1,KM) THINNED GRID FIELDS IF IDIR<0
      KGDS2    - INTEGER (200) GDS PARMS OF FULL GRID IF IDIR>0
      IB2      - INTEGER (KM) FULL BITMAP FLAGS IF IDIR>0
      L2       - LOGICAL*1 (M1,KM) FULL BITMAP FIELDS IF IDIR>0
      F2       - REAL (M2,KM) FULL GRID FIELDS IF IDIR>0
      IRET     - INTEGER RETURN CODE
                 0    SUCCESSFUL TRANSFORMATION
                 1    IMPROPER GRID SPECIFICATION
 
 
iplib.tar Library contains subroutines to be used for interpolating almost any grids used at NCEP. (Fortran90)
Date posted: 2/22/2007