subroutine getlalo(rlats,rlons,wgts,jcap,nlon,nlath,
     *    del2out,trigs,ifax,pln,qln,rln)
c$$$  subprogram documentation block
c                .      .    .                                       .
c subprogram:    getlalo    return gaussian lats and lons and trancons
c   prgmmr: parrish          org: w/nmc22    date: 90-09-21
c
c abstract: return gaussian lats, lons, and transform stuff.
c
c program history log:
c   90-09-21  parrish
c
c   input argument list:
c     jcap     - triangular truncation
c     nlon     - number of longitudes
c     nlath    - number of gaussian lats in one hemisphere
c
c   output argument list:
c     rlats,rlons    - grid latitudes, longitudes (radians)
c     wgts     - gaussian integration weights
c     pln,qln,rln - spherical harmonics.
c     del2     - n*(n+1)/(a**2)
c     trigs,ifax - used by fft
c
c attributes:
c   language: cft77
c   machine:  cray ymp
c
c$$$
c
C-CRA          dimension rlats(2*nlath),wgts(2*nlath),rlons(nlon)
C-CRA          dimension ap(0:jcap,0:jcap),bp(0:jcap,0:jcap)
C-CRA          dimension aqr(0:jcap,0:jcap),bqr(0:jcap,0:jcap)
C-CRA          dimension gr(0:jcap,0:jcap),del2(0:jcap,0:jcap)
C-CRA          dimension slat(nlath),clat(nlath),trigs(nlon*2),ifax(10)
C-CRA          dimension pe0(nlath,0:jcap),qe0(nlath,0:jcap),ro0(nlath,0:jcap)
C-CRA          dimension pln((jcap+1)*(jcap+2),nlath)
C-CRA          dimension qln((jcap+1)*(jcap+2),nlath)
C-CRA          dimension rln((jcap+1)*(jcap+2),nlath)
C-CRA          dimension del2out((jcap+1)*(jcap+2))
C-CRA          dimension w1(nlath-1),w2(nlath-1),w3(nlath-1),w4(nlath-1)
 
          dimension rlats(2*48),wgts(2*48),rlons(192)
          dimension ap(0:62,0:62),bp(0:62,0:62)
          dimension aqr(0:62,0:62),bqr(0:62,0:62)
          dimension gr(0:62,0:62),del2(0:62,0:62)
          dimension slat(48),clat(48),trigs(192*2),ifax(10)
          dimension pe0(48,0:62),qe0(48,0:62)
          dimension ro0(48,0:62)
          dimension pln((62+1)*(62+2),48)
          dimension qln((62+1)*(62+2),48)
          dimension rln((62+1)*(62+2),48)
          dimension del2out((62+1)*(62+2))
          dimension w1(48-1),w2(48-1),w3(48-1),w4(48-1)
c--------
c--------
c-------- get stuff for fft*s
c-------
C-CRA             call fftfax(nlon,ifax,trigs)
             call fax   (ifax,nlon,3)
             call fftrig(trigs,nlon,3)

c--------
c-------- compute recursion constants and initialize legendre functions
c--------
      call m1rcons(ap,bp,aqr,bqr,gr,del2,jcap)
c--------
c-------- get latitudes, longitudes, and integration weights.
c-------
      pih=2.*atan(1.)
      call m1glat(nlath-1,w1,w2,w3,w4)
      rlats(1)=-pih
      wgts(1)=0.
ccdir$ ivdep
      do 100 i=1,nlath-1
        rlats(i+1)=w1(i)-pih
        wgts(i+1)=w2(i)
100   continue
ccdir$ ivdep
      do 200 i=1,nlath
        rlats(nlath*2+1-i)=-rlats(i)
        wgts(nlath*2+1-i)=wgts(i)
        slat(i)=sin(rlats(i))
        clat(i)=cos(rlats(i))
200   continue
      dlon=8.*atan(1.)/nlon
      do 300 i=1,nlon
        rlons(i)=(i-1.)*dlon
300   continue
c--------
c-------- next get initial p,q,r.
c--------
      call m1ipqr(pe0,qe0,ro0,slat,clat,nlath,jcap)
c-------- finally get pln,qln,rln
      call getpln(pln,qln,rln,jcap,nlath,ap,bp,slat,pe0,
     *          qe0,ro0,aqr,bqr,gr,clat,del2,del2out)
      return
      end