[Commits] [svn:einsteintoolkit] IDLinearWaves/ (Rev. 123)

[rhaas at tapir.caltech.edu]

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