[cig-commits] r14494 - in seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS: . station_lists vtk
carltape at geodynamics.org
carltape at geodynamics.org
Fri Mar 27 13:39:59 PDT 2009
Author: carltape
Date: 2009-03-27 13:39:58 -0700 (Fri, 27 Mar 2009)
New Revision: 14494
Added:
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.tcl
Modified:
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in
seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/make_station_lists.m
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_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.tcl
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:
Updated vtk plotting scripts for viewing tomographic models and kernels.
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/EIDS_in 2009-03-27 20:39:58 UTC (rev 14494)
@@ -1,9 +1,119 @@
-9105672 -116.8417 34.3228
+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
+9968977 -120.8385 35.5487
+14018508 -121.1372 35.7032
+9982749 -120.8740 35.5918
+9983625 -120.8693 35.5800
+10005209 -120.8028 35.5067
+14096736 -120.8108 35.5473
+14189556 -121.0838 35.6500
+14263252 -120.7510 35.6360
+13965956 -121.1007 36.5565
+14094528 -120.6661 36.1434
+14094992 -120.3667 35.8155
+14094996 -120.3531 35.8029
+14095540 -120.5134 35.9528
+14096196 -120.5403 35.9821
+10061489 -120.4113 35.8570
+10063349 -120.4963 35.9437
+10100053 -120.4792 35.9269
+9700217 -120.2220 36.2250
+14155056 -120.1020 36.0720
+10222697 -120.2430 36.1860
+10222753 -120.2640 36.1950
+9171679 -119.0298 34.8942
+9983429 -119.1412 35.0118
+14138080 -119.1940 34.9987
+10097009 -119.1958 35.0023
+14186612 -119.0247 35.0178
+14186928 -119.0367 35.0198
+14421592 -119.4320 35.3040
+9094270 -118.3055 35.0944
+9151000 -118.3073 35.0942
+9875657 -118.6677 35.3185
+9875665 -118.6636 35.3102
+9882325 -118.6632 35.3152
+9882329 -118.6585 35.3128
+14095628 -118.6292 35.3852
+14187364 -118.4718 35.3782
+9095528 -118.4718 35.7377
+9150885 -118.0484 36.3269
+3324419 -118.0519 36.3249
+9151609 -118.0507 36.3267
+9644345 -118.3256 35.9817
+9653293 -118.0397 35.7982
+9653349 -118.0407 35.7973
+9653493 -118.0423 35.7957
+12887732 -118.0758 35.7057
+9915909 -118.1038 35.8388
+13986104 -118.2692 36.4782
+9994573 -118.1589 36.1815
+14169456 -118.0652 36.1488
+9044494 -117.6405 36.0778
+3298170 -117.6253 36.0923
+9044650 -117.6505 36.0737
+9045109 -117.6200 36.0912
+9045697 -117.6133 36.0827
+9116921 -117.9045 36.4512
+7179710 -117.6053 36.0882
+9141142 -117.6010 36.0850
+7180136 -117.5993 36.0888
+9163702 -117.8712 36.0392
+9642941 -117.8131 36.0145
+9646589 -117.8868 36.0222
+9671933 -117.8722 36.0293
+10964587 -117.8735 36.0223
+9673577 -117.8742 36.0153
+9674049 -117.8682 36.0135
+10970835 -117.8693 36.0216
+9674093 -117.8679 36.0232
+9674097 -117.8660 36.0392
+9674205 -117.8755 36.0247
+9674213 -117.8753 36.0140
+9674653 -117.8808 36.0268
+10992159 -117.8723 35.9915
+11671240 -117.8650 35.9783
+9686565 -117.8632 36.0495
+9688025 -117.8781 36.0394
+9688709 -117.8720 35.9896
+9828889 -117.2913 35.9468
+9829213 -117.2917 35.9490
+10023841 -117.8910 36.0691
+9152038 -117.6532 35.7822
+9165019 -117.6462 35.7796
+9171064 -117.1652 35.5801
+9631385 -117.5195 35.1143
+14007388 -117.5664 35.6352
+14204000 -117.5450 35.1267
+14215316 -117.5402 35.1228
+14219360 -117.5828 35.6232
+14418600 -117.7860 35.4130
+9058934 -117.5510 36.9673
+9059586 -117.5420 36.9400
+9087073 -118.8180 37.5300
+9090617 -118.8040 37.5270
+9098566 -116.9680 37.5300
+3319204 -117.0890 37.3960
+9098867 -117.0650 37.4170
+9106250 -117.0430 37.4350
+9689717 -117.7910 37.2190
+13657604 -116.3193 36.7103
+9966033 -117.8488 37.1433
9108606 -116.2694 34.5965
3320736 -116.2465 34.4368
-9109243 -115.7045 33.2806
+9109131 -116.3112 34.6862
+9109254 -116.2463 34.4245
+9109287 -116.3013 34.7087
+9109442 -116.2793 34.6940
+9109496 -116.3390 34.6786
+9109636 -116.3570 34.7097
9109752 -116.0675 34.2855
9110281 -116.2527 34.5588
+9110685 -116.1375 34.3465
9111353 -116.1448 34.3512
9112735 -116.3442 34.7110
9113909 -116.3952 34.8632
@@ -11,157 +121,144 @@
3321595 -116.3898 34.8648
3321590 -116.3955 34.8735
3320884 -116.3978 34.8633
+9114042 -116.2663 34.5237
+9114612 -116.3597 34.7005
3324595 -116.2639 34.6055
9114763 -116.2085 34.3300
3321597 -116.1944 34.4042
+9114775 -116.4070 34.8566
9114812 -116.4060 34.8620
-9116921 -117.9045 36.4512
+9114858 -116.3820 34.8292
+3320951 -116.3895 34.8297
+3320940 -116.3947 34.8660
+3320954 -116.3943 34.8630
9117942 -116.2707 34.5200
3321426 -116.2888 34.8031
+9119414 -116.3570 34.8470
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
+9775765 -116.2952 34.5177
+9805021 -116.4447 34.5660
+9854597 -116.2650 34.8068
+13945908 -116.1303 34.3582
+9930549 -116.6665 34.6172
+14408052 -116.4190 34.8130
+13978600 -121.1410 34.7150
+13978612 -121.1133 34.7002
+9968525 -120.4817 34.7847
+12659440 -119.3317 33.6678
+10006857 -120.0142 34.4135
+10065241 -119.6955 34.2883
+14139108 -120.0295 33.6922
+14139160 -120.0258 33.6880
+14165408 -119.7527 33.6853
+7210945 -118.0680 32.7260
+9695397 -118.2513 32.7920
+9695549 -118.3180 32.7340
+14065544 -117.8488 32.4238
+10148829 -118.1450 32.4970
+9096972 -118.6105 34.3935
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
+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
+9941081 -118.6509 34.4034
+14000376 -118.7402 34.2910
+14077668 -119.4365 34.3885
+9038699 -117.7178 33.9462
+9064568 -117.6502 34.3737
+9093975 -118.2253 34.0097
+9644101 -117.7091 33.8725
+9703873 -118.3968 34.0527
+9716853 -118.2795 33.9297
+9735129 -117.7483 33.9552
+9818433 -117.7840 33.9133
+10094253 -118.6268 33.8280
+14383980 -117.7610 33.9530
+3298292 -117.2222 34.0355
+9064093 -116.9232 34.1245
+7112721 -116.9158 34.1155
+9069997 -116.8418 34.3208
+9070083 -116.8455 34.3202
+9105672 -116.8417 34.3228
+9128775 -117.0072 34.1048
+9132433 -117.0028 34.1040
+9140050 -117.2432 34.0588
+9151375 -116.8683 34.2906
+9169867 -116.7722 34.2673
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
+9627953 -116.9370 34.2930
+9652545 -117.4415 34.2269
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
+13936596 -116.8445 34.3100
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
+13939856 -116.8428 34.3040
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
+14155260 -117.0072 34.0612
+14158696 -117.0232 34.0615
+10148369 -116.7725 34.0198
+10148421 -116.7715 34.0182
+10187953 -116.7903 33.9217
+14239184 -117.1103 33.8567
+10370141 -117.3040 34.1070
+9085734 -116.3697 34.0725
+9086693 -116.3623 34.0375
+3317364 -116.3582 34.0378
+9096656 -116.7113 33.6370
+9627557 -116.1394 33.8170
+9700049 -116.3986 33.3807
+9718013 -116.5023 33.5112
+9742277 -116.4308 33.3852
+9774569 -116.7117 33.2065
+13813696 -116.7634 33.5023
+9853417 -116.5685 33.5125
+9915709 -116.0215 33.7508
+14073800 -116.0520 33.7152
14118096 -116.3912 33.9578
+14151344 -116.5675 33.5380
+10223765 -116.0448 33.7063
+13303428 -115.9762 33.1627
+9826789 -116.1118 33.2352
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
+14183744 -116.0260 33.1787
+14236768 -116.0220 33.2450
+14255632 -116.0632 33.2663
+10230869 -116.1357 33.2220
+9109243 -115.7045 33.2806
+9148510 -115.6349 33.1579
+9150059 -115.6362 33.1573
+9152745 -115.5988 33.1751
+9154092 -115.5035 32.8898
+9154179 -115.5082 32.8725
+9154233 -115.5060 32.9034
+9722529 -115.7100 33.3070
+9722633 -115.7012 33.3170
+9722669 -115.7045 33.3135
+9817605 -115.6173 33.0366
+13966396 -115.5538 32.9475
+13966672 -115.5409 32.9455
+13970876 -115.5472 32.9443
14178184 -115.6207 33.1544
14178188 -115.6098 33.1639
14178212 -115.6157 33.1548
@@ -170,35 +267,28 @@
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
+9075784 -115.9196 32.7260
+9075803 -115.9248 32.7190
+12456160 -115.7451 32.5553
+14072464 -115.7441 32.5392
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
+9146641 -115.3868 32.7020
+9158503 -115.1060 32.1622
+9164821 -115.1458 32.1727
+9660449 -115.4224 32.7416
+9744905 -115.0875 32.4633
+9504418 -115.0893 32.4723
+9827109 -115.3773 32.3042
+13917260 -115.7883 32.2498
+9944301 -115.2837 32.5620
+14215812 -115.5730 32.5678
+14137160 -116.8122 32.7233
+14181056 -116.8393 32.5112
Modified: 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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/station_lists/make_station_lists.m 2009-03-27 20:39:58 UTC (rev 14494)
@@ -1,6 +1,6 @@
%
% make_station_lists.m
-% CARL TAPE, 01-Oct-2008
+% CARL TAPE, 20-Jan-2009
% printed xxx
%
% This program reads in a list of receivers and a list of events, and it
@@ -22,11 +22,15 @@
%------------------------------------------
% read in files
+% Here are some example commands to make the input files:
+% awk '{print $2,$8,$9}' /net/sierra/raid1/carltape/results/SOURCES/socal_16/EIDs_lonlat_loc > EIDS_in
+% awk '{print $3"."$4,$1,$2}' /home/carltape/gmt/stations/seismic/Matlab_output/STATIONS_CALIFORNIA_TOMO_INNER_gmt > STATIONS_ini
[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);
+disp(sprintf('%i sources and %i receivers',neid,nrec));
%------------------------------------------
% write files
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs.pl 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs.pl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -7,12 +7,12 @@
# This script makes the event kernel plots by calling vtk scripts.
#
# EXAMPLE (unsmoothed kernels):
-# ~/UTILS/tomo_make_figs.pl 1 1 all m13 0 0
-# ~/UTILS/tomo_make_figs.pl 1 0 all m13 0 0 (make vtu files only)
+# ~/UTILS/tomo_make_figs.pl 1 1 all m16 0 0
+# ~/UTILS/tomo_make_figs.pl 1 0 all m16 0 0 (make vtu files only)
#
# EXAMPLE (smoothed kernels):
-# ~/UTILS/tomo_make_figs.pl 1 1 all m13 6 1
-# ~/UTILS/tomo_make_figs.pl 1 0 all m13 6 1 (make vtu files only)
+# ~/UTILS/tomo_make_figs.pl 1 1 all m16 6 1
+# ~/UTILS/tomo_make_figs.pl 1 0 all m16 6 1 (make vtu files only)
#
#-----------------------------------
@@ -29,7 +29,7 @@
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"; # dir to run vtk script
+$dir_vtk = "/net/denali/raid1/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")}
@@ -43,7 +43,7 @@
# list of event IDs
#$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/SOURCES/socal_16/SOCAL_FINAL_CMT_v16_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;
@@ -75,8 +75,8 @@
# loop over all events
$imin = 1; $imax = $nevent; # default
-#$imin = 1; $imax = 40;
-#$imin = 173; $imax = $imin;
+#$imin = 21; $imax = $nevent;
+#$imin = 29; $imax = $imin;
for ($i = $imin; $i <= $imax; $i = $i+1) {
@@ -102,6 +102,7 @@
# directories
$dir_run_eid = "${dir_run}/${eid}/${smodel}";
+ #$dir_run_eid = "${dir_run}/${eid}/${smodel}_BDK_SAVE"; # TEMPORARY
$dir_mesh = "${dir_run_eid}/MESH_${ftag}";
$dir_output = "${dir_run_eid}/OUTPUT_${ftag}";
@@ -167,7 +168,7 @@
$srvtk = "${dir_output}/sr.vtk";
$nrec = `wc $srvtk | awk '{print \$1}'` - 6;
- #if (1==1) {
+ #if (0==1) {
if ($nfig >= 1) {
print "--> PDF files are already done\n";
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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/tomo_make_figs_pmax.pl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -8,7 +8,7 @@
# by calling vtk scripts.
#
# EXAMPLE:
-# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm13 0.80 0.0 mu_kernel_smooth beta_window
+# ~/UTILS/tomo_make_figs_pmax.pl 1 0 dm15 0.80 0.0 mu_kernel_smooth beta_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
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel.pl 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel.pl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -39,9 +39,12 @@
# directories
$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
-$dir_output = "$dir0/${eid}/${smodel}/OUTPUT_${ftag}";
-$dir_mesh = "$dir0/${eid}/${smodel}/MESH_${ftag}";
+$dir_smodel = "$dir0/${eid}/${smodel}";
+#$dir_smodel = "$dir0/${eid}/${smodel}_BDK_SAVE"; # TEMPORARY
+$dir_output = "${dir_smodel}/OUTPUT_${ftag}";
+$dir_mesh = "${dir_smodel}/MESH_${ftag}";
$dir_ker_lab = "\\/net\\/sierra\\/raid1\\/carltape\\/socal\\/socal_3D\\/RUNS\\/${eid}\\/${smodel}\\/MESH_${ftag}";
+#$dir_ker_lab = "\\/net\\/sierra\\/raid1\\/carltape\\/socal\\/socal_3D\\/RUNS\\/${eid}\\/${smodel}_BDK_SAVE\\/MESH_${ftag}";
# check for vtu files
$nvtu0 = 6;
@@ -76,7 +79,7 @@
$cmax = 5e-11; # defalt
}
-#$cmax = 1e-11;
+$cmax = 3e-13;
# corners of the UTM mesh
$utm_xmin = 0.06623919273678*1e6;
@@ -144,9 +147,9 @@
@dlayers = (-0.25,0.001,2,4,6,8,10,15,20,25,30,35,40);
$Nz = @dlayers;
-$imin = 1; $imax = $Nz; # default
-$imin = 1; $imax = 10; # testing
-$imin = 4; $imax = $imin; # testing
+#$imin = 1; $imax = $Nz; # default
+$imin = 2; $imax = $Nz; # 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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_local.tcl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -40,13 +40,13 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9716853/m13/MESH_all/mu_kernel_low_1.vtu
+ kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9983429/m16/MESH_all/mu_kernel_low_1.vtu
vtkXMLUnstructuredGridReader kReader2
- kReader2 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9716853/m13/MESH_all/mu_kernel_low_2.vtu
+ kReader2 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9983429/m16/MESH_all/mu_kernel_low_2.vtu
vtkXMLUnstructuredGridReader kReader3
- kReader3 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9716853/m13/MESH_all/mu_kernel_low_3.vtu
+ kReader3 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9983429/m16/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 -2.5000e-11 2.5000e-11
+ kMapper1 SetScalarRange -3.0000e-13 3.0000e-13
vtkDataSetMapper kMapper2
kMapper2 SetInput [kReader2 GetOutput]
- kMapper2 SetScalarRange -2.5000e-11 2.5000e-11
+ kMapper2 SetScalarRange -3.0000e-13 3.0000e-13
vtkDataSetMapper kMapper3
kMapper3 SetInput [kReader3 GetOutput]
- kMapper3 SetScalarRange -2.5000e-11 2.5000e-11
+ kMapper3 SetScalarRange -3.0000e-13 3.0000e-13
# data actors
vtkActor kActor1
@@ -85,7 +85,7 @@
# === generate plane cut ====
vtkPlane hrPlane1
- hrPlane1 SetOrigin 385719.8 3823329.2 -4000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -40000.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 -2.5000e-11 2.5000e-11
+ hrMapper1 SetScalarRange -3.0000e-13 3.0000e-13
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
vtkPlane hrPlane2
- hrPlane2 SetOrigin 385719.8 3823329.2 -4000.0
+ hrPlane2 SetOrigin 385719.8 3823329.2 -40000.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 -2.5000e-11 2.5000e-11
+ hrMapper2 SetScalarRange -3.0000e-13 3.0000e-13
hrMapper2 SetLookupTable lut
vtkActor hrActor2
hrActor2 SetMapper hrMapper2
vtkPlane hrPlane3
- hrPlane3 SetOrigin 385719.8 3823329.2 -4000.0
+ hrPlane3 SetOrigin 385719.8 3823329.2 -40000.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 -2.5000e-11 2.5000e-11
+ hrMapper3 SetScalarRange -3.0000e-13 3.0000e-13
hrMapper3 SetLookupTable lut
vtkActor hrActor3
hrActor3 SetMapper hrMapper3
@@ -267,7 +267,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "9716853 : m13 kernel for SHEAR MODULUS -- 85 stations -- Cut at z = 4.0 km"
+ titleActor SetInput "9983429 : m16 kernel for SHEAR MODULUS -- 1 stations -- Cut at z = 40.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 46000.0
- cam1 SetFocalPoint 385719.8 3823329.2 -4000.0
+ cam1 SetPosition 385719.8 3823329.2 10000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -40000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -315,7 +315,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "132_9716853_mu_all_kernel_04.ps"
+ writer SetFileName "029_9983429_mu_all_kernel_13.ps"
writer Write
wm withdraw .
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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_smooth_local.tcl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/SMOOTH_EVENT_KERNELS/m13/9716853/mu_kernel_smooth_h006km_v001km.vtu
+ kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/SMOOTH_EVENT_KERNELS/m16/9944301/mu_kernel_smooth_h006km_v001km.vtu
# vtkXMLUnstructuredGridReader kReader2
# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
@@ -267,7 +267,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "9716853 : m13 kernel for SHEAR MODULUS -- 85 stations -- Cut at z = 4.0 km"
+ titleActor SetInput "9944301 : m16 kernel for SHEAR MODULUS -- 80 stations -- Cut at z = 4.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 "132_9716853_mu_all_kernel_h006km_v001km_04.ps"
+ writer SetFileName "291_9944301_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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.pl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -5,65 +5,60 @@
# Carl Tape, after Qinya Liu
#
# This script creates a PDF of cross-sections of a kernel.
-# Execute from /net/denali/scratch1/carltape/vtk/carl_new/
#
+# CALLED BY: tomo_make_figs.pl
+#
# 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 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)
+# view_kernel_sum.pl m16 sum_abs kappa_kernel_smooth 1 5e-10
+# view_kernel_sum.pl m16 sum_abs mu_kernel_smooth 2 5e-10
+# view_kernel_sum.pl m16 sum_abs_nwin kappa_kernel_smooth 1 2e-12
+# view_kernel_sum.pl m16 sum_abs_nwin mu_kernel_smooth 2 2e-12
#
-# view_kernel_sum.pl 0 1 m00 4e-9 # un-smoothed kappa
-# view_kernel_sum.pl 6 1 m00 4e-9 # smoothed kappa (6 km)
-# view_kernel_sum.pl 10 1 m00 4e-9 # smoothed kappa (10 km)
-# view_kernel_sum.pl 15 1 m00 4e-9 # smoothed kappa (15 km)
-# view_kernel_sum.pl 20 1 m00 4e-9 # smoothed kappa (20 km)
-#
#---------------------------------
use File::Basename;
-if (@ARGV < 4) {die("Usage: view_kernel_sum.pl ismooth iker smodel cmax\n")}
-($gsmooth,$iker,$smodel,$cmax) = @ARGV;
+if (@ARGV < 5) {die("Usage: view_kernel_sum.pl XXX\n")}
+($smodel,$dtag,$ktag,$iker,$cmax) = @ARGV;
$tcl_tag = "view_kernel_sum";
$tcl_file = "${tcl_tag}.tcl";
$tcl_file_local ="${tcl_tag}_local.tcl";
if (not -f $tcl_file) {die("Check if sample file $tcl_file exists or not\n");}
+# remove figures and vtk files in local directory
+`rm *.pdf *.ps`;
+
+# KEY COMMANDS
+#$eid = 9818433;
+#$smodel = "m0";
+#$iker = 2; # kappa, mu, rho, alpha, beta, rhop
+#$cmax = 6e-11; # color scale
+#$cmax = 3e-10; # color scale
+#$clabel = "s\\^2 m\\^-3";
$clabel = "m\\^-3";
-$stg = sprintf("%2.2i",$gsmooth);
# directories
-$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS/GRADIENT/${smodel}";
-$dir_output = "$dir0/smooth_${stg}km/FIGURES";
-$dir_mesh = "$dir0/smooth_${stg}km";
-$dir_ker_lab = "\\/net\\/sierra\\/raid1\\/carltape\\/socal\\/socal_3D\\/RUNS\\/GRADIENT\\/${smodel}\\/smooth_${stg}km";
+$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/COVERAGE";
+$dir_smodel = "$dir0/${dtag}";
+$dir_mesh = $dir_smodel;
+$dir_output = "$dir0/${dtag}/FIGURES";
+$dir_ker_lab = "\\/net\\/sierra\\/raid1\\/carltape\\/socal\\/socal_3D\\/RUNS\\/MODELS\\/COVERAGE\\/${dtag}";
-# remove figures and vtk files in local directory
-`rm *.pdf *.ps`;
-
# check for vtu files
$nvtu0 = 2;
$nvtu = `ls -1 ${dir_mesh}/*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"}
-
# kernel options
@klabs = ("kappa","mu","rho","alpha","beta","rhop");
@ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
$klab = $klabs[$iker-1];
$ktitle = $ktitles[$iker-1];
-$file_tag = "${dir_ker_lab}\\/${klab}_kernel${stag}";
-$file_tag = "${dir_ker_lab}\\/${klab}_kernel_smooth${stag}";
-$tlab = "Summed $smodel kernel for $ktitle";
-$ftag = "${klab}_kernel${stag}";
+$file_tag = "${dir_ker_lab}\\/${ktag}";
+$tlab = "$smodel $ktitle coverage kernel";
+$ftag = "${smodel}_${ktag}_${dtag}_coverage";
-#print "\n ${file_tag} "; die("testing");
-
$pwd = basename($ENV{PWD});
# corners of the UTM mesh
@@ -105,9 +100,6 @@
# volumetric VTU files
$file1 = "${file_tag}.vtu";
-#$file1 = "${file_tag}_1.vtu";
-#$file2 = "${file_tag}_2.vtu";
-#$file3 = "${file_tag}_3.vtu";
$xcen = $x_center;
$ycen = $y_center;
@@ -129,11 +121,13 @@
# 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 = 2; $imax = 9; # testing
+#$imin = 2; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
#$zcen = $utm_zmax - ($i-1)*$dzinc;
@@ -142,22 +136,19 @@
printf ("%2.2i : %.1f, %.1f, %.1f\n",$i,$xcen,$ycen,$zcen);
# title for plot
- $title = sprintf("%s, Depth = %.1f km",$tlab,-$zcen_km);
+ $title = sprintf("%s -- Cut at z = %.1f km",$tlab,-$zcen_km);
# open file for replacing the cross-section line
open(SED,">sed.txt");
# file names
print SED "/kReader1 SetFileName/s/SetFileName.*\$/SetFileName $file1/ \n";
- #print SED "/kReader2 SetFileName/s/SetFileName.*\$/SetFileName $file2/ \n";
- #print SED "/kReader3 SetFileName/s/SetFileName.*\$/SetFileName $file3/ \n";
# color scale
print SED "/scalarBar SetTitle/s/SetTitle.*\$/SetTitle \"$clabel\"/ \n";
$scalar_low = sprintf("%.4e",-$cmax);
$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";
@@ -165,25 +156,6 @@
printf SED ("/hrPlane1 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
print SED "/hrPlane1 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
-# print SED "/kMapper1 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
-# print SED "/hrMapper1 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
-
-# print SED "/kMapper2 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
-# print SED "/hrMapper2 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
-
-# print SED "/kMapper3 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
-# print SED "/hrMapper3 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
-
-# # cross-section
-# printf SED ("/hrPlane1 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
-# print SED "/hrPlane1 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
-
-# printf SED ("/hrPlane2 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
-# print SED "/hrPlane2 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
-
-# printf SED ("/hrPlane3 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
-# print SED "/hrPlane3 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
-
# # color scale
# $scalar_low = $cmean[$i-1] * 1000 * (1 - $pert/100);
# $scalar_high = $cmean[$i-1] * 1000 * (1 + $pert/100);
@@ -223,9 +195,8 @@
#---------------------------
-$ofile = "${dir_output}/${ftag}_set.pdf";
+$ofile = "${ftag}_set.pdf";
#system("\\rm ${ofile}");
system("/home/carltape/bin/pdcat -r *.pdf ${ofile}");
-system("mv *.ps ${dir_output}");
#---------------------------
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.tcl 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum.tcl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -42,12 +42,6 @@
vtkXMLUnstructuredGridReader kReader1
kReader1 SetFileName vtu_files/lin_model/vp_low_1.vtu
-# vtkXMLUnstructuredGridReader kReader2
-# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
-
-# vtkXMLUnstructuredGridReader kReader3
-# kReader3 SetFileName vtu_files/lin_model/vp_low_3.vtu
-
## When the files are merged, everything crashes.
#
#vtkAppendFilter kReader
@@ -68,20 +62,10 @@
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
kMapper1 SetScalarRange 4687.47 5729.13
-# vtkDataSetMapper kMapper2
-# kMapper2 SetInput [kReader2 GetOutput]
-# kMapper2 SetScalarRange 4687.47 5729.13
-# vtkDataSetMapper kMapper3
-# kMapper3 SetInput [kReader3 GetOutput]
-# kMapper3 SetScalarRange 4687.47 5729.13
# data actors
vtkActor kActor1
kActor1 SetMapper kMapper1
-# vtkActor kActor2
-# kActor2 SetMapper kMapper2
-# vtkActor kActor3
-# kActor3 SetMapper kMapper3
# === generate plane cut ====
vtkPlane hrPlane1
@@ -98,34 +82,25 @@
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
-# vtkPlane hrPlane2
-# hrPlane2 SetOrigin 385719.8 3823329.2 0.0
-# hrPlane2 SetNormal 0 0 1
-# vtkCutter hrCut2
-# hrCut2 SetInput [kReader2 GetOutput]
-# hrCut2 SetCutFunction hrPlane2
-# vtkDataSetMapper hrMapper2
-# hrMapper2 SetInput [hrCut2 GetOutput]
-# hrMapper2 InterpolateScalarsBeforeMappingOn
-# hrMapper2 SetScalarRange 4687.47 5729.13
-# hrMapper2 SetLookupTable lut
-# vtkActor hrActor2
-# hrActor2 SetMapper hrMapper2
+# === source (sphere->mapper->actor) ====
+# vtkSphereSource sourceSphere
+# sourceSphere SetCenter 427523.33 3752613.39 0.0
+# sourceSphere SetRadius 15000.0
+# sourceSphere SetThetaResolution 20
+# sourceSphere SetPhiResolution 20
-# vtkPlane hrPlane3
-# hrPlane3 SetOrigin 385719.8 3823329.2 0.0
-# hrPlane3 SetNormal 0 0 1
-# vtkCutter hrCut3
-# hrCut3 SetInput [kReader3 GetOutput]
-# hrCut3 SetCutFunction hrPlane3
-# vtkDataSetMapper hrMapper3
-# hrMapper3 SetInput [hrCut3 GetOutput]
-# hrMapper3 InterpolateScalarsBeforeMappingOn
-# hrMapper3 SetScalarRange 4687.47 5729.13
-# hrMapper3 SetLookupTable lut
-# vtkActor hrActor3
-# hrActor3 SetMapper hrMapper3
+# vtkPolyDataMapper sourceMapper
+# sourceMapper SetInput [sourceSphere GetOutput]
+# vtkActor sourceActor
+# sourceActor SetMapper sourceMapper
+# eval [sourceActor GetProperty] SetColor $hot_pink
+# [sourceActor GetProperty] SetSpecularColor 1 1 1
+# [sourceActor GetProperty] SetSpecular 0.3
+# [sourceActor GetProperty] SetSpecularPower 20
+# [sourceActor GetProperty] SetAmbient 0.2
+# [sourceActor GetProperty] SetDiffuse 0.8
+
# === load the coastline vtk file ===
vtkPolyDataReader cReader1
cReader1 SetFileName vtk_files/coastfile_mod_utm.vtk
@@ -194,13 +169,13 @@
#=== text title ===
vtkTextActor titleActor
- titleActor SetDisplayPosition 250 60
+ titleActor SetDisplayPosition 280 60
titleActor SetInput "Lin model with Harvard model -- Depth = -0.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
$tprop SetFontFamilyToTimes
- $tprop SetFontSize 16
+ $tprop SetFontSize 15
#=== render window and camera positions ===
vtkRenderer ren1
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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_local.tcl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -40,14 +40,8 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/GRADIENT/m03/smooth_00km/mu_kernel_smooth.vtu
+ kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/COVERAGE/sum_abs/mu_kernel_smooth.vtu
-# vtkXMLUnstructuredGridReader kReader2
-# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
-
-# vtkXMLUnstructuredGridReader kReader3
-# kReader3 SetFileName vtu_files/lin_model/vp_low_3.vtu
-
## When the files are merged, everything crashes.
#
#vtkAppendFilter kReader
@@ -67,25 +61,15 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange -2.0000e-09 2.0000e-09
-# vtkDataSetMapper kMapper2
-# kMapper2 SetInput [kReader2 GetOutput]
-# kMapper2 SetScalarRange 4687.47 5729.13
-# vtkDataSetMapper kMapper3
-# kMapper3 SetInput [kReader3 GetOutput]
-# kMapper3 SetScalarRange 4687.47 5729.13
+ kMapper1 SetScalarRange -5.0000e-10 5.0000e-10
# data actors
vtkActor kActor1
kActor1 SetMapper kMapper1
-# vtkActor kActor2
-# kActor2 SetMapper kMapper2
-# vtkActor kActor3
-# kActor3 SetMapper kMapper3
# === generate plane cut ====
vtkPlane hrPlane1
- hrPlane1 SetOrigin 385719.8 3823329.2 -40000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -20000.0
hrPlane1 SetNormal 0 0 1
vtkCutter hrCut1
hrCut1 SetInput [kReader1 GetOutput]
@@ -93,39 +77,30 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -2.0000e-09 2.0000e-09
+ hrMapper1 SetScalarRange -5.0000e-10 5.0000e-10
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
-# vtkPlane hrPlane2
-# hrPlane2 SetOrigin 385719.8 3823329.2 0.0
-# hrPlane2 SetNormal 0 0 1
-# vtkCutter hrCut2
-# hrCut2 SetInput [kReader2 GetOutput]
-# hrCut2 SetCutFunction hrPlane2
-# vtkDataSetMapper hrMapper2
-# hrMapper2 SetInput [hrCut2 GetOutput]
-# hrMapper2 InterpolateScalarsBeforeMappingOn
-# hrMapper2 SetScalarRange 4687.47 5729.13
-# hrMapper2 SetLookupTable lut
-# vtkActor hrActor2
-# hrActor2 SetMapper hrMapper2
+# === source (sphere->mapper->actor) ====
+# vtkSphereSource sourceSphere
+# sourceSphere SetCenter 427523.33 3752613.39 0.0
+# sourceSphere SetRadius 15000.0
+# sourceSphere SetThetaResolution 20
+# sourceSphere SetPhiResolution 20
-# vtkPlane hrPlane3
-# hrPlane3 SetOrigin 385719.8 3823329.2 0.0
-# hrPlane3 SetNormal 0 0 1
-# vtkCutter hrCut3
-# hrCut3 SetInput [kReader3 GetOutput]
-# hrCut3 SetCutFunction hrPlane3
-# vtkDataSetMapper hrMapper3
-# hrMapper3 SetInput [hrCut3 GetOutput]
-# hrMapper3 InterpolateScalarsBeforeMappingOn
-# hrMapper3 SetScalarRange 4687.47 5729.13
-# hrMapper3 SetLookupTable lut
-# vtkActor hrActor3
-# hrActor3 SetMapper hrMapper3
+# vtkPolyDataMapper sourceMapper
+# sourceMapper SetInput [sourceSphere GetOutput]
+# vtkActor sourceActor
+# sourceActor SetMapper sourceMapper
+# eval [sourceActor GetProperty] SetColor $hot_pink
+# [sourceActor GetProperty] SetSpecularColor 1 1 1
+# [sourceActor GetProperty] SetSpecular 0.3
+# [sourceActor GetProperty] SetSpecularPower 20
+# [sourceActor GetProperty] SetAmbient 0.2
+# [sourceActor GetProperty] SetDiffuse 0.8
+
# === load the coastline vtk file ===
vtkPolyDataReader cReader1
cReader1 SetFileName vtk_files/coastfile_mod_utm.vtk
@@ -194,13 +169,13 @@
#=== text title ===
vtkTextActor titleActor
- titleActor SetDisplayPosition 250 60
- titleActor SetInput "Summed m03 kernel for SHEAR MODULUS, Depth = 40.0 km"
+ titleActor SetDisplayPosition 280 60
+ titleActor SetInput "m16 SHEAR MODULUS coverage kernel -- Cut at z = 20.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
$tprop SetFontFamilyToTimes
- $tprop SetFontSize 16
+ $tprop SetFontSize 15
#=== render window and camera positions ===
vtkRenderer ren1
@@ -212,8 +187,8 @@
iren SetRenderWindow renWin
vtkCamera cam1
- cam1 SetPosition 385719.8 3823329.2 10000.0
- cam1 SetFocalPoint 385719.8 3823329.2 -40000.0
+ cam1 SetPosition 385719.8 3823329.2 30000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -20000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -243,7 +218,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "mu_kernel_10.ps"
+ writer SetFileName "m16_mu_kernel_smooth_sum_abs_coverage_09.ps"
writer Write
wm withdraw .
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.pl (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.pl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -0,0 +1,231 @@
+#!/usr/bin/perl -w
+
+#---------------------------------
+# view_kernel_sum.pl
+# Carl Tape, after Qinya Liu
+#
+# This script creates a PDF of cross-sections of a kernel.
+# Execute from /net/denali/scratch1/carltape/vtk/carl_new/
+#
+# 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 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)
+#
+# view_kernel_sum.pl 0 1 m00 4e-9 # un-smoothed kappa
+# view_kernel_sum.pl 6 1 m00 4e-9 # smoothed kappa (6 km)
+# view_kernel_sum.pl 10 1 m00 4e-9 # smoothed kappa (10 km)
+# view_kernel_sum.pl 15 1 m00 4e-9 # smoothed kappa (15 km)
+# view_kernel_sum.pl 20 1 m00 4e-9 # smoothed kappa (20 km)
+#
+#---------------------------------
+
+use File::Basename;
+
+if (@ARGV < 4) {die("Usage: view_kernel_sum.pl ismooth iker smodel cmax\n")}
+($gsmooth,$iker,$smodel,$cmax) = @ARGV;
+
+$tcl_tag = "view_kernel_sum";
+$tcl_file = "${tcl_tag}.tcl";
+$tcl_file_local ="${tcl_tag}_local.tcl";
+if (not -f $tcl_file) {die("Check if sample file $tcl_file exists or not\n");}
+
+$clabel = "m\\^-3";
+$stg = sprintf("%2.2i",$gsmooth);
+
+# directories
+$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS/GRADIENT/${smodel}";
+$dir_output = "$dir0/smooth_${stg}km/FIGURES";
+$dir_mesh = "$dir0/smooth_${stg}km";
+$dir_ker_lab = "\\/net\\/sierra\\/raid1\\/carltape\\/socal\\/socal_3D\\/RUNS\\/GRADIENT\\/${smodel}\\/smooth_${stg}km";
+
+# remove figures and vtk files in local directory
+`rm *.pdf *.ps`;
+
+# check for vtu files
+$nvtu0 = 2;
+$nvtu = `ls -1 ${dir_mesh}/*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"}
+
+# kernel options
+ at klabs = ("kappa","mu","rho","alpha","beta","rhop");
+ at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
+$klab = $klabs[$iker-1];
+$ktitle = $ktitles[$iker-1];
+$file_tag = "${dir_ker_lab}\\/${klab}_kernel${stag}";
+$file_tag = "${dir_ker_lab}\\/${klab}_kernel_smooth${stag}";
+$tlab = "Summed $smodel kernel for $ktitle";
+$ftag = "${klab}_kernel${stag}";
+
+#print "\n ${file_tag} "; die("testing");
+
+$pwd = basename($ENV{PWD});
+
+# corners of the UTM mesh
+$utm_xmin = 0.06623919273678*1e6;
+$utm_xmax = 0.70520037473958*1e6;
+$utm_ymin = 3.57170704093370*1e6;
+$utm_ymax = 4.07495136244577*1e6;
+$utm_zmax = 0.;
+$utm_zmin = -60000.;
+$xran = $utm_xmax - $utm_xmin;
+$yran = $utm_ymax - $utm_ymin;
+$zran = $utm_zmax - $utm_zmin;
+
+# center of the mesh
+$x_center = ($utm_xmin + $utm_xmax)/2;
+$y_center = ($utm_ymin + $utm_ymax)/2;
+$z_center = ($utm_zmin + $utm_zmax)/2;
+
+# increment between cross sections
+$dxinc = 50000; # meters
+$dyinc = 50000; # meters
+$dzinc = 5000; # meters
+
+# number of cuts
+$Nx = int($xran/$dxinc) + 1;
+$Ny = int($yran/$dyinc) + 1;
+$Nz = int($zran/$dzinc) + 1;
+
+print "\n UTM-bounds for mesh:\n";
+print "$utm_xmin, $utm_xmax, $utm_ymin, $utm_ymax, $utm_zmin, $utm_zmax\n";
+print "DX increment between cross-sections: $dxinc\n";
+print "DY increment between cross-sections: $dyinc\n";
+print "DZ increment between cross-sections: $dzinc\n";
+print "Number of cross-sections: $Ny\n";
+
+if ($Nx < 1 || $Ny < 1 || $Nz < 1) {
+ die("Number of cross-sections must be at least 1\n");
+}
+
+# volumetric VTU files
+$file1 = "${file_tag}.vtu";
+#$file1 = "${file_tag}_1.vtu";
+#$file2 = "${file_tag}_2.vtu";
+#$file3 = "${file_tag}_3.vtu";
+
+$xcen = $x_center;
+$ycen = $y_center;
+
+# normal for plane
+$nx = 0; $ny = 0; $nz = 1;
+
+# 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);
+$Nz = @dlayers;
+
+$imin = 1; $imax = $Nz; # default
+#$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);
+
+ # title for plot
+ $title = sprintf("%s, Depth = %.1f km",$tlab,-$zcen_km);
+
+ # open file for replacing the cross-section line
+ open(SED,">sed.txt");
+
+ # file names
+ print SED "/kReader1 SetFileName/s/SetFileName.*\$/SetFileName $file1/ \n";
+ #print SED "/kReader2 SetFileName/s/SetFileName.*\$/SetFileName $file2/ \n";
+ #print SED "/kReader3 SetFileName/s/SetFileName.*\$/SetFileName $file3/ \n";
+
+ # color scale
+ print SED "/scalarBar SetTitle/s/SetTitle.*\$/SetTitle \"$clabel\"/ \n";
+
+ $scalar_low = sprintf("%.4e",-$cmax);
+ $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";
+
+ # cross-section
+ printf SED ("/hrPlane1 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
+ print SED "/hrPlane1 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
+
+# print SED "/kMapper1 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+# print SED "/hrMapper1 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+
+# print SED "/kMapper2 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+# print SED "/hrMapper2 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+
+# print SED "/kMapper3 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+# print SED "/hrMapper3 SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+
+# # cross-section
+# printf SED ("/hrPlane1 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
+# print SED "/hrPlane1 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
+
+# printf SED ("/hrPlane2 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
+# print SED "/hrPlane2 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
+
+# printf SED ("/hrPlane3 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
+# print SED "/hrPlane3 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
+
+ # # color scale
+ # $scalar_low = $cmean[$i-1] * 1000 * (1 - $pert/100);
+ # $scalar_high = $cmean[$i-1] * 1000 * (1 + $pert/100);
+ # print SED "/kMapper SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+ # print SED "/hrMapper SetScalarRange/s/SetScalarRange.*\$/SetScalarRange $scalar_low $scalar_high/ \n";
+
+ # # horizontal cross-section
+ # printf SED ("/hrPlane1 SetOrigin/s/SetOrigin.*\$/SetOrigin %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
+ # print SED "/hrPlane1 SetNormal/s/SetNormal.*\$/SetNormal $nx $ny $nz / \n";
+
+ # SetPosition
+ $cam_dist = 50000;
+ $nxd = $xcen;
+ $nyd = $ycen;
+ $nzd = $zcen + $cam_dist;
+
+ # orientation
+ printf SED ("/cam1 SetFocalPoint/s/SetFocalPoint.*\$/SetFocalPoint %.1f %.1f %.1f / \n",$xcen,$ycen,$zcen);
+ printf SED ("/cam1 SetPosition/s/SetPosition.*\$/SetPosition %.1f %.1f %.1f / \n",$nxd,$nyd,$nzd);
+ print SED "/cam1 SetViewUp/s/SetViewUp.*\$/SetViewUp $ux $uy $uz/ \n";
+ print SED "/titleActor SetInput/s/SetInput.*\$/SetInput \"$title\"/ \n";
+
+ # file name
+ $filename = sprintf("${ftag}_%2.2i.ps",$i);
+ print SED "/writer SetFileName/s/SetFileName.*\$/SetFileName \"$filename\"/ \n";
+
+ # close SED file
+ close(SED);
+
+ # make a local executable file and run it to generate a PS figure
+ #print "\n ${tcl_file_local} \n";
+ system("sed -f sed.txt ${tcl_file} > ${tcl_file_local}");
+ system("vtk ${tcl_file_local}");
+ system("ps2pdf $filename");
+
+}
+
+#---------------------------
+
+$ofile = "${dir_output}/${ftag}_set.pdf";
+#system("\\rm ${ofile}");
+system("/home/carltape/bin/pdcat -r *.pdf ${ofile}");
+system("mv *.ps ${dir_output}");
+
+#---------------------------
Property changes on: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.pl
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.tcl (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.tcl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -0,0 +1,257 @@
+#
+# This program automates the cross-section cutting for Banana-Doughnut kernels
+#
+# We start off by loading some Tcl modules. One is the basic VTK library;
+# the second is a package for rendering, and the last includes a set
+# of color definitions.
+#
+package require vtk
+package require vtkinteraction
+package require vtktesting
+
+# === lookup table ===
+vtkLookupTable lut
+lut SetNumberOfTableValues 25
+lut SetTableValue 0 0.72549 0 0 1
+lut SetTableValue 1 0.847059 0 0 1
+lut SetTableValue 2 0.964706 0 0 1
+lut SetTableValue 3 1 0.0862745 0 1
+lut SetTableValue 4 1 0.203922 0 1
+lut SetTableValue 5 1 0.32549 0 1
+lut SetTableValue 6 1 0.447059 0 1
+lut SetTableValue 7 1 0.564706 0 1
+lut SetTableValue 8 1 0.686275 0 1
+lut SetTableValue 9 1 0.807843 0 1
+lut SetTableValue 10 1 0.92549 0 1
+lut SetTableValue 11 1 1 0 1
+lut SetTableValue 12 1 1 0 1
+lut SetTableValue 13 1 1 0 1
+lut SetTableValue 14 0.858824 1 0.027451 1
+lut SetTableValue 15 0.623529 1 0.0666667 1
+lut SetTableValue 16 0.392157 1 0.109804 1
+lut SetTableValue 17 0.247059 0.980392 0.211765 1
+lut SetTableValue 18 0.117647 0.960784 0.32549 1
+lut SetTableValue 19 0 0.92549 0.443137 1
+lut SetTableValue 20 0 0.701961 0.65098 1
+lut SetTableValue 21 0 0.47451 0.854902 1
+lut SetTableValue 22 0 0.282353 0.980392 1
+lut SetTableValue 23 0 0.168627 0.909804 1
+lut SetTableValue 24 0 0.054902 0.839216 1
+
+# === load the unstructured grid data ===
+vtkXMLUnstructuredGridReader kReader1
+ kReader1 SetFileName vtu_files/lin_model/vp_low_1.vtu
+
+# vtkXMLUnstructuredGridReader kReader2
+# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
+
+# vtkXMLUnstructuredGridReader kReader3
+# kReader3 SetFileName vtu_files/lin_model/vp_low_3.vtu
+
+## When the files are merged, everything crashes.
+#
+#vtkAppendFilter kReader
+# kReader AddInput [kReader1 GetOutput]
+# kReader AddInput [kReader2 GetOutput]
+# kReader AddInput [kReader3 GetOutput]
+#
+## data mappers
+#vtkDataSetMapper kMapper
+# kMapper SetInput [kReader GetOutput]
+# kMapper SetScalarRange 4687.47 5729.13
+#
+## data Actor
+#vtkActor kActor
+# kActor SetMapper kMapper
+
+# data mappers
+vtkDataSetMapper kMapper1
+ kMapper1 SetInput [kReader1 GetOutput]
+ kMapper1 SetScalarRange 4687.47 5729.13
+# vtkDataSetMapper kMapper2
+# kMapper2 SetInput [kReader2 GetOutput]
+# kMapper2 SetScalarRange 4687.47 5729.13
+# vtkDataSetMapper kMapper3
+# kMapper3 SetInput [kReader3 GetOutput]
+# kMapper3 SetScalarRange 4687.47 5729.13
+
+# data actors
+vtkActor kActor1
+ kActor1 SetMapper kMapper1
+# vtkActor kActor2
+# kActor2 SetMapper kMapper2
+# vtkActor kActor3
+# kActor3 SetMapper kMapper3
+
+# === generate plane cut ====
+vtkPlane hrPlane1
+ hrPlane1 SetOrigin 385719.8 3823329.2 0.0
+ hrPlane1 SetNormal 0 0 1
+vtkCutter hrCut1
+ hrCut1 SetInput [kReader1 GetOutput]
+ hrCut1 SetCutFunction hrPlane1
+vtkDataSetMapper hrMapper1
+ hrMapper1 SetInput [hrCut1 GetOutput]
+ hrMapper1 InterpolateScalarsBeforeMappingOn
+ hrMapper1 SetScalarRange 4687.47 5729.13
+ hrMapper1 SetLookupTable lut
+vtkActor hrActor1
+ hrActor1 SetMapper hrMapper1
+
+# vtkPlane hrPlane2
+# hrPlane2 SetOrigin 385719.8 3823329.2 0.0
+# hrPlane2 SetNormal 0 0 1
+# vtkCutter hrCut2
+# hrCut2 SetInput [kReader2 GetOutput]
+# hrCut2 SetCutFunction hrPlane2
+# vtkDataSetMapper hrMapper2
+# hrMapper2 SetInput [hrCut2 GetOutput]
+# hrMapper2 InterpolateScalarsBeforeMappingOn
+# hrMapper2 SetScalarRange 4687.47 5729.13
+# hrMapper2 SetLookupTable lut
+# vtkActor hrActor2
+# hrActor2 SetMapper hrMapper2
+
+# vtkPlane hrPlane3
+# hrPlane3 SetOrigin 385719.8 3823329.2 0.0
+# hrPlane3 SetNormal 0 0 1
+# vtkCutter hrCut3
+# hrCut3 SetInput [kReader3 GetOutput]
+# hrCut3 SetCutFunction hrPlane3
+# vtkDataSetMapper hrMapper3
+# hrMapper3 SetInput [hrCut3 GetOutput]
+# hrMapper3 InterpolateScalarsBeforeMappingOn
+# hrMapper3 SetScalarRange 4687.47 5729.13
+# hrMapper3 SetLookupTable lut
+# vtkActor hrActor3
+# hrActor3 SetMapper hrMapper3
+
+# === load the coastline vtk file ===
+vtkPolyDataReader cReader1
+ cReader1 SetFileName vtk_files/coastfile_mod_utm.vtk
+
+vtkPolyDataMapper cMapper1
+ cMapper1 SetInput [cReader1 GetOutput]
+
+vtkActor cActor1
+ cActor1 SetMapper cMapper1
+ eval [cActor1 GetProperty] SetColor $navy
+ eval [cActor1 GetProperty] SetLineWidth 2
+
+# === load the borders vtk file ===
+vtkPolyDataReader cReader2
+ cReader2 SetFileName vtk_files/borderfile_mod_utm.vtk
+
+vtkPolyDataMapper cMapper2
+ cMapper2 SetInput [cReader2 GetOutput]
+
+vtkActor cActor2
+ cActor2 SetMapper cMapper2
+ eval [cActor2 GetProperty] SetColor $navy
+ eval [cActor2 GetProperty] SetLineWidth 1
+
+# === load the plate boundary vtk file ===
+vtkPolyDataReader cReader3
+ cReader3 SetFileName vtk_files/NA_PA_boundary_utm.vtk
+
+vtkPolyDataMapper cMapper3
+ cMapper3 SetInput [cReader3 GetOutput]
+
+vtkActor cActor3
+ cActor3 SetMapper cMapper3
+ eval [cActor3 GetProperty] SetColor $navy
+ eval [cActor3 GetProperty] SetLineWidth 2
+
+# === load the faults vtk file ===
+vtkPolyDataReader cReader4
+ cReader4 SetFileName vtk_files/jennings_more_utm_air.vtk
+
+vtkPolyDataMapper cMapper4
+ cMapper4 SetInput [cReader4 GetOutput]
+
+vtkActor cActor4
+ cActor4 SetMapper cMapper4
+ eval [cActor4 GetProperty] SetColor $navy
+ eval [cActor4 GetProperty] SetLineWidth 1.5
+
+#=== color bar ====
+# NOTE hrMapper variable here
+vtkScalarBarActor scalarBar
+ scalarBar SetLookupTable lut
+ #scalarBar SetLookupTable [hrMapper GetLookupTable]
+scalarBar SetTitle "P-wave (m\\/s)"
+ [scalarBar GetPositionCoordinate] SetCoordinateSystemToNormalizedViewport
+ [scalarBar GetPositionCoordinate] SetValue 0.1 0.0
+ scalarBar SetOrientationToHorizontal
+ scalarBar SetWidth 0.5
+ scalarBar SetHeight 0.1
+ scalarBar SetPosition 0.2 0.
+ scalarBar SetLabelFormat "%-#6.3g"
+ [scalarBar GetLabelTextProperty] SetColor 0 0 0
+ [scalarBar GetLabelTextProperty] SetFontFamilyToTimes
+ [scalarBar GetTitleTextProperty] SetColor 0 0 0
+ [scalarBar GetTitleTextProperty] SetFontFamilyToTimes
+
+#=== text title ===
+vtkTextActor titleActor
+ titleActor SetDisplayPosition 250 60
+ titleActor SetInput "Lin model with Harvard model -- Depth = -0.0 km"
+set tprop [titleActor GetTextProperty]
+ $tprop SetJustificationToCentered
+ $tprop SetColor 0 0 1
+ $tprop SetFontFamilyToTimes
+ $tprop SetFontSize 16
+
+#=== render window and camera positions ===
+vtkRenderer ren1
+ ren1 SetBackground 1 1 1
+vtkRenderWindow renWin
+ renWin AddRenderer ren1
+ renWin SetSize 600 600
+vtkRenderWindowInteractor iren
+ iren SetRenderWindow renWin
+
+vtkCamera cam1
+ cam1 SetPosition 385719.8 3823329.2 50000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 0.0
+ cam1 SetViewUp 0 1 0
+
+#=== scene ===============
+ren1 SetActiveCamera cam1
+
+# colorbar only works with --> nsActor
+# kActor is the volume
+
+#ren1 AddActor hrActor
+ren1 AddActor hrActor1
+#ren1 AddActor hrActor2
+#ren1 AddActor hrActor3
+ren1 AddActor cActor1
+ren1 AddActor cActor2
+#ren1 AddActor cActor3
+ren1 AddActor cActor4
+#ren1 AddActor sourceActor
+#ren1 AddActor receiverActor
+ren1 AddActor2D scalarBar
+ren1 AddActor2D titleActor
+
+# I don't understand why, but ResetCamera is very important
+ren1 ResetCamera 300000 550000 3572000 4075000 -60000 200000
+renWin Render
+
+vtkWindowToImageFilter w2i
+ w2i SetInput renWin
+vtkPostScriptWriter writer
+ writer SetInputConnection [w2i GetOutputPort]
+ writer SetFileName "xc_hr_02.ps"
+ writer Write
+
+wm withdraw .
+
+exit
+
+#set cam2 [ren1 GetActiveCamera]
+#puts stdout [ $cam2 GetPosition]
+#puts stdout [$cam2 GetFocalPoint]
+#puts stdout [$cam2 GetViewUp]
+
Property changes on: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_kernel_sum_orig.tcl
___________________________________________________________________
Name: svn:executable
+ *
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model.pl 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model.pl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -11,19 +11,19 @@
#
# EXAMPLES:
#
-# view_model.pl 4 vp_m14
+# view_model.pl 4 vp_m16
# view_model.pl 4 vp_m00
-# view_model.pl 5 vs_m14
+# view_model.pl 5 vs_m16
# view_model.pl 5 vs_m00
-# view_model.pl 6 vb_m14
+# view_model.pl 6 vb_m16
# view_model.pl 6 vb_m00
-# view_model.pl 7 poisson_m14
+# view_model.pl 7 poisson_m16
# view_model.pl 7 poisson_m00
-# view_model.pl 8 vp_m14_m13
-# view_model.pl 8 vp_m14_m00
-# view_model.pl 9 vs_m14_m13
-# view_model.pl 9 vs_m14_m00
-# view_model.pl 10 vb_m14_m00
+# view_model.pl 8 vp_m16_m15
+# view_model.pl 8 vp_m16_m00
+# view_model.pl 9 vs_m16_m15
+# view_model.pl 9 vs_m16_m00
+# view_model.pl 10 vb_m16_m00
#
#---------------------------------
@@ -48,11 +48,13 @@
#`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";
-#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m13/*vtu .`; $title_tag = "Model m13 for";
-`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m14/*vtu .`; $title_tag = "Model m14 for";
+`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m13/*vtu .`; $title_tag = "Model m13 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m14/*vtu .`; $title_tag = "Model m14 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m15/*vtu .`; $title_tag = "Model m15 for";
+#`cp /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/m16/*vtu .`; $title_tag = "Model m16 for";
-$title_tag = "LN(m14\\/m13) for"; $cmax = 0.05;
-$title_tag = "LN(m14\\/m00) for"; $cmax = 0.20;
+#$title_tag = "LN(m13\\/m12) for"; $cmax = 0.05;
+#$title_tag = "LN(m16\\/m00) for"; $cmax = 0.20;
$tcl_tag = "view_model";
$tcl_file = "${tcl_tag}.tcl";
@@ -199,7 +201,7 @@
$Nz = @dlayers;
$imin = 1; $imax = $Nz; # default
#$imin = 9; $imax = 13; # testing
-#$imin = 2; $imax = $imin; # testing
+$imin = 2; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_model_local.tcl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName vs_m13_m00.vtu
+ kReader1 SetFileName vs_m13.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-01 2.0000e-01
+ kMapper1 SetScalarRange 2507.5 3392.5
# data actors
vtkActor kActor1
@@ -75,7 +75,7 @@
# === generate plane cut ====
vtkPlane hrPlane1
- hrPlane1 SetOrigin 385719.8 3823329.2 -40000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -1.0
hrPlane1 SetNormal 0 0 1
vtkCutter hrCut1
hrCut1 SetInput [kReader1 GetOutput]
@@ -83,7 +83,7 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -2.0000e-01 2.0000e-01
+ hrMapper1 SetScalarRange 2507.5 3392.5
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
@@ -160,7 +160,7 @@
vtkScalarBarActor scalarBar
scalarBar SetLookupTable lut
#scalarBar SetLookupTable [hrMapper GetLookupTable]
-scalarBar SetTitle "Percent Change"
+scalarBar SetTitle "2950 m/s +- 15 percent"
[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 "LN(m13/m00) for S WAVE-SPEED -- Cut at z = 40.0 km"
+ titleActor SetInput "Model m13 for S WAVE-SPEED -- Cut at z = 0.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
@@ -193,8 +193,8 @@
iren SetRenderWindow renWin
vtkCamera cam1
- cam1 SetPosition 385719.8 3823329.2 10000.0
- cam1 SetFocalPoint 385719.8 3823329.2 -40000.0
+ cam1 SetPosition 385719.8 3823329.2 49999.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -1.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -220,7 +220,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "vs_m13_m00_13.ps"
+ writer SetFileName "vs_m13_02.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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update.pl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -48,7 +48,15 @@
# view_update.pl 2 dm13 82 1.00 0 dm_mu_kernel_smooth_p082 beta_window
# view_update.pl 2 dm13 74 2.80 0 dm_mu_kernel_smooth_p074 beta_window
# view_update.pl 2 dm13 78 2.40 0 dm_mu_kernel_smooth_p074 beta_window
+
+# view_update.pl 2 dm14 60 3.40 0 dm_mu_kernel_smooth_p060 beta_window # 123
+# view_update.pl 2 dm14 80 2.00 0 dm_mu_kernel_smooth_p080 beta_window # 47
+# view_update.pl 2 dm14 72 2.00 0 dm_mu_kernel_smooth_p072 beta_window # 80
+# view_update.pl 2 dm14 80 1.80 0 dm_mu_kernel_smooth_p080 beta_window # 40
#
+# view_update.pl 2 dm15 80 1.80 0 dm_mu_kernel_smooth_p080 beta_window
+# view_update.pl 2 dm15 90 0.80 0 dm_mu_kernel_smooth_p090 beta_window
+#
# 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
#
@@ -152,8 +160,8 @@
@dlayers = (-0.25,0.001,2,4,6,8,10,15,20,25,30,35,40);
$Nz = @dlayers;
-#$imin = 1; $imax = $Nz; # default
-$imin = 2; $imax = $imin; # testing
+$imin = 1; $imax = $Nz; # default
+#$imin = 2; $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 2009-03-27 20:32:32 UTC (rev 14493)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/UTILS/vtk/view_update_local.tcl 2009-03-27 20:39:58 UTC (rev 14494)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/dm13/beta_window/dm_mu_kernel_smooth_p074.vtu
+ kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/dm15/beta_window/dm_mu_kernel_smooth_p090.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.4000e+00 2.4000e+00
+ kMapper1 SetScalarRange -8.0000e-01 8.0000e-01
# 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 -1.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -40000.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 -2.4000e+00 2.4000e+00
+ hrMapper1 SetScalarRange -8.0000e-01 8.0000e-01
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 78 -- Cut at z = 0.0 km"
+ titleActor SetInput "Subspace update for SHEAR-WAVE-SPEED -- model 90 -- Cut at z = 40.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 49999.0
- cam1 SetFocalPoint 385719.8 3823329.2 -1.0
+ cam1 SetPosition 385719.8 3823329.2 10000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -40000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -315,7 +315,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "dm_mu_kernel_smooth_p074_02.ps"
+ writer SetFileName "dm_mu_kernel_smooth_p090_13.ps"
writer Write
wm withdraw .
More information about the CIG-COMMITS
mailing list