[cig-commits] r18972 - seismo/3D/ADJOINT_TOMO/flexwin

[carltape at geodynamics.org]

Author: carltape
Date: 2011-09-26 14:36:50 -0700 (Mon, 26 Sep 2011)
New Revision: 18972

Modified:
   seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90
Log:
important change in indexing for the signal part in a signal-to-noise calculation; old version was is1 = j, new version is is1 = j-1 (index corresponding to signal_start). this was causing some seismograms to exit without picking any windows


Modified: seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90	2011-09-26 20:19:26 UTC (rev 18971)
+++ seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90	2011-09-26 21:36:50 UTC (rev 18972)
@@ -312,12 +312,13 @@
        time_obs_signal = b+(j-1)*dt
        if (time_obs_signal > signal_start) exit
     enddo
-    is1 = j-1     ! index corresponding to signal_start
+    is1 = j     ! index corresponding to signal_start
     do k = 1, npts
        time_obs_signal = b+(k-1)*dt
        if (time_obs_signal > signal_end) exit
     enddo
     is2 = k-1     ! index corresponding to signal_end
+    print *, 'is1 = ', is1, ', is2 = ', is2
     signal_int = sum( (obs_lp(is1:is2))**2 )/(is2-is1)
     signal_amp = maxval( abs(obs_lp(is1:is2)) )