SPLAPLAC

The SPLAPLAC routine computes the Laplacian or the inverse Laplacian
             of a scalar field in spectral space.
             Subprogram speps should be called already.
             The Laplacian of q(l,n) is simply -n*(n+1)/a**2*q(l,n)

USAGE:    CALL SPLAPLAC(I,M,ENN1,Q,QD2,IDIR)

   Input argument list:
      I        - INTEGER SPECTRAL DOMAIN SHAPE
                 (0 FOR TRIANGULAR, 1 FOR RHOMBOIDAL)
      M        - INTEGER SPECTRAL TRUNCATION
      ENN1     - REAL ((M+1)*((I+1)*M+2)/2) N*(N+1)/A**2
      Q        - IF IDIR > 0, REAL ((M+1)*((I+1)*M+2)) SCALAR FIELD
      QD2      - IF IDIR < 0, REAL ((M+1)*((I+1)*M+2)) LAPLACIAN
      IDIR     - INTEGER FLAG
                 IDIR > 0 TO TAKE LAPLACIAN
                 IDIR < 0 TO TAKE INVERSE LAPLACIAN
 
   Output argument list:
      Q        - IF IDIR < 0, REAL ((M+1)*((I+1)*M+2)) SCALAR FIELD
                 (Q(0,0) IS NOT COMPUTED)
      QD2      - IF IDIR > 0, REAL ((M+1)*((I+1)*M+2)) LAPLACIAN
 
 


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

   Date posted: 2/23/2007