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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/10/26 09:38 AM

Modified:
 /trunk/
  param.ccl
 /trunk/src/
  GRHydro_CylindricalExplosionM.F90

Log:
 GRHydro_InitData: add cylindrical explosion ID type to initial_Bvec
 
 this avoids having to rely on GRHydro_InitData running after HydroBase's
 ZeroBvec routine in order to have non-zero magnetic field

File Changes:

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

File [modified]: GRHydro_CylindricalExplosionM.F90
Delta lines: +6 -0
===================================================================
--- trunk/src/GRHydro_CylindricalExplosionM.F90	2012-09-17 16:27:39 UTC (rev 172)
+++ trunk/src/GRHydro_CylindricalExplosionM.F90	2012-10-26 14:38:03 UTC (rev 173)
@@ -59,6 +59,12 @@
   CCTK_REAL :: tmp,tmp2,gam,r_inner,r_outer
 
   CCTK_REAL :: cyl_fr
+
+!!$ Check that user selected proper magnetic field ID. Warn if wrong one is
+!!$ selected and proceed to ignore the value
+  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

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

File [modified]: param.ccl
Delta lines: +1 -0
===================================================================
--- trunk/param.ccl	2012-09-17 16:27:39 UTC (rev 172)
+++ trunk/param.ccl	2012-10-26 14:38:03 UTC (rev 173)
@@ -29,6 +29,7 @@
 EXTENDS KEYWORD initial_Bvec
 {
   "shocktube"  :: "Shocktube type"
+  "cylexp"     :: "Poloidal Magnetic Field"
   "poloidalmagfield" :: "Poloidal Magnetic Field"
 }
 shares:ADMBase