From dbf5b7f2ee9d7a8dec2301cb8ecd554e552831f7 Mon Sep 17 00:00:00 2001 From: rhaas Date: Thu, 28 Mar 2013 14:26:34 +0000 Subject: GRHydro_InitData: fix turning cylexp options into parameters From: Roland Haas git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@209 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_CylindricalExplosionM.F90 | 14 ++++++++++++++ 1 file changed, 14 insertions(+) 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 -- cgit v1.2.3