diff options
Diffstat (limited to 'param.ccl')
-rw-r--r-- | param.ccl | 71 |
1 files changed, 70 insertions, 1 deletions
@@ -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 |