User: rhaas Date: 2013/04/10 12:54 PM
Modified: /trunk/ param.ccl /trunk/src/ GRHydro_CylindricalExplosionM.F90
Log: GRHydro: revert meaning of cyl_r_inner and cyl_r_outer so that they now
make sense
From: Roland Haas rhaas@caltech.edu
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_CylindricalExplosionM.F90 Delta lines: +16 -16 =================================================================== --- trunk/src/GRHydro_CylindricalExplosionM.F90 2013-03-28 14:26:34 UTC (rev 209) +++ trunk/src/GRHydro_CylindricalExplosionM.F90 2013-04-10 17:54:27 UTC (rev 210) @@ -65,28 +65,28 @@ 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
Directory: /trunk/ ==================
File [modified]: param.ccl Delta lines: +4 -4 =================================================================== --- trunk/param.ccl 2013-03-28 14:26:34 UTC (rev 209) +++ trunk/param.ccl 2013-04-10 17:54:27 UTC (rev 210) @@ -284,22 +284,22 @@ CCTK_REAL cyl_rho_inner "density in inner core" { (0:* :: "any positive number" -} 1.d-4 +} 1.d-2
CCTK_REAL cyl_press_inner "pressure in inner core" { (0:* :: "any positive number" -} 3.d-5 +} 1.d0
CCTK_REAL cyl_rho_outer "density in outer region" { (0:* :: "any positive number" -} 1.d-2 +} 1.d-4
CCTK_REAL cyl_press_outer "pressure in outer region" { (0:* :: "any positive number" -} 1.d0 +} 3.d-5
# For magnetic advected loop test: