speps.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	speps (I, M, EPS, EPSTOP, ENN1, ELONN1, EON, EONTOP)
	Computes constant fields indexed in the spectral domain in "IBM ORDER" (Zonal wavenumber is the slower index). More...

Detailed Description

Compute utility spectral fields

Author

IREDELL

Date

92-10-31

Definition in file speps.f.

Function/Subroutine Documentation

speps()

subroutine speps	(		I,
			M,
		real, dimension((m+1)*((i+1)*m+2)/2)	EPS,
		real, dimension(m+1)	EPSTOP,
		real, dimension((m+1)*((i+1)*m+2)/2)	ENN1,
		real, dimension((m+1)*((i+1)*m+2)/2)	ELONN1,
		real, dimension((m+1)*((i+1)*m+2)/2)	EON,
		real, dimension(m+1)	EONTOP
	)

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

Parameters

I	- INTEGER SPECTRAL DOMAIN SHAPE (0 FOR TRIANGULAR, 1 FOR RHOMBOIDAL)
M	- INTEGER SPECTRAL TRUNCATION
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

Definition at line 26 of file speps.f.