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

[rhaas at tapir.caltech.edu]

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 at 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: