w3fa11.f

Go to the documentation of this file.

C> @file
C> @brief Computes coefficients for use in w3fa12.
C> @author Joe Sela @date 1980-10-28
 
C> Subroutine computes double precision  coefficients
C> used in generating  legendre polynomials in subr. w3fa12.
C> on a cray double precision is changed to real, dsqrt to sqrt.
C>
C> Program history log:
C> - Joe Sela 1980-10-28
C> - Ralph Jones 1984-06-01 0change to ibm vs fortran.
C> - Ralph Jones 1993-04-12 0changes for cray, double precision to real.
C>
C> @param[out] EPS Real coefficients used in computing legendre polynomials.
C> dimension of eps is (jcap+2)*(jcap+1)
C> @param[in] JCAP Zonal wave number thirty, etc.
C>
C> @author Joe Sela @date 1980-10-28
 
       SUBROUTINE w3fa11 (EPS,JCAP)
C
       REAL  EPS(*)
       REAL  A
C
       SAVE
C
         jcap1 = jcap + 1
         jcap2 = jcap + 2
C
         DO 100 ll = 1,jcap1
           l = ll - 1
           jle = (ll-1) * jcap2
C
         DO 100 inde = 2,jcap2
           n = l + inde - 1
           a=(n*n-l*l)/(4.0*n*n-1.0)
           eps(jle+inde) = sqrt(a)
 100     CONTINUE
C
         DO 200 ll = 1,jcap1
           jle = (ll-1) * jcap2
           eps(jle+1) = 0.0
 200     CONTINUE
C
         RETURN
       END