IPXWAFS
The IPXWAFS 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.
USAGE: CALL IPXWAFS(IDIR,M1,M2,KM,KGDS1,F1,KGDS2,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).
F1 - REAL (M1,KM) THINNED GRID FIELDS IF IDIR>0
KGDS2 - INTEGER (200) GDS PARMS OF FULL GRID 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
F1 - REAL (M1,KM) THINNED GRID FIELDS IF IDIR<0
KGDS2 - INTEGER (200) GDS PARMS OF FULL GRID 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 |