diff options
Diffstat (limited to 'src/GRHydro_PPMMReconstruct_drv.F90')
| -rw-r--r-- | src/GRHydro_PPMMReconstruct_drv.F90 | 56 |
1 files changed, 6 insertions, 50 deletions
diff --git a/src/GRHydro_PPMMReconstruct_drv.F90 b/src/GRHydro_PPMMReconstruct_drv.F90 index 370d764..f9ceaf7 100644 --- a/src/GRHydro_PPMMReconstruct_drv.F90 +++ b/src/GRHydro_PPMMReconstruct_drv.F90 @@ -135,60 +135,16 @@ subroutine GRHydro_PPMMReconstruct_drv(CCTK_ARGUMENTS) ny = cctk_lsh(2) nz = cctk_lsh(3) -!!$ Initialize variables that store reconstructed quantities: - - !$OMP PARALLEL DO PRIVATE(i,j,k) - do k=1,cctk_lsh(3) - do j=1,cctk_lsh(2) - do i=1,cctk_lsh(1) + !initialize trivial_rp to false + !$OMP PARALLEL DO + do k=1,nz + do j=1,ny + do i=1,nx trivial_rp(i,j,k) = .false. - rhoplus(i,j,k) = 0.0d0 - rhominus(i,j,k)= 0.0d0 - epsplus(i,j,k) = 0.0d0 - epsminus(i,j,k) = 0.0d0 - velxplus(i,j,k) = 0.0d0 - velxminus(i,j,k) = 0.0d0 - velyplus(i,j,k) = 0.0d0 - velyminus(i,j,k) = 0.0d0 - velzplus(i,j,k) = 0.0d0 - velzminus(i,j,k) = 0.0d0 - - Bvecxplus(i,j,k) = 0.0d0 - Bvecxminus(i,j,k) = 0.0d0 - Bvecyplus(i,j,k) = 0.0d0 - Bvecyminus(i,j,k) = 0.0d0 - Bveczplus(i,j,k) = 0.0d0 - Bveczminus(i,j,k) = 0.0d0 - if(clean_divergence.ne.0) then - psidcplus(i,j,k) = 0.0d0 - psidcminus(i,j,k) = 0.0d0 - endif - - if (evolve_entropy .ne. 0) then - entropyplus(i,j,k) = 0.0d0 - entropyminus(i,j,k) = 0.0d0 - endif - - if (evolve_tracer .ne. 0) then - tracerplus(i,j,k,:) = 0.0d0 - tracerminus(i,j,k,:) = 0.0d0 - endif - - if (evolve_Y_e .ne. 0) then - Y_e_plus(i,j,k) = Y_e(i,j,k) - Y_e_minus(i,j,k) = Y_e(i,j,k) - endif - - if(evolve_temper .ne. 0) then - tempplus(i,j,k) = temperature(i,j,k) - tempminus(i,j,k) = temperature(i,j,k) - endif - - enddo enddo enddo - !$OMP END PARALLEL DO + !$OMP END PARALLEL DO !!$ PPM starts: if (flux_direction == 1) then |