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

[alessia at geodynamics.org]

Author: alessia
Date: 2008-10-01 05:53:54 -0700 (Wed, 01 Oct 2008)
New Revision: 12979

Modified:
   seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90
Log:
Fixed incorrect SNR_power calculation for complete seismogram.  New version should reject more traces for a given SNR limit.  Changed variable name for SNR_window parameter to WINDOW_S2N_BASE.

Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90	2008-10-01 12:53:17 UTC (rev 12978)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90	2008-10-01 12:53:54 UTC (rev 12979)
@@ -41,7 +41,7 @@
     CC_LIMIT(i)=CC_BASE
     TSHIFT_LIMIT(i)=TSHIFT_BASE
     STALTA_W_LEVEL(i)=STALTA_BASE
-    S2N_LIMIT(i)=WINDOW_AMP_BASE
+    S2N_LIMIT(i)=WINDOW_S2N_BASE
   enddo
 
   ! these values will be used for signal2noise calculations
@@ -74,7 +74,7 @@
 ! do i = 1, npts
 !   time=b+(i-1)*dt
 !   if (time.gt.R_time) then
-!     S2N_LIMIT(i)=2*WINDOW_AMP_BASE
+!     S2N_LIMIT(i)=2*WINDOW_S2N_BASE
 !   endif
 ! enddo
 !
@@ -101,7 +101,7 @@
    ! if we are beyond the Rayleigh wave, then make the all criteria stronger
    ! ratio criterion stronger
    if (time.gt.R_time) then
-     S2N_LIMIT(i)=10*WINDOW_AMP_BASE    ! only pick big signals
+     S2N_LIMIT(i)=10*WINDOW_S2N_BASE    ! only pick big signals
      CC_LIMIT(i)= 0.95                  ! only pick very similar signals
      TSHIFT_LIMIT(i)=TSHIFT_BASE/3.0    ! only pick small timeshifts
      DLNA_LIMIT(i)=DLNA_BASE/3.0        ! only pick small amplitude anomalies