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