User: rhaas Date: 2014/01/13 12:11 PM
Modified: /branches/ET_2013_11/ /branches/ET_2013_11/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.
Directory Changes:
Directory: /branches/ET_2013_11/ ================================
Property changes on: branches/ET_2013_11 ___________________________________________________________________ Directory: /svn:mergeinfo/ ==========================
+ /trunk:122
File Changes:
Directory: /branches/ET_2013_11/src/ ====================================
File [modified]: teukwaves.F77 Delta lines: +2 -1 =================================================================== --- branches/ET_2013_11/src/teukwaves.F77 2014-01-13 17:16:03 UTC (rev 122) +++ branches/ET_2013_11/src/teukwaves.F77 2014-01-13 18:11:37 UTC (rev 123) @@ -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