diff options
| author | bmundim <bmundim@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2011-04-27 22:25:23 +0000 |
|---|---|---|
| committer | bmundim <bmundim@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2011-04-27 22:25:23 +0000 |
| commit | ed1c9f732192bdc8d4f6d5b3237a621855f1b23a (patch) | |
| tree | c1c980799095f90e7964bf291bffa49ea5295343 /src/GRHydro_ShockTube.F90 | |
| parent | cef4cc9bdbe84dca35260f6f02c2249b4cc99856 (diff) | |
RIT MHD dev:
Balsara's and Komissarov's shock tube tests.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@123 ac85fae7-cede-4708-beff-ae01c7fa1c26
Diffstat (limited to 'src/GRHydro_ShockTube.F90')
| -rw-r--r-- | src/GRHydro_ShockTube.F90 | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/src/GRHydro_ShockTube.F90 b/src/GRHydro_ShockTube.F90 index baf891d..8b9105b 100644 --- a/src/GRHydro_ShockTube.F90 +++ b/src/GRHydro_ShockTube.F90 @@ -93,6 +93,17 @@ subroutine GRHydro_shocktube(CCTK_ARGUMENTS) velzr = 0.d0 epsl = 1.5d0 epsr = 1.2d0 + else if (CCTK_EQUALS(shock_case,"Balsaralike1")) then + rhol = 1.0d0 + rhor = 0.125d0 + velxl = 0.0d0 + velxr = 0.0d0 + velyl = 0.0d0 + velyr = 0.0d0 + velzl = 0.0d0 + velzr = 0.0d0 + epsl = 1.0d0/rhol + epsr = 0.1d0/rhor !!$This line only for polytrope, k=1 !!$ epsr = 0.375d0 else if (CCTK_EQUALS(shock_case,"Blast")) then |