updates coupling_points array format

Author: danielpeter

src/specfem3D/couple_with_injection.f90

@@ -2259,7 +2259,7 @@ subroutine locate_coupling_points(num_coupling_points_total,coupling_points)
   implicit none
 
   integer, intent(in) :: num_coupling_points_total
-  real(kind=CUSTOM_REAL),dimension(7,num_coupling_points_total),intent(in) :: coupling_points
+  real(kind=CUSTOM_REAL),dimension(6,num_coupling_points_total),intent(in) :: coupling_points
 
   ! local parameters
   integer :: ipoin,islice,ispec,ier
@@ -2573,7 +2573,7 @@ subroutine determine_local_coupling_points(num_coupling_points_total,coupling_po
   implicit none
 
   integer, intent(in) :: num_coupling_points_total
-  real(kind=CUSTOM_REAL),dimension(7,num_coupling_points_total),intent(in) :: coupling_points
+  real(kind=CUSTOM_REAL),dimension(6,num_coupling_points_total),intent(in) :: coupling_points
 
   double precision, dimension(num_coupling_points_total), intent(in) :: xi_point,eta_point,gamma_point
   double precision, dimension(num_coupling_points_total), intent(in) :: x_found,y_found,z_found

src/specfem3D/setup_sources_receivers.f90

@@ -2212,13 +2212,13 @@ subroutine setup_coupling_boundary_points()
     call flush_IMAIN()
 
     ! allocate points array
-    allocate(coupling_points(7,num_coupling_points_total),stat=ier)
+    allocate(coupling_points(6,num_coupling_points_total),stat=ier)
     if (ier /= 0) stop 'Error allocating coupling points array'
     coupling_points(:,:) = 0.0_CUSTOM_REAL
 
     ! reads in point infos
     do ipoin = 1,num_coupling_points_total
-      ! format: #x #y #z #nx #ny #nz #iproc
+      ! format: #x #y #z #nx #ny #nz
       read(IIN) coupling_points(:,ipoin)
     enddo
 
@@ -2234,13 +2234,13 @@ subroutine setup_coupling_boundary_points()
 
   ! allocate arrays on other processes
   if (.not. allocated(coupling_points)) then
-    allocate(coupling_points(7,num_coupling_points_total),stat=ier)
+    allocate(coupling_points(6,num_coupling_points_total),stat=ier)
     if (ier /= 0) stop 'Error allocating coupling points array'
     coupling_points(:,:) = 0.0_CUSTOM_REAL
   endif
 
   ! broadcast point infos
-  call bcast_all_cr(coupling_points,7*num_coupling_points_total)
+  call bcast_all_cr(coupling_points,6*num_coupling_points_total)
 
   ! locate coupling points
   call locate_coupling_points(num_coupling_points_total,coupling_points)