diff options
Diffstat (limited to 'src/GRHydro_ReconstructPolyM.F90')
-rw-r--r-- | src/GRHydro_ReconstructPolyM.F90 | 26 |
1 files changed, 13 insertions, 13 deletions
diff --git a/src/GRHydro_ReconstructPolyM.F90 b/src/GRHydro_ReconstructPolyM.F90 index 70d2b12..9b4bcd6 100644 --- a/src/GRHydro_ReconstructPolyM.F90 +++ b/src/GRHydro_ReconstructPolyM.F90 @@ -203,7 +203,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) - !$OMP PARALLEL DO PRIVATE(i, j) + !$OMP PARALLEL DO PRIVATE(i, j, k) do k = 1, nz do j = 1, ny do i = 1, nx @@ -232,7 +232,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) else if (CCTK_EQUALS(recon_method,"ppm")) then if (flux_direction == 1) then - !$OMP PARALLEL DO PRIVATE(i, j) + !$OMP PARALLEL DO PRIVATE(i, j, k) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 do j = GRHydro_stencil, ny - GRHydro_stencil + 1 call SimplePPM_1dM(GRHydro_eos_handle,1,nx,CCTK_DELTA_SPACE(1),& @@ -260,7 +260,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) end do !$OMP END PARALLEL DO if(evolve_tracer.ne.0) then - !$OMP PARALLEL DO PRIVATE(j) + !$OMP PARALLEL DO PRIVATE(j, k) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 do j = GRHydro_stencil, ny - GRHydro_stencil + 1 call SimplePPM_tracer_1d(nx,CCTK_DELTA_SPACE(1),rho(:,j,k), & @@ -272,7 +272,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) !$OMP END PARALLEL DO end if if(evolve_Y_e.ne.0) then - !$OMP PARALLEL DO PRIVATE(j) + !$OMP PARALLEL DO PRIVATE(j, k) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 do j = GRHydro_stencil, ny - GRHydro_stencil + 1 call SimplePPM_tracer_1d(nx,CCTK_DELTA_SPACE(1),rho(:,j,k), & @@ -285,7 +285,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) end if else if (flux_direction == 2) then - !$OMP PARALLEL DO PRIVATE(i, j) + !$OMP PARALLEL DO PRIVATE(i, j, k) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_1dM(GRHydro_eos_handle,1,ny,CCTK_DELTA_SPACE(2),& @@ -313,7 +313,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) end do !$OMP END PARALLEL DO if(evolve_tracer.ne.0) then - !$OMP PARALLEL DO PRIVATE(j) + !$OMP PARALLEL DO PRIVATE(j, k) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_tracer_1d(ny,CCTK_DELTA_SPACE(2),rho(j,:,k), & @@ -325,7 +325,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) !$OMP END PARALLEL DO end if if(evolve_Y_e.ne.0) then - !$OMP PARALLEL DO PRIVATE(j) + !$OMP PARALLEL DO PRIVATE(j, k) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_tracer_1d(ny,CCTK_DELTA_SPACE(2),rho(j,:,k), & @@ -338,7 +338,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) end if else if (flux_direction == 3) then - !$OMP PARALLEL DO PRIVATE(i, j) + !$OMP PARALLEL DO PRIVATE(i, j, k) do k = GRHydro_stencil, ny - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_1dM(GRHydro_eos_handle,1,nz,CCTK_DELTA_SPACE(3),& @@ -366,7 +366,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) end do !$OMP END PARALLEL DO if(evolve_tracer.ne.0) then - !$OMP PARALLEL DO PRIVATE(j) + !$OMP PARALLEL DO PRIVATE(j, k) do k = GRHydro_stencil, ny - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 @@ -379,7 +379,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) !$OMP END PARALLEL DO end if if(evolve_Y_e.ne.0) then - !$OMP PARALLEL DO PRIVATE(j) + !$OMP PARALLEL DO PRIVATE(j, k) do k = GRHydro_stencil, ny - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_tracer_1d(nz,CCTK_DELTA_SPACE(3),rho(j,k,:), & @@ -414,7 +414,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) endif if (flux_direction == 1) then - !$OMP PARALLEL DO PRIVATE(i, j) + !$OMP PARALLEL DO PRIVATE(i, j, k) do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + 1 do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + 1 if (CCTK_EQUALS(recon_vars,"primitive")) then @@ -476,7 +476,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) end do !$OMP END PARALLEL DO else if (flux_direction == 2) then - !$OMP PARALLEL DO PRIVATE(i, j) + !$OMP PARALLEL DO PRIVATE(i, j, k) do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + 1 do j = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + 1 if (CCTK_EQUALS(recon_vars,"primitive")) then @@ -537,7 +537,7 @@ subroutine ReconstructionPolytypeM(CCTK_ARGUMENTS) end do !$OMP END PARALLEL DO else if (flux_direction == 3) then - !$OMP PARALLEL DO PRIVATE(i, j) + !$OMP PARALLEL DO PRIVATE(i, j, k) do k = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + 1 do j = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + 1 if (CCTK_EQUALS(recon_vars,"primitive")) then |