aboutsummaryrefslogtreecommitdiff
path: root/src/qlm_killing_axial.F90
diff options
context:
space:
mode:
Diffstat (limited to 'src/qlm_killing_axial.F90')
-rw-r--r--src/qlm_killing_axial.F9035
1 files changed, 35 insertions, 0 deletions
diff --git a/src/qlm_killing_axial.F90 b/src/qlm_killing_axial.F90
new file mode 100644
index 0000000..feda999
--- /dev/null
+++ b/src/qlm_killing_axial.F90
@@ -0,0 +1,35 @@
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+
+
+subroutine qlm_killing_axial (CCTK_ARGUMENTS, hn)
+ use cctk
+ implicit none
+ DECLARE_CCTK_ARGUMENTS
+ DECLARE_CCTK_FUNCTIONS
+ DECLARE_CCTK_PARAMETERS
+ integer :: hn
+
+ CCTK_REAL :: xi(2), chi
+ integer :: i, j
+
+ do j = 1, qlm_nphi(hn)
+ do i = 1, qlm_ntheta(hn)
+
+ xi(1) = 0
+ xi(2) = 1
+
+ chi = 0
+
+ qlm_xi_t(i,j,hn) = xi(1)
+ qlm_xi_p(i,j,hn) = xi(2)
+
+ qlm_chi(i,j,hn) = chi
+
+ end do
+ end do
+
+end subroutine qlm_killing_axial