SPEPS

The SPEPS routine computes constant fields indexed in the 
          spectral domain in "IBM order" (zonal wavenumber is the
          slower index).
          If l is the zonal wavenumber and n is the total wavenumber
          and a is the earth radius, then the fields returned are:
          (1) Normalizing factor epsilon=sqrt((n**2-l**2)/(4*n**2-1))
          (2) Laplacian factor n*(n+1)/a**2
          (3) Zonal derivative/laplacian factor l/(n*(n+1))*a
          (4) Meridional derivative/laplacian factor epsilon/n*a

USAGE:    CALL SPEPS(I,M,EPS,EPSTOP,ENN1,ELONN1,EON,EONTOP)

   Input argument list:
      I        - INTEGER SPECTRAL DOMAIN SHAPE
                 (0 FOR TRIANGULAR, 1 FOR RHOMBOIDAL)
      M        - INTEGER SPECTRAL TRUNCATION
 
   Output argument list:
      EPS      - REAL ((M+1)*((I+1)*M+2)/2) SQRT((N**2-L**2)/(4*N**2-1))
      EPSTOP   - REAL (M+1) SQRT((N**2-L**2)/(4*N**2-1)) OVER TOP
      ENN1     - REAL ((M+1)*((I+1)*M+2)/2) N*(N+1)/A**2
      ELONN1   - REAL ((M+1)*((I+1)*M+2)/2) L/(N*(N+1))*A
      EON      - REAL ((M+1)*((I+1)*M+2)/2) EPSILON/N*A
      EONTOP   - REAL (M+1) EPSILON/N*A OVER TOP

 


  
SPLIB.tar
  
Library contains routines to be
   be used for a variety of spectral transform functions. (Fortran90) 


   Date posted: 2/23/2007