From 6670e935f9a18108f5a9a16a69de3ccaf8996dbe Mon Sep 17 00:00:00 2001 From: rhaas Date: Wed, 10 Apr 2013 17:54:27 +0000 Subject: GRHydro: revert meaning of cyl_r_inner and cyl_r_outer so that they now make sense From: Roland Haas git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@210 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_CylindricalExplosionM.F90 | 32 ++++++++++++++++---------------- 1 file changed, 16 insertions(+), 16 deletions(-) (limited to 'src') diff --git a/src/GRHydro_CylindricalExplosionM.F90 b/src/GRHydro_CylindricalExplosionM.F90 index d038c18..f6feeac 100644 --- a/src/GRHydro_CylindricalExplosionM.F90 +++ b/src/GRHydro_CylindricalExplosionM.F90 @@ -65,27 +65,27 @@ subroutine GRHydro_cylindricalexplosionM(CCTK_ARGUMENTS) if (.not. CCTK_EQUALS(initial_Bvec, "cylexp")) then call CCTK_WARN(1, "When using the cyclindrical explosion initial data, please also select 'cylexp' for initial_Bvec. I will proceed as if this was set") end if - + !!$Uses Bx_init, By_init and Bz_init to set magnetic field strength !!$Original tests had Bx_init = 0.1, 1.0 with By_init=Bz_init = 0 - + bvcxl = Bx_init 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 + +!!$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 + rhor = cyl_rho_inner + pressr = cyl_press_inner + +!!$Outer values + rhol = cyl_rho_outer + pressl = cyl_press_outer !!$Adiabatic index for test gam = gl_gamma -- cgit v1.2.3