[cig-commits] r13109 - in seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS: . station_lists vtk
carltape at geodynamics.org
carltape at geodynamics.org
Mon Oct 20 22:14:27 PDT 2008
Author: carltape
Date: 2008-10-20 22:14:27 -0700 (Mon, 20 Oct 2008)
New Revision: 13109
Added:
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/README
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/STATIONS_in
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/make_station_lists.m
Modified:
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/kernel_combine.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs_pmax.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update_local.tcl
Log:
Added script to read in STATIONS file and list of event IDs, and then output sets of files (1) per event, with stations sorted by azimuth; (2) per station, with events sorted by back-azimuths. These are very handy for making composite plots and for quick checks to see nearby stations or events.
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/kernel_combine.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/kernel_combine.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/kernel_combine.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -4,23 +4,27 @@
#
# kernel_combine.pl
# Carl Tape
-# 26-Nov-2007
+# 19-June-2008
#
# EXAMPLE:
-# ~/UTILS/kernel_combine.pl 3 mu_all_kernel_03
+# ~/UTILS/kernel_combine.pl 5 7 mu_all_kernel_set
#
#==========================================================
-if (@ARGV < 2) {die("Usage: kernel_combine.pl imodel suffix\n")}
-($imodel,$ksuffix) = @ARGV;
+if (@ARGV < 3) {die("Usage: kernel_combine.pl imodel_min imodel_max suffix\n")}
+($imodel_min,$imodel_max,$ksuffix) = @ARGV;
-$smodel = sprintf("m%2.2i",$imodel);
-
# base directory
$dir0 = "/net/sierra/raid1/carltape/results/KERNELS";
+$smodel = sprintf("m%2.2i",$imodel_max);
$kdirmax = "$dir0/kernel_${smodel}";
if (not -e $kdirmax) {die("check if kdirmax $kdirmax exist or not\n")}
+# list of event IDs
+#$file_eids = "/net/sierra/raid1/carltape/results/EID_LISTS/kernels_use_m05";
+#if (not -f $file_eids) {die("\n check if $file_eids exists\n")}
+#open(IN,$file_eids); @eids = <IN>; $numk = @eids;
+
@kfiles = glob("$kdirmax/*${ksuffix}.pdf");
$numk = @kfiles;
print "$numk kernel files in $kdirmax\n";
@@ -37,15 +41,18 @@
for ($k = 1; $k <= $numk; $k++) {
- $kfilemax = @kfiles[$k-1]; chomp($kfilemax);
+ $kfilemax = $kfiles[$k-1]; chomp($kfilemax);
$ktag = `basename $kfilemax`; chomp($ktag);
- print "$kfilemax -- $ktag\n";
+ print "$kfilemax -- $ktag --\n";
- for ($m = 0; $m <= $imodel; $m++) {
+ for ($m = $imodel_min; $m <= $imodel_max; $m++) {
$smod = sprintf("m%2.2i",$m);
$kfilem = "$dir0/kernel_${smod}/$ktag";
- if (not -f $kfilem) {die("check if kfilem $kfilem exist or not\n")}
- $pdcat[$p] = $kfilem; $p = $p+1;
+ if (not -f $kfilem) {
+ print "check if kfilem $kfilem exist or not\n";
+ } else {
+ $pdcat[$p] = $kfilem; $p = $p+1;
+ }
}
}
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in 2008-10-21 05:14:27 UTC (rev 13109)
@@ -0,0 +1,204 @@
+9105672 -116.8417 34.3228
+9108606 -116.2694 34.5965
+3320736 -116.2465 34.4368
+9109243 -115.7045 33.2806
+9109752 -116.0675 34.2855
+9110281 -116.2527 34.5588
+9111353 -116.1448 34.3512
+9112735 -116.3442 34.7110
+9113909 -116.3952 34.8632
+3320883 -116.3982 34.8668
+3321595 -116.3898 34.8648
+3321590 -116.3955 34.8735
+3320884 -116.3978 34.8633
+3324595 -116.2639 34.6055
+9114763 -116.2085 34.3300
+3321597 -116.1944 34.4042
+9114812 -116.4060 34.8620
+9116921 -117.9045 36.4512
+9117942 -116.2707 34.5200
+3321426 -116.2888 34.8031
+9120741 -116.2980 34.7903
+9122706 -116.4105 34.8343
+9128775 -117.0072 34.1048
+9130422 -116.2635 34.5927
+9132433 -117.0028 34.1040
+7177729 -116.3710 34.8038
+9140050 -117.2432 34.0588
+7179710 -117.6053 36.0882
+9141142 -117.6010 36.0850
+9146641 -115.3868 32.7020
+9147453 -116.2704 34.8046
+9150885 -118.0484 36.3269
+9151000 -118.3073 35.0942
+3324419 -118.0519 36.3249
+9151609 -118.0507 36.3267
+9153800 -116.3622 34.6845
+9154092 -115.5035 32.8898
+9154233 -115.5060 32.9034
+9155518 -116.2985 34.7878
+9165761 -118.9170 34.5555
+9169867 -116.7722 34.2673
+9171064 -117.1652 35.5801
+9171679 -119.0298 34.8942
+9173365 -118.4178 34.2833
+9173374 -118.4173 34.2872
+9627557 -116.1394 33.8170
+9627721 -116.9397 34.2910
+9631385 -117.5195 35.1143
+9644101 -117.7091 33.8725
+9644345 -118.3256 35.9817
+9653493 -118.0423 35.7957
+9655209 -116.7523 34.0299
+9660449 -115.4224 32.7416
+9666905 -116.7605 34.2594
+10964587 -117.8735 36.0223
+9674049 -117.8682 36.0135
+9674213 -117.8753 36.0140
+10972299 -117.4629 34.2696
+10992159 -117.8723 35.9915
+11671240 -117.8650 35.9783
+9686565 -117.8632 36.0495
+9688709 -117.8720 35.9896
+9695397 -118.2513 32.7920
+9695549 -118.3180 32.7340
+9700049 -116.3986 33.3807
+9700217 -120.2220 36.2250
+9703873 -118.3968 34.0527
+9716853 -118.2795 33.9297
+9718013 -116.5023 33.5112
+9722529 -115.7100 33.3070
+9722633 -115.7012 33.3170
+9722669 -115.7045 33.3135
+9734033 -116.7048 34.1178
+9735129 -117.7483 33.9552
+9742277 -116.4308 33.3852
+9744905 -115.0875 32.4633
+9504418 -115.0893 32.4723
+9753485 -118.6642 34.3638
+9753489 -118.6642 34.3670
+9753497 -118.6645 34.3655
+9753949 -118.6669 34.3631
+9755013 -118.6670 34.3647
+12456160 -115.7451 32.5553
+12659440 -119.3317 33.6678
+9772973 -116.7693 34.2619
+9774569 -116.7117 33.2065
+9775765 -116.2952 34.5177
+12887732 -118.0758 35.7057
+13303428 -115.9762 33.1627
+13692644 -117.4322 34.1653
+9818433 -117.7840 33.9133
+9826789 -116.1118 33.2352
+9828889 -117.2913 35.9468
+9829213 -117.2917 35.9490
+9853417 -116.5685 33.5125
+9854597 -116.2650 34.8068
+9875657 -118.6677 35.3185
+9875665 -118.6636 35.3102
+9882325 -118.6632 35.3152
+9882329 -118.6585 35.3128
+13935988 -116.8460 34.3103
+13936432 -116.8547 34.3208
+13936812 -116.8482 34.3097
+13938812 -116.8407 34.3137
+13945908 -116.1303 34.3582
+13965956 -121.1007 36.5565
+13966396 -115.5538 32.9475
+13966672 -115.5409 32.9455
+13970876 -115.5472 32.9443
+9930549 -116.6665 34.6172
+13978600 -121.1410 34.7150
+13978612 -121.1133 34.7002
+13986104 -118.2692 36.4782
+9941081 -118.6509 34.4034
+9944301 -115.2837 32.5620
+14007388 -117.5664 35.6352
+9966869 -120.8890 35.5432
+9967025 -121.0993 35.6872
+9967137 -121.1113 35.6933
+9967249 -121.0983 35.6665
+9967541 -120.9095 35.5330
+9967901 -121.0428 35.6493
+9968525 -120.4817 34.7847
+9968977 -120.8385 35.5487
+14018508 -121.1372 35.7032
+9982749 -120.8740 35.5918
+9983429 -119.1412 35.0118
+9983625 -120.8693 35.5800
+9994573 -118.1589 36.1815
+10005209 -120.8028 35.5067
+10006857 -120.0142 34.4135
+14072464 -115.7441 32.5392
+14073800 -116.0520 33.7152
+14077668 -119.4365 34.3885
+14079184 -117.4478 34.1358
+14094528 -120.6661 36.1434
+14094992 -120.3667 35.8155
+14094996 -120.3531 35.8029
+14095540 -120.5134 35.9528
+14095628 -118.6292 35.3852
+14096196 -120.5403 35.9821
+14096736 -120.8108 35.5473
+10059745 -116.8413 34.3533
+10061489 -120.4113 35.8570
+10063349 -120.4963 35.9437
+10065241 -119.6955 34.2883
+10066389 -116.9755 34.1535
+14116920 -117.4420 34.1225
+14116972 -117.4438 34.1272
+14118096 -116.3912 33.9578
+14133048 -116.2515 33.2884
+14137160 -116.8122 32.7233
+14138080 -119.1940 34.9987
+14139108 -120.0295 33.6922
+14139160 -120.0258 33.6880
+10094253 -118.6268 33.8280
+10097009 -119.1958 35.0023
+10100053 -120.4792 35.9269
+14151344 -116.5675 33.5380
+14155056 -120.1020 36.0720
+14155260 -117.0072 34.0612
+14158696 -117.0232 34.0615
+14165408 -119.7527 33.6853
+14169456 -118.0652 36.1488
+14178184 -115.6207 33.1544
+14178188 -115.6098 33.1639
+14178212 -115.6157 33.1548
+14178236 -115.5924 33.1748
+14178248 -115.5969 33.1712
+14179288 -115.6168 33.1538
+14179292 -115.6064 33.1643
+14179736 -115.6295 33.1479
+14181056 -116.8393 32.5112
+14183744 -116.0260 33.1787
+14186612 -119.0247 35.0178
+14186928 -119.0367 35.0198
+14187364 -118.4718 35.3782
+14189556 -121.0838 35.6500
+10148369 -116.7725 34.0198
+10148421 -116.7715 34.0182
+10148829 -118.1450 32.4970
+14204000 -117.5450 35.1267
+14215316 -117.5402 35.1228
+14215812 -115.5730 32.5678
+14219360 -117.5828 35.6232
+10186185 -115.8518 32.7050
+10187953 -116.7903 33.9217
+14236768 -116.0220 33.2450
+14239184 -117.1103 33.8567
+10207681 -116.0402 32.7333
+14255632 -116.0632 33.2663
+10215753 -116.0520 32.7165
+10217433 -116.0312 32.7170
+14263252 -120.7510 35.6360
+14263544 -115.9628 32.8423
+14263712 -115.9672 32.8385
+14263716 -115.9672 32.8377
+14263768 -115.9632 32.8318
+10222697 -120.2430 36.1860
+10222753 -120.2640 36.1950
+10223765 -116.0448 33.7063
+10226877 -116.2947 32.9945
+10230869 -116.1357 33.2220
+14383980 -117.7630 33.9560
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/README
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/README (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/README 2008-10-21 05:14:27 UTC (rev 13109)
@@ -0,0 +1,15 @@
+Carl Tape, 01-Oct-2008
+
+/net/denali/scratch1/carltape/svn
+ /cig/seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/
+
+First, get two files, STATIONS_in and EIDS_in, each with the three columns
+ name longitude latitude
+
+This is what I did to get the files:
+awk '{print $2,$8,$9}' /net/sierra/raid1/carltape/results/SOURCES/socal_10/EIDs_lonlat_otime > EIDS_in
+awk '{print $3"."$4,$1,$2}' /net/denali/home2/carltape/gmt/stations/seismic/Matlab_output/STATIONS_CALIFORNIA_TOMO_INNER_gmt > STATIONS_in
+
+Then run make_station_lists.m from Matlab to generate the output lists.
+
+-----------------------
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/STATIONS_in
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/STATIONS_in (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/STATIONS_in 2008-10-21 05:14:27 UTC (rev 13109)
@@ -0,0 +1,369 @@
+BVA1.AZ -116.365500 33.326500
+BVA2.AZ -116.365500 33.326500
+BVDA1.AZ -116.365500 33.326500
+BVDA2.AZ -116.365500 33.326500
+BZN.AZ -116.667000 33.491500
+CPE.AZ -117.105100 32.888900
+CRY.AZ -116.737300 33.565400
+CSLB.AZ -118.112600 33.779300
+DHL.AZ -115.787500 33.388900
+FLV.AZ -117.061900 32.911400
+FLV2.AZ -117.061900 32.911400
+FRD.AZ -116.602200 33.494700
+FSE.AZ -116.592600 33.005700
+HWB.AZ -116.959600 33.026200
+KNW.AZ -116.711900 33.714100
+LVA2.AZ -116.561500 33.351600
+MONP.AZ -116.422500 32.892700
+MTRP.AZ -117.056700 32.820300
+PFO.AZ -116.459400 33.611700
+RDM.AZ -116.847800 33.630000
+SCI2.AZ -118.487900 32.915000
+SCV1.AZ -116.459700 33.611800
+SCV2.AZ -116.459700 33.611800
+SMER.AZ -117.170800 33.457700
+SMER2.AZ -117.170800 33.457700
+SMTC.AZ -115.799900 32.944900
+SND.AZ -116.612900 33.551900
+SOL.AZ -117.248000 32.841000
+THSB.AZ -117.226900 32.878800
+TRO.AZ -116.425700 33.523400
+WMC.AZ -116.674700 33.573600
+YAQ.AZ -116.353900 33.166600
+HAST.BK -121.551400 36.388700
+PKD.BK -120.541600 35.945170
+PKD1.BK -120.426100 35.889400
+PKD2.BK -120.430500 35.888200
+RAMR.BK -120.869800 35.636000
+CCRB.BP -120.551610 35.957160
+EADB.BP -120.422600 35.895220
+FROB.BP -120.486880 35.910950
+GHIB.BP -120.347280 35.832250
+JCNB.BP -120.431120 35.938970
+JCSB.BP -120.434000 35.921200
+LCCB.BP -120.514230 35.980007
+MMNB.BP -120.496000 35.956500
+RMNB.BP -120.477720 36.000870
+SCYB.BP -120.536610 36.009420
+SMNB.BP -120.579880 35.972970
+VARB.BP -120.447050 35.926080
+VCAB.BP -120.533920 35.921620
+14405.CE -118.356000 33.787000
+14901.CE -118.117000 33.562000
+23732.CE -117.330000 34.190000
+24400.CE -118.178000 34.037000
+24851.CE -118.346000 34.070000
+24853.CE -118.286000 34.077000
+25131.CE -119.198000 34.151000
+26903.CE -120.731000 34.611000
+G405.CE -118.356000 33.787000
+J732.CE -117.330000 34.190000
+K851.CE -118.346000 34.070000
+K853.CE -118.286000 34.077000
+ADO.CI -117.433910 34.550460
+AGA.CI -116.401100 33.638380
+AGO.CI -118.766990 34.146470
+ALP.CI -118.299460 34.687080
+ARG.CI -117.480600 36.083700
+ARV.CI -118.830090 35.126900
+BAK.CI -119.104450 35.344440
+BAR.CI -116.672150 32.680050
+BBA.CI -118.353400 34.195510
+BBB.CI -115.733250 33.352620
+BBR.CI -116.920750 34.262300
+BBS.CI -116.980580 33.921390
+BC3.CI -115.453660 33.655150
+BCC.CI -117.261190 33.575080
+BEL.CI -115.998200 34.000600
+BFS.CI -117.658530 34.238830
+BKR.CI -116.070300 35.269300
+BLA.CI -116.388960 34.069480
+BOR.CI -116.417160 33.268200
+BR2.CI -115.549360 32.978990
+BRE.CI -117.981160 33.807760
+BTC.CI -115.219870 33.012130
+BTP.CI -118.573980 34.682240
+BVH.CI -118.395900 34.076180
+CAB.CI -118.640930 34.155730
+CAC.CI -118.121990 34.136780
+CAP.CI -117.194950 33.388220
+CAR.CI -119.845830 35.308190
+CBC.CI -118.127500 34.140110
+CCC.CI -117.364530 35.524950
+CFS.CI -117.282870 34.105390
+CGO.CI -117.802950 36.550400
+CHF.CI -118.025850 34.333410
+CHN.CI -117.680440 33.998790
+CIA.CI -118.413720 33.401860
+CLC.CI -117.597510 35.815740
+CLT.CI -117.316870 34.092840
+CPP.CI -117.809000 34.060200
+CRN.CI -117.561060 33.875740
+CRP.CI -118.127050 34.136240
+CTC.CI -115.990060 33.655160
+CTD.CI -118.444360 33.354860
+CWC.CI -118.080160 36.439880
+DAN.CI -115.381150 34.637450
+DEC.CI -118.333830 34.253530
+DEV.CI -116.577940 33.935970
+DGR.CI -117.009470 33.650010
+DJJ.CI -118.455050 34.106180
+DJJB.CI -118.455050 34.106180
+DLA.CI -118.096240 33.848220
+DNR.CI -116.630560 33.566670
+DPP.CI -116.941530 32.998620
+DRC.CI -115.446540 32.805400
+DRE.CI -115.446790 32.805320
+DSC.CI -116.103950 35.142550
+DVT.CI -116.100610 32.659150
+EDW.CI -117.991060 34.883030
+EDW2.CI -117.993880 34.881100
+ELC.CI -115.535600 32.781400
+EML.CI -116.845660 32.890830
+ERR.CI -115.822670 33.116750
+EW2.CI -116.408390 33.940590
+FIG.CI -119.988030 34.728320
+FLL.CI -118.918070 34.397270
+FLS.CI -117.039090 34.970360
+FMP.CI -118.293810 33.712640
+FON.CI -117.438760 34.099570
+FPC.CI -117.582670 35.082000
+FUL.CI -117.922510 33.871700
+FUR.CI -116.863220 36.467030
+GMR.CI -115.659940 34.784570
+GOR.CI -117.229210 33.153700
+GR2.CI -118.299400 34.118300
+GRF.CI -118.300400 34.119200
+GRH.CI -118.559540 34.308030
+GSA.CI -118.128300 34.136770
+GSC.CI -116.805740 35.301770
+GVR.CI -118.119970 34.049700
+HEC.CI -116.335000 34.829400
+HLL.CI -118.359670 34.176430
+HLN.CI -117.218970 34.121280
+IRM.CI -115.145130 34.157380
+ISA.CI -118.474030 35.662780
+JCS.CI -116.595900 33.085900
+JRC.CI -117.807600 35.982300
+JRC2.CI -117.808850 35.982490
+JVA.CI -116.612660 34.366220
+KIK.CI -118.101560 34.150370
+KML.CI -117.879090 33.597380
+LAF.CI -118.331430 33.868890
+LAX.CI -118.413900 33.943820
+LBW1.CI -118.088370 33.798000
+LBW2.CI -118.088370 33.798000
+LCG.CI -118.377940 34.000330
+LCP.CI -120.279960 34.735510
+LDF.CI -115.184160 35.130660
+LDR.CI -118.341560 34.990600
+LEV.CI -118.291040 34.614620
+LFP.CI -118.488050 34.305290
+LGB.CI -118.149180 33.975300
+LGU.CI -119.065870 34.108190
+LJR.CI -118.867750 34.807620
+LKL.CI -117.824930 34.615940
+LLS.CI -117.943040 33.684470
+LMR.CI -117.696610 34.934260
+LMR2.CI -117.696290 34.934390
+LRL.CI -117.682120 35.479540
+LTP.CI -118.175680 33.881100
+LTR.CI -117.990300 34.521100
+LUG.CI -117.366830 34.365600
+MAG.CI -118.922820 35.361650
+MCT.CI -116.040730 34.226450
+MGE.CI -116.368740 33.818410
+MIK.CI -118.126080 34.137070
+MIS.CI -118.335170 33.740300
+MLS.CI -117.561620 34.004600
+MOP.CI -118.904900 34.280850
+MPI.CI -119.145300 34.812580
+MPM.CI -117.489010 36.057990
+MPP.CI -119.813620 34.888480
+MSJ.CI -116.967890 33.808010
+MTL.CI -118.237400 34.268800
+MTP.CI -115.553200 35.484340
+MUR.CI -117.195430 33.599990
+MWC.CI -118.058320 34.223620
+NBS.CI -116.557980 34.780350
+NJQ.CI -120.177370 34.534120
+NOT.CI -118.558290 34.228690
+NSS.CI -115.946500 33.556100
+NSS2.CI -115.945860 33.555530
+OGC.CI -117.844000 33.788160
+OKV.CI -119.298740 34.396970
+OLI.CI -117.923720 33.945390
+OLP.CI -116.930360 32.607830
+OSI.CI -118.723500 34.614500
+PAS.CI -118.171170 34.148440
+PASC.CI -118.185230 34.171410
+PDE.CI -118.582150 34.441990
+PDR.CI -118.437020 33.962730
+PDU.CI -117.638080 34.120700
+PER.CI -117.205290 33.861640
+PHL.CI -120.545560 35.407730
+PLC.CI -116.511950 33.824360
+PLM.CI -116.862650 33.353610
+PLS.CI -117.609060 33.795300
+PMD.CI -116.377690 33.647850
+PSD.CI -116.550260 33.823930
+PSR.CI -117.807090 34.091810
+QUG.CI -118.498370 34.395580
+RCT.CI -119.243840 36.305230
+RIN.CI -118.479250 34.281960
+RINB.CI -118.479250 34.281960
+RIO.CI -117.979560 34.104730
+RPV.CI -118.404120 33.743460
+RRS.CI -117.366460 33.882170
+RRX.CI -116.996840 34.875330
+RSB.CI -117.327200 33.973100
+RSS.CI -117.326740 33.973270
+RUS.CI -118.080850 34.050730
+RVR.CI -117.375450 33.993510
+RXH.CI -115.622570 33.183130
+SAL.CI -115.985850 33.280100
+SAN.CI -117.885780 33.704320
+SBB2.CI -117.824480 34.688500
+SBC.CI -119.714920 34.440760
+SBI.CI -119.029860 33.480460
+SBP.CI -117.234880 34.232360
+SBPX.CI -117.234840 34.232400
+SCI.CI -118.546970 32.979900
+SCI2.CI -118.546970 32.979900
+SCZ.CI -119.634350 33.995320
+SCZ2.CI -119.635100 33.995430
+SDD.CI -117.661710 33.552590
+SDG.CI -117.138050 32.784000
+SDP.CI -120.501370 34.565470
+SDR.CI -116.942410 32.735610
+SES.CI -119.137500 34.436920
+SHO.CI -116.275300 35.899530
+SIO.CI -119.164600 34.293000
+SJU.CI -117.681140 33.487250
+SLA.CI -117.283320 35.890950
+SLR.CI -116.797370 33.833590
+SMB.CI -120.447090 34.902180
+SMF.CI -118.446750 34.021590
+SMM.CI -119.995810 35.314200
+SMS.CI -118.456170 34.014380
+SMV.CI -118.744070 34.270910
+SNCC.CI -119.524000 33.248000
+SOT.CI -118.449000 34.416000
+SPF.CI -118.646140 34.059330
+SPG.CI -118.810990 36.135500
+SPG2.CI -118.766240 36.200570
+SRN.CI -117.789380 33.828540
+SSW.CI -115.615640 33.177470
+STC.CI -119.186760 34.303020
+STG.CI -117.768560 33.664000
+STS.CI -118.198780 33.790330
+SVD.CI -117.098220 34.106470
+SWS.CI -115.795800 32.940800
+SYP.CI -119.978340 34.527750
+TA2.CI -117.678220 34.382030
+TCF.CI -118.599000 34.083770
+TEH.CI -118.420790 35.291300
+TFT.CI -119.419460 35.145920
+THX.CI -116.164020 33.634950
+TOR.CI -116.225840 33.575260
+TOV.CI -118.820390 34.156070
+TUQ.CI -115.923890 35.435840
+USB.CI -119.842700 34.413000
+USC.CI -118.286310 34.019190
+VCS.CI -118.117830 34.483640
+VES.CI -119.084690 35.840890
+VRD.CI -118.279640 34.214590
+VTV.CI -117.329600 34.560650
+WBS.CI -118.140350 35.536640
+WER.CI -119.027110 35.060530
+WES.CI -115.731610 32.759040
+WGR.CI -119.274070 34.510850
+WLT.CI -117.950770 34.009480
+WNS.CI -118.380000 34.115000
+WSS.CI -118.649710 34.171700
+WTT.CI -118.255470 33.948690
+FABS.FA -118.571113 34.066670
+DAC.LB -117.593700 36.277000
+PAGB.NC -120.249860 35.730650
+PHOB.NC -120.479603 35.866610
+PHSB.NC -121.053930 35.824040
+PMPB.NC -120.801263 36.215905
+PSM.NC -120.596180 36.068807
+CC12.NN -115.170900 36.133700
+CC16.NN -115.061800 36.161600
+CC22.NN -115.238200 36.113400
+CHYN.NN -115.105700 36.222300
+CLRK.NN -115.204700 36.084600
+CWC.NN -118.080200 36.439900
+DOEL.NN -115.139300 36.209700
+FMT.NN -116.782600 36.637500
+GMR.NN -115.659900 34.784600
+GSLV.NN -115.044200 36.037900
+GWY.NN -116.668900 36.186000
+HF05.NN -115.089400 36.030700
+HYDE.NN -115.192300 36.160900
+LV04.NN -115.131700 36.160200
+LV07.NN -115.314600 36.181600
+LV09.NN -115.241900 36.247500
+MOHS.NN -115.138400 36.254700
+NELL.NN -115.034700 36.220300
+PAN.NN -117.099000 36.393200
+QSM.NN -116.868200 35.965000
+RRK.NN -115.466000 36.065700
+SHP.NN -115.160200 36.505500
+UNLV.NN -115.141600 36.108200
+VOTK.NN -115.069900 36.080900
+WTNK.NN -115.009400 36.191600
+0655.NP -118.499180 34.313130
+286.NP -115.824303 32.954768
+5030.NP -117.991000 34.521333
+5062.NP -115.615640 33.177470
+5081.NP -118.600000 34.084000
+5271.NP -115.733000 33.352667
+5398.NP -118.353667 34.194833
+5399.NP -118.411000 33.943000
+5401.NP -118.918000 34.397000
+5402.NP -118.395000 34.076000
+5403.NP -119.298667 34.396833
+5404.NP -118.641000 34.155000
+5442.NP -116.940323 34.049048
+5443.NP -116.310790 33.844470
+5444.NP -116.095340 33.634610
+5445.NP -119.027940 35.061310
+5470.NP -119.026440 35.060810
+5471.NP -118.362010 34.199880
+BBA.NP -118.353400 34.195510
+BBB.NP -115.733250 33.352620
+BVH.NP -118.395900 34.076180
+FLL.NP -118.918070 34.397270
+GRF.NP -118.300400 34.119200
+JAB.NP -118.496700 34.310920
+JGB.NP -118.499180 34.313130
+LAX.NP -118.413900 33.943820
+LT2.NP -117.990990 34.521490
+LTR.NP -117.990300 34.521100
+OKV.NP -119.298740 34.396970
+TCF.NP -118.599000 34.083770
+B081.PB -116.714170 33.711170
+B082.PB -116.596010 33.598180
+B084.PB -116.456370 33.611570
+B086.PB -116.530000 33.557500
+B087.PB -116.602670 33.495500
+DCD.PG -120.840833 35.212167
+DPD.PG -120.781667 35.233000
+SHD.PG -120.761333 35.168167
+WRD.PG -120.881667 35.458833
+GVDA.SB -116.673798 33.668900
+WLA.SB -115.530602 33.097598
+AMD.SN -116.280900 36.452600
+LEC.SN -116.689600 36.562700
+109C.TA -117.105100 32.888900
+HAST.TA -121.551400 36.388700
+U04C.TA -120.782800 36.362700
+U05C.TA -120.120500 36.335600
+U10A.TA -116.329700 36.419300
+U11A.TA -115.383500 36.423000
+V03C.TA -121.235600 36.021400
+V04C.TA -120.869900 35.636000
+V05C.TA -119.902800 35.866700
+V11A.TA -115.430500 35.838400
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/make_station_lists.m
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/make_station_lists.m (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/make_station_lists.m 2008-10-21 05:14:27 UTC (rev 13109)
@@ -0,0 +1,90 @@
+%
+% make_station_lists.m
+% CARL TAPE, 01-Oct-2008
+% printed xxx
+%
+% This program reads in a list of receivers and a list of events, and it
+% outputs a set of lists sorted by azimuth and distance. These are useful
+% for sorting sets of data or plots, or also for a quick reference to find
+% nearby stations for a particular event.
+%
+% calls xxx
+% called by xxx
+%
+
+clc, clear, close all, format short, format compact
+
+% output directory to dump the files into
+odir = '/net/sierra/raid1/carltape/results/SOURCES/EID_STATION_LISTS/';
+
+% format statement for output files
+stfmt = '%12s%10.4f%10.4f%12s%10.4f%10.4f%10.4f%10.4f\n';
+
+%------------------------------------------
+% read in files
+
+[eid_name,elon,elat] = textread('EIDS_in','%s%f%f');
+[rec_name,rlon,rlat] = textread('STATIONS_in','%s%f%f');
+neid = length(elon);
+nrec = length(rlon);
+
+%------------------------------------------
+% write files
+
+for ii = 1:neid % loop over events
+ ename = eid_name{ii};
+ elon0 = elon(ii);
+ elat0 = elat(ii);
+ disp(sprintf('Event %s : %i out of %i',ename,ii,neid));
+
+ % distances and azimuths to all receivers
+ [dists, azis] = distance(repmat(elat0,nrec,1),repmat(elon0,nrec,1),rlat,rlon);
+ dists_km = dists*pi/180*6371;
+
+ [junk1, idist] = sort(dists);
+ [junk2, iaz] = sort(azis);
+
+ % write files
+ fid1 = fopen([odir 'STATIONS_by_dist_from_' ename],'w');
+ fid2 = fopen([odir 'STATIONS_by_az_from_' ename],'w');
+ for jj = 1:nrec
+ k1 = idist(jj);
+ k2 = iaz(jj);
+ fprintf(fid1,stfmt,...
+ rec_name{k1},rlon(k1),rlat(k1),ename,elon0,elat0,dists_km(k1),azis(k1));
+ fprintf(fid2,stfmt,...
+ rec_name{k2},rlon(k2),rlat(k2),ename,elon0,elat0,dists_km(k2),azis(k2));
+ end
+ fclose(fid1);
+ fclose(fid2);
+end
+
+for jj = 1:nrec % loop over receivers
+ rname = rec_name{jj};
+ rlon0 = rlon(jj);
+ rlat0 = rlat(jj);
+ disp(sprintf('Receiver %s : %i out of %i',rname,jj,nrec));
+
+ % distances and (back-)azimuths to all events
+ [dists, azis] = distance(repmat(rlat0,neid,1),repmat(rlon0,neid,1),elat,elon);
+ dists_km = dists*pi/180*6371;
+
+ [junk1, idist] = sort(dists);
+ [junk2, iaz] = sort(azis);
+
+ % write files
+ fid1 = fopen([odir 'EIDS_by_dist_from_' rname],'w');
+ fid2 = fopen([odir 'EIDS_by_az_from_' rname],'w');
+ for ii = 1:neid
+ k1 = idist(ii);
+ k2 = iaz(ii);
+ fprintf(fid1,stfmt,...
+ eid_name{k1},elon(k1),elat(k1),rname,rlon0,rlat0,dists_km(k1),azis(k1));
+ fprintf(fid2,stfmt,...
+ eid_name{k2},elon(k2),elat(k2),rname,rlon0,rlat0,dists_km(k2),azis(k2));
+ end
+ fclose(fid1);
+ fclose(fid2);
+end
+
+%====================================================================
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -1,33 +1,36 @@
#!/usr/bin/perl -w
#-----------------------------------
-# Carl Tape, 15-Dec-2007
+# Carl Tape, 04-Oct-2008
# tomo_make_figs.pl
#
-# This script makes the event kernels by calling vtk scripts.
+# This script makes the event kernel plots by calling vtk scripts.
#
-# EXAMPLE:
-# ~/UTILS/tomo_make_figs.pl 1 1 all m03
-# ~/UTILS/tomo_make_figs.pl 1 0 all m03 (make mesh files only)
+# EXAMPLE (unsmoothed kernels):
+# ~/UTILS/tomo_make_figs.pl 1 1 all m10 0 0
+# ~/UTILS/tomo_make_figs.pl 1 0 all m10 0 0 (make vtu files only)
#
+# EXAMPLE (smoothed kernels):
+# ~/UTILS/tomo_make_figs.pl 1 1 all m10 6 1
+# ~/UTILS/tomo_make_figs.pl 1 0 all m10 6 1 (make vtu files only)
+#
#-----------------------------------
-if (@ARGV < 4) {die("Usage: tomo_make_figs.pl imesh ifig ftag smodel\n")}
-($imesh,$ifig,$ftag,$smodel) = @ARGV;
+if (@ARGV < 6) {die("Usage: tomo_make_figs.pl imesh ifig ftag smodel hsmooth vsmooth\n")}
+($imesh,$ifig,$ftag,$smodel,$hsmooth,$vsmooth) = @ARGV;
-#$ismooth = 0;
-$ismooth = 1; $gsmooth = 6;
-#$ismooth = 1; $gsmooth = 3;
+# to smooth or not to smooth
+$ismooth = 1;
+if($hsmooth==0 && $vsmooth==0) {$ismooth = 0};
+$stg = sprintf("h%3.3ikm_v%3.3ikm",$hsmooth,$vsmooth);
-$stg = sprintf("%2.2i",$gsmooth);
-
$pwd = $ENV{PWD};
$dir_run = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
if (not -e ${dir_run}) {die("check if dir_run ${dir_run} exist or not\n")}
# directories
-$dir_vtk = "/net/denali/scratch1/carltape/vtk/carl_new"; # to run vtk script
-$dir_ker_out = "/net/sierra/raid1/carltape/results/KERNELS/kernel_${smodel}"; # collect copies of the figures
+$dir_vtk = "/net/denali/scratch1/carltape/vtk/carl_new"; # dir to run vtk script
+$dir_ker_out = "/net/sierra/raid1/carltape/results/KERNELS/kernel_${smodel}"; # dir for output figures
if (not -e ${dir_vtk}) {die("check if dir_vtk ${dir_vtk} exist or not\n")}
# directory containing smoothed event kernels
@@ -39,18 +42,23 @@
if (not -e $dir_ker_out) {die("check if dir_ker_out $dir_ker_out exist or not\n")}
# list of event IDs
-#$file_eids = "/net/sierra/raid1/carltape/results/SOURCES/socal_6/SOCAL_FINAL_CMT_v6_eid";
-$file_eids = "/net/sierra/raid1/carltape/results/KERNELS/kernel_${smodel}/kernels_${smodel}";
+#$file_eids = "/net/sierra/raid1/carltape/results/EID_LISTS/kernels_run_${smodel}";
+$file_eids = "/net/sierra/raid1/carltape/results/SOURCES/socal_09/SOCAL_FINAL_CMT_v09_eid";
#$file_eids = "/net/sierra/raid1/carltape/results/WINDOWS/EIDs_pass_10_plus";
-
if (not -f $file_eids) {die("\n check if $file_eids exists\n")}
open(IN,$file_eids); @eids = <IN>; $nevent = @eids;
-# read in list of kernels TO EXCLUDE
-$kexclude = "/net/sierra/raid1/carltape/results/KERNELS/kernel_${smodel}/kernels_exclude_${smodel}";
-open(IN,"$kexclude"); @kex = <IN>; close(IN);
-$nex = @kex;
+# COLOR SCALE for kernels
+$cmax0 = 2.5e-11; # default value 5.0e-11
+#$cmax0 = 15.0e-11;
+#========================================================
+
+# # read in list of kernels TO EXCLUDE
+# $kexclude = "/net/sierra/raid1/carltape/results/KERNELS/kernel_${smodel}/kernels_exclude_${smodel}";
+# open(IN,"$kexclude"); @kex = <IN>; close(IN);
+# $nex = @kex;
+
if (0==1) {
for ($i = 1; $i <= $nevent; $i = $i+1) {
$eid = $eids[$i-1]; chomp($eid);
@@ -68,7 +76,7 @@
# loop over all events
$imin = 1; $imax = $nevent; # default
#$imin = 1; $imax = 40;
-#$imin = 133; $imax = $imin;
+#$imin = 173; $imax = $imin;
for ($i = $imin; $i <= $imax; $i = $i+1) {
@@ -78,14 +86,14 @@
print "------------------------------\n";
print "$imin to $imax, $i, $eid --\n";
- # check if event has been excluded
- $imatch = 1;
- for ($j = 1; $j <= $nex; $j++) {
- $keid = $kex[$j-1]; chomp($keid);
- if ($eid == $keid) {
- $imatch = $imatch*0.0;
- }
- }
+# # check if event has been excluded
+ $imatch = 1;
+# for ($j = 1; $j <= $nex; $j++) {
+# $keid = $kex[$j-1]; chomp($keid);
+# if ($eid == $keid) {
+# $imatch = $imatch*0.0;
+# }
+# }
if($imatch == 0) {
print "--> event has been excluded\n";
@@ -109,8 +117,8 @@
} else {
# smoothed event kernels
- $mesh_files = "$sdir/$eid/*${stg}km.mesh";
- $vtu_files = "$sdir/$eid/*${stg}km.vtu";
+ $mesh_files = "$sdir/$eid/*${stg}.mesh";
+ $vtu_files = "$sdir/$eid/*${stg}.vtu";
$nmesh0 = 2;
$nvtu0 = 2;
$nmesh = `ls -1 ${mesh_files} | wc | awk '{print \$1}'`; chomp($nmesh);
@@ -160,7 +168,7 @@
$nrec = `wc $srvtk | awk '{print \$1}'` - 6;
#if (1==1) {
- if ($nfig == 1) {
+ if ($nfig >= 1) {
print "--> PDF files are already done\n";
} else {
@@ -175,7 +183,6 @@
# color scale value
$cfile = "${dir_mesh}/cmax_kernel";
- $cmax0 = 5.0e-11; # default value
`echo $cmax0 > $cfile`;
#if(not -f $cfile) {
# $cmax0 = 1.0e-11; # default value
@@ -210,10 +217,10 @@
`rm ${dir_vtk}/*.vtk`;
`cp ${dir_output}/*.vtk ${dir_vtk}`;
- #`${dir_vtk}/view_kernel_smooth.pl $eid 1 $smodel $ftag $elab $gsmooth`; # kappa kernel
+ #`${dir_vtk}/view_kernel_smooth.pl $eid 1 $smodel $ftag $elab $hsmooth $vsmooth`; # kappa kernel
#`cp ${dir_vtk}/*ps ${dir_vtk}/*pdf ${dir_fig1}`; # kappa kernel
- `${dir_vtk}/view_kernel_smooth.pl $eid 2 $smodel $ftag $elab $gsmooth`; # mu kernel
+ `${dir_vtk}/view_kernel_smooth.pl $eid 2 $smodel $ftag $elab $hsmooth $vsmooth`; # mu kernel
`cp ${dir_vtk}/*ps ${dir_vtk}/*pdf ${dir_fig1}`; # mu kernel
`cp ${dir_vtk}/view_kernel_smooth*pl ${dir_vtk}/view_kernel_smooth*tcl ${dir_fig1}`;
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs_pmax.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs_pmax.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs_pmax.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -1,15 +1,17 @@
#!/usr/bin/perl -w
#-----------------------------------
-# Carl Tape, 15-Dec-2007
+# Carl Tape, 15-Oct-2008
# tomo_make_figs_pmax.pl
#
-# This script makes the event kernels by calling vtk scripts.
+# This script plots model updates from the subspace method
+# by calling vtk scripts.
#
# EXAMPLE:
+# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm10 0.80 0.0 mu_kernel_smooth beta
#
-# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm03 0.80 0.0 mu_kernel_smooth beta_dcov_window
-# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm02 0.80 0.0 mu_kernel_smooth beta_dcov_window
+# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm08 0.80 0.0 mu_kernel_smooth beta
+# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm08 0.40 0.0 kappa_kernel_smooth bulk
#
# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm00 0.1 0.0 mu_kernel_smooth beta
# ~/UTILS/tomo_make_figs_pmax.pl 1 1 dm00 0.1 0.0 mu_kernel_smooth beta
@@ -45,8 +47,8 @@
#$npmax = @clines;
$imin = 1; $imax = $npmax;
-#$imin = 13; $imax = $npmax;
-#$imin = 7; $imax = $imin;
+#$imin = 10; $imax = $npmax;
+#$imin = 5; $imax = $imin;
for ($i = $imin; $i <= $imax; $i++) {
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -9,9 +9,9 @@
# CALLED BY: tomo_make_figs.pl
#
# EXAMPLE:
-# view_kernel.pl 9818433 2 m0 all 133_9818433
-# view_kernel.pl 9818433 2 m0 T02 133_9818433
-# view_kernel.pl 9818433 2 m0 T06 133_9818433
+# view_kernel.pl 9818433 2 m00 all 133_9818433
+# view_kernel.pl 9818433 2 m00 T02 133_9818433
+# view_kernel.pl 9818433 2 m00 T06 133_9818433
#
#---------------------------------
@@ -140,11 +140,12 @@
# depth layers
# Why does depth z = 0 not work, so we must use z = 0.001 ?
- at dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+#@dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+ at dlayers = (-0.25,0.001,2,4,6,8,10,15,20,25,30,35,40);
$Nz = @dlayers;
#$imin = 1; $imax = $Nz; # default
-$imin = 3; $imax = $imin; # testing
+$imin = 4; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
#$zcen = $utm_zmax - ($i-1)*$dzinc;
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_local.tcl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_local.tcl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -40,13 +40,13 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/14263768/m03/MESH_all/kappa_kernel_low_1.vtu
+ kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/14079184/m11/MESH_all/mu_kernel_low_1.vtu
vtkXMLUnstructuredGridReader kReader2
- kReader2 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/14263768/m03/MESH_all/kappa_kernel_low_2.vtu
+ kReader2 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/14079184/m11/MESH_all/mu_kernel_low_2.vtu
vtkXMLUnstructuredGridReader kReader3
- kReader3 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/14263768/m03/MESH_all/kappa_kernel_low_3.vtu
+ kReader3 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/14079184/m11/MESH_all/mu_kernel_low_3.vtu
## When the files are merged, everything crashes.
#
@@ -67,13 +67,13 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange -5.0000e-11 5.0000e-11
+ kMapper1 SetScalarRange -2.5000e-11 2.5000e-11
vtkDataSetMapper kMapper2
kMapper2 SetInput [kReader2 GetOutput]
- kMapper2 SetScalarRange -5.0000e-11 5.0000e-11
+ kMapper2 SetScalarRange -2.5000e-11 2.5000e-11
vtkDataSetMapper kMapper3
kMapper3 SetInput [kReader3 GetOutput]
- kMapper3 SetScalarRange -5.0000e-11 5.0000e-11
+ kMapper3 SetScalarRange -2.5000e-11 2.5000e-11
# data actors
vtkActor kActor1
@@ -85,7 +85,7 @@
# === generate plane cut ====
vtkPlane hrPlane1
- hrPlane1 SetOrigin 385719.8 3823329.2 -5000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -4000.0
hrPlane1 SetNormal 0 0 1
vtkCutter hrCut1
hrCut1 SetInput [kReader1 GetOutput]
@@ -93,13 +93,13 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -5.0000e-11 5.0000e-11
+ hrMapper1 SetScalarRange -2.5000e-11 2.5000e-11
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
vtkPlane hrPlane2
- hrPlane2 SetOrigin 385719.8 3823329.2 -5000.0
+ hrPlane2 SetOrigin 385719.8 3823329.2 -4000.0
hrPlane2 SetNormal 0 0 1
vtkCutter hrCut2
hrCut2 SetInput [kReader2 GetOutput]
@@ -107,13 +107,13 @@
vtkDataSetMapper hrMapper2
hrMapper2 SetInput [hrCut2 GetOutput]
hrMapper2 InterpolateScalarsBeforeMappingOn
- hrMapper2 SetScalarRange -5.0000e-11 5.0000e-11
+ hrMapper2 SetScalarRange -2.5000e-11 2.5000e-11
hrMapper2 SetLookupTable lut
vtkActor hrActor2
hrActor2 SetMapper hrMapper2
vtkPlane hrPlane3
- hrPlane3 SetOrigin 385719.8 3823329.2 -5000.0
+ hrPlane3 SetOrigin 385719.8 3823329.2 -4000.0
hrPlane3 SetNormal 0 0 1
vtkCutter hrCut3
hrCut3 SetInput [kReader3 GetOutput]
@@ -121,7 +121,7 @@
vtkDataSetMapper hrMapper3
hrMapper3 SetInput [hrCut3 GetOutput]
hrMapper3 InterpolateScalarsBeforeMappingOn
- hrMapper3 SetScalarRange -5.0000e-11 5.0000e-11
+ hrMapper3 SetScalarRange -2.5000e-11 2.5000e-11
hrMapper3 SetLookupTable lut
vtkActor hrActor3
hrActor3 SetMapper hrMapper3
@@ -267,7 +267,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "14263768 : m03 kernel for BULK MODULUS -- 35 stations -- Cut at z = 5.0 km"
+ titleActor SetInput "14079184 : m11 kernel for SHEAR MODULUS -- 105 stations -- Cut at z = 4.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
@@ -284,8 +284,8 @@
iren SetRenderWindow renWin
vtkCamera cam1
- cam1 SetPosition 385719.8 3823329.2 45000.0
- cam1 SetFocalPoint 385719.8 3823329.2 -5000.0
+ cam1 SetPosition 385719.8 3823329.2 46000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -4000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -315,7 +315,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "177_14263768_kappa_all_kernel_03.ps"
+ writer SetFileName "117_14079184_mu_all_kernel_04.ps"
writer Write
wm withdraw .
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -9,15 +9,17 @@
# CALLED BY: tomo_make_figs.pl
#
# EXAMPLE:
-# view_kernel_smooth.pl 9818433 2 m0 all 133_9818433 3
+# view_kernel_smooth.pl 14095628 2 m07 all 049_14095628 6 1
#
#---------------------------------
use File::Basename;
-if (@ARGV < 5) {die("Usage: view_kernel_smooth.pl eid iker smodel elab\n")}
-($eid,$iker,$smodel,$ftag,$elab,$gsmooth) = @ARGV;
+if (@ARGV < 7) {die("Usage: view_kernel_smooth.pl eid iker smodel elab\n")}
+($eid,$iker,$smodel,$ftag,$elab,$hsmooth,$vsmooth) = @ARGV;
+$stg = sprintf("h%3.3ikm_v%3.3ikm",$hsmooth,$vsmooth);
+
$tcl_tag = "view_kernel_smooth";
$tcl_file = "${tcl_tag}.tcl";
$tcl_file_local ="${tcl_tag}_local.tcl";
@@ -27,7 +29,6 @@
`rm *.pdf *.ps`;
$clabel = "m\\^-3";
-$stg = sprintf("%2.2i",$gsmooth);
# directories
$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
@@ -37,11 +38,11 @@
# check for vtu files
$nvtu0 = 2;
-$nvtu = `ls -1 ${dir_mesh}/*${stg}km.vtu | wc | awk '{print \$1}'`; chomp($nvtu);
+$nvtu = `ls -1 ${dir_mesh}/*${stg}.vtu | wc | awk '{print \$1}'`; chomp($nvtu);
if ( $nvtu != $nvtu0 ) {die("You have $nvtu files, not $nvtu0")}
# smoothing options
-if ( $gsmooth == 0 ) {$stag = ""} else {$stag = "_smooth_${stg}km"}
+if ( $hsmooth == 0 ) {$stag = ""} else {$stag = "_smooth_${stg}"}
# kernel options
@klabs = ("kappa","mu","rho","alpha","beta","rhop");
@@ -50,7 +51,7 @@
$ktitle = $ktitles[$iker-1];
$file_tag = "${dir_ker_lab}\\/${klab}_kernel${stag}";
$tlab = "$eid : $smodel kernel for $ktitle";
-$ftag = "${elab}_${klab}_${ftag}_kernel";
+$filetag = "${elab}_${klab}_${ftag}_kernel_${stg}";
$pwd = basename($ENV{PWD});
@@ -59,8 +60,9 @@
#print "\n -- $nrec -- \n"; die("testing");
# copy cmax value here
-$cmax = 5e-11; # default
-#$cmax = 1e-11; # default
+#$cmax = 15e-11; # default
+#$cmax = 5e-11; # default
+$cmax = 2.5e-11; # default
# corners of the UTM mesh
$utm_xmin = 0.06623919273678*1e6;
@@ -125,11 +127,12 @@
# depth layers
# Why does depth z = 0 not work, so we must use z = 0.001 ?
- at dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+#@dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+ at dlayers = (-0.25,0.001,2,4,6,8,10,15,20,25,30,35,40);
$Nz = @dlayers;
#$imin = 1; $imax = $Nz; # default
-$imin = 3; $imax = $imin; # testing
+$imin = 4; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
#$zcen = $utm_zmax - ($i-1)*$dzinc;
@@ -197,7 +200,7 @@
print SED "/titleActor SetInput/s/SetInput.*\$/SetInput \"$title\"/ \n";
# file name
- $filename = sprintf("${ftag}_%2.2i.ps",$i);
+ $filename = sprintf("${filetag}_%2.2i.ps",$i);
print SED "/writer SetFileName/s/SetFileName.*\$/SetFileName \"$filename\"/ \n";
# close SED file
@@ -213,7 +216,7 @@
#---------------------------
-$ofile = "${ftag}_set.pdf";
+$ofile = "${filetag}_set.pdf";
#system("\\rm ${ofile}");
system("/home/carltape/bin/pdcat -r *.pdf ${ofile}");
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth_local.tcl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth_local.tcl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/SMOOTH_EVENT_KERNELS/m03/9105672/mu_kernel_smooth_06km.vtu
+ kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/SMOOTH_EVENT_KERNELS/m11/14079184/mu_kernel_smooth_h006km_v001km.vtu
# vtkXMLUnstructuredGridReader kReader2
# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
@@ -67,7 +67,7 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange -5.0000e-11 5.0000e-11
+ kMapper1 SetScalarRange -2.5000e-11 2.5000e-11
# vtkDataSetMapper kMapper2
# kMapper2 SetInput [kReader2 GetOutput]
# kMapper2 SetScalarRange 4687.47 5729.13
@@ -85,7 +85,7 @@
# === generate plane cut ====
vtkPlane hrPlane1
- hrPlane1 SetOrigin 385719.8 3823329.2 -5000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -4000.0
hrPlane1 SetNormal 0 0 1
vtkCutter hrCut1
hrCut1 SetInput [kReader1 GetOutput]
@@ -93,7 +93,7 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -5.0000e-11 5.0000e-11
+ hrMapper1 SetScalarRange -2.5000e-11 2.5000e-11
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
@@ -267,7 +267,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "9105672 : m03 kernel for SHEAR MODULUS -- 45 stations -- Cut at z = 5.0 km"
+ titleActor SetInput "14079184 : m11 kernel for SHEAR MODULUS -- 105 stations -- Cut at z = 4.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
@@ -284,8 +284,8 @@
iren SetRenderWindow renWin
vtkCamera cam1
- cam1 SetPosition 385719.8 3823329.2 45000.0
- cam1 SetFocalPoint 385719.8 3823329.2 -5000.0
+ cam1 SetPosition 385719.8 3823329.2 46000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -4000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -315,7 +315,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "062_9105672_mu_all_kernel_03.ps"
+ writer SetFileName "117_14079184_mu_all_kernel_h006km_v001km_04.ps"
writer Write
wm withdraw .
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -10,7 +10,7 @@
# EXAMPLE:
# view_kernel_sum.pl 0 2 m00 4e-9 # un-smoothed mu
# view_kernel_sum.pl 6 2 m00 4e-9 # smoothed mu (6 km)
-# view_kernel_sum.pl 10 2 m004e-9 # smoothed mu (10 km)
+# view_kernel_sum.pl 10 2 m00 4e-9 # smoothed mu (10 km)
# view_kernel_sum.pl 15 2 m00 4e-9 # smoothed mu (15 km)
# view_kernel_sum.pl 20 2 m00 4e-9 # smoothed mu (20 km)
#
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_local.tcl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_local.tcl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -67,7 +67,7 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange -4.0000e-09 4.0000e-09
+ kMapper1 SetScalarRange -2.0000e-09 2.0000e-09
# vtkDataSetMapper kMapper2
# kMapper2 SetInput [kReader2 GetOutput]
# kMapper2 SetScalarRange 4687.47 5729.13
@@ -93,7 +93,7 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -4.0000e-09 4.0000e-09
+ hrMapper1 SetScalarRange -2.0000e-09 2.0000e-09
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -10,10 +10,17 @@
# CALLED BY: tomo_make_figs.pl
#
# EXAMPLE:
-# view_model.pl 5 vs_new
+# view_model.pl 5 vs_m12
+# view_model.pl 5 vs
+# view_model.pl 7 vs_m12_m00
+# view_model.pl 7 vs_m12_m11
# view_model.pl 5 vs_low
# view_model.pl 7 Beta
#
+# view_model.pl 4 vp_m12
+# view_model.pl 4 vp
+# view_model.pl 6 vp_m12_m00
+#
#---------------------------------
use File::Basename;
@@ -24,12 +31,23 @@
# copy files into local directory
#`cp vtu_files/lin_model_2/*vtu .`;
#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m00/vtu_old/*vtu .`;
-#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m00/*vtu .`;
-#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m01/*vtu .`;
-#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m02/*vtu .`;
-#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m03/*vtu .`;
-`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m04/*vtu .`;
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m00/*vtu .`; $title_tag = "Model m00 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m01/*vtu .`; $title_tag = "Model m01 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m02/*vtu .`; $title_tag = "Model m02 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m03/*vtu .`; $title_tag = "Model m03 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m04/*vtu .`; $title_tag = "Model m04 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m05/*vtu .`; $title_tag = "Model m05 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m06/*vtu .`; $title_tag = "Model m06 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m07/*vtu .`; $title_tag = "Model m07 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m08/*vtu .`; $title_tag = "Model m08 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m09/*vtu .`; $title_tag = "Model m09 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m10/*vtu .`; $title_tag = "Model m10 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m11/*vtu .`; $title_tag = "Model m11 for";
+`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m12/*vtu .`; $title_tag = "Model m12 for";
+$title_tag = "ln(m12\\/m11) for"; $cmax = 0.05;
+#$title_tag = "ln(m12\\/m00) for"; $cmax = 0.30;
+
$tcl_tag = "view_model";
$tcl_file = "${tcl_tag}.tcl";
$tcl_file_local ="${tcl_tag}_local.tcl";
@@ -44,13 +62,13 @@
#$dir_ker_lab = "\\/net\\/sierra\\/raid1\\/carltape\\/socal\\/socal_3D\\/RUNS\\/${eid}\\/${smodel}\\/MESH_${ftag}";
# kernel options
- at klabs = ("kappa","mu","rho","vp","vs","Bulk","Beta");
- at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","BULK SOUND","S-WAVE SPEED");
+ at klabs = ("kappa","mu","rho","vp","vs","vp","vs");
+ at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P WAVE-SPEED","S WAVE-SPEED","P WAVE-SPEED","S WAVE-SPEED");
$klab = $klabs[$iker-1];
$ktitle = $ktitles[$iker-1];
#$file_tag = "${klab}_low"; $ftag = "${klab}_low";
$file_tag = $kfile; $ftag = $kfile;
-$tlab = "Model for $ktitle";
+$tlab = "$title_tag $ktitle";
$pwd = basename($ENV{PWD});
@@ -68,28 +86,37 @@
# depth layers
# Why does depth z = 0 not work, so we must use z = 0.001 ?
- at dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+
+#@dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+ at dlayers = (-0.25,0.001,2,4,6,8,10,15,20,25,30,35,40);
#@dlayers = (-0.25,0.001,1,2,3,4,5,6,7,8,9,10); # density plots
#@dlayers = (-0.25,0.001,2.75,10.75,24,46);
$Nz = @dlayers;
+# percent perturbation for each layer
+ at players = (15,15,15,15,15,15,10,10,10,10,10,10,5);
+
$pert = 10;
if($iker == 1) {
- $cmin = 0e10; $cmax = 10e10; $clabel = "Pa";
+ $cmin = 0e10; $cmax = 10e10; $cunits = "Pa";
} elsif($iker == 2) {
- $cmin = 0e10; $cmax = 6e10; $clabel = "Pa";
+ $cmin = 0e10; $cmax = 6e10; $cunits = "Pa";
} elsif($iker == 3) {
- $cmin = 1800; $cmax = 3200; $clabel = "kg m\\^-3";
+ $cmin = 1800; $cmax = 3200; $cunits = "kg m\\^-3";
+
} elsif($iker == 4) {
- $cmin = 1000; $cmax = 6000; $clabel = "m s\\^-1";
- @cmean = (5.2083,5.2083,6.1101,6.2208,6.4695,6.5064,6.5516,6.6604,6.6604,6.6604);
+ $cmin = 1000; $cmax = 6000; $cunits = "m\\/s";
+ @cmean = (5.20,5.20,5.50,5.9,6.11,6.11,6.22,6.46,6.50,6.55,6.7,7.30,7.80);
+
} elsif($iker == 5) {
- $cmin = 1000; $cmax = 5000; $clabel = "m s\\^-1";
- @cmean = (2.9194,2.9194,3.35,3.5586,3.6642,3.6656,3.6844,3.741,3.741,3.741);
+ $cmin = 1000; $cmax = 5000; $cunits = "m\\/s";
+ @cmean = (2.95,2.95,3.1,3.25,3.4,3.45,3.55,3.65,3.7,3.7,3.75,4.2,4.5);
+
} elsif($iker == 6) {
- $cmin = -0.5; $cmax = 0.5; $clabel = " ";
+ $cmin = -$cmax; $cunits = " ";
+
} elsif($iker == 7) {
- $cmin = -0.5; $cmax = 0.5; $clabel = " ";
+ $cmin = -$cmax; $cunits = " ";
}
# corners of the UTM mesh
@@ -143,17 +170,32 @@
$Nz = @dlayers;
$imin = 1; $imax = $Nz; # default
+#$imin = 9; $imax = 13; # testing
#$imin = 3; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
+
#$zcen = $utm_zmax - ($i-1)*$dzinc;
$zcen = -1000 * $dlayers[$i-1];
$zcen_km = $zcen/1000;
printf ("%2.2i : %.1f, %.1f, %.1f\n",$i,$xcen,$ycen,$zcen);
+ $pert = $players[$i-1];
+ $cmid = $cmean[$i-1];
# title for plot
$title = sprintf("%s -- Cut at z = %.1f km",$tlab,-$zcen_km);
+ if($iker <= 5) {
+ $clabel = sprintf("%.0f $cunits +- %.0f percent",$cmid*1000,$pert);
+ $scalar_low = $cmid* 1000 * (1 - $pert/100);
+ $scalar_high = $cmid * 1000 * (1 + $pert/100);
+
+ } else {
+ $clabel = "Percent Change";
+ $scalar_low = sprintf("%.4e",$cmin);
+ $scalar_high = sprintf("%.4e",$cmax);
+ }
+
# open file for replacing the cross-section line
open(SED,">sed.txt");
@@ -163,10 +205,6 @@
# color scale
print SED "/scalarBar SetTitle/s/SetTitle.*\$/SetTitle \"$clabel\"/ \n";
- $scalar_low = $cmean[$i-1] * 1000 * (1 - $pert/100);
- $scalar_high = $cmean[$i-1] * 1000 * (1 + $pert/100);
- #$scalar_low = sprintf("%.4e",$cmin);
- #$scalar_high = sprintf("%.4e",$cmax);
print SED "/kMapper1 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
print SED "/hrMapper1 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model_local.tcl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model_local.tcl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName vs_new.vtu
+ kReader1 SetFileName vs_m12_m11.vtu
vtkXMLUnstructuredGridReader kReader2
kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
@@ -67,7 +67,7 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange 3366.9 4115.1
+ kMapper1 SetScalarRange -5.0000e-02 5.0000e-02
# data actors
vtkActor kActor1
@@ -83,7 +83,7 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange 3366.9 4115.1
+ hrMapper1 SetScalarRange -5.0000e-02 5.0000e-02
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
@@ -160,7 +160,7 @@
vtkScalarBarActor scalarBar
scalarBar SetLookupTable lut
#scalarBar SetLookupTable [hrMapper GetLookupTable]
-scalarBar SetTitle "m s^-1"
+scalarBar SetTitle "Percent Change"
[scalarBar GetPositionCoordinate] SetCoordinateSystemToNormalizedViewport
[scalarBar GetPositionCoordinate] SetValue 0.1 0.0
scalarBar SetOrientationToHorizontal
@@ -176,7 +176,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "Model for S-WAVE-SPEED -- Cut at z = 40.0 km"
+ titleActor SetInput "ln(m12/m11) for S WAVE-SPEED -- Cut at z = 40.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
@@ -220,7 +220,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "vs_new_10.ps"
+ writer SetFileName "vs_m12_m11_13.ps"
writer Write
wm withdraw .
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update.pl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update.pl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -3,6 +3,7 @@
#---------------------------------
# view_update.pl
# Carl Tape and Qinya Liu
+# July 29, 2008
#
# This script creates a PDF of cross-sections of a kernel.
#
@@ -18,6 +19,11 @@
# view_update.pl 1 dm02 31 0.40 16 dm_kappa_kernel_smooth_p031_smooth_16km bulk_smooth_16km
# view_update.pl 1 dm03 25 0.80 0 dm_kappa_kernel_smooth_p025 bulk
# view_update.pl 1 dm03 25 0.80 10 dm_kappa_kernel_smooth_p025_smooth_10km bulk_smooth_10km
+# view_update.pl 1 dm04 40 0.80 0 dm_kappa_kernel_smooth_p040 bulk
+# view_update.pl 1 dm05 44 0.50 0 dm_kappa_kernel_smooth_p044 bulk
+# view_update.pl 1 dm05 44 0.50 10 dm_kappa_kernel_smooth_p044_smooth_10km bulk_smooth_10km
+# view_update.pl 1 dm06 48 0.50 0 dm_kappa_kernel_smooth_p048 bulk
+# view_update.pl 1 dm07 44 0.50 10 dm_kappa_kernel_smooth_p044_smooth_h006km_v001km bulk_smooth_h006km_v001km
#
# view_update.pl 2 dm00 68 0.10 0 dm_mu_kernel_smooth_p068 beta
# view_update.pl 2 dm00 68 0.10 6 dm_mu_kernel_smooth_p068_smooth_06km beta_smooth_06km
@@ -29,6 +35,15 @@
# view_update.pl 2 dm03 25 0.80 0 dm_mu_kernel_smooth_p025 beta
# view_update.pl 2 dm03 25 0.80 10 dm_mu_kernel_smooth_p025_smooth_10km beta_smooth_10km
# view_update.pl 2 dm03 25 2.00 10 dm_mu_kernel_smooth_p025_smooth_10km beta_smooth_10km
+# view_update.pl 2 dm04 40 1.40 0 dm_mu_kernel_smooth_p040 beta
+# view_update.pl 2 dm04 40 2.00 10 dm_mu_kernel_smooth_p040_smooth_10km beta_smooth_10km
+# view_update.pl 2 dm05 44 1.00 0 dm_mu_kernel_smooth_p044 beta
+# view_update.pl 2 dm05 44 0.80 10 dm_mu_kernel_smooth_p044_smooth_10km beta_smooth_10km
+# view_update.pl 2 dm06 48 1.00 0 dm_mu_kernel_smooth_p048 beta
+# view_update.pl 2 dm06 48 0.80 10 dm_mu_kernel_smooth_p048_smooth_h006km_v001km beta_smooth_h006km_v001km
+# view_update.pl 2 dm07 44 1.00 0 dm_mu_kernel_smooth_p044 beta
+# view_update.pl 2 dm10 28 1.60 0 dm_mu_kernel_smooth_p028 beta
+# view_update.pl 2 dm11 40 1.60 0 dm_mu_kernel_smooth_p040 beta
#
# view_update.pl 2 dm00 1 0.05 6 mu_kernel_smooth_dm beta_cg_smooth_06km
# view_update.pl 1 dm00 1 0.05 6 kappa_kernel_smooth_dm bulk_cg_smooth_06km
@@ -127,22 +142,14 @@
# SetViewUp
$ux = 0; $uy = 1; $uz = 0;
-# 0 2.8109 5.2083 6.2289 1.4284 2.9194 3.7731
-# 3.0000 3.7365 5.6315 6.7879 1.8581 3.1706 3.8945
-# 6.0000 4.4321 6.1101 7.3782 2.5197 3.4845 4.1260
-# 10.0000 5.3617 6.2208 6.8249 3.0108 3.5586 4.3172
-# 15.0000 6.0794 6.4695 7.1689 3.3740 3.6642 4.1778
-# 17.0000 6.2422 6.5064 6.8832 3.5313 3.6656 4.1101
-# 22.0000 6.4608 6.5516 6.8584 3.5813 3.6844 3.9630
-# 31.0000 6.6050 6.6604 6.7609 3.7213 3.7419 3.8484
-
# depth layers
# Why does depth z = 0 not work, so we must use z = 0.001 ?
- at dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+#@dlayers = (-0.25,0.001,5,10,15,20,25,30,35,40);
+ at dlayers = (-0.25,0.001,2,4,6,8,10,15,20,25,30,35,40);
$Nz = @dlayers;
$imin = 1; $imax = $Nz; # default
-#$imin = 3; $imax = $imin; # testing (3)
+#$imin = 4; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
#$zcen = $utm_zmax - ($i-1)*$dzinc;
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update_local.tcl 2008-10-21 02:56:11 UTC (rev 13108)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update_local.tcl 2008-10-21 05:14:27 UTC (rev 13109)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/dm03/beta_smooth_10km/dm_mu_kernel_smooth_p025_smooth_10km.vtu
+ kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/dm11/beta/dm_mu_kernel_smooth_p040.vtu
# vtkXMLUnstructuredGridReader kReader2
# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
@@ -67,7 +67,7 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange -2.0000e+00 2.0000e+00
+ kMapper1 SetScalarRange -1.6000e+00 1.6000e+00
# vtkDataSetMapper kMapper2
# kMapper2 SetInput [kReader2 GetOutput]
# kMapper2 SetScalarRange 4687.47 5729.13
@@ -93,7 +93,7 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -2.0000e+00 2.0000e+00
+ hrMapper1 SetScalarRange -1.6000e+00 1.6000e+00
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
@@ -267,7 +267,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "Subspace update for SHEAR-WAVE-SPEED -- model 25 -- Cut at z = 40.0 km"
+ titleActor SetInput "Subspace update for SHEAR-WAVE-SPEED -- model 40 -- Cut at z = 40.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
@@ -315,7 +315,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "dm_mu_kernel_smooth_p025_smooth_10km_10.ps"
+ writer SetFileName "dm_mu_kernel_smooth_p040_13.ps"
writer Write
wm withdraw .
More information about the CIG-COMMITS
mailing list