1***********************************************************************
2*     INITIALIZE BASE LEGENDRE FUNCTIONS
3************************************************************************
4*     Copyright (C) 2000-2004 GFD Dennou Club. All rights reserved.
5*-----------------------------------------------------------------------
6      SUBROUTINE SHINIP(MM,JM,YS,YC,PY)
7
8      REAL YS(0:JM),YC(0:JM)
9      REAL PY(2,0:JM,0:MM)
10
11*     M=0: P^0_0 and P^0_1,
12*     M>0: P^M_M/\cos\phi and P^M_{M+1}/\cos\phi
13
14      DO 10 J=0,JM
15        PY(1,J,0)=1
16        PY(2,J,0)=SQRT(3.0)*YS(J)
17   10 CONTINUE
18
19      M=1
20      A=SQRT(1.0*(2*M+1)/(2*M))
21      B=SQRT(2*M+3.0)
22      DO 20 J=0,JM
23        PY(1,J,M)=A*PY(1,J,M-1)
24        PY(2,J,M)=B*YS(J)*PY(1,J,M)
25   20 CONTINUE
26
27      DO 40 M=2,MM
28        A=SQRT(1.0*(2*M+1)/(2*M))
29        B=SQRT(2*M+3.0)
30        DO 30 J=0,JM
31          PY(1,J,M)=A*YC(J)*PY(1,J,M-1)
32          PY(2,J,M)=B*YS(J)*PY(1,J,M)
33   30   CONTINUE
34   40 CONTINUE
35
36      END
37