[cig-commits] r15866 - in seismo/3D/ADJOINT_TOMO/measure_adj: . PLOTS USER_MANUAL USER_MANUAL/figures scripts_tomo
carltape at geodynamics.org
carltape at geodynamics.org
Wed Oct 21 20:45:36 PDT 2009
Author: carltape
Date: 2009-10-21 20:45:35 -0700 (Wed, 21 Oct 2009)
New Revision: 15866
Added:
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker01_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker02_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker03_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker04_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker05_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker06_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker07_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker08_win_adj.eps
seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_win.eps
Removed:
seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/9818433_T006_T030_MPM_CI_m16_ikerXX_win_adj.pdf
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_sub.f90
seismo/3D/ADJOINT_TOMO/measure_adj/run_mt_measure_adj.csh
seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_mt_measure_adj.pl
seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl
seismo/3D/ADJOINT_TOMO/measure_adj/write_par_file.pl
Log:
Updates to manual and to run scripts.
Deleted: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/9818433_T006_T030_MPM_CI_m16_ikerXX_win_adj.pdf
===================================================================
(Binary files differ)
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker01_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker01_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker02_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker02_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker03_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker03_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker04_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker04_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker05_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker05_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker06_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker06_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker07_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker07_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker08_win_adj.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_iker08_win_adj.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_win.eps
===================================================================
(Binary files differ)
Property changes on: seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/figures/9818433_T006_T030_MPM_CI_m16_win.eps
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
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-22 00:43:10 UTC (rev 15865)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/USER_MANUAL/measure_adj_manual.tex 2009-10-22 03:45:35 UTC (rev 15866)
@@ -142,19 +142,103 @@
\section{Running MEAS--ADJ}
-\subsection{Input files}
+%-------------
+\subsection{Input file 1: {\tt MEASUREMENT.WINDOWS}}
+
+\begin{verbatim}
+3
+DATA/9818433.CI.MPM.BHR.sac.d.T006_T030
+SYN/MPM.CI.BHR.semd.sac.m16.T006_T030
+ 3
+ 11.2770 58.4770
+ 58.4770 81.3770
+ 81.3770 101.6770
+DATA/9818433.CI.MPM.BHT.sac.d.T006_T030
+SYN/MPM.CI.BHT.semd.sac.m16.T006_T030
+ 1
+ 49.3270 95.3270
+DATA/9818433.CI.MPM.BHZ.sac.d.T006_T030
+SYN/MPM.CI.BHZ.semd.sac.m16.T006_T030
+ 1
+ 60.9770 99.3770
+\end{verbatim}
+
+
+
+%-------------
+
+\pagebreak
+\subsection{Input file 2: {\tt MEASUREMENT.PAR}}
+
+\begin{verbatim}
+ -0.585 0.0110 18200 # tstart, DT, npts: time vector for simulations
+ 7 # imeas (1-8; see manual)
+ BH # channel: BH or LH
+ 30.000 6.000 # TLONG and TSHORT: band-pass periods for records
+ .false. # RUN_BANDPASS: use band-pass on records
+ .true. # DISPLAY_DETAILS
+ .true. # OUTPUT_MEASUREMENT_FILES
+ .true. # COMPUTE_ADJOINT_SOURCE
+ -4.5000 4.5000 # TSHIFT_MIN; TSHIFT_MAX
+ -1.5000 1.5000 # DLNA_MIN; DLNA_MAX
+ 0.690 # CC_MIN
+ 1 # ERROR_TYPE -- 0 none; 1 CC, MT-CC; 2 MT-jack-knife
+ 1.000 # DT_SIGMA_MIN
+ 0.500 # DLNA_SIGMA_MIN
+ 1 # itaper -- taper type: 1 multi-taper; 2 cosine; 3 boxcar
+ 0.020 2.50 # wtr, npi (ntaper = 2*npi)
+ 2.000 # DT_FAC
+ 2.500 # ERR_FAC
+ 3.500 # DT_MAX_SCALE
+ 1.500 # NCYCLE_IN_WINDOW
+\end{verbatim}
+
+\begin{itemize}
+\item Time vector: $t_0$ (start time), $dt$ (time step), and $n$ (number of points).
+\item Option for measurement: 1--8 (see \refSec{sec:meas}).
+\item Channel: BH or LH (based on options for SEM synthetics).
+\item Bandpass range: long-period, short-period
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\end{itemize}
+
+%-------------
+
\subsection{Output files}
+
+
+%-------------
+
\subsection{Scripts}
%-------------------------------------------------------------------
\section{Measurement options}
+\label{sec:meas}
\begin{enumerate}
-\item No adjoint source\verb+iker = 0+
+\item
+\item
+\item
+\item
+\item
+\item
+\item
+\item
\end{enumerate}
%-------------------------------------------------------------------
@@ -198,16 +282,96 @@
%\normalsize
%-------------------------------------------------------------------
-%\begin{figure}
-%\hspace{-1cm}
-%\includegraphics[width=18cm]{Hollywood_CMT.eps}
-%\caption[]
-%{{
-%CMT example figure for the {\tt SPECFEM3D$\_$BASIN} manual.
-%\label{fig:CMT_basin}
-%}}
-%\end{figure}
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_win.eps}
+\caption[]
+{{
+Test case seismograms with {\tt COMPUTE$\_$ADJOINT$\_$SOURCE = .false.}. Plot shows a three-component seismogram: Z = vertical, R = radial, T = transverse. Black is the observed records, red is the synthetics, and blue is the reconstructed synthetics, made by applying the $\Delta T$ and $\Delta \ln A$ measurements within each window to the synthetics. No adjoint sources are plotted.
+\label{fig:noadj}
+}}
+\end{figure}
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker01_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources constructed from the data alone ({\tt imeas = 1}), $-d(t)/\|d(t)\|^2$.
+\label{fig:iker01}
+}}
+\end{figure}
+
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker02_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources for a waveform difference measurement ({\tt imeas = 2}), $\bs - \bd$.
+\label{fig:iker02}
+}}
+\end{figure}
+
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker03_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources for a sensitivity kernel based on a cross-correlation traveltime difference, $\Delta T$ ({\tt imeas = 3}).
+\label{fig:iker03}
+}}
+\end{figure}
+
+
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker04_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources for a sensitivity kernel based on an amplitude difference, $\Delta \ln A$ ({\tt imeas = 4}).
+\label{fig:iker04}
+}}
+\end{figure}
+
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker05_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources for a cross-correlation traveltime difference, $\Delta T$ ({\tt imeas = 5}), .
+\label{fig:iker05}
+}}
+\end{figure}
+
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker06_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources for an amplitude difference, $\Delta \ln A$ ({\tt imeas = 6}), .
+\label{fig:iker06}
+}}
+\end{figure}
+
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker07_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources for a multitaper traveltime difference ({\tt imeas = 7}), .
+\label{fig:iker07}
+}}
+\end{figure}
+
+\begin{figure}
+\includegraphics[width=17cm]{9818433_T006_T030_MPM_CI_m16_iker08_win_adj.eps}
+\caption[]
+{{
+{\em Left:} Data (black), synthetics (red), and measurement windows.
+{\em Right:} Adjoint sources for a multitaper amplitude difference ({\tt imeas = 8}), .
+\label{fig:iker08}
+}}
+\end{figure}
+
%-------------------------------------------------------------------
\end{document}
%-------------------------------------------------------------------
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90 2009-10-22 00:43:10 UTC (rev 15865)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90 2009-10-22 03:45:35 UTC (rev 15866)
@@ -733,14 +733,13 @@
fa_bar_t = 0.
Mnorm = dt * sum( syn_dtw(1:nlen) * syn_dtw(1:nlen) )
fa_bar_t(1:nlen) = syn_dtw(1:nlen) / Mnorm
-
endif
! -------------------------------
! MULTITAPER ADJOINT SOURCES
! -------------------------------
- if (is_mtm == 1) then
+ if ( (is_mtm == 1).and.COMPUTE_ADJOINT_SOURCE ) then
! allocate MT variables
ntaper = int(NPI * 2.0)
@@ -834,9 +833,7 @@
if (COMPUTE_ADJOINT_SOURCE) then
- ! compute traveltime and amplitude adjoint sources
do i = 1,nlen
-
i1 = istart + i
! waveform
@@ -845,21 +842,20 @@
am_adj_src(i1) = ( syn_dtw(i) - dat_dtw(i) ) * time_window(i)
! consider normalizing this by waveform_d2
- ! banana-doughnut kernel adjoint source (no measurement)
+ ! banana-doughnut kernel adjoint source (no measurement)
elseif(imeas==3 .or. imeas==4) then
tr_adj_src(i1) = ft_bar_t(i) * time_window(i)
am_adj_src(i1) = fa_bar_t(i) * time_window(i)
- ! cross-correlation
+ ! cross-correlation
elseif(imeas==5 .or. imeas==6) then
tr_adj_src(i1) = -(tshift / sigma_dt_cc**2 ) * ft_bar_t(i) * time_window(i)
am_adj_src(i1) = -(dlnA / sigma_dlnA_cc**2 ) * fa_bar_t(i) * time_window(i)
- ! multitaper
+ ! multitaper
elseif(imeas==7 .or. imeas==8) then
tr_adj_src(i1) = fp(i) * time_window(i)
am_adj_src(i1) = fq(i) * time_window(i)
-
endif
enddo
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/run_mt_measure_adj.csh
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/run_mt_measure_adj.csh 2009-10-22 00:43:10 UTC (rev 15865)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/run_mt_measure_adj.csh 2009-10-22 03:45:35 UTC (rev 15866)
@@ -15,5 +15,5 @@
cp STATIONS_ADJOINT ADJOINT_SOURCES
\mv STATIONS_ADJOINT PLOTS
cd PLOTS
-plot_win_adj_all.pl -l -10/200 -m ../CMTSOLUTION_9818433 -n BH -b 0 -k 7 -a STATIONS_ADJOINT -d DATA -s SYN -c RECON -w MEASUREMENT.WINDOWS -i m16 -j 6/30
+plot_win_adj_all.pl -l -10/200 -m ../CMTSOLUTION_9818433 -n BH -b 0 -k 7/1 -a STATIONS_ADJOINT -d DATA -s SYN -c RECON -w MEASUREMENT.WINDOWS -i m16 -j 6/30
cd /data1/cig/seismo/3D/ADJOINT_TOMO/measure_adj_work
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_mt_measure_adj.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_mt_measure_adj.pl 2009-10-22 00:43:10 UTC (rev 15865)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_mt_measure_adj.pl 2009-10-22 03:45:35 UTC (rev 15866)
@@ -34,6 +34,7 @@
#$iboth = 0; # option to plot both cross-correlation and multitaper adjoint sources
($Tmin,$Tmax) = split("/",$Ts);
+(undef,undef,undef,$iadj) = split("/",$iparbools);
#$cmtfile = "CMTSOLUTION_9818433";
#$stafile1 = "EXAMPLE_2/STATIONS_CHT";
@@ -118,28 +119,20 @@
#print CSH "ascii2sac.csh ${plot_adj_dir}*.adj\n";
}
-# create adjoint sources and STATIONS_ADJOINT file for SPECFEM3D
-# prepare_adj_src.pl dumps the ZEN adjoint sources into $adj_dir
-$adj_dir = "ADJOINT_SOURCES";
-print CSH "\\rm -rf ${adj_dir}\n mkdir ${adj_dir}\n";
-print CSH "prepare_adj_src.pl -m $cmtfile -s PLOTS/$stafile1 -o ${adj_dir} OUTPUT_FILES/*adj\n";
-print CSH "cp STATIONS_ADJOINT ${adj_dir}\n";
+if ($iadj == 1) {
+ # create adjoint sources and STATIONS_ADJOINT file for SPECFEM3D
+ # prepare_adj_src.pl dumps the ZEN adjoint sources into $adj_dir
+ $adj_dir = "ADJOINT_SOURCES";
+ print CSH "\\rm -rf ${adj_dir}\n mkdir ${adj_dir}\n";
+ print CSH "prepare_adj_src.pl -m $cmtfile -s PLOTS/$stafile1 -o ${adj_dir} OUTPUT_FILES/*adj\n";
+ print CSH "cp STATIONS_ADJOINT ${adj_dir}\n";
+}
print CSH "\\mv STATIONS_ADJOINT PLOTS\n";
# make plots of (filtered) data, synthetics, windows, and adjoint sources
if ($iplot == 1) {
- # remove any plot files
- #print CSH "\\rm PLOTS/*pdf PLOTS/*jpg PLOTS/*ps\n";
-
print CSH "cd PLOTS\n";
- print CSH "plot_win_adj_all.pl -l $lcut -m ../$cmtfile -n $chan -b $iboth -k $imeas -a $stafile2 -d $dir_data -s $dir_syn -c $dir_recon -w MEASUREMENT.WINDOWS -i $smodel -j $Ts\n";
-
-# # make a single PDF file
-# $isort = 3;
-# $otag = "mt_cc_all";
-# print CSH "make_pdf_by_event.pl $isort $otag\n";
-
- #print CSH "/home/carltape/bin/pdcat -r *.pdf mt_cc_all.pdf\n"; # replace with make_pdf_by_event.pl
+ print CSH "plot_win_adj_all.pl -l $lcut -m ../$cmtfile -n $chan -b $iboth -k $imeas/$iadj -a $stafile2 -d $dir_data -s $dir_syn -c $dir_recon -w MEASUREMENT.WINDOWS -i $smodel -j $Ts\n";
print CSH "cd $pwd\n";
}
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl 2009-10-22 00:43:10 UTC (rev 15865)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl 2009-10-22 03:45:35 UTC (rev 15866)
@@ -12,8 +12,15 @@
# NOTE: The command run_mt_measure_adj.pl below contains additional input parameters.
#
# EXAMPLES (19-Oct-2009):
-# run_tomo.pl m16 0.011 1.0 0 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 5 BH 6/30 0/1/1/0 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 1 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 2 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 3 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 4 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 5 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 6 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 3/0.02/2.5/2.0/2.5/3.5/1.5
# run_tomo.pl m16 0.011 1.0 7 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 1/0.02/2.5/2.0/2.5/3.5/1.5
+# run_tomo.pl m16 0.011 1.0 8 BH 6/30 0/1/1/1 -5.0/5.0/-1.5/1.5/0.7 1/1.0/0.5 1/0.02/2.5/2.0/2.5/3.5/1.5
# run_tomo.pl m16 0.011 1.0 7 BH 3/30 0/1/1/1 -3.0/3.0/-1.5/1.5/0.7 1/1.0/0.5 1/0.02/2.5/0.25/0.5/3.5/1.5
# run_tomo.pl m16 0.011 1.0 7 BH 2/30 0/1/1/1 -2.0/2.0/-1.5/1.5/0.7 1/1.0/0.5 1/0.02/2.5/0.25/0.5/3.5/1.5
#
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/write_par_file.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/write_par_file.pl 2009-10-22 00:43:10 UTC (rev 15865)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/write_par_file.pl 2009-10-22 03:45:35 UTC (rev 15866)
@@ -35,6 +35,7 @@
if($ibool1==1) {$RUN_BANDPASS = ".true."} else {$RUN_BANDPASS = ".false."}
if($ibool2==1) {$DISPLAY_DETAILS = ".true."} else {$DISPLAY_DETAILS = ".false."}
if($ibool3==1) {$OUTPUT_MEASUREMENT_FILES = ".true."} else {$OUTPUT_MEASUREMENT_FILES = ".false."}
+if($ibool4==1) {$COMPUTE_ADJOINT_SOURCE = ".true."} else {$COMPUTE_ADJOINT_SOURCE = ".false."}
($TSHIFT_MIN,$TSHIFT_MAX,$DLNA_MIN,$DLNA_MAX,$CC_MIN) = split("/",$par1);
($ERROR_TYPE,$DT_SIGMA_MIN,$DLNA_SIGMA_MIN) = split("/",$par2);
($itaper,$wtr,$npi,$DT_FAC,$ERR_FAC,$DT_MAX_SCALE,$NCYCLE_IN_WINDOW) = split("/",$par3);
@@ -42,7 +43,7 @@
# comments for the PAR_FILE
#$line00 = "output directory"; # cannot have a string here
$line01 = "tstart, DT, npts: time vector for simulations";
-$line02 = "imeas (0-8) -- 0 no adj src; 1-8 see manual;";
+$line02 = "imeas (1-8; see manual)";
$line03 = "channel: BH or LH";
$line04 = "TLONG and TSHORT: band-pass periods for records";
$line05 = "RUN_BANDPASS: use band-pass on records";
More information about the CIG-COMMITS
mailing list