SPTRANF0
The SPTRANF0 routine performs an initialization for
Subprogram SPTRANF. Use this subprogram outside
the SPTRANF family context at your own risk.
USAGE: CALL SPTRANF0(IROMB,MAXWV,IDRT,IMAX,JMAX,JB,JE,
EPS,EPSTOP,ENN1,ELONN1,EON,EONTOP,
AFFT,CLAT,SLAT,WLAT,PLN,PLNTOP)
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
Output arguments:
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)
AFFT - REAL(8) (50000+4*IMAX) AUXILIARY ARRAY IF IDIR=0
CLAT - REAL (JB:JE) COSINES OF LATITUDE
SLAT - REAL (JB:JE) SINES OF LATITUDE
WLAT - REAL (JB:JE) GAUSSIAN WEIGHTS
PLN - REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2,JB:JE)
LEGENDRE POLYNOMIALS
PLNTOP - REAL (MAXWV+1,JB:JE) LEGENDRE POLYNOMIAL OVER TOP
SPLIB.tar |
Library contains routines to be
be used for a variety of spectral transform functions. (Fortran90)
Date posted: 2/23/2007 |