diff options
Diffstat (limited to 'src/GRHydro_ReconstructM.F90')
| -rw-r--r-- | src/GRHydro_ReconstructM.F90 | 26 |
1 files changed, 13 insertions, 13 deletions
diff --git a/src/GRHydro_ReconstructM.F90 b/src/GRHydro_ReconstructM.F90 index 654d9a7..a48d67f 100644 --- a/src/GRHydro_ReconstructM.F90 +++ b/src/GRHydro_ReconstructM.F90 @@ -205,7 +205,7 @@ subroutine ReconstructionM(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 @@ -234,7 +234,7 @@ subroutine ReconstructionM(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,0,nx,CCTK_DELTA_SPACE(1),& @@ -262,7 +262,7 @@ subroutine ReconstructionM(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), & @@ -274,7 +274,7 @@ subroutine ReconstructionM(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), & @@ -288,7 +288,7 @@ subroutine ReconstructionM(CCTK_ARGUMENTS) 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,0,ny,CCTK_DELTA_SPACE(2),& @@ -316,7 +316,7 @@ subroutine ReconstructionM(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), & @@ -328,7 +328,7 @@ subroutine ReconstructionM(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), & @@ -341,7 +341,7 @@ subroutine ReconstructionM(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,0,nz,CCTK_DELTA_SPACE(3),& @@ -369,7 +369,7 @@ subroutine ReconstructionM(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 @@ -382,7 +382,7 @@ subroutine ReconstructionM(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,:), & @@ -417,7 +417,7 @@ subroutine ReconstructionM(CCTK_ARGUMENTS) 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 @@ -486,7 +486,7 @@ subroutine ReconstructionM(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 @@ -553,7 +553,7 @@ subroutine ReconstructionM(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 |