[cig-commits] r22965 - seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-10-19 15:12:33 -0700 (Sat, 19 Oct 2013)
New Revision: 22965

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90
Log:
fixed the test about the total number of receivers found: it did not take into account the fact that some stations located outside of the mesh can be excluded when the mesh consists of a single chunk


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90	2013-10-19 01:10:22 UTC (rev 22964)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90	2013-10-19 22:12:33 UTC (rev 22965)
@@ -445,20 +445,16 @@
   if(myrank == 0) then
     write(IMAIN,*)
     write(IMAIN,*) 'found a total of ',nrec_tot_found,' receivers in all slices'
-    if(nrec_tot_found /= nrec_simulation) then
-      call exit_MPI(myrank,'problem when dispatching the receivers')
-    else
-      write(IMAIN,*) 'this total is okay'
+    if(nrec_tot_found <= 1) then
+      call exit_MPI(myrank,'error: no receiver found in the mesh')
+!! DK DK we do not check if the total found is equal to nrec_simulation here because some stations in the STATIONS file
+!! DK DK can be located outside of the mesh, thus the only thing we can check is that the total found is smaller or equal
+    else if(nrec_tot_found > nrec_simulation) then
+      call exit_MPI(myrank,'error: too many receivers found, problem when dispatching the receivers')
     endif
     call flush_IMAIN()
   endif
 
-  ! check that the sum of the number of receivers in each slice is nrec
-  if( SIMULATION_TYPE == 1 .or. SIMULATION_TYPE == 3 ) then
-    if(myrank == 0 .and. nrec_tot_found /= nrec) &
-      call exit_MPI(myrank,'total number of receivers is incorrect')
-  endif
-
   end subroutine setup_receivers
 
 !