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

Last change on this file since 379 was 379, checked in by vondreele, 11 years ago

finish spherical harmonics preferred orientation

File size: 1.1 KB
Line 
1      SUBROUTINE PLMPSI(L,M,PSI,P,DPDPS)
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          DPDPS = DPDPS-RS*APR*SIND(RS*PSI)
34        END DO
35      ELSE                                           
36        DO S=2,L,2       
37          CALL APLMS(L,M,S,APR)     
38          RS = S
39          P = P+APR*SIND(RS*PSI)       
40          DPDPS = DPDPS+RS*APR*COSD(RS*PSI)
41        END DO     
42      END IF
43      RETURN
44      END
45
Note: See TracBrowser for help on using the repository browser.