SPTRANF1
The SPTRANF1 routine performs an single latitude transform for
subprogram SPTRANF. Use this subprogram outside
the SPTRANF family context at your own risk.
USAGE: CALL SPTRANF1(IROMB,MAXWV,IDRT,IMAX,JMAX,JB,JE,
EPS,EPSTOP,ENN1,ELONN1,EON,EONTOP,
AFFT,CLAT,SLAT,WLAT,PLN,PLNTOP,MP,
W,WTOP,G,IDIR)
Input arguments:
IROMB - INTEGER SPECTRAL DOMAIN SHAPE
(0 FOR TRIANGULAR, 1 FOR RHOMBOIDAL)
MAXWV - INTEGER SPECTRAL TRUNCATION
IDRT - INTEGER GRID IDENTIFIER
(IDRT=4 FOR GAUSSIAN GRID,
IDRT=0 FOR EQUALLY-SPACED GRID INCLUDING POLES,
IDRT=256 FOR EQUALLY-SPACED GRID EXCLUDING POLES)
IMAX - INTEGER EVEN NUMBER OF LONGITUDES
JMAX - INTEGER NUMBER OF LATITUDES
JB - INTEGER LATITUDE INDEX (FROM POLE) TO BEGIN TRANSFORM
JE - INTEGER LATITUDE INDEX (FROM POLE) TO END TRANSFORM
EPS - REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
EPSTOP - REAL (MAXWV+1)
ENN1 - REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
ELONN1 - REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
EON - REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
EONTOP - REAL (MAXWV+1)
CLAT - REAL (JB:JE) COSINES OF LATITUDE
SLAT - REAL (JB:JE) SINES OF LATITUDE
WLAT - REAL (JB:JE) GAUSSIAN WEIGHTS
AFFT - REAL(8) (50000+4*IMAX) AUXILIARY ARRAY IF IDIR=0
PLN - REAL ((M+1)*((I+1)*M+2)/2,JB:JE) LEGENDRE POLYNOMIALS
PLNTOP - REAL (M+1,JB:JE) LEGENDRE POLYNOMIAL OVER TOP
MP - INTEGER IDENTIFIER (0 FOR SCALAR, 1 FOR VECTOR)
W - REAL (*) WAVE FIELD IF IDIR>0
WTOP - REAL (*) WAVE FIELD OVER TOP IF IDIR>0
G - REAL (IMAX,2,JB:JE) GRID FIELD IF IDIR<0
IDIR - INTEGER TRANSFORM FLAG
(IDIR>0 FOR WAVE TO GRID, IDIR<0 FOR GRID TO WAVE)
Output arguments:
W - REAL (*) WAVE FIELD IF IDIR<0
WTOP - REAL (*) WAVE FIELD OVER TOP IF IDIR<0
G - REAL (IMAX,2,JB:JE) GRID FIELD IF IDIR>0
SPLIB.tar |
Library contains routines to be
be used for a variety of spectral transform functions. (Fortran90)
Date posted: 2/23/2007 |