User: rhaas Date: 2012/11/08 07:59 PM

Modified: /trunk/ param.ccl /trunk/src/ CheckParam.c

Log: GRHydro_Init_Data: add "magnetized Bondi" to hydrobase::initial_bvec

this is so that there is a proper initial_bvec setting rather than just "zero".

From: Roland Haas roland.haas@physics.gatech.edu

File Changes:

Directory: /trunk/src/ ======================

File [modified]: CheckParam.c Delta lines: +8 -2 =================================================================== --- trunk/src/CheckParam.c 2012-11-09 01:59:16 UTC (rev 182) +++ trunk/src/CheckParam.c 2012-11-09 01:59:19 UTC (rev 183) @@ -50,15 +50,21 @@ if(CCTK_Equals(Bvec_evolution_method,"GRHydro") && CCTK_Equals(initial_hydro,"hydro_bondi_solution") ) { - CCTK_PARAMWARN("Please use set initial_hydro=magnetized_bondi_solution instead to initialize the magnetic field correctly!"); + CCTK_PARAMWARN("Please set initial_hydro='magnetized_bondi_solution' instead to initialize the magnetic field correctly!"); }

if(!CCTK_Equals(Bvec_evolution_method,"GRHydro") && CCTK_Equals(initial_hydro,"magnetized_bondi_solution") ) { - CCTK_PARAMWARN("Please use set initial_hydro=hydro_bondi_solution instead to NOT initialize the magnetic field!"); + CCTK_PARAMWARN("Please set initial_hydro='hydro_bondi_solution' instead to NOT initialize the magnetic field!"); }

+ if((CCTK_Equals(initial_hydro,"magnetized_bondi_solution") || CCTK_Equals(initial_hydro,"magnetized_bondi_solution_iso")) && + !CCTK_Equals(initial_Bvec, "magnetized Bondi") ) + { + CCTK_PARAMWARN("Please set initial_Bec='magnetized Bondi' to properly initialize the magnetic field for Bondi solutions!"); + } + if( CCTK_Equals(initial_hydro,"magnetized_bondi_solution") || CCTK_Equals(initial_hydro,"hydro_bondi_solution") ) { if( num_bondi_sols > 1 ) { CCTK_PARAMWARN("Currently only one Bondi solution is supported, please change [num_bondi_sols] ");

Directory: /trunk/ ==================

File [modified]: param.ccl Delta lines: +1 -0 =================================================================== --- trunk/param.ccl 2012-11-09 01:59:16 UTC (rev 182) +++ trunk/param.ccl 2012-11-09 01:59:19 UTC (rev 183) @@ -31,6 +31,7 @@ "shocktube" :: "Shocktube type" "cylexp" :: "Poloidal Magnetic Field" "poloidalmagfield" :: "Poloidal Magnetic Field" + "magnetized Bondi" :: "radial magnetic field appropriate for Bondi test" } shares:ADMBase