diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:11:04 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:11:04 +0000 |
| commit | fcbaf698714e9594ef35fbba18530a4a1ab164ba (patch) | |
| tree | 9581535efd0b04a3ba9dfcaeabc43a876f214a6a | |
| parent | c4270e1b76a25d22cbb8a20d23235ddd4bb70822 (diff) | |
GRHydro: pass correct temperature o prim2con call on coarse levels
Fix a bug in prim2con_hot that was causing temp to go negative near the
shock in 3D CCSN simulations with WENO reconstruction.
From: abdik <abdik@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@554 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_Prim2Con.F90 | 1 | ||||
| -rw-r--r-- | src/GRHydro_Prim2ConM.F90 | 1 |
2 files changed, 2 insertions, 0 deletions
diff --git a/src/GRHydro_Prim2Con.F90 b/src/GRHydro_Prim2Con.F90 index 80ce321..680ecfc 100644 --- a/src/GRHydro_Prim2Con.F90 +++ b/src/GRHydro_Prim2Con.F90 @@ -378,6 +378,7 @@ subroutine prim2con_hot(handle, keytemp, GRHydro_reflevel, cctk_iteration, ii, j ! away), we use the average temperature between cells and call ! the EOS with keytemp=1 lkeytemp=1 + temp = temp0 call EOS_Omni_press(handle,lkeytemp,GRHydro_eos_rf_prec,n,& drho,deps,temp,ye,dpress,keyerr,anyerr) lkeytemp=0 diff --git a/src/GRHydro_Prim2ConM.F90 b/src/GRHydro_Prim2ConM.F90 index 7c6118f..cb459ee 100644 --- a/src/GRHydro_Prim2ConM.F90 +++ b/src/GRHydro_Prim2ConM.F90 @@ -294,6 +294,7 @@ subroutine prim2conM_hot(handle, GRHydro_reflevel, ii, jj, kk, x, y, z, gxx, gxy ! away), we use the average temperature between cells and call ! the EOS with keytemp=1 keytemp=1 + temp = temp0 call EOS_Omni_press(handle,keytemp,GRHydro_eos_rf_prec,n,& drho,deps,temp,ye,dpress,keyerr,anyerr) keytemp=0 |