Changeset 379 for trunk/fsource/texturesubs

Timestamp:

Sep 22, 2011 10:11:18 AM (14 years ago)

Author:

vondreele

Message:

finish spherical harmonics preferred orientation

File:

• trunk/fsource/texturesubs/plmpsi.for (modified) (3 diffs)

trunk/fsource/texturesubs/plmpsi.for ¶

@@ -1 @@
-      SUBROUTINE PLMPSI(L,M,PSI,P)
+      SUBROUTINE PLMPSI(L,M,PSI,P,DPDPS)
 
 !PURPOSE: Compute P(l,m,psi)
@@ -31 +31 @@
           RS = S
           P = P+APR*COSD(RS*PSI)
+          DPDPS = DPDPS-RS*APR*SIND(RS*PSI)
         END DO
       ELSE                                            
@@ -37 +38 @@
           RS = S
           P = P+APR*SIND(RS*PSI)        
+          DPDPS = DPDPS+RS*APR*COSD(RS*PSI)
         END DO      
       END IF