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