speps.f File Reference

Compute utility spectral fields. More...

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).

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	spectral domain shape (0 for triangular, 1 for rhomboidal)
M	spectral truncation
EPS	((M+1)*((I+1)*M+2)/2) SQRT((N**2-L**2)/(4*N**2-1))
EPSTOP	(M+1) SQRT((N**2-L**2)/(4*N**2-1)) OVER TOP
ENN1	((M+1)*((I+1)*M+2)/2) N*(N+1)/A**2
ELONN1	((M+1)*((I+1)*M+2)/2) L/(N*(N+1))*A
EON	((M+1)*((I+1)*M+2)/2) EPSILON/N*A
EONTOP	(M+1) EPSILON/N*A OVER TOP

Author

Iredell

Date

92-10-31

Definition at line 25 of file speps.f.

Referenced by spwget().