diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-08-27 19:19:35 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-08-27 19:19:35 +0000 |
| commit | e114f22b1a50d5ad237ef7a94c75ff20591b8c4c (patch) | |
| tree | 3a12cea7029bd180efe2c06c97779c96dfd2700f /src/GRHydro_Con2Prim.F90 | |
| parent | df593a5d92af8f243eaade956d17d6b8791fb36c (diff) | |
GRHydro: Changes fixing issues with hot EOS treatment:
1) Set correct keytemp in C2P
2) Fix OpenMP private variable declaration in HLLEM
3) Fix velocity arguments for Y_e 1D PPM reconstruction
4) If in atmosphere and if evolve a Y_e also reset Y_e with cell averaged value in Reconstruction
5) Fix index permutation bug in ReconstructPoly routines for divergence cleaning field
6) Remove extra comment character in UpdateMaskM
7) Fix a missing velocity permutation when calling 1D reconstruction for Y_e.
8) Insert missing sqrt(det) factor for pressure term in resetting tau when temperature got too cold with a hot EOS.
From: Philipp Moesta <pmoesta@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@415 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Con2Prim.F90')
| -rw-r--r-- | src/GRHydro_Con2Prim.F90 | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/GRHydro_Con2Prim.F90 b/src/GRHydro_Con2Prim.F90 index 22afcd3..bb53c0a 100644 --- a/src/GRHydro_Con2Prim.F90 +++ b/src/GRHydro_Con2Prim.F90 @@ -327,7 +327,7 @@ subroutine Conservative2Primitive(CCTK_ARGUMENTS) call EOS_Omni_press(GRHydro_eos_handle,keytemp,GRHydro_eos_rf_prec,n,& rho(i,j,k),eps(i,j,k),temperature(i,j,k),y_e(i,j,k),press(i,j,k),keyerr,anyerr) tau(i,j,k) = sqrt(det) * (rho(i,j,k)+rho(i,j,k)*eps(i,j,k) +press(i,j,k))& - * w_lorentz(i,j,k)**2 - dens(i,j,k) - press(i,j,k) + * w_lorentz(i,j,k)**2 - dens(i,j,k) - sqrt(det)*press(i,j,k) keytemp = 0 endif |