[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 209)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/03/28 09:26 AM

Modified:
 /trunk/src/
  GRHydro_CylindricalExplosionM.F90

Log:
 GRHydro_InitData: fix turning cylexp options into parameters
 
 From: Roland Haas <rhaas at tapir.caltech.edu>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_CylindricalExplosionM.F90
Delta lines: +14 -0
===================================================================
--- trunk/src/GRHydro_CylindricalExplosionM.F90	2013-03-28 01:49:11 UTC (rev 208)
+++ trunk/src/GRHydro_CylindricalExplosionM.F90	2013-03-28 14:26:34 UTC (rev 209)
@@ -73,6 +73,20 @@
   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                                                                                       
+
 !!$Adiabatic index for test
   gam = gl_gamma