NCEPLIBS-ip 5.2.0
|
Sptranf spectral transform. More...
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(). | |
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.
IROMB | spectral domain shape (0 for triangular, 1 for rhomboidal) | |
MAXWV | spectral truncation | |
IDRT | grid identifier
| |
IMAX | even number of longitudes | |
JMAX | number of latitudes | |
JB | latitude index (from pole) to begin transform | |
JE | latitude index (from pole) to end transform | |
EPS | ||
EPSTOP | ||
ENN1 | ||
ELONN1 | ||
EON | ||
EONTOP | ||
CLAT | cosines of latitude | |
SLAT | sines of latitude | |
WLAT | Gaussian weights | |
AFFT | auxiliary array if IDIR=0 | |
PLN | Legendre polynomials | |
PLNTOP | Legendre polynomial over top | |
MP | identifier (0 for scalar, 1 for vector) | |
[out] | W | wave field if IDIR>0 |
[out] | WTOP | wave field over top if IDIR>0 |
[out] | G | grid field if IDIR<0 |
IDIR | transform flag (IDIR>0 for wave to grid, IDIR<0 for grid to wave) |
Definition at line 40 of file sptranf1.f.
References spanaly(), spffte(), and spsynth().
Referenced by sptranf(), and sptranfv().