diff options
Diffstat (limited to 'src/GRHydro_TVDReconstruct_drv.F90')
| -rw-r--r-- | src/GRHydro_TVDReconstruct_drv.F90 | 12 |
1 files changed, 12 insertions, 0 deletions
diff --git a/src/GRHydro_TVDReconstruct_drv.F90 b/src/GRHydro_TVDReconstruct_drv.F90 index a51da93..8ab09ba 100644 --- a/src/GRHydro_TVDReconstruct_drv.F90 +++ b/src/GRHydro_TVDReconstruct_drv.F90 @@ -149,6 +149,14 @@ subroutine GRHydro_TVDReconstruct_drv(CCTK_ARGUMENTS) Y_e_plus(i,j,k) = 0.0d0 Y_e_minus(i,j,k) = 0.0d0 endif + + if(evolve_temper .ne. 0) then + ! set this to cell center value to have + ! good initial guesses at interfaces + ! in case we don't reconstruct temp + tempplus(i,j,k) = temperature(i,j,k) + tempminus(i,j,k) = temperature(i,j,k) + endif enddo enddo @@ -196,6 +204,10 @@ subroutine GRHydro_TVDReconstruct_drv(CCTK_ARGUMENTS) call tvdreconstruct(nx, ny, nz, xoffset, yoffset, zoffset, & entropy, entropyplus, entropyminus, trivial_rp, hydro_excision_mask) endif + if (evolve_temper .ne. 0 .and. reconstruct_temper .ne. 0) then + call tvdreconstruct(nx, ny, nz, xoffset, yoffset, zoffset, & + temperature, tempplus, tempminus, trivial_rp, hydro_excision_mask) + endif else if (CCTK_EQUALS(recon_vars,"conservative")) then call tvdreconstruct(nx, ny, nz, xoffset, yoffset, zoffset, & |