[Commits] [svn:einsteintoolkit] NullSHRExtract/branches/tapir/src/ (Rev. 21)

[bela at caltech.edu]

User: szilagyi
Date: 2013/11/07 03:56 PM

Modified:
 /branches/tapir/src/
  NullSHRE_modRadius.F90

Log:
 work on debugging code

File Changes:

Directory: /branches/tapir/src/
===============================

File [modified]: NullSHRE_modRadius.F90
Delta lines: +26 -1
===================================================================
--- branches/tapir/src/NullSHRE_modRadius.F90	2013-11-07 17:52:13 UTC (rev 20)
+++ branches/tapir/src/NullSHRE_modRadius.F90	2013-11-07 21:56:12 UTC (rev 21)
@@ -25,10 +25,35 @@
        end do
     end do
 
-    r0%d = cr * sqrt(sqrt(detg%d*temp%d))
+    r0%d = cr * (detg%d*temp%d)**0.25d0
 
   end subroutine wt_r_long
 
+  ! This is debugging code.  given the metric and derivatives,
+  ! will compute drBondi/d(cartesian R)
+  subroutine wt_dr(g3, detg3, g3inv, dg3_dR, j0inv, r0, dr_dR)
+    use NullSHRE_modGFdef
+    implicit none
+
+    type (gf2d), dimension (3,3),intent (in)    :: g3, g3inv ! g_ij, g^ij
+    type (gf2d),                 intent (in)    :: detg3     ! det(g_ij)
+    type (gf2d), dimension (3,3),intent (in)    :: dg3_dR    ! g_ij{,r}(cart)
+    ! x,y,z derivatives of (r, q, p) coordinates :
+    type (gf2d), dimension (3,3),intent (in)    :: j0inv
+    type (gf2d),                 intent (in)    :: r0
+    type (gf2d),                 intent (inout) :: dr_dR
+
+
+    ! radial deriv of j0inv:
+    !   dj0inv(a,:)/dr = 0 for a=1
+    !   dj0inv(a,:)/dr = -j0inv(a,:)/cr for a=2,3
+
+    ! rB = R ( detg g^rr)^(1/4)
+
+
+
+  end subroutine wt_dr
+
   subroutine wt_r(nn1, nn2, pp, eta0, r0)
 ! the bondi radius eq. (53) 
     use NullSHRE_modGFdef