SPTGPTD
The SPTGPTD routine performs a spherical transform
from spectral coefficients of scalar fields
to specified sets of station point gradients on the globe.
The wave-space can be either Triangular or Rhomboidal.
The wave and point fields may have general indexing,
but each wave field is in sequential 'IBM order',
i.e. with zonal wavenumber as the slower index.
The transforms are all multiprocessed over stations.
transform several fields at a time to improve vectorization.
Subprogram can be called from a multiprocessing environment.
USAGE: CALL SPTGPTD(IROMB,MAXWV,KMAX,NMAX,
KWSKIP,KGSKIP,NRSKIP,NGSKIP,
RLAT,RLON,WAVE,XP,YP)
Input arguments:
IROMB - INTEGER SPECTRAL DOMAIN SHAPE
(0 FOR TRIANGULAR, 1 FOR RHOMBOIDAL)
MAXWV - INTEGER SPECTRAL TRUNCATION
KMAX - INTEGER NUMBER OF FIELDS TO TRANSFORM.
NMAX - INTEGER NUMBER OF STATION POINTS TO RETURN
KWSKIP - INTEGER SKIP NUMBER BETWEEN WAVE FIELDS
(DEFAULTS TO (MAXWV+1)*((IROMB+1)*MAXWV+2) IF KWSKIP=0)
KGSKIP - INTEGER SKIP NUMBER BETWEEN STATION POINT SETS
(DEFAULTS TO NMAX IF KGSKIP=0)
NRSKIP - INTEGER SKIP NUMBER BETWEEN STATION LATS AND LONS
(DEFAULTS TO 1 IF NRSKIP=0)
NGSKIP - INTEGER SKIP NUMBER BETWEEN STATION POINTS
(DEFAULTS TO 1 IF NGSKIP=0)
RLAT - REAL (*) STATION LATITUDES IN DEGREES
RLON - REAL (*) STATION LONGITUDES IN DEGREES
WAVE - REAL (*) WAVE FIELDS
Output arguments:
XP - REAL (*) STATION POINT X-GRADIENT SETS
YP - REAL (*) STATION POINT Y-GRADIENT SETS
SPLIB.tar |
Library contains routines to be
be used for a variety of spectral transform functions. (Fortran90)
Date posted: 2/23/2007 |