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 --- src/GRHydro_CylindricalExplosionM.F90 | 6 ++++++ 1 file changed, 6 insertions(+) (limited to 'src') 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