GDSAWT
The GDSAWT routine decodes the grib grid description section
and returns the earth's area covered by each grid point.
The grid and earth coordinates are also returned. The
edge points and corner points are given full weight now.
USAGE: CALL GDSAWT(KGDS,KB,KA,XPTS,YPTS,RLAT,RLON,AWT)
Input argument list:
KGDS - INTEGER (200) GDS PARAMETERS AS DECODED BY W3FI63
KB - INTEGER MAXIMUM NUMBER OF COORDINATES
Output argument list:
KA - INTEGER NUMBER OF VALID POINTS COMPUTED
XPTS - REAL (KB) GRID X POINT COORDINATES
YPTS - REAL (KB) GRID Y POINT COORDINATES
RLON - REAL (KB) EARTH LONGITUDES IN DEGREES
RLAT - REAL (KB) EARTH LATITUDES IN DEGREES N
AWT - REAL (KB) AREA WEIGHTS IN M**2
iplib.tar |
Library contains subroutines to be used for
interpolating almost any grids used at NCEP. (Fortran90)
Date posted: 2/22/2007 |