aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorrhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26>2013-03-28 14:26:34 +0000
committerrhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26>2013-03-28 14:26:34 +0000
commitdbf5b7f2ee9d7a8dec2301cb8ecd554e552831f7 (patch)
tree7535034cf503174bddc0dc339235c172ee510cbd
parent37086b8a20d91e71f8a6e260aaa25ed490cdee10 (diff)
GRHydro_InitData: fix turning cylexp options into parameters
From: Roland Haas <rhaas@tapir.caltech.edu> git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@209 ac85fae7-cede-4708-beff-ae01c7fa1c26
-rw-r--r--src/GRHydro_CylindricalExplosionM.F9014
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