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_PPMReconstruct_drv.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_PPMReconstruct_drv.F90')
| -rw-r--r-- | src/GRHydro_PPMReconstruct_drv.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/GRHydro_PPMReconstruct_drv.F90 b/src/GRHydro_PPMReconstruct_drv.F90 index 02486d0..0a6bcb6 100644 --- a/src/GRHydro_PPMReconstruct_drv.F90 +++ b/src/GRHydro_PPMReconstruct_drv.F90 @@ -246,7 +246,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_ye_1d(ny,CCTK_DELTA_SPACE(2),rho(j,:,k), & - velx(j,:,k),vely(j,:,k),velz(j,:,k), & + vely(j,:,k),velz(j,:,k),velx(j,:,k), & Y_e(j,:,k),Y_e_minus(j,:,k),Y_e_plus(j,:,k), & press(j,:,k)) end do @@ -298,7 +298,7 @@ subroutine GRHydro_PPMReconstruct_drv(CCTK_ARGUMENTS) do k = GRHydro_stencil, ny - GRHydro_stencil + 1 do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_ye_1d(nz,CCTK_DELTA_SPACE(3),rho(j,k,:), & - velx(j,k,:),vely(j,k,:),velz(j,k,:), & + velz(j,k,:),velx(j,k,:),vely(j,k,:), & Y_e(j,k,:),Y_e_minus(j,k,:),Y_e_plus(j,k,:), & press(j,k,:)) end do |