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