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
|
# Schedule definitions for thorn GRHydro_Init_Data
schedule GRHydro_InitData_CheckParameters AT CCTK_PARAMCHECK
{
LANG: C
} "Check parameters"
if (CCTK_Equals(initial_hydro,"shocktube")) {
if(CCTK_Equals(Bvec_evolution_method,"GRHydro"))
{
schedule GRHydro_shocktubeM in HydroBase_Initial
{
LANG: Fortran
} "Shocktube initial data - MHD version"
} else {
schedule GRHydro_shocktube in HydroBase_Initial
{
LANG: Fortran
} "Shocktube initial data"
}
}
if (CCTK_Equals(initial_data,"con2primtest")) {
STORAGE:GRHydro_init_data_reflevel
schedule GRHydro_Init_Data_RefinementLevel IN HydroBase_Initial BEFORE GRHydro_con2primtest
{
LANG: Fortran
} "Calculate current refinement level"
schedule GRHydro_con2primtest in HydroBase_Initial
{
LANG: Fortran
} "Testing the conservative to primitive solver"
}
if (CCTK_Equals(initial_data,"con2prim2con_test")) {
STORAGE:GRHydro_init_data_reflevel
schedule GRHydro_Init_Data_RefinementLevel IN HydroBase_Initial BEFORE c2p2c_call
{
LANG: Fortran
} "Calculate current refinement level"
if(CCTK_Equals(Bvec_evolution_method,"GRHydro"))
{
schedule c2p2cM in HydroBase_Initial AS c2p2c_call
{
LANG: Fortran
} "Testing conservative to primitive to conservative - MHD version"
} else {
schedule c2p2c in HydroBase_Initial AS c2p2c_call
{
LANG: Fortran
} "Testing conservative to primitive to conservative"
}
}
if (CCTK_Equals(initial_data,"prim2con2prim_test")) {
STORAGE:GRHydro_init_data_reflevel
schedule GRHydro_Init_Data_RefinementLevel IN HydroBase_Initial BEFORE p2c2p_call
{
LANG: Fortran
} "Calculate current refinement level"
if(CCTK_Equals(Bvec_evolution_method,"GRHydro"))
{
schedule p2c2pM in HydroBase_Initial AS p2c2p_call
{
LANG: Fortran
} "Testing primitive to conservative to primitive - MHD version"
} else {
schedule p2c2p in HydroBase_Initial AS p2c2p_call
{
LANG: Fortran
} "Testing primitive to conservative to primitive"
}
}
if (CCTK_Equals(initial_data,"prim2con2prim_polytype_test")) {
STORAGE:GRHydro_init_data_reflevel
schedule GRHydro_Init_Data_RefinementLevel IN HydroBase_Initial BEFORE p2c2p_call
{
LANG: Fortran
} "Calculate current refinement level"
if(CCTK_Equals(Bvec_evolution_method,"GRHydro"))
{
schedule p2c2pM_polytype in HydroBase_Initial AS p2c2p_call
{
LANG: Fortran
} "Testing primitive to conservative to primitive - MHD polytype version"
}
}
if (CCTK_Equals(initial_data,"reconstruction_test")) {
schedule GRHydro_reconstruction_test in HydroBase_Initial
{
LANG: Fortran
STORAGE: GRHydro_prim_bext
STORAGE: GRHydro_scalars
OPTIONS: global loop-local
} "Testing the reconstruction"
}
if (CCTK_Equals(initial_hydro,"only_atmo")) {
schedule GRHydro_Only_Atmo in HydroBase_Initial
{
LANG: Fortran
} "Only atmosphere as initial data"
}
if (CCTK_Equals(initial_hydro,"read_conformal")) {
schedule GRHydro_ReadConformalData in HydroBase_Initial
{
LANG: Fortran
} "Set the missing quantities, after reading in from file initial data from conformally-flat codes (Garching)"
}
if (CCTK_Equals(initial_hydro,"simple_wave")) {
STORAGE: simple_wave_grid_functions
STORAGE: simple_wave_scalars
schedule GRHydro_SimpleWave in HydroBase_Initial
{
LANG: Fortran
} "Set initial data from Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)"
STORAGE: simple_wave_output
schedule GRHydro_SimpleWave_Analysis AT CCTK_ANALYSIS AFTER GRHydro_Entropy
{
LANG: Fortran
} "Compute some output variables for the Simple Wave"
}
|
