diff options
Diffstat (limited to 'src/Bin_NS.cc')
-rw-r--r-- | src/Bin_NS.cc | 10 |
1 files changed, 9 insertions, 1 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]; |