diff options
author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-10-26 14:38:03 +0000 |
---|---|---|
committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-10-26 14:38:03 +0000 |
commit | aee279db11e1a21e1d04da5e4a78e31140fe9a23 (patch) | |
tree | 95e89beff8cd3fffcb71e5bd0dd6585d311130fc | |
parent | 91a2ee778dc93e45a759810d790d04e9dfcbdd04 (diff) |
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
-rw-r--r-- | param.ccl | 1 | ||||
-rw-r--r-- | src/GRHydro_CylindricalExplosionM.F90 | 6 |
2 files changed, 7 insertions, 0 deletions
@@ -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 |