diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-29 16:35:12 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-29 16:35:12 +0000 |
commit | f1d2591498f08a6912fa6cb81466683552173874 (patch) | |
tree | 91561fd830d90e11d22a93c44396d641a98a954d /src/GRHydro_PPMMReconstruct_drv.F90 | |
parent | c86286abf77e6ccc59f6a09aabe1080e4f79c0b4 (diff) |
GRHydro: initialy set y_e_plus/minus to cell center value in reconstruct
This ensures that all EOS calls have some sensible y_e to work with.
From: Christian Ott <cott@bethe.tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@531 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_PPMMReconstruct_drv.F90')
-rw-r--r-- | src/GRHydro_PPMMReconstruct_drv.F90 | 10 |
1 files changed, 8 insertions, 2 deletions
diff --git a/src/GRHydro_PPMMReconstruct_drv.F90 b/src/GRHydro_PPMMReconstruct_drv.F90 index 096d056..370d764 100644 --- a/src/GRHydro_PPMMReconstruct_drv.F90 +++ b/src/GRHydro_PPMMReconstruct_drv.F90 @@ -175,10 +175,16 @@ subroutine GRHydro_PPMMReconstruct_drv(CCTK_ARGUMENTS) endif if (evolve_Y_e .ne. 0) then - Y_e_plus(i,j,k) = 0.0d0 - Y_e_minus(i,j,k) = 0.0d0 + 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 |