[cig-commits] r18825 - in seismo/3D/ADJOINT_TOMO/measure_adj: . scripts_tomo
carltape at geodynamics.org
carltape at geodynamics.org
Fri Aug 12 16:21:10 PDT 2011
Author: carltape
Date: 2011-08-12 16:21:10 -0700 (Fri, 12 Aug 2011)
New Revision: 18825
Modified:
seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90
seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90
seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/prepare_mt_measure_adj.pl
seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_mt_measure_adj.pl
seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl
Log:
added output check statement for COMPUTE_ADJOINT_SOURCE that was missing
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90 2011-08-12 20:39:42 UTC (rev 18824)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/mt_measure_adj.f90 2011-08-12 23:21:10 UTC (rev 18825)
@@ -4,7 +4,7 @@
! and compute the corresponding adjoint sources
! input parameter:
- ! 1. imeas = 1, normalized waveform difference. Adjoint source is constructed from the data
+ ! 1. imeas = 1, normalized waveform difference. Adjoint source is constructed from the data
! only, with the form −d(t)/ || d(t) || 2
! 2. imeas = 2, waveform difference, s(t) − d(t).
! 3. imeas = 3, cross-correlation traveltime difference for a (banana-doughtnut) sensitivity ker-
@@ -21,7 +21,6 @@
! 8. imeas = 8, multitaper ampltidue difference for an event kernel. The measurement between
! data and synthetics is used in constructing the adjoint source. See multitaper_notes.pdf.
-
use mt_variables
use mt_constants
use ascii_rw ! dwascii()
@@ -111,7 +110,6 @@
print *, ' min/max: ',minval(syn(:)),maxval(syn(:))
!print *, ' ',syn(1:5)
-
if (max(npt1,npt2) > NDIM) then
print *, 'Error: Too many number of points in data or syn'
stop
@@ -122,7 +120,6 @@
dt = dt1
npts = min(npt1,npt2)
-
! check
if (abs(t01-t02) > dt) then
print*,'data t0: ',t01
@@ -134,7 +131,6 @@
print*,' npts: ',npts
endif
-
! apply bandpass filter to data and synthetics, if desired
! http://www.iris.washington.edu/pipermail/sac-help/2008-March/000376.html
! Access to the kidate, xapiir, and getfil is not simple and not
@@ -208,7 +204,6 @@
if( cmp(1:1) == "T" ) ipick_type = Love_T
endif
-
! write values to output file
nwin = nwin + 1 ! overall window counter
write(12,'(a3,a8,a5,a5,3i5,2f12.3)') net,sta,chan_syn,chan_dat,nwin,ipair,j,tstart,tend
@@ -244,7 +239,7 @@
!call mt_measure_select_0(nlen,dt,i_left,i_right,fstart,fend,use_trace)
! make measurements
- ! CHT: also compute reconstructed synthetics for CC (and MT, if specified) measurements
+ ! also compute reconstructed synthetics for CC (and MT, if specified) measurements
call mt_measure(datafile,measure_file_prefix,data,syn,t0,dt,npts,tstart,tend,&
istart,data_dtw,syn_dtw,nlen,tshift,sigma_dt_cc,dlnA,sigma_dlnA_cc,cc_max,syn_dtw_cc,&
i_pmax_dat,i_pmax_syn,i_right0,trans_mtm,dtau_w,dlnA_w,sigma_dt,sigma_dlnA,err_dt,err_dlnA)
@@ -264,7 +259,7 @@
print *, 'fstart/fend :', fstart, fend
print *, 'Tpmax :', T_pmax_dat, T_pmax_syn
- ! CHT: If MT measurement window is rejected by mt_measure_select, then use a CC measurement.
+ ! if MT measurement window is rejected by mt_measure_select, then use a CC measurement
if(.not. use_trace) then
!stop 'Check why this MT measurement was rejected'
print *, 'Reverting from multitaper measurement to cross-correlation measurement'
@@ -412,7 +407,6 @@
! but this seems to be the best and fairly easy what we can do here...
call bandpass(adj_syn_all,npts,dt,fstart0,fend0)
-
! cut and interpolate to match time-stepping for SEM
! NOTE: This can leave a non-zero value to start the record,
! which is NOT GOOD for the SEM simulation.
@@ -424,7 +418,8 @@
itmax = int(TSHORT/dtt)
call taper_start(adj_syn_all,nn,itmax)
- ! we can have choices of outputing the adjoint source as ASCII or SAC format
+ ! output the adjoint source (or ray density) as ASCII or SAC format
+ print *, 'writing adjoint source to file for the full seismogram'
if( DO_RAY_DENSITY_SOURCE ) then
call dwascii(trim(adj_file_prefix)//'.density.adj',adj_syn_all,nn,tt,dtt)
else
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90 2011-08-12 20:39:42 UTC (rev 18824)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/mt_sub.f90 2011-08-12 23:21:10 UTC (rev 18825)
@@ -1862,14 +1862,15 @@
! check the read-in values
print *, 'INPUTS FROM MEASUREMENT.PAR :'
print *, ' tt, dtt, nn : ',tt,dtt,nn
+ print *, ' imeas : ',imeas
print *, ' chan : ',chan
- print *, ' imeas : ',imeas
print *, ' TLONG, TSHORT : ',TLONG, TSHORT
fstart0 = 1./TLONG ; fend0 = 1./TSHORT
print *, ' fstart, fend : ', fstart0, fend0
print *, ' RUN_BANDPASS : ',RUN_BANDPASS
print *, ' DISPLAY_DETAILS : ',DISPLAY_DETAILS
print *, ' OUTPUT_MEASUREMENT_FILES : ',OUTPUT_MEASUREMENT_FILES
+ print *, ' COMPUTE_ADJOINT_SOURCE : ',COMPUTE_ADJOINT_SOURCE
print *, ' TSHIFT_MIN, TSHIFT_MAX : ',TSHIFT_MIN, TSHIFT_MAX
print *, ' DLNA_MIN, DLNA_MAX : ',DLNA_MIN, DLNA_MAX
print *, ' CC_MIN : ',CC_MIN
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/prepare_mt_measure_adj.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/prepare_mt_measure_adj.pl 2011-08-12 20:39:42 UTC (rev 18824)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/prepare_mt_measure_adj.pl 2011-08-12 23:21:10 UTC (rev 18825)
@@ -52,7 +52,7 @@
# directories for data, synthetics, and window files : MUST BE CHANGED FOR EACH USER
if ($iwindow==1) {
# windowing code run directory
- $dir_win_run = "/data1/cig/seismo/3D/flexwin_run";
+ $dir_win_run = "/data2/SVN/seismo/3D/flexwin_run";
$dir_win_run_syn = "${dir_win_run}/${dir_syn}";
$dir_win_run_data = "${dir_win_run}/${dir_data}";
$win_in = "${dir_win_run}/${dir_meas}/MEASUREMENT_WINDOWS";
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 2011-08-12 20:39:42 UTC (rev 18824)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_mt_measure_adj.pl 2011-08-12 23:21:10 UTC (rev 18825)
@@ -15,7 +15,7 @@
# INPUT:
# smodel index label for model
# lcut cut times for plotting records
-# itest DATA_TEST/SYN_TEST or DATA/SYN
+# itest use DATA_TEST/SYN_TEST or DATA/SYN
# iplot make plots of data, syn, windows, and adjoint sources
# iboth plot both MT and CC adjoint sources (=1) or not (=0)
# ... see write_par_file.pl
Modified: seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl 2011-08-12 20:39:42 UTC (rev 18824)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/scripts_tomo/run_tomo.pl 2011-08-12 23:21:10 UTC (rev 18825)
@@ -12,12 +12,12 @@
# 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 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 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 # CC: no adjoint source
# 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 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 # CC: adjoint source
# 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
More information about the CIG-COMMITS
mailing list