#include <misc.h>#include <params.h>module comspe!----------------------------------------------------------------------- ! ! Purpose: Spectral space arrays! ! Method: ! ! Author: CCM Core Group! $Author: rosinski $! $Id: comspe.F90,v 1.3 2001/09/16 22:13:32 rosinski Exp $! !-----------------------------------------------------------------------  use precision  use infnan  use pmgrid, only: plev, plat  use pspect  implicit none  real(r8) :: vz(psp,plev)   = inf     ! Vorticity spectral coefficients  real(r8) :: d(psp,plev)    = inf     ! Divergence spectral coefficients  real(r8) :: t(psp,plev)    = inf     ! Temperature spectral coefficients  real(r8) :: q(psp,plev)    = inf     ! Moisture     spectral coefficients  real(r8) :: alps(psp)      = inf     ! Log-pressure spectral coefficients  real(r8) :: hs(psp,plev)   = inf     ! hydrostatic matrix for "real" atmosphere  real(r8) :: hsnm(psp,plev) = inf   ! vertical normal modes of "hs"  real(r8) :: dsnm(psp,plev) = inf   ! vertical normal modes of "ds"  real(r8) :: dnm(psp,plev)  = inf   ! vertical normal modes of "d"  real(r8) :: vznm(psp,plev) = inf   ! vertical normal modes of "vz"  real(r8) :: lnpstar(psp)   = inf   ! ln (Ps*) (SLD term; Ritchie & Tanguay, 1995)  real(r8) :: a0nm(psp)      = inf   ! wave # coefs (use in vert normal mode space)  real(r8) :: bmnm(psp)      = inf   ! wave # coefs (use in vert normal mode space)  real(r8) :: bpnm(psp)      = inf   ! wave # coefs (use in vert normal mode space)  real(r8) :: atri(psp)      = inf   ! wave # coefs (use in vert normal mode space)  real(r8) :: btri(psp)      = inf   ! wave # coefs (use in vert normal mode space)  real(r8) :: ctri(psp)      = inf   ! wave # coefs (use in vert normal mode space)  integer :: ncutoff         = bigint   ! Break-even point for vector lengths in GRCALC  integer :: nalp(pmax)      = bigint   ! Pointer into polynomial arrays#if ( defined SPMD )  integer :: begm(0:plat-1)  = bigint  integer :: endm(0:plat-1)  = bigint#else  integer :: begm(0:0)       = 1  integer :: endm(0:0)       = pmmax#endif  integer :: ncoefi(pmaxp)   = bigint   ! Pointer to start of coefficient diagonals  integer :: nm(pmax)        = bigint   ! Number of coeffs stored on a given diagonal  integer :: nco2(pmax)      = bigint   ! Complex form of ncoefi  integer :: nstart(pmmax)   = bigint   ! Starting indices for spectral arrays (real)  integer :: nlen(pmmax)     = bigint   ! Length vectors for spectral arrays  real(r8) :: alp(pspt,plat/2)  = inf   ! Legendre polynomials  real(r8) :: dalp(pspt,plat/2) = inf   ! Legendre polynomial derivativesend module comspe