From 80d71d7dab11773fa2ede82ba14db59920d40c44 Mon Sep 17 00:00:00 2001 From: reisswig Date: Sat, 26 Nov 2011 07:52:10 +0000 Subject: Implemented enhanced PPM scheme by Colella & Sekora 2008, McCorquodale & Colella 2011. Can be actived by setting use_enhanced_ppm = yes. Also fixed a bug in GRHydro_Reconstruct.F90 when using MP. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@302 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_PPMReconstruct_drv.F90 | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) (limited to 'src/GRHydro_PPMReconstruct_drv.F90') diff --git a/src/GRHydro_PPMReconstruct_drv.F90 b/src/GRHydro_PPMReconstruct_drv.F90 index 02ab359..b074bec 100644 --- a/src/GRHydro_PPMReconstruct_drv.F90 +++ b/src/GRHydro_PPMReconstruct_drv.F90 @@ -168,7 +168,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) w_lorentz(:,j,k), & 1, j, k, nx, ny, nz, GRHydro_mppm_eigenvalue_x_left, & GRHydro_mppm_eigenvalue_x_right, & - GRHydro_mppm_xwind) + GRHydro_mppm_xwind, GRHydro_rho_min) do i = 1, nx if (trivial_rp(i,j,k)) then SpaceMask_SetStateBitsF90(space_mask, i, j, k, type_bitsx, trivialx) @@ -198,7 +198,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) call SimplePPM_ye_1d(nx,CCTK_DELTA_SPACE(1),rho(:,j,k), & velx(:,j,k),vely(:,j,k),velz(:,j,k), & Y_e(:,j,k),Y_e_minus(:,j,k),Y_e_plus(:,j,k), & - press(:,j,k)) + press(:,j,k), GRHydro_rho_min) end do end do !$OMP END PARALLEL DO @@ -219,7 +219,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) w_lorentz(j,:,k), & 2, j, k, nx, ny, nz, GRHydro_mppm_eigenvalue_y_left, & GRHydro_mppm_eigenvalue_y_right, & - GRHydro_mppm_xwind) + GRHydro_mppm_xwind, GRHydro_rho_min) do i = 1, ny if (trivial_rp(j,i,k)) then SpaceMask_SetStateBitsF90(space_mask, j, i, k, type_bitsy, trivialy) @@ -249,7 +249,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) call SimplePPM_ye_1d(ny,CCTK_DELTA_SPACE(2),rho(j,:,k), & velx(j,:,k),vely(j,:,k),velz(j,:,k), & Y_e(j,:,k),Y_e_minus(j,:,k),Y_e_plus(j,:,k), & - press(j,:,k)) + press(j,:,k), GRHydro_rho_min) end do end do !$OMP END PARALLEL DO @@ -270,7 +270,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) w_lorentz(j,k,:), & 3, j, k, nx, ny, nz, GRHydro_mppm_eigenvalue_z_left, & GRHydro_mppm_eigenvalue_z_right, & - GRHydro_mppm_xwind) + GRHydro_mppm_xwind, GRHydro_rho_min) do i = 1, nz if (trivial_rp(j,k,i)) then SpaceMask_SetStateBitsF90(space_mask, j, k, i, type_bitsz, trivialz) @@ -301,7 +301,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) call SimplePPM_ye_1d(nz,CCTK_DELTA_SPACE(3),rho(j,k,:), & velx(j,k,:),vely(j,k,:),velz(j,k,:), & Y_e(j,k,:),Y_e_minus(j,k,:),Y_e_plus(j,k,:), & - press(j,k,:)) + press(j,k,:), GRHydro_rho_min) end do end do !$OMP END PARALLEL DO -- cgit v1.2.3