[Commits] [svn:einsteintoolkit] IDLinearWaves/trunk/src/ (Rev. 122)

[rhaas at tapir.caltech.edu]

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