diff options
Diffstat (limited to 'src')
-rw-r--r-- | src/Bin_NS.cc | 10 | ||||
-rw-r--r-- | src/check_parameters.cc | 12 |
2 files changed, 16 insertions, 6 deletions
diff --git a/src/Bin_NS.cc b/src/Bin_NS.cc index 4e2696f..c04bf66 100644 --- a/src/Bin_NS.cc +++ b/src/Bin_NS.cc @@ -1,5 +1,5 @@ -#include <cassert> #include <cstdio> +#include <cassert> #include <vector> #include <cctk.h> @@ -94,6 +94,14 @@ void ID_Bin_NS_initialise (CCTK_ARGUMENTS) betay[i] = bin_ns.beta_y[i]; betaz[i] = bin_ns.beta_z[i]; + // These initial data assume a helical Killing vector field + // TODO: calculate spatial derivatives to set this correctly + dtalp[i] = 0.0; + + dtbetax[i] = 0.0; + dtbetay[i] = 0.0; + dtbetaz[i] = 0.0; + CCTK_REAL g[3][3]; g[0][0] = bin_ns.g_xx[i]; g[0][1] = bin_ns.g_xy[i]; diff --git a/src/check_parameters.cc b/src/check_parameters.cc index bef3136..40a760f 100644 --- a/src/check_parameters.cc +++ b/src/check_parameters.cc @@ -10,11 +10,13 @@ void ID_Bin_NS_check_parameters (CCTK_ARGUMENTS) DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; - if (not CCTK_EQUALS (initial_data , "ID_Bin_NS") or - not CCTK_EQUALS (initial_lapse, "ID_Bin_NS") or - not CCTK_EQUALS (initial_shift, "ID_Bin_NS") or - not CCTK_EQUALS (initial_hydro, "ID_Bin_NS")) + if (not CCTK_EQUALS (initial_data, "ID_Bin_NS") or + not CCTK_EQUALS (initial_lapse, "ID_Bin_NS") or + not CCTK_EQUALS (initial_shift, "ID_Bin_NS") or + not CCTK_EQUALS (initial_dtlapse, "ID_Bin_NS") or + not CCTK_EQUALS (initial_dtshift, "ID_Bin_NS") or + not CCTK_EQUALS (initial_hydro, "ID_Bin_NS")) { - CCTK_PARAMWARN ("The parameters ADMBase::initial_data, ADMBase::initial_lapse, ADMBase::initial_shift, and HydroBase::initial_hydro must all be set to the value \"ID_Bin_NS\""); + CCTK_PARAMWARN ("The parameters ADMBase::initial_data, ADMBase::initial_lapse, ADMBase::initial_shift, ADMBase::initial_dtlapse, ADMBase::initial_dtshift, and HydroBase::initial_hydro must all be set to the value \"ID_Bin_NS\""); } } |