[cig-commits] r15024 - in seismo/3D/ADJOINT_TOMO/flexwin: . scripts scripts/prepare_scripts/socal
carltape at geodynamics.org
carltape at geodynamics.org
Thu May 21 04:51:50 PDT 2009
Author: carltape
Date: 2009-05-21 04:51:50 -0700 (Thu, 21 May 2009)
New Revision: 15024
Modified:
seismo/3D/ADJOINT_TOMO/flexwin/README
seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh
seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input
seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test
Log:
Successful compilation and run using 64-bit processing, including a beta version of 64-bit SAC (v101.3). gfortran 4.3.3 on Linux Ubuntu. This involved removing the -m32 flag in the Makefiles.
Modified: seismo/3D/ADJOINT_TOMO/flexwin/README
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/README 2009-05-20 21:07:03 UTC (rev 15023)
+++ seismo/3D/ADJOINT_TOMO/flexwin/README 2009-05-21 11:51:50 UTC (rev 15024)
@@ -44,6 +44,8 @@
SAC. Future versions of the SAC distribution may no longer require this
compatibility flag.
+21-MAY-2009: CHT successfully compiled and ran FLEXWIN WITHOUT the -m32 option using a beta version of 64-bit SAC (101.3) with gfortran 4.3.3 on Ubuntu Linux.
+
Steps to compile the flexwin package:
(1) Compile libtau.a and create iasp91.hed and iasp91.tbl. In the
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh 2009-05-20 21:07:03 UTC (rev 15023)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh 2009-05-21 11:51:50 UTC (rev 15024)
@@ -28,7 +28,7 @@
#CHT
rectxt=rectext
-gmtset PAPER_MEDIA letter MEASURE_UNIT inch
+gmtset PAPER_MEDIA letter MEASURE_UNIT inch BASEMAP_TYPE plain PLOT_DEGREE_FORMAT D LABEL_FONT_SIZE 12 ANOT_FONT_SIZE 10 HEADER_FONT_SIZE 12
red=255/0/0
black=0/0/0
grey=220/220/220
@@ -193,9 +193,9 @@
proj=X1.75/1.5
## plot window number statistics
-awk '{print $6}' $t0 | pshistogram -J$proj -C -B:"CC":/WS -W0.02 -Lthin -G$grey -O -K -X-4.5 -Y-2.0 >> $out
-awk '{print $5}' $t0 | pshistogram -J$proj -C -B:"Tshift":/WS -W1 -Lthin -G$grey -O -K -X2.75 >> $out
-awk '{print $7}' $t0 | pshistogram -J$proj -C -B:"dlnA":/WS -W0.1 -Lthin -G$grey -O -X2.75 -K >> $out
+awk '{print $6}' $t0 | pshistogram -J$proj -F -B:"CC":/WS -W0.02 -Lthin -G$grey -O -K -X-4.5 -Y-2.0 >> $out
+awk '{print $5}' $t0 | pshistogram -J$proj -F -B:"Tshift":/WS -W1 -Lthin -G$grey -O -K -X2.75 >> $out
+awk '{print $7}' $t0 | pshistogram -J$proj -F -B:"dlnA":/WS -W0.1 -Lthin -G$grey -O -X2.75 -K >> $out
# set projection for stats plotting
# plot Tshift vs CC
@@ -228,9 +228,9 @@
## plot window number statistics
## NOTE: IF THERE IS A CONSTANT NUMBER OF WINDOWS FOR A COMPONENT, THEN PSHISTOGRAM CRASHES
-grep Z $t1 | awk '{print $2}' | pshistogram -J$proj -C -B:"Windows per Z trace":/WS -W1 -Lthin -G$red -O -K -X-5.5 -Y-2.25 >> $out
-grep R $t1 | awk '{print $2}' | pshistogram -J$proj -C -B:"Windows per R trace":/WS -W1 -Lthin -G$green -O -K -X2.75 >> $out
-grep T $t1 | awk '{print $2}' | pshistogram -J$proj -C -B:"Windows per T trace":/WS -W1 -Lthin -G$blue -O -X2.75 >> $out # FINISH (no -K)
+grep Z $t1 | awk '{print $2}' | pshistogram -J$proj -F -B:"Windows per Z trace":/WS -W1 -Lthin -G$red -O -K -X-5.5 -Y-2.25 >> $out
+grep R $t1 | awk '{print $2}' | pshistogram -J$proj -F -B:"Windows per R trace":/WS -W1 -Lthin -G$green -O -K -X2.75 >> $out
+grep T $t1 | awk '{print $2}' | pshistogram -J$proj -F -B:"Windows per T trace":/WS -W1 -Lthin -G$blue -O -X2.75 >> $out # FINISH (no -K)
# CHT comment out
#echo "%!PS-Adobe-3.0" >t1
@@ -240,7 +240,7 @@
echo $out
-#ggv $out
+#gv $out
#-------------------------------
@@ -249,7 +249,7 @@
out2=${basename}/event_recordsection.eps
out2_ps=${basename}/event_recordsection.ps
out2_pdf=${basename}/event_recordsection.pdf
-gmtset PAPER_MEDIA letter MEASURE_UNIT inch
+#gmtset PAPER_MEDIA letter MEASURE_UNIT inch
region=$min_time/$max_time/$min_ddg/$max_ddg
proj=X3.0/-6.8
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2009-05-20 21:07:03 UTC (rev 15023)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2009-05-21 11:51:50 UTC (rev 15024)
@@ -75,9 +75,9 @@
# directories for data and synthetics
if($idataset == 1) {
- $dir_data = "/net/sierra/raid1/carltape/socal/socal_3D/DATA/FINAL"; # Socal data (Carl)
- $dir_syn = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_${smodel}"; # Socal syn
- #$dir_syn = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_pre_${smodel}"; # Socal syn
+ $dir_data = "/home/carltape/SOCAL_ADJOINT/DATA/FINAL"; # Socal data (Carl)
+ $dir_syn = "/home/carltape/SOCAL_ADJOINT/SYN/model_${smodel}"; # Socal syn
+ #$dir_syn = "/home/carltape/SOCAL_ADJOINT/SYN/model_pre_${smodel}"; # Socal syn
$sdataset = "socal";
} elsif ($idataset == 2) {
@@ -92,7 +92,7 @@
#-------------------------------------
# specify various directories (MUST BE MODIFIED FOR EACH USER)
-$dir0 = "/net/denali/raid1/carltape/svn/cig/seismo/3D";
+$dir0 = "/data1/cig/seismo/3D"; # SVN path
$dir_win_code = "$dir0/ADJOINT_TOMO/flexwin_work";
$dir_win_run = "$dir0/flexwin_run";
#$dir_win_code = "$dir0/ADJOINT_TOMO/flexwin_work_copy2";
@@ -106,10 +106,10 @@
$dir_user = "${dir_win_code}/user_files";
# directory to collect COPIES of various output files
-$odir = "/net/sierra/raid1/carltape/results/WINDOWS/model_${smodel}";
+$odir = "/home/carltape/results/WINDOWS/model_${smodel}";
# run directory for windows, measurements, adjoint sources, and kernels
-$rundir = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
+$rundir = "/home/carltape/SOCAL_ADJOINT/RUNS";
#-------------------------------------
@@ -144,8 +144,8 @@
$dir_win_run_meas = "${dir_win_run}/MEASURE";
# EVENT LIST
-$eid_list = "/net/sierra/raid1/carltape/results/EID_LISTS/syn_run_${smodel}";
-#$eid_list = "/net/sierra/raid1/carltape/results/EID_LISTS/syn_run_iterate";
+$eid_list = "/home/carltape/results/EID_LISTS/syn_run_${smodel}";
+#$eid_list = "/home/carltape/results/EID_LISTS/syn_run_iterate";
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;
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input 2009-05-20 21:07:03 UTC (rev 15023)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input 2009-05-21 11:51:50 UTC (rev 15024)
@@ -15,14 +15,14 @@
comp=`echo $chan | awk '{print substr($1,3,1)}'`
#echo $net $sta $chan $comp
- # synthetics are always BH_
+ # synthetics are always BH_ (for socal problem)
#synt=SYN/${sta}.${net}.BH${comp}.semd.sac.d
#synt=`ls SYN/${sta}.${net}.BH${comp}*`
nsynt=`ls -1 SYN/${sta}.${net}.BH${comp}* | wc | awk '{print $1}'`
out=MEASURE/${sta}.${net}.${chan}
- if [ $nsynt == 1 ] ; then # data with unique syn match
+ if [ $nsynt -eq 1 ] ; then # data with unique syn match
synt=`ls SYN/${sta}.${net}.BH${comp}*`
b1=`saclst b f $file | awk '{print $2}'`
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test 2009-05-20 21:07:03 UTC (rev 15023)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test 2009-05-21 11:51:50 UTC (rev 15024)
@@ -15,14 +15,14 @@
comp=`echo $chan | awk '{print substr($1,3,1)}'`
#echo $net $sta $chan $comp
- # synthetics are always BH_
+ # synthetics are always BH_ (for socal problem)
#synt=SYN_TEST/${sta}.${net}.BH${comp}.semd.sac.d
#synt=`ls SYN_TEST/${sta}.${net}.BH${comp}*`
nsynt=`ls -1 SYN_TEST/${sta}.${net}.BH${comp}* | wc | awk '{print $1}'`
out=MEASURE_TEST/${sta}.${net}.${chan}
- if [ $nsynt == 1 ] ; then # data with unique syn match
+ if [ $nsynt -eq 1 ] ; then # data with unique syn match
synt=`ls SYN_TEST/${sta}.${net}.BH${comp}*`
b1=`saclst b f $file | awk '{print $2}'`
More information about the CIG-COMMITS
mailing list