diff options
-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 |