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

[alessia at geodynamics.org]

Author: alessia
Date: 2008-09-24 00:39:00 -0700 (Wed, 24 Sep 2008)
New Revision: 12942

Modified:
   seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran
   seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90
   seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90
Log:
Added modify_T0_T1_on_condition subroutine that permits modification of frequency limits on user-given conditions in the user_functions file.  Useful for scenarios with heterogeneous instrumentation or a large range of earthquake magnitudes.  The stadard implementation of the subroutine does nothing.  NOTE: YOU MUST UPDATE YOUR USER_FUNCTIONS FILE TO CONTAIN AT LEAST THE DO-NOTHING VERSION OF THIS ROUTINE.

Modified: seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran	2008-09-24 03:05:57 UTC (rev 12941)
+++ seismo/3D/ADJOINT_TOMO/flexwin/make_gfortran	2008-09-24 07:39:00 UTC (rev 12942)
@@ -37,5 +37,8 @@
 
 clean:
 	rm -f  ${LIB2} ${PROGS} *.o *.mod
+
+cleanall:
 	rm -f  iasp91.*
+	cd $TAULIBDUR ; make -f make_gfortran clean
 

Modified: seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90	2008-09-24 03:05:57 UTC (rev 12941)
+++ seismo/3D/ADJOINT_TOMO/flexwin/seismo_subs.f90	2008-09-24 07:39:00 UTC (rev 12942)
@@ -341,6 +341,10 @@
   ! calculate distances and azimuths
   call distaz(evla,evlo,stla,stlo,azimuth,backazimuth,dist_deg,dist_km)
 
+  ! Frequency limits may be conditional on station or event information
+  ! so call user function to modify them if required
+  call modify_T0_T1_on_condition
+
   if(RUN_BANDPASS) then
      ! clean up the seismograms
      call detrend(obs,npts)

Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90	2008-09-24 03:05:57 UTC (rev 12941)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_functions.f90	2008-09-24 07:39:00 UTC (rev 12942)
@@ -1,9 +1,28 @@
-  ! -------------------------------------------------------------
-  ! edit here to change the time dependent properties of the 
-  ! selection criteria
-  ! Note, this function is called AFTER the seismogram has been 
-  ! read.
-  ! -------------------------------------------------------------
+
+! -------------------------------------------------------------
+! edit here to change T0 and T1 on some condition 
+! Note, this function is called AFTER the seismogram has been 
+! read but before it is filtered.
+! -------------------------------------------------------------
+
+  subroutine modify_T0_T1_on_condition
+  use seismo_variables
+
+  ! do nothing
+
+  ! adjust fstart and fend accordingly
+  FSTART=1./WIN_MAX_PERIOD
+  FEND=1./WIN_MIN_PERIOD
+
+  end subroutine
+
+! -------------------------------------------------------------
+! edit here to change the time dependent properties of the 
+! selection criteria
+! Note, this function is called AFTER the seismogram has been 
+! read and filtered.
+! -------------------------------------------------------------
+
   subroutine set_up_criteria_arrays
   use seismo_variables