User: rhaas Date: 2013/05/13 12:30 AM
Modified: /branches/ET_2013_05/ /branches/ET_2013_05/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
Directory Changes:
Directory: /branches/ET_2013_05/ ================================
Property changes on: branches/ET_2013_05 ___________________________________________________________________ Directory: /svn:mergeinfo/ ==========================
- /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 /trunk:510 + /branches/divclean:231-243 /branches/hot_and_MHD_temp_dev:166-185 /trunk:510,512
File Changes:
Directory: /branches/ET_2013_05/src/ ====================================
File [modified]: GRHydro_Prim2ConM.F90 Delta lines: +9 -9 =================================================================== --- branches/ET_2013_05/src/GRHydro_Prim2ConM.F90 2013-05-13 05:27:57 UTC (rev 513) +++ branches/ET_2013_05/src/GRHydro_Prim2ConM.F90 2013-05-13 05:29:59 UTC (rev 514) @@ -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))