diff options
| -rw-r--r-- | src/GRHydro_CylindricalExplosionM.F90 | 14 |
1 files changed, 14 insertions, 0 deletions
diff --git a/src/GRHydro_CylindricalExplosionM.F90 b/src/GRHydro_CylindricalExplosionM.F90 index f8a0a7a..d038c18 100644 --- a/src/GRHydro_CylindricalExplosionM.F90 +++ b/src/GRHydro_CylindricalExplosionM.F90 @@ -73,6 +73,20 @@ subroutine GRHydro_cylindricalexplosionM(CCTK_ARGUMENTS) bvcyl = By_init bvczl = Bz_init +!!$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 = cyl_rho_inner + pressl = cyl_press_inner + +!!$Outer values + rhor = cyl_rho_outer + pressr = cyl_press_outer + !!$Adiabatic index for test gam = gl_gamma |