From aee279db11e1a21e1d04da5e4a78e31140fe9a23 Mon Sep 17 00:00:00 2001 From: rhaas Date: Fri, 26 Oct 2012 14:38:03 +0000 Subject: 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 git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@173 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- param.ccl | 1 + src/GRHydro_CylindricalExplosionM.F90 | 6 ++++++ 2 files changed, 7 insertions(+) diff --git a/param.ccl b/param.ccl index cf0a7ba..e36a67a 100644 --- a/param.ccl +++ b/param.ccl @@ -29,6 +29,7 @@ EXTENDS KEYWORD initial_hydro "" EXTENDS KEYWORD initial_Bvec { "shocktube" :: "Shocktube type" + "cylexp" :: "Poloidal Magnetic Field" "poloidalmagfield" :: "Poloidal Magnetic Field" } shares:ADMBase diff --git a/src/GRHydro_CylindricalExplosionM.F90 b/src/GRHydro_CylindricalExplosionM.F90 index 61a8ab0..6c76542 100644 --- a/src/GRHydro_CylindricalExplosionM.F90 +++ b/src/GRHydro_CylindricalExplosionM.F90 @@ -59,6 +59,12 @@ subroutine GRHydro_cylindricalexplosionM(CCTK_ARGUMENTS) 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 -- cgit v1.2.3