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
|
# 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,"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"
}
|