1 files changed, 46 insertions, 0 deletions
diff --git a/src/qlm_set_coordinates.F90 b/src/qlm_set_coordinates.F90
new file mode 100644
index 0000000..309a361
--- /dev/null
+++ b/src/qlm_set_coordinates.F90
@@ -0,0 +1,46 @@
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+
+
+subroutine qlm_set_coordinates (CCTK_ARGUMENTS, hn)
+  use cctk
+  implicit none
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_FUNCTIONS
+  DECLARE_CCTK_PARAMETERS
+  integer :: hn
+  
+  CCTK_REAL :: theta, phi
+  CCTK_REAL :: sin_theta, cos_theta, sin_phi, cos_phi
+  
+  integer :: i, j
+  
+  ! Calculate the coordinates
+  do j = 1, qlm_nphi(hn)
+     do i = 1, qlm_ntheta(hn)
+        theta = qlm_origin_theta(hn) + (i-1)*qlm_delta_theta(hn)
+        phi   = qlm_origin_phi(hn)   + (j-1)*qlm_delta_phi(hn)
+        
+        sin_theta = sin(theta)
+        cos_theta = cos(theta)
+        sin_phi   = sin(phi)
+        cos_phi   = cos(phi)
+        
+        qlm_x(i,j,hn) = qlm_origin_x(hn) + qlm_shape(i,j,hn) * sin_theta * cos_phi
+        qlm_y(i,j,hn) = qlm_origin_y(hn) + qlm_shape(i,j,hn) * sin_theta * sin_phi
+        qlm_z(i,j,hn) = qlm_origin_z(hn) + qlm_shape(i,j,hn) * cos_theta
+        
+        qlm_x_p(i,j,hn) = qlm_origin_x_p(hn) + qlm_shape_p(i,j,hn) * sin_theta * cos_phi
+        qlm_y_p(i,j,hn) = qlm_origin_y_p(hn) + qlm_shape_p(i,j,hn) * sin_theta * sin_phi
+        qlm_z_p(i,j,hn) = qlm_origin_z_p(hn) + qlm_shape_p(i,j,hn) * cos_theta
+        
+        qlm_x_p_p(i,j,hn) = qlm_origin_x_p_p(hn) + qlm_shape_p_p(i,j,hn) * sin_theta * cos_phi
+        qlm_y_p_p(i,j,hn) = qlm_origin_y_p_p(hn) + qlm_shape_p_p(i,j,hn) * sin_theta * sin_phi
+        qlm_z_p_p(i,j,hn) = qlm_origin_z_p_p(hn) + qlm_shape_p_p(i,j,hn) * cos_theta
+     end do
+  end do
+  
+end subroutine qlm_set_coordinates