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

[alessia at geodynamics.org]

Author: alessia
Date: 2008-10-01 00:49:48 -0700 (Wed, 01 Oct 2008)
New Revision: 12977

Modified:
   seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90
Log:
Improved calculation of early STA:LTA

Modified: seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90	2008-10-01 03:38:35 UTC (rev 12976)
+++ seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90	2008-10-01 07:49:48 UTC (rev 12977)
@@ -538,17 +538,20 @@
   TOL=1e-9
 
 ! set pre-extension for synthetic data and allocate extended_syn
-  n_extend=5*12*WIN_MIN_PERIOD/dt
+!  n_extend=5*12*WIN_MIN_PERIOD/dt
+  n_extend=1000*WIN_MIN_PERIOD/dt
   allocate(extended_syn(npts+n_extend))
 
 ! set noise level
-  n_sample=12*WIN_MIN_PERIOD/dt
-  noise=sum(env_synt_lp(1:n_sample))/n_sample
+!  n_sample=12*WIN_MIN_PERIOD/dt
+!  noise=sum(env_synt_lp(1:n_sample))/n_sample
+  noise=maxval(env_synt_lp)/10**5
 
 ! copy the original synthetic into the extended array, right justified  
-  call random_number(extended_syn)
-  extended_syn=noise*extended_syn
-  extended_syn(n_extend+1:n_extend+npts)=env_synt_lp(1:npts)
+!  call random_number(extended_syn)
+!  extended_syn=noise*extended_syn
+  extended_syn=noise
+  extended_syn(n_extend+1:n_extend+npts)=env_synt_lp(1:npts)+noise
 
   if (DEBUG) write(*,*) 'DEBUG : Cs, Cl = ', Cs, Cl
   if (DEBUG) write(*,*) 'Number of points used to pre-extend synthetics ', n_extend
@@ -795,7 +798,8 @@
 
   call prepare_sta_lta
 
-! set up the selection criteria arrays
+! set up the selection criteria arrays (subroutine is in user_functions.f90
+! file)
   call set_up_criteria_arrays
 
   end subroutine