[cig-commits] r15889 - in seismo/3D/ADJOINT_TOMO/measure_adj: . USER_MANUAL
carltape at geodynamics.org
carltape at geodynamics.org
Tue Oct 27 17:58:46 PDT 2009
Author: carltape
Date: 2009-10-27 17:58:45 -0700 (Tue, 27 Oct 2009)
New Revision: 15889
Modified:
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/measure_adj_manual.pdf
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/measure_adj_manual.tex
seismo/3D/ADJOINT_TOMO/measure_adj/mt_constants.f90
Log:
Cleaned up mt_constants.f90 and updated the manual.
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/measure_adj_manual.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/measure_adj_manual.tex
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/measure_adj_manual.tex 2009-10-27 23:17:35 UTC (rev 15888)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/measure_adj_manual.tex 2009-10-28 00:58:45 UTC (rev 15889)
@@ -207,27 +207,75 @@
%
\begin{itemize}
\item Time vector: \verb+tt+ (start time), \verb+dtt+ (time step), \verb+nn+ (number of points). These values are particularly important if adjoint sources are requested. Based on the convention in \verb+SPECFEM3D+, for example, it is best to select the start time to be $t_0 = -1.5 \tau_{\rm h}$, where $\tau_{\rm h}$ is the half duration of the CMT source. The time step and the number of time steps should exactly match the time file of the \verb+SPECFEM3D+ synthetics.
+
\item Option for measurement: \verb+imeas0+ = 1--8 (see \refSec{sec:meas}).
+
\item Channel: \verb+chan+ is BH or LH (based on options for SEM synthetics).
+
\item Bandpass range: \verb+TLONG+ (long-period), \verb+TSHORT+ short-period
+
\item \verb+RUN_BANDPASS+: whether to bandpass the seismograms.
+
\item \verb+OUTPUT_MEASUREMENT_FILES+: whether to output measurement files for each window.
+
\item Limits on cross-correlation measurement: \verb+TSHIFT_MIN+; \verb+TSHIFT_MAX+. These can optionally be read in from the \verb+PAR_FILE+ used in the FLEXWIN algorithm.
+
\item Limits on amplitude measurement: \verb+DLNA_MIN+; \verb+DLNA_MAX+. These can optionally be read in from the \verb+PAR_FILE+ used in the FLEXWIN algorithm.
+
\item Limits on the cross-correlation measurement: \verb+CC_MIN+. This can optionally be read in from the \verb+PAR_FILE+ used in the FLEXWIN algorithm.
+
\item \verb+ERROR_TYPE+ = 0--2: Type of estimated errors associated with each measurement (and adjoint source). No error (0), cross-correlation error estimate (1), multitaper jack-knife error estimate (2). The cross-correlation error estimate is presented in \citet[][Appendix~A]{Tape2010}.
+
\item \verb+DT_SIGMA_MIN+: water-level minimum cross-correlation traveltime difference uncertainty estimate.
+
\item \verb+DLNA_SIGMA_MIN+: water-level minimum amplitude difference uncertainty estimate.
-\item Multitaper variable, \verb+ITAPER+: Type of taper to use in constructing the transfer function between synthetics and data. Taper options are multitaper (1), cosine taper (2), and boxcar taper (3). For the single-taper options (2--3) the transfer function is not used, as the adjoint sources are constructed directly from the synthetic seismograms. For the multitaper option, the number of tapers is fixed to be twice \verb+NPI+ (see \verb+multitaper_vala.pdf+).
-\item Multitaper variables, \verb+WTR+, \verb+NPI+ (see \verb+multitaper_vala.pdf+).
-\item Multitaper variable, \verb+DT_FAC+, used for assessing whether to make a multitaper measurement (or revert to a cross-correlation measurement); see subroutine \verb+mt_measure_select.f90+.
-\item Multitaper variable, \verb+ERR_FAC+, used for assessing whether to make a multitaper measurement (or revert to a cross-correlation measurement); see subroutine \verb+mt_measure_select.f90+.
-\item Multitaper variable, \verb+DT_MAX_SCALE+,used for assessing whether to make a multitaper measurement (or revert to a cross-correlation measurement); see subroutine \verb+mt_measure_select.f90+.
-\item Multitaper variable, \verb+NCYCLE_IN_WINDOW+, used for assessing whether to make a multitaper measurement (or revert to a cross-correlation measurement); see subroutine \verb+mt_measure_select.f90+.
+
+\item Multitaper parameter, \verb+ITAPER+: Type of taper to use in constructing the transfer function between synthetics and data. Taper options are multitaper (1), cosine taper (2), and boxcar taper (3). For the single-taper options (2--3) the transfer function is not used, as the adjoint sources are constructed directly from the synthetic seismograms. For the multitaper option, the number of tapers is fixed to be twice \verb+NPI+ (see \verb+multitaper_vala.pdf+).
+
+\item Multitaper parameters, \verb+WTR+, \verb+NPI+ (see \verb+multitaper_vala.pdf+).
+
+\item Multitaper parameter, \verb+DT_FAC+ (see \refSec{sec:MTparm}).
+
+\item Multitaper parameter, \verb+ERR_FAC+ (see \refSec{sec:MTparm}).
+
+\item Multitaper parameter, \verb+DT_MAX_SCALE+ (see \refSec{sec:MTparm}).
+
+\item Multitaper parameter, \verb+NCYCLE_IN_WINDOW+, (see \refSec{sec:MTparm}).
\end{itemize}
%-------------
+\subsubsection{Multitaper parameters}
+\label{sec:MTparm}
+
+The selection of the multitaper parameters \verb+DT_FAC+, \verb+ERR_FAC+, \verb+DT_MAX_SCALE+, and \verb+NCYCLE_IN_WINDOW+ are not easy to determine for each particular dataset. These parameters determine whether a multitaper measurement is reasonable enough to retain (otherwise revert to a cross-correlation measurement). The key subroutine is \verb+mt_measure_select.f90+.
+
+This information used to be specified in \verb+mt_constants.f90+, and two example sets of parameters are:
+%
+\begin{verbatim}
+! ********* adjust following for global traces ************************
+ ! ratio of current period with respect to dt and err_dt measurements
+ double precision, parameter :: DT_FAC = 1.0
+ double precision, parameter :: ERR_FAC = 8.0
+ ! max time shift allowed at all freq should be DT_MAX_SCALE * Tshift_xc
+ double precision, parameter :: DT_MAX_SCALE = 5.0
+ ! use 3 cycles for surfaces waves
+ double precision, parameter :: NCYCLE_IN_WINDOW = 2.0
+!***********************************************************************
+
+! ********* adjust following for socal 3-30s traces ********************
+ ! ratio of current period with respect to dt and err_dt measurements
+ double precision, parameter :: DT_FAC = 2.0
+ double precision, parameter :: ERR_FAC = 2.5
+ ! max time shift allowed at all freq should be DT_MAX_SCALE * Tshift_xc
+ double precision, parameter :: DT_MAX_SCALE = 3.5
+ ! use 3 cycles for surfaces waves
+ double precision, parameter :: NCYCLE_IN_WINDOW = 1.5
+!***********************************************************************
+\end{verbatim}
+
+%-------------
+
\subsection{Output files}
With \verb+OUTPUT_MEASUREMENT_FILES = .true.+ in \verb+MEASUREMENT.PAR+, you should find numerous output files in \verb+OUTPUT_FILES+.
@@ -279,7 +327,7 @@
%-------------------------------------------------------------------
-\pagebreak
+%\pagebreak
\section{Measurement options}
\label{sec:meas}
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/mt_constants.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/mt_constants.f90 2009-10-27 23:17:35 UTC (rev 15888)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/mt_constants.f90 2009-10-28 00:58:45 UTC (rev 15889)
@@ -1,73 +1,26 @@
module mt_constants
-!!$! ********* adjust following for global traces ************************
-!!$ ! ratio of current period with respect to dt and err_dt measurements
-!!$ double precision, parameter :: DT_FAC = 1.0
-!!$ double precision, parameter :: ERR_FAC = 8.0
-!!$ ! max time shift allowed at all freq should be DT_MAX_SCALE * Tshift_xc
-!!$ double precision, parameter :: DT_MAX_SCALE = 5.0
-!!$ ! use 3 cycles for surfaces waves
-!!$ double precision, parameter :: NCYCLE_IN_WINDOW = 2.0
-!!$
-!!$ ! criteria to determine if a trace is used for adjiont source or not
-!!$ double precision, parameter :: BEFORE_QUALITY = 0.6
-!!$ double precision, parameter :: AFTER_QUALITY = 0.8
-!!$ double precision, parameter :: BEFORE_TSHIFT = 30.0
-!!$ double precision, parameter :: AFTER_TSHIFT = 3.0
-!!$!***********************************************************************
-
-!!$! ********* adjust following for basin traces ************************
-!!$ ! ratio of current period with respect to dt and err_dt measurements
-!!$ double precision, parameter :: DT_FAC = 2.5
-!!$! double precision, parameter :: DT_FAC = 0.25 ! big bear 2s only
-!!$ double precision, parameter :: ERR_FAC = 2.5
-!!$! double precision, parameter :: ERR_FAC = 0.5 ! big bear 2s only
-!!$ ! max time shift allowed at all freq should be DT_MAX_SCALE * Tshift_xc
-!!$ double precision, parameter :: DT_MAX_SCALE = 3.0
-!!$ ! use 3 cycles for surfaces waves
-!!$ double precision, parameter :: NCYCLE_IN_WINDOW = 3.0
-!!$
-!!$ ! criteria to determine if a trace is used for adjoint source or not
-!!$ ! --> BEFORE_QUALITY and AFTER_QUALITY should be <= (threshold used in the windowing code)
-!!$ double precision, parameter :: BEFORE_QUALITY = 0.7
-!!$ double precision, parameter :: AFTER_QUALITY = 0.7
-!!$ double precision, parameter :: BEFORE_TSHIFT = 8.0
-!!$ double precision, parameter :: AFTER_TSHIFT = 1.0
-!!$!***********************************************************************
-
-! total number of possible measuement types (waveform, MT-tt, MT-amp, XC-tt, XC-amp)
+ ! number of entries in window_chi output file
integer, parameter :: N_MEASUREMENT = 5
-
-! number of elements in window_chi
integer, parameter :: NCHI = 3*(N_MEASUREMENT-1) + 8
- integer, parameter :: IKER_WV = 0
- integer, parameter :: IKER_MT = 1
- integer, parameter :: IKER_CC = 2
-
-! constants
-
+ ! constants
double precision, parameter :: PI = 3.141592653589793d+00
double precision, parameter :: TWOPI = 2.0 * PI
- complex*16 ,parameter :: CCI = cmplx(0.,1.)
- double precision, parameter :: LARGE_VAL = 1.0d8
-! TOLERRANCE CONTROL
- double precision, parameter :: TOL2=1e-2
+ complex*16, parameter :: CCI = cmplx(0.,1.)
+ double precision, parameter :: LARGE_VAL = 1.0d8
-! FFT parameters
-
+ ! FFT parameters
integer, parameter :: LNPT = 13, NPT = 2**LNPT, NDIM = 40000
double precision, parameter :: FORWARD_FFT = 1.0
double precision, parameter :: REVERSE_FFT = -1.0
-! phase correction control parameters
-! set this between (PI, 2PI), use a higher value for
-! conservative phase wrapping
+ ! phase correction control parameters, set this between (PI, 2PI),
+ ! use a higher value for conservative phase wrapping
double precision, parameter :: PHASE_STEP = 1.5 * PI
-! filter parameters for xapiir subroutine (filter type is BP)
-! (These should EXACTLY match the filter used in the windowing code.)
-
+ ! filter parameters for xapiir bandpass subroutine (filter type is BP)
+ ! (These should match the filter used in pre-processing.)
double precision, parameter :: TRBDNDW = 0.3
double precision, parameter :: APARM = 30.
integer, parameter :: IORD = 4
More information about the CIG-COMMITS
mailing list