blob: 57de555021a2b20876afc9529ca79cbc5d591587 (
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
|
# Schedule definitions for thorn HydroBase
if (timelevels == 1)
{
STORAGE: rho[1]
STORAGE: press[1]
STORAGE: eps[1]
STORAGE: vel[1]
if (!CCTK_EQUALS(initial_Y_e, "none"))
{
STORAGE: Y_e[1]
}
}
else if (timelevels == 2)
{
STORAGE: rho[2]
STORAGE: press[2]
STORAGE: eps[2]
STORAGE: vel[2]
if (!CCTK_EQUALS(initial_Y_e, "none"))
{
STORAGE: Y_e[2]
}
}
else if (timelevels == 3)
{
STORAGE: rho[3]
STORAGE: press[3]
STORAGE: eps[3]
STORAGE: vel[3]
if (!CCTK_EQUALS(initial_Y_e, "none"))
{
STORAGE: Y_e[3]
}
}
schedule group HydroBase_Initial at CCTK_INITIAL \
after (ADMBase_InitialData ADMBase_InitialGauge \
IOUtil_RecoverIDFromDatafiles) \
before (ADMBase_PostInitial SetTmunu)
{
} "HydroBase initial data group"
schedule HydroBase_StartUp at CCTK_STARTUP
{
LANG: C
} "Startup banner"
schedule group HydroBase_RHS in MoL_CalcRHS
{
} "Groups for scheduling tasks for calculating RHS of hydro variables"
schedule group HydroBase_PostStep in MoL_PostStep before SetTmunu
{
} "Post step tasks for hydro thorns"
schedule group HydroBase_PostStep at CCTK_POSTRESTRICTINITIAL
{
} "Post step tasks for hydro thorns"
# Hydro thorns should not schedule anything in group
# HydroBase_Boundaries; they should schedule in
# HydroBase_Select_Boundaries instead. However, they may refer to
# HydroBase_Boundaries in order to schedule own routines before or
# after the boundary treatment, or in order to schedule the whole
# group.
schedule group HydroBase_Boundaries in HydroBase_PostStep \
before HydroBase_Con2Prim
{
} "HydroBase-internal Boundary conditions group"
# This is where hydro thorns schedule functions that select which
# boundary conditions should be applied to which variables.
schedule group HydroBase_Select_Boundaries in HydroBase_Boundaries
{
} "Group to schedule the boundary condition functions"
# This group actually applies all selected boundary conditions. Hydro
# thorns should not need to change anything here.
schedule group ApplyBCs as HydroBase_ApplyBCs in HydroBase_Boundaries \
after HydroBase_Select_Boundaries
{
} "Apply the boundary conditions of HydroBase"
# Hydro thorns should schedule their Con2Prim routines in this group.
# This group is then scheduled at different times during an evolution.
schedule group HydroBase_Con2Prim in HydroBase_PostStep
{
} "Convert from conservative to primitive variables"
schedule group HydroBase_Con2Prim at CCTK_PostPostInitial as Con2Prim \
before ADMConstraintsGroup
{
} "Convert from conservative to primitive variables (might be redundant)"
# Hydro thorns should schedule their Prim2Con routines in this group.
# This group is then scheduled at different times during an evolution.
schedule group HydroBase_Prim2ConInitial at CCTK_INITIAL \
after HydroBase_Initial
{
} "Recover the conservative variables from the primitive variables"
if (CCTK_EQUALS(initial_hydro, "zero"))
{
schedule HydroBase_Zero in HydroBase_Initial
{
LANG: C
} "Set up zero (vacuum without atmosphere) hydro initial data"
}
if (CCTK_Equals(initial_Y_e, "one"))
{
SCHEDULE HydroBase_Y_e_one in HydroBase_Initial
{
LANG: C
} "Set electron fraction to 1"
}
|