User: bmundim Date: 2011/09/26 10:50 AM

Modified: /trunk/src/ GRHydro_SourceM.F90

Log: RIT dev: fix source term in divergence cleaning. Thanks Tanja for catching this up too.

File Changes:

Directory: /trunk/src/ ======================

File [modified]: GRHydro_SourceM.F90 Delta lines: +10 -13 =================================================================== --- trunk/src/GRHydro_SourceM.F90 2011-09-25 22:59:03 UTC (rev 279) +++ trunk/src/GRHydro_SourceM.F90 2011-09-26 15:50:04 UTC (rev 280) @@ -457,27 +457,24 @@

bvcx_source = -1.d0 * ( Bconsx(i,j,k)*dx_betax + & Bconsy(i,j,k)*dy_betax + Bconsz(i,j,k)*dz_betax ) + & - psidc(i,j,k)*sqrtdet*( uxx*dx_alp+uxy*dy_alp+uxz*dz_alp ) + & - psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxx*dx_det_bydet + & - uxy*dy_det_bydet + uxz*dz_det_bydet ) - & - psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxx*gdg_x + uxy*gdg_y + & - uxz*gdg_z ) + psidc(i,j,k)*sqrtdet*(( uxx*dx_alp+uxy*dy_alp+uxz*dz_alp ) + & + alp(i,j,k)*(0.5d0*( uxx*dx_det_bydet + & + uxy*dy_det_bydet + uxz*dz_det_bydet) - & + ( uxx*gdg_x + uxy*gdg_y + uxz*gdg_z )))

bvcy_source = -1.d0 * ( Bconsx(i,j,k)*dx_betay + & Bconsy(i,j,k)*dy_betay + Bconsz(i,j,k)*dz_betay ) + & - psidc(i,j,k)*sqrtdet*( uxy*dx_alp+uyy*dy_alp+uyz*dz_alp ) + & - psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxy*dx_det_bydet + & + psidc(i,j,k)*sqrtdet*(( uxy*dx_alp+uyy*dy_alp+uyz*dz_alp ) + & + alp(i,j,k)*(0.5d0*( uxy*dx_det_bydet + & uyy*dy_det_bydet + uyz*dz_det_bydet ) - & - psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxy*gdg_x + uyy*gdg_y + & - uyz*gdg_z ) + ( uxy*gdg_x + uyy*gdg_y + uyz*gdg_z )))

bvcz_source = -1.d0 * ( Bconsx(i,j,k)*dx_betaz + & Bconsy(i,j,k)*dy_betaz + Bconsz(i,j,k)*dz_betaz ) + & - psidc(i,j,k)*sqrtdet*( uxz*dx_alp+uyz*dy_alp+uzz*dz_alp ) + & - psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxz*dx_det_bydet + & + psidc(i,j,k)*sqrtdet*(( uxz*dx_alp+uyz*dy_alp+uzz*dz_alp ) + & + alp(i,j,k)*(0.5d0*( uxz*dx_det_bydet + & uyz*dy_det_bydet + uzz*dz_det_bydet ) - & - psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxz*gdg_x + uyz*gdg_y + & - uzz*gdg_z ) + ( uxz*gdg_x + uyz*gdg_y + uzz*gdg_z )))

Bconsrhsx(i,j,k) = bvcx_source Bconsrhsy(i,j,k) = bvcy_source