SPTGPTV
The SPTGPTV routine performs a spherical transform
from spectral coefficients of divergences and curls
to specified sets of station point vectors 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 SPTGPTV(IROMB,MAXWV,KMAX,NMAX,
KWSKIP,KGSKIP,NRSKIP,NGSKIP,
RLAT,RLON,WAVED,WAVEZ,UP,VP)
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
WAVED - REAL (*) WAVE DIVERGENCE FIELDS
WAVEZ - REAL (*) WAVE VORTICITY FIELDS
Output arguments:
UP - REAL (*) STATION POINT U-WIND SETS
VP - REAL (*) STATION POINT V-WIND SETS
SPLIB.tar |
Library contains routines to be
be used for a variety of spectral transform functions. (Fortran90)
Date posted: 2/23/2007 |