diff options
author | reisswig <reisswig@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-26 07:52:10 +0000 |
---|---|---|
committer | reisswig <reisswig@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-26 07:52:10 +0000 |
commit | 80d71d7dab11773fa2ede82ba14db59920d40c44 (patch) | |
tree | 777f60434f5e62d6f48d5459c349900548ba747b /src/GRHydro_PPMReconstruct_drv.F90 | |
parent | 1849ffd65e2fe0ea4210cd92102d07061fce28d4 (diff) |
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
Diffstat (limited to 'src/GRHydro_PPMReconstruct_drv.F90')
-rw-r--r-- | src/GRHydro_PPMReconstruct_drv.F90 | 12 |
1 files changed, 6 insertions, 6 deletions
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 |