From b62c3c52a5c5c421cfa2d17626ceab0ba5a90ea7 Mon Sep 17 00:00:00 2001
From: rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26>
Date: Thu, 9 Aug 2012 06:26:39 +0000
Subject: GRHydro_InitData: make rotor parameters Cactus parameters

* cite paper with description of system
* reverse sense of rotation
* add openmp statements

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

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@139 ac85fae7-cede-4708-beff-ae01c7fa1c26
---
 param.ccl | 71 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++-
 1 file changed, 70 insertions(+), 1 deletion(-)

(limited to 'param.ccl')

diff --git a/param.ccl b/param.ccl
index 09c6d3a..950be6d 100644
--- a/param.ccl
+++ b/param.ccl
@@ -14,7 +14,7 @@ EXTENDS KEYWORD initial_hydro ""
   "simple_wave" :: "Set initial data from Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)"
   "monopole"    :: "Monopole at the center"
   "cylexp"      :: "Cylindrical Explosion"
-  "rotor"       :: "Magnetic Rotor test"
+  "rotor"       :: "Magnetic Rotor test from DelZanna,Bucciantini, and Londrillo A&A 400, 397–413 (2003)"
   "advectedloop":: "Magnetic advected loop test"
   "alfvenwave"  :: "Circularly polarized Alfven wave"
   "hydro_bondi_solution" :: "Spherical single black hole Bondi solution"
@@ -275,6 +275,75 @@ CCTK_REAL poloidal_rho_max  "Maximum initial density"
   (0:*	:: "Anything positive."
 } 1.0e-3
 
+# for the Magnetic Rotor test:
+# default values are from DelZanna,Bucciantini, and Londrillo A&A 400, 397–413 (2003) though notation differs
+
+CCTK_REAL rotor_xc "center of rotation"
+{
+  *:* :: "Any location"
+} 0.5
+
+CCTK_REAL rotor_yc "center of rotation"
+{
+  *:* :: "Any location"
+} 0.5
+
+CCTK_REAL rotor_bvcxl "intial component of Bvec[0]"
+{
+  *:* :: "any real number"
+} 1.0
+
+CCTK_REAL rotor_bvcyl "intial component of Bvec[1]"
+{
+  *:* :: "any real number"
+} 0.0
+
+CCTK_REAL rotor_bvczl "intial component of Bvec[2]"
+{
+  *:* :: "any real number"
+} 0.0
+  
+CCTK_REAL rotor_r_rot "radius of rotor"
+{
+  (0:* :: "any positive number"
+} 0.1
+  
+CCTK_REAL rotor_v_max "Maximum velocity"
+{
+  (-1:1) :: "any subluminal speed (negative is clockwise)"
+} 0.995
+
+
+CCTK_REAL rotor_rhoin "initial density inside rotor"
+{
+  (0:* :: "any positive number"
+} 10.d0
+
+CCTK_REAL rotor_pressin "initial pressure inside rotor"
+{
+  (0:* :: "any positive number"
+} 1.d0
+
+CCTK_REAL rotor_rhoout "initial density outside rotor"
+{
+  (0:* :: "any positive number"
+} 1.d0
+
+CCTK_REAL rotor_pressout "initial pressure outside rotor"
+{
+  (0:* :: "any positive number"
+} 1.d0
+
+CCTK_BOOLEAN rotor_use_smoothing "Smooth the edge?"
+{
+} yes
+
+CCTK_REAL rotor_rsmooth_rel "Define the radius in relative terms if so"
+{
+  (0:* :: "any positive number"
+} 0.05
+
+
 shares:GRHydro
 
 USES real GRHydro_eos_rf_prec
-- 
cgit v1.2.3