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