[cig-commits] r12975 - in seismo/3D/ADJOINT_TOMO/flexwin: . scripts user_files/socal_3D
carltape at geodynamics.org
carltape at geodynamics.org
Tue Sep 30 16:47:09 PDT 2008
Author: carltape
Date: 2008-09-30 16:47:09 -0700 (Tue, 30 Sep 2008)
New Revision: 12975
Modified:
seismo/3D/ADJOINT_TOMO/flexwin/README
seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00_orig.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m01.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m02.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m03.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m04.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m05.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m06.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90
Log:
Modified some DEBUG comments in select_windows_stalta2.f90. Fixed typo in socal user_functions files. Added a couple lines to README.
Modified: seismo/3D/ADJOINT_TOMO/flexwin/README
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/README 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/README 2008-09-30 23:47:09 UTC (rev 12975)
@@ -59,12 +59,18 @@
FLEXWIN requires the following libraries external to the package in order to
compile and run: libsacio.a and libsac.a. Both libraries
-are distributed by IRIS as part of the SAC package (version 101.2 and above).
+are distributed by IRIS as part of the SAC package (version 101.2 and above).
+The IRIS download site (as of 30-Sept-2008) is here:
+ http://www.iris.edu/software/sac/sac.request.htm
+(To check your version, type XXXXX from within sac.)
If your compiler of choice is gfortran, then you should be able to use the
make_gfortran makefiles with only minor modifications (notably you may need to
change the search path for the libsacio.a library). If you prefer another
-compiler, you should modify the OPT and FC lines in the makefiles accordingly.
+compiler, you should modify the OPT and FC lines in the makefiles accordingly.
+We tested using gfortran version 4.1.2
+(To check your version, type "gfortran --version".)
+
IMPORTANT NOTE: All the code is compiled with the -m32 option, which makes
32bit binaries. This option is currently required to enable compatibility with
SAC. Future versions of the SAC distribution may no longer require this
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2008-09-30 23:47:09 UTC (rev 12975)
@@ -142,22 +142,30 @@
$dir_win_run_data = "${dir_win_run}/DATA";
$dir_win_run_meas = "${dir_win_run}/MEASURE";
-# obtain list of events
-#@datadirs = `ls -1 -d ${dir_data}/[1-9]* | sort -g`;
- at datadirs = glob("${dir_data}/[1-9]*");
-$nevent0 = @datadirs;
+# NEW: EVENT LIST
+$eid_list = "/net/sierra/raid1/carltape/results/EID_LISTS/syn_run_m07";
+if (not -f $eid_list) {die("check if eid_list ${eid_list} exist or not\n")}
+open(IN,$eid_list); @eids = <IN>; close(IN);
+$nevent0 = @eids;
+## obtain list of events
+##@datadirs = `ls -1 -d ${dir_data}/[1-9]* | sort -g`;
+#@datadirs = glob("${dir_data}/[1-9]*");
+#$nevent0 = @datadirs;
+
# adjust indices for events
print "\n imin = $imin ; imax = $imax ; nevent = $nevent0";
if($imin < 1) {$imin = 1}
+if($imin > $nevent0) {$imin = $nevent0}
if($imax > $nevent0) {$imax = $nevent0}
$nevent = $imax - $imin + 1;
print "\n imin = $imin ; imax = $imax ; nevent = $nevent\n";
print "\n $nevent events to use in the windowing file\n";
for ($ievent = $imin; $ievent <= $imax; $ievent++) {
- $datadir1 = $datadirs[$ievent-1]; chomp($datadir1);
- $eid = `basename $datadir1`; chomp($eid); # event ID
+ #$datadir1 = $datadirs[$ievent-1]; chomp($datadir1);
+ #$eid = `basename $datadir1`; chomp($eid); # event ID
+ $eid = $eids[$ievent-1]; chomp($eid);
print "\n Event $ievent : $eid";
}
#die("testing");
@@ -172,7 +180,6 @@
# copy PAR_FILE into local directory
# directories for data and synthetics
if($idataset == 1) {
- #$par_file_in = "${dir_user}/socal_3D/PAR_FILE_T06_14095540";
$par_file_in = "${dir_user}/socal_3D/PAR_FILE_${sTmin}_${smodel}";
$userfun_file_in = "${dir_user}/socal_3D/user_functions_${smodel}.f90";
@@ -187,11 +194,20 @@
`cp ${par_file_in} ${par_file}`;
`cp ${userfun_file_in} ${userfun_file}`;
+# copy IASP files into run directory
+$iasp_in1 = "${dir_win_code}/ttimes_mod/iasp91.tbl";
+$iasp_in2 = "${dir_win_code}/ttimes_mod/iasp91.hed";
+$iasp_out1 = "${dir_win_run}/iasp91.tbl";
+$iasp_out2 = "${dir_win_run}/iasp91.hed";
+`cp $iasp_in1 $iasp_out1`;
+`cp $iasp_in2 $iasp_out2`;
+
# name of executable windowing code (in $dir_win_code)
$win_execute = "flexwin";
# make command for windowing code
-$make = "make -f make_intel_caltech";
+$make = "make -f make_gfortran"; # change as of 9-30-08
+#$make = "make -f make_intel_caltech";
# location of plot_windows_all.pl
$plot_windows_perl = "${dir_scripts}/plot_windows_all.pl";
@@ -259,8 +275,10 @@
for ($ievent = $imin; $ievent <= $imax; $ievent++) {
# directories
- $datadir1 = $datadirs[$ievent-1]; chomp($datadir1);
- $eid = `basename $datadir1`; chomp($eid); # event ID
+ $eid = $eids[$ievent-1]; chomp($eid);
+ $datadir1 = "${dir_data}/$eid";
+ #$datadir1 = $datadirs[$ievent-1]; chomp($datadir1);
+ #$eid = `basename $datadir1`; chomp($eid); # event ID
$kk = $kk + 1;
print "Event ID is $eid -- $kk out of $nevent\n";
Modified: seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/select_windows_stalta2.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -40,7 +40,6 @@
call check_data_quality(data_ok,signal_int,noise_int,snr_int,time_obs_signal,time_obs_noise,signal_amp,noise_amp,snr_amp)
! if data quality is not ok, then return without selecting windows
if (.not. data_ok) then
- write(*,*) 'DATA QUALITY NOT OK'
num_win = 0
return
endif
@@ -257,6 +256,7 @@
snr_amp = 0.0
if(DEBUG) write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
+ if(DEBUG) write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
! Integrating the signal and noise. We will take the RMS of the signal
! and noise, and take the ratio. Noise is integrated from start time to
@@ -290,8 +290,13 @@
snr_int = signal_int / noise_int
snr_amp = signal_amp / noise_amp
- ! If snr_int and snr_amp are less than a certain value, data_ok is .false.
+ ! If snr_int or snr_amp is less than a certain value, then data_ok is .false.
if (snr_int < SNR_INTEGRATE_BASE .or. snr_amp < SNR_MAX_BASE) then
+ write(*,*) 'DATA QUALITY NOT OK'
+ if(DEBUG) write(*,*) 'snr_int < SNR_INTEGRATE_BASE: ', &
+ sngl(snr_int), sngl(SNR_INTEGRATE_BASE), snr_int < SNR_INTEGRATE_BASE
+ if(DEBUG) write(*,*) 'snr_amp < SNR_MAX_BASE: ', &
+ sngl(snr_amp), sngl(SNR_MAX_BASE), snr_amp < SNR_MAX_BASE
data_ok = .false.
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00 2008-09-30 23:47:09 UTC (rev 12975)
@@ -26,7 +26,7 @@
# -------------------------------------------------------------
# E(t) water level (0.23)
-STALTA_BASE = 0.22
+STALTA_BASE = 0.18
# -------------------------------------------------------------
# TSHIFT
@@ -34,7 +34,7 @@
# -------------------------------------------------------------
# limit on CC for window acceptance
-CC_BASE = 0.74
+CC_BASE = 0.71
# -------------------------------------------------------------
# limit on dlnA (dA/A) for window acceptance
@@ -48,7 +48,7 @@
# signal-to-noise ratios exceeds these two base levels,
# then the data time series (and syn) is kept
SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
+SNR_MAX_BASE = 2.5
# -------------------------------------------------------------
# limit on signal to noise ratio in a particular window.
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07 2008-09-30 23:47:09 UTC (rev 12975)
@@ -26,7 +26,7 @@
# -------------------------------------------------------------
# E(t) water level (0.23)
-STALTA_BASE = 0.22
+STALTA_BASE = 0.18
# -------------------------------------------------------------
# TSHIFT
@@ -34,7 +34,7 @@
# -------------------------------------------------------------
# limit on CC for window acceptance
-CC_BASE = 0.74
+CC_BASE = 0.71
# -------------------------------------------------------------
# limit on dlnA (dA/A) for window acceptance
@@ -48,7 +48,7 @@
# signal-to-noise ratios exceeds these two base levels,
# then the data time series (and syn) is kept
SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
+SNR_MAX_BASE = 2.5
# -------------------------------------------------------------
# limit on signal to noise ratio in a particular window.
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00_orig.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00_orig.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m00_orig.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m01.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m01.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m01.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m02.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m02.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m02.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m03.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m03.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m03.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m04.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m04.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m04.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m05.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m05.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m05.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m06.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m06.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m06.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90 2008-09-30 22:41:11 UTC (rev 12974)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90 2008-09-30 23:47:09 UTC (rev 12975)
@@ -61,7 +61,7 @@
noise_end = Pnl_start
if(DEBUG) then
write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
More information about the cig-commits
mailing list