Snapshot of ET GRHydro_InitData rev. 132 with a poloidal magnetic field routine

added.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@134 ac85fae7-cede-4708-beff-ae01c7fa1c26

1 files changed, 154 insertions, 1 deletions

diff --git a/param.ccl b/param.ccl
index 8d4da47..09c6d3a 100644
--- a/param.ccl
+++ b/param.ccl
@@ -12,8 +12,15 @@ EXTENDS KEYWORD initial_hydro ""
   "only_atmo"	:: "Set only a low atmosphere"
   "read_conformal":: "After reading in initial alp, rho and gxx from h5 files, sets the other quantities"
   "simple_wave" :: "Set initial data from Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)"
-  "Monopole"    :: "Monopole at the center"
+  "monopole"    :: "Monopole at the center"
   "cylexp"      :: "Cylindrical Explosion"
+  "rotor"       :: "Magnetic Rotor test"
+  "advectedloop":: "Magnetic advected loop test"
+  "alfvenwave"  :: "Circularly polarized Alfven wave"
+  "hydro_bondi_solution" :: "Spherical single black hole Bondi solution"
+  "hydro_bondi_solution_iso" :: "Spherical single black hole Bondi solution - TEST ISO CASE!!!!!!"
+  "magnetized_bondi_solution" :: "Magnetized Spherical single black hole Bondi solution"
+  "magnetized_bondi_solution_iso" :: "Magnetized Spherical single black hole Bondi solution - TEST ISO CASE!!!!!!"
 }
 
 shares:ADMBase
@@ -30,6 +37,12 @@ EXTENDS KEYWORD initial_data ""
 
 private:
 
+KEYWORD initial_Bvec "Initial magnetic field."
+{
+  "none"             :: "No initial magnetic field is set."
+  "poloidalmagfield" :: "Poloidal Magnetic Field"
+} "none"
+
 KEYWORD shocktube_type "Diagonal or parallel shock?"
 {
   "diagshock"	:: "Diagonal across all axes"
@@ -124,6 +137,143 @@ REAL Bz_init "Initial B-field in the z-dir"
   *:*     :: "Anything"
 } 0.0
 
+KEYWORD monopole_type "Which kind of monopole?"
+{
+  "Point" :: "Single point with Bx /= 0"
+  "Gauss" :: "Gaussian w/radius R_Gauss"
+  "1dalt" :: "1-d alternating"
+  "2dalt" :: "2-d alternating"
+} "Point"
+
+CCTK_REAL R_Gauss "Radius for a Gaussian monopole"
+{
+  0:* :: "Any positive number"
+} 1.0
+
+# For cylindrical explosion test:
+CCTK_REAL cyl_r_inner "Inner Radius"
+{
+  (0:* :: "Any positive number"
+} 0.8
+
+CCTK_REAL cyl_r_outer "Inner Radius"
+{
+  (0:* :: "Any positive number"
+} 1.0
+
+# For magnetic advected loop test:
+
+KEYWORD advectedloop_type "2-dimensional or 3-dimensional?"
+{
+  "2D"	:: "2-dimensional (B^z=0)"
+  "3D"	:: "3-dimensional (B^3=0, where B^3 || oblique cylinder axis."
+} "2D"
+
+KEYWORD advectedloop_case "V^z=0 or not?"
+{
+  "V^z=0"	:: "Useful to evaluate divB deviations"
+  "V^z/=0"	:: "Useful to evaluate con2prim robustness in keeping V^z const."
+} "V^z=0"
+
+KEYWORD advectedloop_delA "How to calculate B^i field from the potential A^b"
+{
+  "Exact"	:: "Analytic, exact closed formula applied"
+  "Numeric"	:: "Finite difference approximation of the derivatives applied"
+} "Exact"
+
+##################################################################################
+# BONDI PARAMETERS: (black hole mass specified by parameters from the "Exact" thorn)
+##################################################################################
+# Default value is that used in the HARM test
+CCTK_REAL mdot_sonicpt_bondi "Accretion rate at sonic point in hydro units"
+{
+  (0:* :: "positive"
+} 12.566370614359172954
+
+CCTK_REAL  r_sonicpt_bondi  "Radial distance of the sonic point from the black hole in units of mass_bh_bondi"
+{
+  (0:* :: "positive"
+} 8.0
+
+keyword bondi_coordinates  "Which coordinate system to use"
+{
+  "Boyer-Lindquist" :: "Schwarzschild or Boyer-Lindquist (Cartesian) Coordinates"
+  "Kerr-Schild"   :: "Kerr-Schild (Cartesian) Coordinates"
+  "Isotropic"     :: "Isotropic (Cartesian) Coordinates"
+} "Isotropic" 
+
+
+CCTK_INT num_bondi_sols  "Number of central masses about which to calculate Bondi solutions"
+{
+  1:100 :: "positive"
+} 1
+
+CCTK_REAL bh_bondi_pos_x [100] "X-coordinate of black hole in Bondi solution"
+{
+  *:* :: "anything"
+} 0.0
+
+CCTK_REAL bh_bondi_pos_y [100] "Y-coordinate of black hole in Bondi solution"
+{
+  *:* :: "anything"
+} 0.0
+
+CCTK_REAL bh_bondi_pos_z [100] "Z-coordinate of black hole in Bondi solution"
+{
+  *:* :: "anything"
+} 0.0
+
+CCTK_REAL bondi_central_mass [100] "Mass of central object to find Bondi solution about  "
+{
+  (0:* :: "positive"
+} 1.0
+
+CCTK_REAL bondi_central_spin [100] "Dimensionless spin of central object within Bondi solution "
+{
+  (-1.:1.) :: "dimensionless spin so any real number between -1 and 1"
+} 0.0
+
+CCTK_REAL bondi_rmin [100] "Smallest radius in units of central mass at which the solution is found"
+{
+  (0:*  :: "dimensionless inner radius for Bondi solution"
+} 1.e-15
+
+CCTK_REAL bondi_rmax [100] "Largest radius in units of central mass at which the solution is found"
+{
+  (0:*  :: "dimensionless outer radius for Bondi solution"
+} 400. 
+
+CCTK_REAL n_bondi_pts [100] "Number of points to use in determining global Bondi solution"
+{
+  (1:*  :: "number of points in global Bondi solution"
+} 2000
+
+CCTK_REAL bondi_bmag "B_0 parameter for magnetized Bondi"
+{
+  0:* :: "Anything positive"
+} 0.01
+
+# For Poloidal Magnetic field test:
+
+CCTK_REAL poloidal_A_b  "Vector potential strength"
+{
+  *:*	:: "Anything."
+} 0.1
+
+CCTK_INT poloidal_n_p  "Vector potential strength"
+{
+  (0:*	:: "Any positive integer."
+} 3
+
+CCTK_REAL poloidal_P_cut  "Pressure used to confine the B field inside a star"
+{
+  (0:*	:: "Anything positive."
+} 1.0e-8
+
+CCTK_REAL poloidal_rho_max  "Maximum initial density"
+{
+  (0:*	:: "Anything positive."
+} 1.0e-3
 
 shares:GRHydro
 
@@ -141,3 +291,6 @@ USES REAL initial_rho_rel_min
 USES REAL initial_atmosphere_factor
 USES int GRHydro_stencil
 USES BOOLEAN clean_divergence
+
+shares:EOS_Omni
+USES REAL gl_gamma