User: rhaas Date: 2014/01/13 11:16 AM

Modified: /trunk/src/ teukwaves.F77

Log: add m=2 ftp coefficient for teukolsky waves

The theta-phi coefficient was not set at all and contained random data. This commit implements the coefficient as found in Eq. 7 of PRD 26:745 (1982) which is the paper referenced in the docs.

File Changes:

Directory: /trunk/src/ ======================

File [modified]: teukwaves.F77 Delta lines: +2 -1 =================================================================== --- trunk/src/teukwaves.F77 2013-12-28 04:56:11 UTC (rev 121) +++ trunk/src/teukwaves.F77 2014-01-13 17:16:03 UTC (rev 122) @@ -358,7 +358,7 @@ sinp = yp/(xp**2+yp**2)**0.5 cosp = xp/(xp**2+yp**2)**0.5

-c mvalue +c mvalue (Eq. 7 of PRD 26:745 (1982)) if (m.eq.0) then frr = 2-3*sint**2 frt = -3*sint*cost @@ -393,6 +393,7 @@ frp = sint*(1-2*sinp**2) ftt1 = (1+cost**2)*2*sinp*cosp ftt2 = -2*sinp*cosp + ftp = cost*(1-2*cosp**2) fpp1 = -ftt1 fpp2 = cost**2*2*sinp*cosp drt = 8*sint*sinp*cosp