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
|
#include "cctk.h"
#include "cctk_Arguments.h"
#include "cctk_Parameters.h"
void GRHydro_InitData_CheckParameters(CCTK_ARGUMENTS)
{
DECLARE_CCTK_ARGUMENTS
DECLARE_CCTK_PARAMETERS
if (timelevels < 2)
{
CCTK_PARAMWARN("You have to set 'HydroBase::timelevels to at least 2");
}
if(CCTK_Equals(Bvec_evolution_method,"GRHydro") &&
((CCTK_Equals(initial_hydro,"ony_atmo")) ||
(CCTK_Equals(initial_hydro,"read_conformal")) ||
(CCTK_Equals(initial_hydro,"simple_wave")) ||
(CCTK_Equals(initial_data,"con2primtest")) ||
(CCTK_Equals(initial_data,"reconstruction_test")) ||
(CCTK_Equals(shocktube_type,"sphere"))))
{
CCTK_PARAMWARN("That test not yet implemented in MHD!");
}
if(!CCTK_Equals(Bvec_evolution_method,"GRHydro") &&
((CCTK_Equals(shock_case,"Balsara0" )) ||
(CCTK_Equals(shock_case,"Balsara1" )) ||
(CCTK_Equals(shock_case,"Balsara2" )) ||
(CCTK_Equals(shock_case,"Balsara3" )) ||
(CCTK_Equals(shock_case,"Balsara4" )) ||
(CCTK_Equals(shock_case,"Balsara5" )) ||
(CCTK_Equals(shock_case,"Alfven" )) ||
(CCTK_Equals(shock_case,"Komissarov1")) ||
(CCTK_Equals(shock_case,"Komissarov2")) ||
(CCTK_Equals(shock_case,"Komissarov3")) ||
(CCTK_Equals(shock_case,"Komissarov4")) ||
(CCTK_Equals(shock_case,"Komissarov5")) ||
(CCTK_Equals(shock_case,"Komissarov6")) ||
(CCTK_Equals(shock_case,"Komissarov7")) ||
(CCTK_Equals(shock_case,"Komissarov8")) ||
(CCTK_Equals(shock_case,"Komissarov9")) ||
(CCTK_Equals(initial_hydro,"cylexp"))
))
{
CCTK_PARAMWARN("That test requires MHD! Set Bvec_evolution_method=GRHYDRO!");
}
/* Checks for Bondi solution initial data : */
if(CCTK_Equals(Bvec_evolution_method,"GRHydro") &&
CCTK_Equals(initial_hydro,"hydro_bondi_solution") )
{
CCTK_PARAMWARN("Please use set initial_hydro=magnetized_bondi_solution instead to initialize the magnetic field correctly!");
}
if(!CCTK_Equals(Bvec_evolution_method,"GRHydro") &&
CCTK_Equals(initial_hydro,"magnetized_bondi_solution") )
{
CCTK_PARAMWARN("Please use set initial_hydro=hydro_bondi_solution instead to NOT initialize the magnetic field!");
}
if( CCTK_Equals(initial_hydro,"magnetized_bondi_solution") || CCTK_Equals(initial_hydro,"hydro_bondi_solution") ) {
if( num_bondi_sols > 1 ) {
CCTK_PARAMWARN("Currently only one Bondi solution is supported, please change [num_bondi_sols] ");
}
}
}
|