ZNME2PXX.inc

1       SUBROUTINE ZNME2PXX
2      I(MEND1,MNWAV,JEQUH,
3      O PNM  ,HNM  ,
4      W PPNM ,HHNM)
5 C
6 C     THIS SUBROUTINE COMPUTES LEGENDRE FUNCTION VALUES ON HEMISPHERIC
7 C     EQUIDISTANT GRIDS
8 C
9 C    NOTE1: THE ORDER OF WAVE COMPONENTS IS COLUMNWISE
10 C
11       IMPLICIT REAL*8(A-H,O-Z)
12 Crizvi      REAL*4 PNM (MNWAV,JEQUH),HNM (MNWAV,JEQUH)
13       REAL PNM (MNWAV,JEQUH),HNM (MNWAV,JEQUH)
14       REAL*8 PPNM(MNWAV      ),HHNM(MNWAV      )
15 C
16       NLATH=JEQUH-1
17       FLIV =1.E0/DFLOAT(NLATH)
18       FLAT =2.E0*ATAN(1.E0)
19       DLAT =FLAT*FLIV
20 C
21       DO 100 J=1,JEQUH
22 cch      SL=DSIN(FLAT)
23       SL=SIN(FLAT)
24       CALL LGNDR1(SL,MEND1,PPNM,HHNM)
25       DO 120 L=1,MNWAV
26 cch      PNM (L,J)=SNGL(PPNM(L))
27 cch      HNM (L,J)=SNGL(HHNM(L))
28       PNM (L,J)=PPNM(L)
29       HNM (L,J)=HHNM(L)
30   120 CONTINUE
31       FLAT=FLAT-DLAT
32   100 CONTINUE
33 C
34 C
35       RETURN
36       END SUBROUTINE ZNME2PXX
37