sptranf1.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	sptranf1 (IROMB, MAXWV, IDRT, IMAX, JMAX, JB, JE, EPS, EPSTOP, ENN1, ELONN1, EON, EONTOP, AFFT, CLAT, SLAT, WLAT, PLN, PLNTOP, MP, W, WTOP, G, IDIR)
	THIS SUBPROGRAM PERFORMS AN SINGLE LATITUDE TRANSFORM FOR SUBPROGRAM SPTRANF. More...

Detailed Description

Sptranf spectral transform

Author

IREDELL

Date

96-02-29

Definition in file sptranf1.f.

Function/Subroutine Documentation

sptranf1()

subroutine sptranf1	(		IROMB,
			MAXWV,
			IDRT,
			IMAX,
			JMAX,
			JB,
			JE,
		real, dimension((maxwv+1)*((iromb+1)*maxwv+2)/2)	EPS,
		real, dimension(maxwv+1)	EPSTOP,
		real, dimension((maxwv+1)*((iromb+1)*maxwv+2)/2)	ENN1,
		real, dimension((maxwv+1)*((iromb+1)*maxwv+2)/2)	ELONN1,
		real, dimension((maxwv+1)*((iromb+1)*maxwv+2)/2)	EON,
		real, dimension(maxwv+1)	EONTOP,
		real(8), dimension(50000+4*imax)	AFFT,
		real, dimension(jb:je)	CLAT,
		real, dimension(jb:je)	SLAT,
		real, dimension(jb:je)	WLAT,
		real, dimension((maxwv+1)*((iromb+1)*maxwv+2)/2,jb:je)	PLN,
		real, dimension(maxwv+1,jb:je)	PLNTOP,
			MP,
		real, dimension((maxwv+1)*((iromb+1)*maxwv+2))	W,
		real, dimension(2*(maxwv+1))	WTOP,
		real, dimension(imax,2,jb:je)	G,
			IDIR
	)

THIS SUBPROGRAM PERFORMS AN SINGLE LATITUDE TRANSFORM FOR SUBPROGRAM SPTRANF.

USE THIS SUBPROGRAM OUTSIDE THE SPTRANF FAMILY CONTEXT AT YOUR OWN RISK.

Parameters

	IROMB	- INTEGER SPECTRAL DOMAIN SHAPE (0 FOR TRIANGULAR, 1 FOR RHOMBOIDAL)
	MAXWV	- INTEGER SPECTRAL TRUNCATION
	IDRT	- INTEGER GRID IDENTIFIER (IDRT=4 FOR GAUSSIAN GRID, IDRT=0 FOR EQUALLY-SPACED GRID INCLUDING POLES, IDRT=256 FOR EQUALLY-SPACED GRID EXCLUDING POLES)
	IMAX	- INTEGER EVEN NUMBER OF LONGITUDES
	JMAX	- INTEGER NUMBER OF LATITUDES
	JB	- INTEGER LATITUDE INDEX (FROM POLE) TO BEGIN TRANSFORM
	JE	- INTEGER LATITUDE INDEX (FROM POLE) TO END TRANSFORM
	EPS	- REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
	EPSTOP	- REAL (MAXWV+1)
	ENN1	- REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
	ELONN1	- REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
	EON	- REAL ((MAXWV+1)*((IROMB+1)*MAXWV+2)/2)
	EONTOP	- REAL (MAXWV+1)
	CLAT	- REAL (JB:JE) COSINES OF LATITUDE
	SLAT	- REAL (JB:JE) SINES OF LATITUDE
	WLAT	- REAL (JB:JE) GAUSSIAN WEIGHTS
	AFFT	- REAL(8) (50000+4*IMAX) AUXILIARY ARRAY IF IDIR=0
	PLN	- REAL ((M+1)*((I+1)*M+2)/2,JB:JE) LEGENDRE POLYNOMIALS
	PLNTOP	- REAL (M+1,JB:JE) LEGENDRE POLYNOMIAL OVER TOP
	MP	- INTEGER IDENTIFIER (0 FOR SCALAR, 1 FOR VECTOR)
[out]	W	- REAL (*) WAVE FIELD IF IDIR>0
[out]	WTOP	- REAL (*) WAVE FIELD OVER TOP IF IDIR>0
[out]	G	- REAL (IMAX,2,JB:JE) GRID FIELD IF IDIR<0
	IDIR	- INTEGER TRANSFORM FLAG (IDIR>0 FOR WAVE TO GRID, IDIR<0 FOR GRID TO WAVE)

SUBPROGRAMS CALLED:

• SPSYNTH SYNTHESIZE FOURIER FROM SPECTRAL
• SPANALY ANALYZE SPECTRAL FROM FOURIER
• SPFFTE PERFORM FAST FOURIER TRANSFORM

Definition at line 48 of file sptranf1.f.