[cig-commits] r16246 - seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk
carltape at geodynamics.org
carltape at geodynamics.org
Tue Feb 9 16:04:19 PST 2010
Author: carltape
Date: 2010-02-09 16:04:18 -0800 (Tue, 09 Feb 2010)
New Revision: 16246
Added:
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_mesa_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa_local.tcl
Modified:
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/README
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum_local.tcl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_model.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_model_local.tcl
Log:
Updated VTK scripts to allow mesa option for off-screen rendering.
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/README
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/README 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/README 2010-02-10 00:04:18 UTC (rev 16246)
@@ -5,3 +5,10 @@
/net/denali/scratch1/carltape/vtk/carl_new/
This directory needs a lot of work still.
+
+---------------------------------------
+
+link user directory of vtk files to vtk_files
+
+---------------------------------------
+
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.pl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.pl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -3,9 +3,15 @@
#---------------------------------
# view_kernel.pl
# Carl Tape, after Qinya Liu
+# 02-Nov-2009
#
# This script creates a PDF of cross-sections of a kernel.
#
+# Originally, this was designed to use a three-part vtu file,
+# since the high-res files can be too large for the memory of
+# some computers. Now, we assume that there is a single
+# vtu file (note: can be a low-res version).
+#
# CALLED BY: tomo_make_figs.pl
#
# EXAMPLE:
@@ -13,6 +19,12 @@
# view_kernel.pl 9818433 2 m00 T02 133_9818433
# view_kernel.pl 9818433 2 m00 T06 133_9818433
#
+# view_kernel.pl 14236768 7 m16 T003_T010 14236768
+# view_kernel.pl 14236768 8 m16 T003_T010 14236768
+# view_kernel.pl 14236768 9 m16 T003_T010 14236768
+#
+# view_kernel.pl 10992159 9 m16 T003_T010 10992159
+#
#---------------------------------
use File::Basename;
@@ -21,6 +33,7 @@
($eid,$iker,$smodel,$ftag,$elab) = @ARGV;
$tcl_tag = "view_kernel";
+$tcl_tag = "view_kernel_mesa";
$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");}
@@ -38,27 +51,32 @@
$clabel = "m\\^-3";
# directories
-$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
+$dir0 = "/media/raid/carltape/SOCAL_ADJOINT/RUNS";
+
$dir_smodel = "$dir0/${eid}/${smodel}";
$dir_smodel = "$dir0/${eid}/${smodel}_BDK_SAVE"; # TEMPORARY
+$dir_smodel = "$dir0/${eid}/${smodel}/KERNEL_ORIG"; # 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}";
+$dir_ker_lab = "\\/media\\/raid\\/carltape\\/SOCAL_ADJOINT\\/RUNS\\/${eid}\\/${smodel}\\/KERNEL_ORIG\\/MESH_${ftag}";
# check for vtu files
-$nvtu0 = 6;
+$nvtu0 = 1;
$nvtu = `ls -1 ${dir_mesh}/*vtu | wc | awk '{print \$1}'`; chomp($nvtu);
-if ( $nvtu != $nvtu0 ) {die("You have $nvtu files, not $nvtu0")}
+if ( $nvtu < $nvtu0 ) {die("You have $nvtu files, not $nvtu0")}
# kernel options
- at klabs = ("kappa","mu","rho","alpha","beta","rhop");
- at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
+#@klabs = ("kappa","mu","rho","alpha","beta","rhop");
+#@ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
+ at klabs = ("kappa_kmr","mu_kmr","rho_kmr","alpha_abr","beta_abr","rho_abr","c_cbr","beta_cbr","rho_cbr");
+ at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","IMPEDANCE","BULK-SOUND-SPEED","SHEAR-WAVE-SPEED","IMPEDANCE");
$klab = $klabs[$iker-1];
$ktitle = $ktitles[$iker-1];
$file_tag = "${dir_ker_lab}\\/${klab}_kernel_low";
$tlab = "$eid : $smodel kernel for $ktitle";
-$ftag = "${elab}_${klab}_${ftag}_kernel";
+$ftag = "${elab}_${klab}_${ftag}_${smodel}_kernel";
$pwd = basename($ENV{PWD});
@@ -79,7 +97,8 @@
$cmax = 5e-11; # defalt
}
-$cmax = 1e-11;
+$cmax = 3e-13;
+#$cmax = 5e-12; # 9875657 example
# corners of the UTM mesh
$utm_xmin = 0.06623919273678*1e6;
@@ -119,9 +138,10 @@
}
# volumetric VTU files
-$file1 = "${file_tag}_1.vtu";
-$file2 = "${file_tag}_2.vtu";
-$file3 = "${file_tag}_3.vtu";
+$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;
@@ -144,12 +164,13 @@
# depth layers
# Why does depth z = 0 not work, so we must use z = 0.001 ?
#@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,2,4,6,8,10,15,20,25,30,35,40);
+ at dlayers = (0.001,1,2,3,4,6,8,10,15);
$Nz = @dlayers;
-#$imin = 1; $imax = $Nz; # default
+$imin = 1; $imax = $Nz; # default
#$imin = 2; $imax = 9; # testing
-$imin = 9; $imax = $imin; # testing
+$imin = 4; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
#$zcen = $utm_zmax - ($i-1)*$dzinc;
@@ -165,8 +186,8 @@
# 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";
+ #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";
@@ -176,21 +197,21 @@
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 "/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";
+ #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 ("/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";
+ #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);
@@ -225,10 +246,12 @@
#print "\n ${tcl_file_local} \n";
system("sed -f sed.txt ${tcl_file} > ${tcl_file_local}");
system("vtk ${tcl_file_local}");
- system("ps2pdf $filename");
+ #system("ps2pdf $filename");
}
+die("TESTING");
+
#---------------------------
$ofile = "${ftag}_set.pdf";
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.tcl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -9,6 +9,40 @@
package require vtkinteraction
package require vtktesting
+#============================================
+# Create Mesa specific render window and renderer
+# When using Mesa, you should always create Mesa specific renderers,
+# render windows, mappers, actors, lights, cameras etc...
+# See the Mesa classes in the Rendering directory for a list of
+# existing classes.
+
+#vtkXMesaRenderWindow renWin
+##vtkRenderWindow renWin
+## Will render in memory.
+#renWin OffScreenRenderingOn
+
+#vtkMesaRenderer ren1
+##vtkRenderer ren1
+#ren1 SetBackground 1 1 1
+#renWin AddRenderer ren1
+#renWin SetSize 600 600
+#============================================
+
+#=== 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
+#============================================
+
# === lookup table ===
vtkLookupTable lut
lut SetNumberOfTableValues 25
@@ -42,11 +76,11 @@
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 kReader2
+# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
-vtkXMLUnstructuredGridReader kReader3
- kReader3 SetFileName vtu_files/lin_model/vp_low_3.vtu
+# vtkXMLUnstructuredGridReader kReader3
+# kReader3 SetFileName vtu_files/lin_model/vp_low_3.vtu
## When the files are merged, everything crashes.
#
@@ -68,20 +102,20 @@
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
+# 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
+# vtkActor kActor2
+# kActor2 SetMapper kMapper2
+# vtkActor kActor3
+# kActor3 SetMapper kMapper3
# === generate plane cut ====
vtkPlane hrPlane1
@@ -98,33 +132,33 @@
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 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
+# 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
# === source (sphere->mapper->actor) ====
# vtkSphereSource sourceSphere
@@ -274,20 +308,6 @@
$tprop SetFontFamilyToTimes
$tprop SetFontSize 15
-#=== 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
@@ -296,8 +316,8 @@
#ren1 AddActor hrActor
ren1 AddActor hrActor1
-ren1 AddActor hrActor2
-ren1 AddActor hrActor3
+#ren1 AddActor hrActor2
+#ren1 AddActor hrActor3
ren1 AddActor cActor1
ren1 AddActor cActor2
#ren1 AddActor cActor3
@@ -318,9 +338,12 @@
writer SetFileName "xc_hr_02.ps"
writer Write
+# Exit with display
wm withdraw .
-exit
+## Exit without displaying anything
+#vtkCommand DeleteAllObjects
+#exit
#set cam2 [ren1 GetActiveCamera]
#puts stdout [ $cam2 GetPosition]
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_local.tcl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_local.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -9,6 +9,25 @@
package require vtkinteraction
package require vtktesting
+#============================================
+# Create Mesa specific render window and renderer
+# When using Mesa, you should always create Mesa specific renderers,
+# render windows, mappers, actors, lights, cameras etc...
+# See the Mesa classes in the Rendering directory for a list of
+# existing classes.
+
+vtkXMesaRenderWindow renWin
+#vtkRenderWindow renWin
+# Will render in memory.
+renWin OffScreenRenderingOn
+
+vtkMesaRenderer ren
+#vtkRenderer ren
+ren SetBackground 0.1 0.2 0.4
+renWin AddRenderer ren
+renWin SetSize 300 300
+#============================================
+
# === lookup table ===
vtkLookupTable lut
lut SetNumberOfTableValues 25
@@ -40,13 +59,13 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/10006857/m16_BDK_SAVE/MESH_all/mu_kernel_low_1.vtu
+ kReader1 SetFileName /media/raid/carltape/SOCAL_ADJOINT/RUNS/14169456/m61/KERNEL/MESH_T003_T010/rho_cbr_kernel_low.vtu
-vtkXMLUnstructuredGridReader kReader2
- kReader2 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/10006857/m16_BDK_SAVE/MESH_all/mu_kernel_low_2.vtu
+# vtkXMLUnstructuredGridReader kReader2
+# kReader2 SetFileName vtu_files/lin_model/vp_low_2.vtu
-vtkXMLUnstructuredGridReader kReader3
- kReader3 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/10006857/m16_BDK_SAVE/MESH_all/mu_kernel_low_3.vtu
+# vtkXMLUnstructuredGridReader kReader3
+# kReader3 SetFileName vtu_files/lin_model/vp_low_3.vtu
## When the files are merged, everything crashes.
#
@@ -67,25 +86,25 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange -1.0000e-11 1.0000e-11
-vtkDataSetMapper kMapper2
- kMapper2 SetInput [kReader2 GetOutput]
- kMapper2 SetScalarRange -1.0000e-11 1.0000e-11
-vtkDataSetMapper kMapper3
- kMapper3 SetInput [kReader3 GetOutput]
- kMapper3 SetScalarRange -1.0000e-11 1.0000e-11
+ kMapper1 SetScalarRange -3.0000e-13 3.0000e-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
+# vtkActor kActor2
+# kActor2 SetMapper kMapper2
+# vtkActor kActor3
+# kActor3 SetMapper kMapper3
# === generate plane cut ====
vtkPlane hrPlane1
- hrPlane1 SetOrigin 385719.8 3823329.2 -20000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -3000.0
hrPlane1 SetNormal 0 0 1
vtkCutter hrCut1
hrCut1 SetInput [kReader1 GetOutput]
@@ -93,38 +112,38 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -1.0000e-11 1.0000e-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 -20000.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 -1.0000e-11 1.0000e-11
- hrMapper2 SetLookupTable lut
-vtkActor hrActor2
- hrActor2 SetMapper hrMapper2
+# 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 -20000.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 -1.0000e-11 1.0000e-11
- hrMapper3 SetLookupTable lut
-vtkActor hrActor3
- hrActor3 SetMapper hrMapper3
+# 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
# === source (sphere->mapper->actor) ====
# vtkSphereSource sourceSphere
@@ -267,7 +286,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "10006857 : m16 kernel for SHEAR MODULUS -- 133 stations -- Cut at z = 20.0 km"
+ titleActor SetInput "14169456 : m61 kernel for IMPEDANCE -- 46 stations -- Cut at z = 3.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
@@ -284,8 +303,8 @@
iren SetRenderWindow renWin
vtkCamera cam1
- cam1 SetPosition 385719.8 3823329.2 30000.0
- cam1 SetFocalPoint 385719.8 3823329.2 -20000.0
+ cam1 SetPosition 385719.8 3823329.2 47000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -3000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -296,8 +315,8 @@
#ren1 AddActor hrActor
ren1 AddActor hrActor1
-ren1 AddActor hrActor2
-ren1 AddActor hrActor3
+#ren1 AddActor hrActor2
+#ren1 AddActor hrActor3
ren1 AddActor cActor1
ren1 AddActor cActor2
#ren1 AddActor cActor3
@@ -315,11 +334,14 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "155_10006857_mu_all_kernel_09.ps"
+ writer SetFileName "14169456_rho_cbr_T003_T010_kernel_04.ps"
writer Write
-wm withdraw .
+# Exit with display
+#wm withdraw .
+# Exit without displaying anything
+vtkCommand DeleteAllObjects
exit
#set cam2 [ren1 GetActiveCamera]
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_mesa_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_mesa_local.tcl (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_mesa_local.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -0,0 +1,356 @@
+# 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
+
+#============================================
+# Create Mesa specific render window and renderer
+# When using Mesa, you should always create Mesa specific renderers,
+# render windows, mappers, actors, lights, cameras etc...
+# See the Mesa classes in the Rendering directory for a list of
+# existing classes.
+
+vtkXMesaRenderWindow renWin
+#vtkRenderWindow renWin
+# Will render in memory.
+renWin OffScreenRenderingOn
+
+vtkMesaRenderer ren1
+#vtkRenderer ren1
+ren1 SetBackground 1 1 1
+renWin AddRenderer ren1
+renWin SetSize 600 600
+#============================================
+
+#=== 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
+
+vtkMesaCamera cam1
+ cam1 SetPosition 385719.8 3823329.2 47000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -3000.0
+ cam1 SetViewUp 0 1 0
+#============================================
+
+# === 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 /media/raid/carltape/SOCAL_ADJOINT/RUNS/14169456/m16/KERNEL_ORIG/MESH_T003_T010/rho_cbr_kernel_low.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
+#vtkMesaActor kActor
+# kActor SetMapper kMapper
+
+# data mappers
+vtkMesaDataSetMapper kMapper1
+ kMapper1 SetInput [kReader1 GetOutput]
+ kMapper1 SetScalarRange -3.0000e-13 3.0000e-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
+vtkMesaActor kActor1
+ kActor1 SetMapper kMapper1
+# vtkMesaActor kActor2
+# kActor2 SetMapper kMapper2
+# vtkMesaActor kActor3
+# kActor3 SetMapper kMapper3
+
+# === generate plane cut ====
+vtkPlane hrPlane1
+ hrPlane1 SetOrigin 385719.8 3823329.2 -3000.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 -3.0000e-13 3.0000e-13
+ hrMapper1 SetLookupTable lut
+vtkMesaActor 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
+# vtkMesaActor 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
+# vtkMesaActor hrActor3
+# hrActor3 SetMapper hrMapper3
+
+# === source (sphere->mapper->actor) ====
+# vtkSphereSource sourceSphere
+# sourceSphere SetCenter 427523.33 3752613.39 0.0
+# sourceSphere SetRadius 15000.0
+# sourceSphere SetThetaResolution 20
+# sourceSphere SetPhiResolution 20
+
+# vtkMesaPolyDataMapper sourceMapper
+# sourceMapper SetInput [sourceSphere GetOutput]
+
+# vtkMesaActor 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
+
+# === source from file ====
+
+vtkSphereSource sourceSphere
+ sourceSphere SetRadius 15000.0
+ sourceSphere SetThetaResolution 20
+ sourceSphere SetPhiResolution 20
+
+vtkPolyDataReader sourceReader
+ sourceReader SetFileName source.vtk
+
+vtkGlyph3D sourceGlyph
+ sourceGlyph SetInputConnection [sourceReader GetOutputPort]
+ sourceGlyph SetSource [sourceSphere GetOutput]
+
+vtkMesaPolyDataMapper sourceMapper
+ sourceMapper SetInput [sourceGlyph GetOutput]
+
+vtkMesaActor 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
+
+# === receivers from file ===
+
+vtkSphereSource receiverSphere
+ receiverSphere SetRadius 5000.0
+ receiverSphere SetThetaResolution 20
+ receiverSphere SetPhiResolution 20
+
+vtkPolyDataReader receiverReader
+ receiverReader SetFileName receiver.vtk
+
+vtkGlyph3D receiverGlyph
+ receiverGlyph SetInputConnection [receiverReader GetOutputPort]
+ receiverGlyph SetSource [receiverSphere GetOutput]
+
+vtkMesaPolyDataMapper receiverMapper
+ receiverMapper SetInput [receiverGlyph GetOutput]
+ #receiverMapper SetInputConnection [receiverGlyph GetOutputPort]
+
+vtkMesaActor receiverActor
+ receiverActor SetMapper receiverMapper
+ eval [receiverActor GetProperty] SetColor $grey
+ [receiverActor GetProperty] SetSpecularColor 1 1 1
+ [receiverActor GetProperty] SetSpecular 0.3
+ [receiverActor GetProperty] SetSpecularPower 20
+ [receiverActor GetProperty] SetAmbient 0.2
+ [receiverActor GetProperty] SetDiffuse 0.8
+
+# === load the coastline vtk file ===
+vtkPolyDataReader cReader1
+ cReader1 SetFileName vtk_files/coastfile_mod_utm.vtk
+
+vtkMesaPolyDataMapper cMapper1
+ cMapper1 SetInput [cReader1 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper2
+ cMapper2 SetInput [cReader2 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper3
+ cMapper3 SetInput [cReader3 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper4
+ cMapper4 SetInput [cReader4 GetOutput]
+
+vtkMesaActor 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 "m^-3"
+ [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 280 60
+ titleActor SetInput "14169456 : m16 kernel for IMPEDANCE -- 46 stations -- Cut at z = 3.0 km"
+set tprop [titleActor GetTextProperty]
+ $tprop SetJustificationToCentered
+ $tprop SetColor 0 0 1
+ $tprop SetFontFamilyToTimes
+ $tprop SetFontSize 15
+
+#=== 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
+
+# save the window to a PS file
+vtkWindowToImageFilter w2if
+ w2if SetInput renWin
+
+#vtkPNGWriter wr
+#wr SetInput [w2if GetOutput]
+#wr SetFileName "MangledMesaTest.png"
+#wr Write
+#vtkCommand DeleteAllObjects
+#exit
+
+vtkPostScriptWriter wr
+wr SetInput [w2if GetOutput]
+# wr SetInputConnection [w2if GetOutputPort]
+ wr SetFileName "xc_hr_02.ps"
+ wr Write
+
+# Exit without displaying anything
+vtkCommand DeleteAllObjects
+exit
+
+#set cam2 [ren1 GetActiveCamera]
+#puts stdout [ $cam2 GetPosition]
+#puts stdout [$cam2 GetFocalPoint]
+#puts stdout [$cam2 GetViewUp]
+
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth.pl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth.pl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -11,16 +11,22 @@
# EXAMPLE:
# view_kernel_smooth.pl 14095628 2 m07 all 049_14095628 6 1
#
+# view_kernel_smooth.pl 13935988 9 m16 T003_T010 4 1
+# view_kernel_smooth.pl 13935988 9 m61 T003_T010 4 1
+#
#---------------------------------
use File::Basename;
-if (@ARGV < 7) {die("Usage: view_kernel_smooth.pl eid iker smodel elab\n")}
-($eid,$iker,$smodel,$ftag,$elab,$hsmooth,$vsmooth) = @ARGV;
+if (@ARGV < 6) {die("Usage: view_kernel_smooth.pl eid iker smodel elab\n")}
+($eid,$iker,$smodel,$ftag,$hsmooth,$vsmooth) = @ARGV;
+#($eid,$iker,$smodel,$ftag,$elab,$hsmooth,$vsmooth) = @ARGV;
+$elab = $eid; # temporary
$stg = sprintf("h%3.3ikm_v%3.3ikm",$hsmooth,$vsmooth);
-$tcl_tag = "view_kernel_smooth";
+#$tcl_tag = "view_kernel_smooth";
+$tcl_tag = "view_kernel_smooth_mesa"; # offscreen rendering
$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");}
@@ -30,39 +36,58 @@
$clabel = "m\\^-3";
-# directories
-$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
-$dir_output = "$dir0/${eid}/${smodel}/OUTPUT_${ftag}";
-$dir_mesh = "$dir0/SMOOTH_EVENT_KERNELS/${smodel}/${eid}";
-$dir_ker_lab = "\\/net\\/sierra\\/raid1\\/carltape\\/socal\\/socal_3D\\/RUNS\\/SMOOTH_EVENT_KERNELS\\/${smodel}\\/${eid}";
+# directories -- these paths cannot be linked
+#$dir0 = "/home/carltape/RUNS";
+#$dir_mesh = "$dir0/SMOOTH_EVENT_KERNELS/${smodel}/${eid}";
+#$dir_ker_lab = "\\/home\\/carltape\\/RUNS\\/SMOOTH_EVENT_KERNELS\\/${smodel}\\/${eid}";
+#$dir0 = "/media/%export%raid1/carltape/SOCAL_ADJOINT/RUNS";
+#$dir_mesh = "$dir0/SMOOTH_EVENT_KERNELS/${smodel}/${eid}";
+#$dir_ker_lab = "\\/media\\/%export%raid1\\/carltape\\/SOCAL_ADJOINT\\/RUNS\\/SMOOTH_EVENT_KERNELS\\/${smodel}\\/${eid}";
+
+$dir0 = "/media/raid/carltape/SOCAL_ADJOINT/RUNS";
+$dir_smodel = "$dir0/${eid}/${smodel}/KERNEL";
+$dir_output = "${dir_smodel}/OUTPUT_${ftag}";
+$dir_mesh = "${dir_smodel}/MESH_${ftag}"; # smoothed directory
+$dir_ker_lab = "\\/media\\/raid\\/carltape\\/SOCAL_ADJOINT\\/RUNS\\/${eid}\\/${smodel}\\/KERNEL\\/MESH_${ftag}";
+
# check for vtu files
-$nvtu0 = 2;
+$nvtu0 = 1;
$nvtu = `ls -1 ${dir_mesh}/*${stg}.vtu | wc | awk '{print \$1}'`; chomp($nvtu);
-if ( $nvtu != $nvtu0 ) {die("You have $nvtu files, not $nvtu0")}
+#if ( $nvtu < $nvtu0 ) {die("You have $nvtu files, not $nvtu0")}
+if ( $nvtu < $nvtu0 ) {
+ `mesh2vtu_all.pl $dir_mesh`;
+}
+#die("TESTING");
+
# smoothing options
if ( $hsmooth == 0 ) {$stag = ""} else {$stag = "_smooth_${stg}"}
# kernel options
- at klabs = ("kappa","mu","rho","alpha","beta","rhop");
- at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
+#@klabs = ("kappa","mu","rho","alpha","beta","rhop");
+#@ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
+ at klabs = ("kappa_kmr","mu_kmr","rho_kmr","alpha_abr","beta_abr","rho_abr","c_cbr","beta_cbr","rho_cbr");
+ at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","IMPEDANCE","BULK-SOUND-SPEED","SHEAR-WAVE-SPEED","IMPEDANCE");
$klab = $klabs[$iker-1];
$ktitle = $ktitles[$iker-1];
$file_tag = "${dir_ker_lab}\\/${klab}_kernel${stag}";
$tlab = "$eid : $smodel kernel for $ktitle";
-$filetag = "${elab}_${klab}_${ftag}_kernel_${stg}";
+$filetag = "${elab}_${klab}_${ftag}_${smodel}_kernel_${stg}";
+$filetag = "${filetag}_tar1";
$pwd = basename($ENV{PWD});
+# copy vtk files here
+#`split_sr_vtk.pl sr.vtk`;
+`cp ${dir_output}/*vtk .`;
+
# number of receivers
$nrec = `wc sr.vtk | awk '{print \$1}'` - 6;
#print "\n -- $nrec -- \n"; die("testing");
# copy cmax value here
-#$cmax = 15e-11; # default
-#$cmax = 5e-11; # default
-$cmax = 2.5e-11; # default
+$cmax = 5e-14;
# corners of the UTM mesh
$utm_xmin = 0.06623919273678*1e6;
@@ -128,11 +153,14 @@
# depth layers
# Why does depth z = 0 not work, so we must use z = 0.001 ?
#@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,2,4,6,8,10,15,20,25,30,35,40);
+#@dlayers = (0.001,2,4,6,8,10,15,20,25,30);
+ at dlayers = (0.001,1,2,3,4,6,8,10,15);
$Nz = @dlayers;
-#$imin = 1; $imax = $Nz; # default
-$imin = 4; $imax = $imin; # testing
+$imin = 1; $imax = $Nz; # default
+#$imin = 2; $imax = 9; # testing
+$imin = 3; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
#$zcen = $utm_zmax - ($i-1)*$dzinc;
@@ -201,7 +229,7 @@
# file name
$filename = sprintf("${filetag}_%2.2i.ps",$i);
- print SED "/writer SetFileName/s/SetFileName.*\$/SetFileName \"$filename\"/ \n";
+ print SED "/writer SetFileName/s/SetFileName.*\$/SetFileName \"$filename\"/ \n";
# close SED file
close(SED);
@@ -209,6 +237,7 @@
# 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}");
+ print "executing: vtk ${tcl_file_local}\n";
system("vtk ${tcl_file_local}");
system("ps2pdf $filename");
@@ -220,4 +249,10 @@
#system("\\rm ${ofile}");
system("/home/carltape/bin/pdcat -r *.pdf ${ofile}");
+die("TESTING");
+
+# optional
+`cp ${ofile} $dir_mesh`;
+`cp ${ofile} /media/raid/carltape/SOCAL_ADJOINT/results/KERNELS/kernel_m16_imaging/`;
+
#---------------------------
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_local.tcl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_local.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/SMOOTH_EVENT_KERNELS/m16/9944301/mu_kernel_smooth_h006km_v001km.vtu
+ kReader1 SetFileName /media/raid/carltape/SOCAL_ADJOINT/RUNS/9967901/m16/KERNEL/MESH_T003_T010/rho_cbr_kernel_smooth_h004km_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 -2.5000e-11 2.5000e-11
+ kMapper1 SetScalarRange -5.0000e-14 5.0000e-14
# 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 -4000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -15000.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.5000e-11 2.5000e-11
+ hrMapper1 SetScalarRange -5.0000e-14 5.0000e-14
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
@@ -267,7 +267,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "9944301 : m16 kernel for SHEAR MODULUS -- 80 stations -- Cut at z = 4.0 km"
+ titleActor SetInput "9967901 : m16 kernel for IMPEDANCE -- 10 stations -- Cut at z = 15.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 35000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -15000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -315,7 +315,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "291_9944301_mu_all_kernel_h006km_v001km_04.ps"
+ writer SetFileName "9967901_rho_cbr_T003_T010_m16_kernel_h004km_v001km_tar4_09.ps"
writer Write
wm withdraw .
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa.tcl (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -0,0 +1,355 @@
+#
+# 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
+
+#============================================
+# Create Mesa specific render window and renderer
+# When using Mesa, you should always create Mesa specific renderers,
+# render windows, mappers, actors, lights, cameras etc...
+# See the Mesa classes in the Rendering directory for a list of
+# existing classes.
+
+vtkXMesaRenderWindow renWin
+#vtkRenderWindow renWin
+# Will render in memory.
+renWin OffScreenRenderingOn
+
+vtkMesaRenderer ren1
+#vtkRenderer ren1
+ren1 SetBackground 1 1 1
+renWin AddRenderer ren1
+renWin SetSize 600 600
+
+vtkMesaCamera cam1
+ cam1 SetPosition 385719.8 3823329.2 47000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -3000.0
+ cam1 SetViewUp 0 1 0
+
+# === 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
+vtkMesaActor 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
+vtkMesaActor 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
+
+# === source (sphere->mapper->actor) ====
+# vtkSphereSource sourceSphere
+# sourceSphere SetCenter 427523.33 3752613.39 0.0
+# sourceSphere SetRadius 15000.0
+# sourceSphere SetThetaResolution 20
+# sourceSphere SetPhiResolution 20
+
+# 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
+
+# === source from file ====
+
+vtkSphereSource sourceSphere
+ sourceSphere SetRadius 15000.0
+ sourceSphere SetThetaResolution 20
+ sourceSphere SetPhiResolution 20
+
+vtkPolyDataReader sourceReader
+ sourceReader SetFileName source.vtk
+
+vtkGlyph3D sourceGlyph
+ sourceGlyph SetInputConnection [sourceReader GetOutputPort]
+ sourceGlyph SetSource [sourceSphere GetOutput]
+
+vtkMesaPolyDataMapper sourceMapper
+ sourceMapper SetInput [sourceGlyph GetOutput]
+
+vtkMesaActor 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
+
+# === receivers from file ===
+
+vtkSphereSource receiverSphere
+ receiverSphere SetRadius 5000.0
+ receiverSphere SetThetaResolution 20
+ receiverSphere SetPhiResolution 20
+
+vtkPolyDataReader receiverReader
+ receiverReader SetFileName receiver.vtk
+
+vtkGlyph3D receiverGlyph
+ receiverGlyph SetInputConnection [receiverReader GetOutputPort]
+ receiverGlyph SetSource [receiverSphere GetOutput]
+
+vtkMesaPolyDataMapper receiverMapper
+ receiverMapper SetInput [receiverGlyph GetOutput]
+ #receiverMapper SetInputConnection [receiverGlyph GetOutputPort]
+
+vtkMesaActor receiverActor
+ receiverActor SetMapper receiverMapper
+ eval [receiverActor GetProperty] SetColor $grey
+ [receiverActor GetProperty] SetSpecularColor 1 1 1
+ [receiverActor GetProperty] SetSpecular 0.3
+ [receiverActor GetProperty] SetSpecularPower 20
+ [receiverActor GetProperty] SetAmbient 0.2
+ [receiverActor GetProperty] SetDiffuse 0.8
+
+# === load the coastline vtk file ===
+vtkPolyDataReader cReader1
+ cReader1 SetFileName vtk_files/coastfile_mod_utm.vtk
+
+vtkMesaPolyDataMapper cMapper1
+ cMapper1 SetInput [cReader1 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper2
+ cMapper2 SetInput [cReader2 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper3
+ cMapper3 SetInput [cReader3 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper4
+ cMapper4 SetInput [cReader4 GetOutput]
+
+vtkMesaActor 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 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 15
+
+#=== 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
+
+# save the window to a PS file
+vtkWindowToImageFilter w2if
+ w2if SetInput renWin
+
+vtkPostScriptWriter writer
+writer SetInput [w2if GetOutput]
+# writer SetInputConnection [w2if GetOutputPort]
+ writer SetFileName "xc_hr_02.ps"
+ writer Write
+
+# Exit without displaying anything
+vtkCommand DeleteAllObjects
+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/model_plot/vtk/view_kernel_smooth_mesa.tcl
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa_local.tcl (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_smooth_mesa_local.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -0,0 +1,355 @@
+#
+# 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
+
+#============================================
+# Create Mesa specific render window and renderer
+# When using Mesa, you should always create Mesa specific renderers,
+# render windows, mappers, actors, lights, cameras etc...
+# See the Mesa classes in the Rendering directory for a list of
+# existing classes.
+
+vtkXMesaRenderWindow renWin
+#vtkRenderWindow renWin
+# Will render in memory.
+renWin OffScreenRenderingOn
+
+vtkMesaRenderer ren1
+#vtkRenderer ren1
+ren1 SetBackground 1 1 1
+renWin AddRenderer ren1
+renWin SetSize 600 600
+
+vtkMesaCamera cam1
+ cam1 SetPosition 385719.8 3823329.2 48000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -2000.0
+ cam1 SetViewUp 0 1 0
+
+# === 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 /media/raid/carltape/SOCAL_ADJOINT/RUNS/14169456/m61/KERNEL/MESH_T003_T010/rho_cbr_kernel_smooth_h004km_v001km.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 -5.0000e-14 5.0000e-14
+# 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
+vtkMesaActor 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 -2000.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 -5.0000e-14 5.0000e-14
+ hrMapper1 SetLookupTable lut
+vtkMesaActor 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
+
+# === source (sphere->mapper->actor) ====
+# vtkSphereSource sourceSphere
+# sourceSphere SetCenter 427523.33 3752613.39 0.0
+# sourceSphere SetRadius 15000.0
+# sourceSphere SetThetaResolution 20
+# sourceSphere SetPhiResolution 20
+
+# 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
+
+# === source from file ====
+
+vtkSphereSource sourceSphere
+ sourceSphere SetRadius 15000.0
+ sourceSphere SetThetaResolution 20
+ sourceSphere SetPhiResolution 20
+
+vtkPolyDataReader sourceReader
+ sourceReader SetFileName source.vtk
+
+vtkGlyph3D sourceGlyph
+ sourceGlyph SetInputConnection [sourceReader GetOutputPort]
+ sourceGlyph SetSource [sourceSphere GetOutput]
+
+vtkMesaPolyDataMapper sourceMapper
+ sourceMapper SetInput [sourceGlyph GetOutput]
+
+vtkMesaActor 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
+
+# === receivers from file ===
+
+vtkSphereSource receiverSphere
+ receiverSphere SetRadius 5000.0
+ receiverSphere SetThetaResolution 20
+ receiverSphere SetPhiResolution 20
+
+vtkPolyDataReader receiverReader
+ receiverReader SetFileName receiver.vtk
+
+vtkGlyph3D receiverGlyph
+ receiverGlyph SetInputConnection [receiverReader GetOutputPort]
+ receiverGlyph SetSource [receiverSphere GetOutput]
+
+vtkMesaPolyDataMapper receiverMapper
+ receiverMapper SetInput [receiverGlyph GetOutput]
+ #receiverMapper SetInputConnection [receiverGlyph GetOutputPort]
+
+vtkMesaActor receiverActor
+ receiverActor SetMapper receiverMapper
+ eval [receiverActor GetProperty] SetColor $grey
+ [receiverActor GetProperty] SetSpecularColor 1 1 1
+ [receiverActor GetProperty] SetSpecular 0.3
+ [receiverActor GetProperty] SetSpecularPower 20
+ [receiverActor GetProperty] SetAmbient 0.2
+ [receiverActor GetProperty] SetDiffuse 0.8
+
+# === load the coastline vtk file ===
+vtkPolyDataReader cReader1
+ cReader1 SetFileName vtk_files/coastfile_mod_utm.vtk
+
+vtkMesaPolyDataMapper cMapper1
+ cMapper1 SetInput [cReader1 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper2
+ cMapper2 SetInput [cReader2 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper3
+ cMapper3 SetInput [cReader3 GetOutput]
+
+vtkMesaActor 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
+
+vtkMesaPolyDataMapper cMapper4
+ cMapper4 SetInput [cReader4 GetOutput]
+
+vtkMesaActor 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 "m^-3"
+ [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 280 60
+ titleActor SetInput "14169456 : m61 kernel for IMPEDANCE -- 17 stations -- Cut at z = 2.0 km"
+set tprop [titleActor GetTextProperty]
+ $tprop SetJustificationToCentered
+ $tprop SetColor 0 0 1
+ $tprop SetFontFamilyToTimes
+ $tprop SetFontSize 15
+
+#=== 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 48000.0
+# cam1 SetFocalPoint 385719.8 3823329.2 -2000.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
+
+# save the window to a PS file
+vtkWindowToImageFilter w2if
+ w2if SetInput renWin
+
+vtkPostScriptWriter writer
+writer SetInput [w2if GetOutput]
+# writer SetInputConnection [w2if GetOutputPort]
+ writer SetFileName "14169456_rho_cbr_T003_T010_m61_kernel_h004km_v001km_tar1_03.ps"
+ writer Write
+
+# Exit without displaying anything
+vtkCommand DeleteAllObjects
+exit
+
+#set cam2 [ren1 GetActiveCamera]
+#puts stdout [ $cam2 GetPosition]
+#puts stdout [$cam2 GetFocalPoint]
+#puts stdout [$cam2 GetViewUp]
+
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum.pl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum.pl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -14,6 +14,8 @@
# 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 m16 sum rho_cbr_kernel_smooth_h004km_v001km 9 3e-13
+#
#---------------------------------
use File::Basename;
@@ -39,25 +41,30 @@
$clabel = "m\\^-3";
# directories
-$dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/COVERAGE";
+#$dir0 = "/media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/COVERAGE";
+$dir0 = "/media/raid/carltape/SOCAL_ADJOINT/RUNS/GRADIENT/${smodel}";
$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}";
+#$dir_output = "$dir0/${dtag}/FIGURES";
+$dir_ker_lab = "\\/media\\/raid\\/carltape\\/SOCAL_ADJOINT\\/RUNS\\/GRADIENT\\/${smodel}\\/${dtag}";
# 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")}
+#$nvtu0 = 2;
+#$nvtu = `ls -1 ${dir_mesh}/*vtu | wc | awk '{print \$1}'`; chomp($nvtu);
+#if ( $nvtu != $nvtu0 ) {die("You have $nvtu files, not $nvtu0")}
# kernel options
- at klabs = ("kappa","mu","rho","alpha","beta","rhop");
- at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
+#@klabs = ("kappa","mu","rho","alpha","beta","rhop");
+#@ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","DENSITY");
+ at klabs = ("kappa_kmr","mu_kmr","rho_kmr","alpha_abr","beta_abr","rho_abr","c_cbr","beta_cbr","rho_cbr");
+ at ktitles = ("BULK MODULUS","SHEAR MODULUS","DENSITY","P-WAVE-SPEED","S-WAVE-SPEED","IMPEDANCE","BULK-SOUND-SPEED","SHEAR-WAVE-SPEED","IMPEDANCE");
$klab = $klabs[$iker-1];
$ktitle = $ktitles[$iker-1];
$file_tag = "${dir_ker_lab}\\/${ktag}";
-$tlab = "$smodel $ktitle coverage kernel";
-$ftag = "${smodel}_${ktag}_${dtag}_coverage";
+#$tlab = "$smodel $ktitle coverage kernel";
+#$ftag = "${smodel}_${ktag}_${dtag}_coverage";
+$tlab = "$smodel $ktitle summed kernel";
+$ftag = "${smodel}_${ktag}_${dtag}";
$pwd = basename($ENV{PWD});
@@ -100,6 +107,7 @@
# volumetric VTU files
$file1 = "${file_tag}.vtu";
+if (not -f $file1) {print "file1 $file1 does not exist\n";}
$xcen = $x_center;
$ycen = $y_center;
@@ -122,11 +130,12 @@
# depth layers
# Why does depth z = 0 not work, so we must use z = 0.001 ?
#@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,2,4,6,8,10,15,20,25,30,35,40);
+ at dlayers = (0.001,1,2,3,4,6,8,10,15);
$Nz = @dlayers;
$imin = 1; $imax = $Nz; # default
-$imin = 2; $imax = 9; # testing
+#$imin = 2; $imax = 9; # testing
#$imin = 2; $imax = $imin; # testing
for ($i = $imin; $i <= $imax; $i++) {
@@ -199,4 +208,8 @@
#system("\\rm ${ofile}");
system("/home/carltape/bin/pdcat -r *.pdf ${ofile}");
+# optional
+`cp ${ofile} $dir0`;
+`cp ${ofile} /media/raid/carltape/SOCAL_ADJOINT/results/KERNELS/kernel_m16_imaging/`;
+
#---------------------------
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum_local.tcl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_sum_local.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/MODELS/COVERAGE/sum_abs/mu_kernel_smooth.vtu
+ kReader1 SetFileName /media/raid/carltape/SOCAL_ADJOINT/RUNS/GRADIENT/m16/sum/rho_cbr_kernel_smooth_h000km_v000km.vtu
## When the files are merged, everything crashes.
#
@@ -61,7 +61,7 @@
# data mappers
vtkDataSetMapper kMapper1
kMapper1 SetInput [kReader1 GetOutput]
- kMapper1 SetScalarRange -5.0000e-10 5.0000e-10
+ kMapper1 SetScalarRange -3.0000e-13 3.0000e-13
# data actors
vtkActor kActor1
@@ -69,7 +69,7 @@
# === generate plane cut ====
vtkPlane hrPlane1
- hrPlane1 SetOrigin 385719.8 3823329.2 -20000.0
+ hrPlane1 SetOrigin 385719.8 3823329.2 -15000.0
hrPlane1 SetNormal 0 0 1
vtkCutter hrCut1
hrCut1 SetInput [kReader1 GetOutput]
@@ -77,7 +77,7 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange -5.0000e-10 5.0000e-10
+ hrMapper1 SetScalarRange -3.0000e-13 3.0000e-13
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
@@ -170,7 +170,7 @@
#=== text title ===
vtkTextActor titleActor
titleActor SetDisplayPosition 280 60
- titleActor SetInput "m16 SHEAR MODULUS coverage kernel -- Cut at z = 20.0 km"
+ titleActor SetInput "m16 IMPEDANCE summed kernel -- Cut at z = 15.0 km"
set tprop [titleActor GetTextProperty]
$tprop SetJustificationToCentered
$tprop SetColor 0 0 1
@@ -187,8 +187,8 @@
iren SetRenderWindow renWin
vtkCamera cam1
- cam1 SetPosition 385719.8 3823329.2 30000.0
- cam1 SetFocalPoint 385719.8 3823329.2 -20000.0
+ cam1 SetPosition 385719.8 3823329.2 35000.0
+ cam1 SetFocalPoint 385719.8 3823329.2 -15000.0
cam1 SetViewUp 0 1 0
#=== scene ===============
@@ -218,7 +218,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "m16_mu_kernel_smooth_sum_abs_coverage_09.ps"
+ writer SetFileName "m16_rho_cbr_kernel_smooth_h000km_v000km_sum_09.ps"
writer Write
wm withdraw .
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_model.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_model.pl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_model.pl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -10,6 +10,7 @@
# CALLED BY: tomo_make_figs.pl
#
# EXAMPLES:
+# view_model.pl 5 vs_m61
#
# view_model.pl 4 vp_m16
# view_model.pl 4 vp_m00
@@ -34,25 +35,27 @@
# 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 .`; $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";
-`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";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m00/vtu_old/*vtu .`;
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m00/*vtu .`; $title_tag = "Model m00 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m01/*vtu .`; $title_tag = "Model m01 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m02/*vtu .`; $title_tag = "Model m02 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m03/*vtu .`; $title_tag = "Model m03 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m04/*vtu .`; $title_tag = "Model m04 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m05/*vtu .`; $title_tag = "Model m05 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m06/*vtu .`; $title_tag = "Model m06 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m07/*vtu .`; $title_tag = "Model m07 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m08/*vtu .`; $title_tag = "Model m08 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m09/*vtu .`; $title_tag = "Model m09 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m10/*vtu .`; $title_tag = "Model m10 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m11/*vtu .`; $title_tag = "Model m11 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m12/*vtu .`; $title_tag = "Model m12 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m13/*vtu .`; $title_tag = "Model m13 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m14/*vtu .`; $title_tag = "Model m14 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m15/*vtu .`; $title_tag = "Model m15 for";
+#`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m16/*vtu .`; $title_tag = "Model m16 for";
+`cp /media/raid/carltape/SOCAL_ADJOINT/RUNS/MODELS/m61/*vtu .`; $title_tag = "Model m16-smooth for";
+
#$title_tag = "LN(m13\\/m12) for"; $cmax = 0.05;
#$title_tag = "LN(m16\\/m00) for"; $cmax = 0.20;
@@ -201,7 +204,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/model_plot/vtk/view_model_local.tcl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_model_local.tcl 2010-02-09 23:55:51 UTC (rev 16245)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_model_local.tcl 2010-02-10 00:04:18 UTC (rev 16246)
@@ -40,7 +40,7 @@
# === load the unstructured grid data ===
vtkXMLUnstructuredGridReader kReader1
- kReader1 SetFileName vs_m13.vtu
+ kReader1 SetFileName vs_m61.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 2507.5 3392.5
+ kMapper1 SetScalarRange 4275 4725
# data actors
vtkActor kActor1
@@ -75,7 +75,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]
@@ -83,7 +83,7 @@
vtkDataSetMapper hrMapper1
hrMapper1 SetInput [hrCut1 GetOutput]
hrMapper1 InterpolateScalarsBeforeMappingOn
- hrMapper1 SetScalarRange 2507.5 3392.5
+ hrMapper1 SetScalarRange 4275 4725
hrMapper1 SetLookupTable lut
vtkActor hrActor1
hrActor1 SetMapper hrMapper1
@@ -160,7 +160,7 @@
vtkScalarBarActor scalarBar
scalarBar SetLookupTable lut
#scalarBar SetLookupTable [hrMapper GetLookupTable]
-scalarBar SetTitle "2950 m/s +- 15 percent"
+scalarBar SetTitle "4500 m/s +- 5 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 "Model m13 for S WAVE-SPEED -- Cut at z = 0.0 km"
+ titleActor SetInput "Model m16-smooth for S WAVE-SPEED -- Cut at z = 40.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 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 ===============
@@ -220,7 +220,7 @@
w2i SetInput renWin
vtkPostScriptWriter writer
writer SetInputConnection [w2i GetOutputPort]
- writer SetFileName "vs_m13_02.ps"
+ writer SetFileName "vs_m61_13.ps"
writer Write
wm withdraw .
More information about the CIG-COMMITS
mailing list