[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 183)

[rhaas at tapir.caltech.edu]

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 at 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