User: bmundim
Date: 2011/11/01 05:19 PM
Modified:
/trunk/src/
GRHydro_SourceM.F90, GRHydro_TmunuM.F90
Log:
RIT dev: Implement comoving contravariant magnetic 4-vector.
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_SourceM.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_SourceM.F90 2011-10-25 18:45:10 UTC (rev 286)
+++ trunk/src/GRHydro_SourceM.F90 2011-11-01 22:19:55 UTC (rev 287)
@@ -252,7 +252,7 @@
rhohstarW2 = (rho(i,j,k)*(one + eps(i,j,k)) + press(i,j,k)+ b2)*&
w_lorentz(i,j,k)**2
- pstar = press(i,j,k)+b2/2.d0
+ pstar = press(i,j,k)+0.5d0*b2
!!$ For a change, these are T^{ij}
File [modified]: GRHydro_TmunuM.F90
Delta lines: +13 -1
===================================================================
--- trunk/src/GRHydro_TmunuM.F90 2011-10-25 18:45:10 UTC (rev 286)
+++ trunk/src/GRHydro_TmunuM.F90 2011-11-01 22:19:55 UTC (rev 287)
@@ -67,12 +67,16 @@
#define Bvecx(i,j,k) Bvec(i,j,k,1)
#define Bvecy(i,j,k) Bvec(i,j,k,2)
#define Bvecz(i,j,k) Bvec(i,j,k,3)
+#define bcomx(i,j,k) bcom(i,j,k,1)
+#define bcomy(i,j,k) bcom(i,j,k,2)
+#define bcomz(i,j,k) bcom(i,j,k,3)
subroutine GRHydro_TmunuM(CCTK_ARGUMENTS)
implicit none
DECLARE_CCTK_ARGUMENTS
+ DECLARE_CCTK_PARAMETERS
CCTK_REAL :: velxlow, velylow, velzlow
CCTK_REAL :: betaxlow, betaylow, betazlow, beta2
@@ -119,7 +123,7 @@
!!$ square of the lorentz factor.
rhohstarw2 = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0D0 + eps(i,j,k)) + press(i,j,k) + b2)
- pstar = press(i,j,k)+b2/2.d0
+ pstar = press(i,j,k)+0.5d0*b2
!!$ Calculate lower components of 4-velocity (without the Lorent factor).
!!$ uxlow = velxlow
@@ -142,6 +146,14 @@
eTxz(i,j,k) = eTxz(i,j,k) + rhohstarw2*velxlow*velzlow + pstar*gxz(i,j,k) - bxlow*bzlow
eTyz(i,j,k) = eTyz(i,j,k) + rhohstarw2*velylow*velzlow + pstar*gyz(i,j,k) - bylow*bzlow
+ if(calculate_bcom .ne. 0) then
+ bcom_sq(i,j,k) = b2
+ bcom0(i,j,k) = w_lorentz(i,j,k)*bdotv/alp(i,j,k)
+ bcomx(i,j,k) = Bvecx(i,j,k)/w_lorentz(i,j,k) + bcom0(i,j,k)*(alp(i,j,k)*velx(i,j,k)-betax(i,j,k))
+ bcomy(i,j,k) = Bvecy(i,j,k)/w_lorentz(i,j,k) + bcom0(i,j,k)*(alp(i,j,k)*vely(i,j,k)-betay(i,j,k))
+ bcomz(i,j,k) = Bvecz(i,j,k)/w_lorentz(i,j,k) + bcom0(i,j,k)*(alp(i,j,k)*velz(i,j,k)-betaz(i,j,k))
+ endif
+
end do
end do
end do
