[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 208)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Wed Mar 27 20:49:11 CDT 2013
User: rhaas
Date: 2013/03/27 08:49 PM
Modified:
/trunk/
param.ccl
/trunk/src/
GRHydro_CylindricalExplosionM.F90
Log:
GRHydro_InitData: make cylindrical explosion parameters parameters
tie Gamma value used to that of eos_omni's gamma law eos
From: Roland Haas <rhaas at tapir.caltech.edu>
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_CylindricalExplosionM.F90
Delta lines: +3 -15
===================================================================
--- trunk/src/GRHydro_CylindricalExplosionM.F90 2013-03-28 01:49:10 UTC (rev 207)
+++ trunk/src/GRHydro_CylindricalExplosionM.F90 2013-03-28 01:49:11 UTC (rev 208)
@@ -74,22 +74,8 @@
bvczl = Bz_init
!!$Adiabatic index for test
- gam = (4.d0/3.d0)
+ gam = gl_gamma
-!!$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 = 1.d-4
- pressl = 3.d-5
-
-!!$Outer values
- rhor = 1.d-2
- pressr = 1.d0
-
nx = cctk_lsh(1)
ny = cctk_lsh(2)
nz = cctk_lsh(3)
@@ -99,6 +85,8 @@
do k=1,nz
!!$direction represents the cylindrical radius here.
+!!$TODO: maybe switch these over so that shocktube_tube choses the axis of the
+!!$cylinder instead of the apparently random mapping
if (CCTK_EQUALS(shocktube_type,"xshock")) then
direction = sqrt((x(i,j,k)-shock_xpos)**2+&
(y(i,j,k)-shock_ypos)**2)
Directory: /trunk/
==================
File [modified]: param.ccl
Delta lines: +21 -0
===================================================================
--- trunk/param.ccl 2013-03-28 01:49:10 UTC (rev 207)
+++ trunk/param.ccl 2013-03-28 01:49:11 UTC (rev 208)
@@ -280,6 +280,27 @@
(0:* :: "Any positive number"
} 1.0
+# the defaults are Komissarov's values
+CCTK_REAL cyl_rho_inner "density in inner core"
+{
+ (0:* :: "any positive number"
+} 1.d-4
+
+CCTK_REAL cyl_press_inner "pressure in inner core"
+{
+ (0:* :: "any positive number"
+} 3.d-5
+
+CCTK_REAL cyl_rho_outer "density in outer region"
+{
+ (0:* :: "any positive number"
+} 1.d-2
+
+CCTK_REAL cyl_press_outer "pressure in outer region"
+{
+ (0:* :: "any positive number"
+} 1.d0
+
# For magnetic advected loop test:
KEYWORD advectedloop_type "2-dimensional or 3-dimensional?"
More information about the Commits
mailing list