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
|
/*@@
@file GRHydro_ParamCheck.F90
@date Sat Feb 9 23:48:01 2002
@author
@desc
Parameter checking routine.
@enddesc
@@*/
#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"
#include "cctk_Functions.h"
/*@@
@routine GRHydro_ParamCheck
@date Sat Feb 9 23:48:43 2002
@author Ian Hawke
@desc
Checks the parameters.
@enddesc
@calls
@calledby
@history
@endhistory
@@*/
subroutine GRHydro_ParamCheck(CCTK_ARGUMENTS)
implicit none
DECLARE_CCTK_ARGUMENTS
DECLARE_CCTK_PARAMETERS
DECLARE_CCTK_FUNCTIONS
integer :: coordinates_is_active
if (GRHydro_stencil > minval(cctk_nghostzones)) then
call CCTK_PARAMWARN("The stencil is larger than the number of ghost zones. Answer will be dependent on processor number...")
end if
if (CCTK_EQUALS(recon_method,"tvd").and.(GRHydro_stencil < 2)) then
call CCTK_PARAMWARN("The stencil size must be at least 2 to use TVD reconstruction")
end if
if (CCTK_EQUALS(recon_method,"ppm")) then
if (CCTK_EQUALS(ppm_flatten,"stencil_3").and.(GRHydro_stencil < 3)) then
call CCTK_PARAMWARN("The stencil size must be at least 3 to use PPM reconstruction with the stencil-3 variant of the flattening procedure")
else if (CCTK_EQUALS(ppm_flatten,"stencil_4").and.(GRHydro_stencil < 4)) then
call CCTK_PARAMWARN("The stencil size must be at least 4 to use PPM reconstruction with the stencil-4 (original) flattening procedure")
end if
end if
if (CCTK_EQUALS(recon_method,"eno").and.(GRHydro_stencil < eno_order)) then
call CCTK_PARAMWARN("The stencil size must be at least the order of the reconstruction to use ENO reconstruction")
end if
if (CCTK_EQUALS(GRHydro_eos_table,"2D_Polytrope").and.&
(.not.CCTK_EQUALS(GRHydro_eos_type,"Polytype"))) then
call CCTK_PARAMWARN("When using the 2D_Polytrope EOS you need to set eos_type to Polytype")
end if
if (CCTK_EQUALS(GRHydro_eos_table,"Ideal_Fluid").and.&
(.not.CCTK_EQUALS(GRHydro_eos_type,"General"))) then
call CCTK_PARAMWARN("When using the Ideal_Fluid EOS you need to set eos_type to General")
end if
if (.not.(CCTK_EQUALS(metric_type, "Physical").or.&
CCTK_EQUALS(metric_type, "Static Conformal"))) then
call CCTK_WARN(0,"GRHydro only knows how to deal with physical metric type at the moment. Complain to the maintainers...")
end if
if (use_mask .eq. 0) then
call CCTK_PARAMWARN("GRHydro now requires you to set SpaceMask::use_mask = ""yes""")
end if
if (number_of_particles .gt. 0) then
if (number_of_arrays .ne. 3) then
call CCTK_PARAMWARN("If tracking particles then you must have number_of_arrays = 3")
end if
end if
! This check and 'static' as valid keyword option can be removed after the
! year 2010
if (CCTK_EQUALS(bound,"static")) then
call CCTK_PARAMWARN("GRHydro::bound = 'static' is no longer supported, use 'none' instead");
end if
if (CCTK_EQUALS(Bvec_evolution_method,"GRHydro")) then
evolve_MHD = 1
else
evolve_MHD = 0
endif
if (CCTK_EQUALS(Y_e_evolution_method,"GRHydro")) then
evolve_Y_e = 1
else
evolve_Y_e = 0
endif
if (CCTK_EQUALS(temperature_evolution_method,"GRHydro")) then
evolve_temper = 1
else
evolve_temper = 0
endif
call CCTK_IsImplementationActive(coordinates_is_active, "Coordinates")
! this test is somewhat overzealous but we cannot access
! coordinates::general_coordinates yet since it is only set in BaseGrid
if (CCTK_Equals(riemann_solver,"HLLE").eq.0.and.coordinates_is_active.ne.0) then
call CCTK_WARN(0, "There is currently no Riemann solver other than HLLE compatible with multipatch!")
end if
if (CCTK_EQUALS(riemann_solver,"Roe").and.CCTK_EQUALS(Bvec_evolution_method,"GRHydro")) then
call CCTK_PARAMWARN("Roe solver is not implemented yet for MHD")
end if
if (CCTK_EQUALS(riemann_solver,"Marquina").and.CCTK_EQUALS(Bvec_evolution_method,"GRHydro")) then
call CCTK_PARAMWARN("Marquina solver is not implemented yet for MHD")
end if
if(clean_divergence.ne.0) then
if (CCTK_EQUALS(psidcspeed,"char speed")) then
whichpsidcspeed = 0
else if (CCTK_EQUALS(psidcspeed,"light speed")) then
whichpsidcspeed = 1
else if (CCTK_EQUALS(psidcspeed,"set speed")) then
whichpsidcspeed = 2
else
call CCTK_PARAMWARN("Unknown type of speed setting for psidc (psidcspeed)")
end if
end if
if (Tmunu_damping_radius_min .gt. Tmunu_damping_radius_max) then
call CCTK_WARN(0, "Minimum damping radius is greater than maximum damping radius!")
end if
end subroutine GRHydro_ParamCheck
|
