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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-09-30 05:27:11 -0700 (Mon, 30 Sep 2013)
New Revision: 22903

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90
Log:
fix the value of NSTEP when using UNDO_ATTENUATION in setup_sources_receivers.f90


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90	2013-09-30 05:26:15 UTC (rev 22902)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90	2013-09-30 12:27:11 UTC (rev 22903)
@@ -289,6 +289,12 @@
   ! time steps needs to be doubled, due to +/- branches
   if ( NOISE_TOMOGRAPHY /= 0 )   NSTEP = 2*NSTEP-1
 
+!! DK DK make sure NSTEP is a multiple of NT_DUMP_ATTENUATION
+  if(UNDO_ATTENUATION .and. mod(NSTEP,NT_DUMP_ATTENUATION) /= 0) then
+    NSTEP = (NSTEP/NT_DUMP_ATTENUATION + 1)*NT_DUMP_ATTENUATION
+  endif
+  it_end = NSTEP
+
   ! subsets used to save seismograms must not be larger than the whole time series,
   ! otherwise we waste memory
   if(NTSTEP_BETWEEN_OUTPUT_SEISMOS > NSTEP .or. is_initial_guess) NTSTEP_BETWEEN_OUTPUT_SEISMOS = NSTEP