blob: 8d4da47f8c85727db1e1f804028a986ff7db723b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
|
# Parameter definitions for thorn GRHydro_Init_Data
# $Header$
shares:HydroBase
USES CCTK_INT timelevels
USES KEYWORD Bvec_evolution_method
EXTENDS KEYWORD initial_hydro ""
{
"shocktube" :: "Shocktube type"
"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"
"cylexp" :: "Cylindrical Explosion"
}
shares:ADMBase
EXTENDS KEYWORD initial_data ""
{
# "shocktube" :: "Shock tube initial data for GRHydro"
"con2primtest" :: "Testing the con -> prim conversion"
"con2prim2con_test" :: "Testing the con -> prim -> con conversion"
"prim2con2prim_test" :: "Testing the prim -> con -> prim conversion"
"prim2con2prim_polytype_test" :: "Testing the prim -> con -> prim conversion - polytype version"
"reconstruction_test" :: "Testing reconstruction"
}
private:
KEYWORD shocktube_type "Diagonal or parallel shock?"
{
"diagshock" :: "Diagonal across all axes"
"diagshock2d" :: "Diagonal across x-y axes"
"xshock" :: "Parallel to x axis"
"yshock" :: "Parallel to y axis"
"zshock" :: "Parallel to z axis"
"sphere" :: "spherically symmetric shock"
} "xshock"
KEYWORD shock_case "Simple, Sod's problem or other?"
{
"Simple" :: "GRAstro_Hydro test case"
"Sod" :: "Sod's problem"
"Blast" :: "Strong blast wave"
"Balsaralike1" :: "Hydro version of Balsara Test #1"
"Balsara0" :: "Balsara Test #1, but unmagnetized"
"Balsara1" :: "Balsara Test #1"
"Balsara2" :: "Balsara Test #2"
"Balsara3" :: "Balsara Test #3"
"Balsara4" :: "Balsara Test #4"
"Balsara5" :: "Balsara Test #5"
"Alfven" :: "Generical Alfven Test"
"Komissarov1" :: "Komissarov Test #1"
"Komissarov2" :: "Komissarov Test #2"
"Komissarov3" :: "Komissarov Test #3"
"Komissarov4" :: "Komissarov Test #4"
"Komissarov5" :: "Komissarov Test #5"
"Komissarov6" :: "Komissarov Test #6"
"Komissarov7" :: "Komissarov Test #7"
"Komissarov8" :: "Komissarov Test #8"
"Komissarov9" :: "Komissarov Test #9"
} "Sod"
REAL shock_xpos "Position of shock plane: x"
{
*:* :: "Anything"
} 0.0
REAL shock_ypos "Position of shock plane: y"
{
*:* :: "Anything"
} 0.0
REAL shock_zpos "Position of shock plane: z"
{
*:* :: "Anything"
} 0.0
REAL shock_radius "Radius of sperical shock"
{
0.0:* :: "Anything positive"
} 1.0
BOOLEAN change_shock_direction "Change the shock direction"
{
} "no"
REAL simple_wave_constant_c_0 "The c_0 constant in Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)"
{
0:1 :: "It is the sound speed where the fluid velocity is zero"
} 0.3
REAL simple_wave_v_max "The v_max constant in Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)"
{
0:1 :: "It is the maximum velocity in the initial configuration (see p. 1457, bottom of first column)"
} 0.7
# For the "atmosphere"
REAL atmosphere_vel[3] "Velocity of the atmosphere if non-trivial"
{
*:* :: "Anything"
} 0.0
BOOLEAN attenuate_atmosphere "Attenuate the velocity in the atmosphere"
{
} "no"
REAL Bx_init "Initial B-field in the x-dir"
{
*:* :: "Anything"
} 0.0
REAL By_init "Initial B-field in the y-dir"
{
*:* :: "Anything"
} 0.0
REAL Bz_init "Initial B-field in the z-dir"
{
*:* :: "Anything"
} 0.0
shares:GRHydro
USES real GRHydro_eos_rf_prec
USES real GRHydro_rho_central
USES real GRHydro_eps_min ""
USES real GRHydro_perc_ptol ""
USES real GRHydro_del_ptol ""
USES string GRHydro_eos_type ""
USES string GRHydro_eos_table ""
USES real rho_abs_min
USES real rho_rel_min
USES REAL initial_rho_abs_min
USES REAL initial_rho_rel_min
USES REAL initial_atmosphere_factor
USES int GRHydro_stencil
USES BOOLEAN clean_divergence
|
