diff options
author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-03-28 01:49:11 +0000 |
---|---|---|
committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-03-28 01:49:11 +0000 |
commit | 37086b8a20d91e71f8a6e260aaa25ed490cdee10 (patch) | |
tree | 28ab526c910d14933bf82be793a52d044a58af91 /src | |
parent | 0d4b0b4617f07e601074368d5b4b01a873750c38 (diff) |
GRHydro_InitData: make cylindrical explosion parameters parameters
tie Gamma value used to that of eos_omni's gamma law eos
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@208 ac85fae7-cede-4708-beff-ae01c7fa1c26
Diffstat (limited to 'src')
-rw-r--r-- | src/GRHydro_CylindricalExplosionM.F90 | 18 |
1 files changed, 3 insertions, 15 deletions
diff --git a/src/GRHydro_CylindricalExplosionM.F90 b/src/GRHydro_CylindricalExplosionM.F90 index 6c76542..f8a0a7a 100644 --- a/src/GRHydro_CylindricalExplosionM.F90 +++ b/src/GRHydro_CylindricalExplosionM.F90 @@ -74,21 +74,7 @@ subroutine GRHydro_cylindricalexplosionM(CCTK_ARGUMENTS) bvczl = Bz_init !!$Adiabatic index for test - gam = (4.d0/3.d0) - -!!$Inner radius and outer radius (Komissarov's defaults) -! r_inner = 8.d-1 -! r_outer = 1.d0 - r_inner = cyl_r_inner - r_outer = cyl_r_outer - -!!$Inner values - rhol = 1.d-4 - pressl = 3.d-5 - -!!$Outer values - rhor = 1.d-2 - pressr = 1.d0 + gam = gl_gamma nx = cctk_lsh(1) ny = cctk_lsh(2) @@ -99,6 +85,8 @@ subroutine GRHydro_cylindricalexplosionM(CCTK_ARGUMENTS) do k=1,nz !!$direction represents the cylindrical radius here. +!!$TODO: maybe switch these over so that shocktube_tube choses the axis of the +!!$cylinder instead of the apparently random mapping if (CCTK_EQUALS(shocktube_type,"xshock")) then direction = sqrt((x(i,j,k)-shock_xpos)**2+& (y(i,j,k)-shock_ypos)**2) |