User: rhaas Date: 2013/05/11 01:29 AM

Modified: /trunk/src/ GRHydro_Prim2ConM.F90

Log: GRHydro: Fix a problem with ghost zones in Prim2ConCells calls in MHD case.

From: Christian Reisswig reisswig@tapir.caltech.edu

File Changes:

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

File [modified]: GRHydro_Prim2ConM.F90 Delta lines: +9 -9 =================================================================== --- trunk/src/GRHydro_Prim2ConM.F90 2013-05-11 06:29:34 UTC (rev 511) +++ trunk/src/GRHydro_Prim2ConM.F90 2013-05-11 06:29:37 UTC (rev 512) @@ -464,9 +464,9 @@ if(evolve_temper.ne.1) then

!$OMP PARALLEL DO PRIVATE(k,j,i,det), REDUCTION(MAX:maxtau0) - do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1 - do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1 - do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1 + do k = 1,cctk_lsh(3) + do j = 1,cctk_lsh(2) + do i = 1,cctk_lsh(1)

det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k), \ gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) @@ -505,9 +505,9 @@

else !$OMP PARALLEL DO PRIVATE(k,j,i,det) - do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1 - do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1 - do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1 + do k = 1,cctk_lsh(3) + do j = 1,cctk_lsh(2) + do i = 1,cctk_lsh(1)

det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k), \ gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) @@ -747,9 +747,9 @@ CCTK_REAL :: det

!$OMP PARALLEL DO PRIVATE(k,j,i,det) - do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1 - do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1 - do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1 + do k = 1,cctk_lsh(3) + do j = 1,cctk_lsh(2) + do i = 1,cctk_lsh(1)

det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\ gyy(i,j,k),gyz(i,j,k),gzz(i,j,k))