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@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?"