diff options
Diffstat (limited to 'src/GRHydro_CylindricalExplosionM.F90')
-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) |