source: trunk/fsource/texturesubs/plmpsi.for @ 305

Last change on this file since 305 was 305, checked in by vondreele, 12 years ago

remove pypowder & powsubs - no longer needed
add texturesubs & pytexture

File size: 1.0 KB
Line 
1      SUBROUTINE PLMPSI(L,M,PSI,P)
2
3!PURPOSE: Compute P(l,m,psi)
4
5      INCLUDE       '../INCLDS/COPYRIGT.FOR' 
6
7!CALLING ARGUMENTS:
8
9      INTEGER*4     L,M                 !Order & index
10      REAL*4        PSI                 !Angle (in deg)
11      REAL*4        P                   !Value returned
12
13!INCLUDE STATEMENTS:
14
15!LOCAL VARIABLES:
16                                                 
17      INTEGER*4     S                   
18      REAL*4        APR,RS             
19                                                 
20!FUNCTION DEFINITIONS:
21
22!DATA STATEMENTS:
23
24!CODE:                 
25         
26      P = 0.0
27      DPDPS = 0.0
28      IF ( MOD(ABS(M),2).EQ.0 ) THEN
29        DO S=0,L,2
30          CALL APLMS(L,M,S,APR)
31          RS = S
32          P = P+APR*COSD(RS*PSI)
33        END DO
34      ELSE                                           
35        DO S=2,L,2       
36          CALL APLMS(L,M,S,APR)     
37          RS = S
38          P = P+APR*SIND(RS*PSI)       
39        END DO     
40      END IF
41      RETURN
42      END
43
Note: See TracBrowser for help on using the repository browser.