[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