[cig-commits] r15034 - seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk

carltape at geodynamics.org carltape at geodynamics.org
Fri May 22 06:38:20 PDT 2009 

Author: carltape
Date: 2009-05-22 06:38:20 -0700 (Fri, 22 May 2009)
New Revision: 15034

Modified:
   seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.pl
   seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_local.tcl
Log:
Minor update to vtk-perl scripts.


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	2009-05-22 13:36:23 UTC (rev 15033)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel.pl	2009-05-22 13:38:20 UTC (rev 15034)
@@ -40,11 +40,11 @@
 # directories
 $dir0 = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
 $dir_smodel = "$dir0/${eid}/${smodel}";
-#$dir_smodel = "$dir0/${eid}/${smodel}_BDK_SAVE";    # TEMPORARY
+$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}";
+#$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;
@@ -79,7 +79,7 @@
   $cmax = 5e-11;   # defalt
 }
 
-$cmax = 3e-13;
+$cmax = 1e-11;
 
 # corners of the UTM mesh
 $utm_xmin = 0.06623919273678*1e6;
@@ -148,8 +148,8 @@
 $Nz = @dlayers;
 
 #$imin = 1; $imax = $Nz;   # default
-$imin = 2; $imax = $Nz;  # testing
-#$imin = 4; $imax = $imin;  # testing
+#$imin = 2; $imax = 9;  # testing
+$imin = 9; $imax = $imin;  # testing
 
 for ($i = $imin; $i <= $imax; $i++) {
   #$zcen = $utm_zmax - ($i-1)*$dzinc;

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	2009-05-22 13:36:23 UTC (rev 15033)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/vtk/view_kernel_local.tcl	2009-05-22 13:38:20 UTC (rev 15034)
@@ -40,13 +40,13 @@
 
 # === load the unstructured grid data ===
 vtkXMLUnstructuredGridReader kReader1
-    kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9983429/m16/MESH_all/mu_kernel_low_1.vtu
+    kReader1 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/10006857/m16_BDK_SAVE/MESH_all/mu_kernel_low_1.vtu
 
 vtkXMLUnstructuredGridReader kReader2
-    kReader2 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9983429/m16/MESH_all/mu_kernel_low_2.vtu
+    kReader2 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/10006857/m16_BDK_SAVE/MESH_all/mu_kernel_low_2.vtu
 
 vtkXMLUnstructuredGridReader kReader3
-    kReader3 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/9983429/m16/MESH_all/mu_kernel_low_3.vtu
+    kReader3 SetFileName /net/sierra/raid1/carltape/socal/socal_3D/RUNS/10006857/m16_BDK_SAVE/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 -3.0000e-13 3.0000e-13
+    kMapper1 SetScalarRange -1.0000e-11 1.0000e-11
 vtkDataSetMapper kMapper2
     kMapper2 SetInput [kReader2 GetOutput]
-    kMapper2 SetScalarRange -3.0000e-13 3.0000e-13
+    kMapper2 SetScalarRange -1.0000e-11 1.0000e-11
 vtkDataSetMapper kMapper3
     kMapper3 SetInput [kReader3 GetOutput]
-    kMapper3 SetScalarRange -3.0000e-13 3.0000e-13
+    kMapper3 SetScalarRange -1.0000e-11 1.0000e-11
 
 # data actors
 vtkActor kActor1
@@ -85,7 +85,7 @@
 
 # === 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,13 +93,13 @@
 vtkDataSetMapper hrMapper1
    hrMapper1 SetInput [hrCut1 GetOutput]
    hrMapper1 InterpolateScalarsBeforeMappingOn
-   hrMapper1 SetScalarRange -3.0000e-13 3.0000e-13
+   hrMapper1 SetScalarRange -1.0000e-11 1.0000e-11
    hrMapper1 SetLookupTable lut
 vtkActor hrActor1
     hrActor1 SetMapper hrMapper1
 
 vtkPlane hrPlane2
-    hrPlane2 SetOrigin 385719.8 3823329.2 -40000.0 
+    hrPlane2 SetOrigin 385719.8 3823329.2 -20000.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 -3.0000e-13 3.0000e-13
+   hrMapper2 SetScalarRange -1.0000e-11 1.0000e-11
    hrMapper2 SetLookupTable lut
 vtkActor hrActor2
     hrActor2 SetMapper hrMapper2
 
 vtkPlane hrPlane3
-    hrPlane3 SetOrigin 385719.8 3823329.2 -40000.0 
+    hrPlane3 SetOrigin 385719.8 3823329.2 -20000.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 -3.0000e-13 3.0000e-13
+   hrMapper3 SetScalarRange -1.0000e-11 1.0000e-11
    hrMapper3 SetLookupTable lut
 vtkActor hrActor3
     hrActor3 SetMapper hrMapper3
@@ -267,7 +267,7 @@
 #=== text title ===
 vtkTextActor titleActor
    titleActor SetDisplayPosition 280 60
-   titleActor SetInput "9983429 : m16 kernel for SHEAR MODULUS -- 1 stations -- Cut at z = 40.0 km"
+   titleActor SetInput "10006857 : m16 kernel for SHEAR MODULUS -- 133 stations -- Cut at z = 20.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 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 ===============
@@ -315,7 +315,7 @@
   w2i SetInput renWin
 vtkPostScriptWriter writer
   writer SetInputConnection [w2i GetOutputPort]
-  writer SetFileName "029_9983429_mu_all_kernel_13.ps"
+  writer SetFileName "155_10006857_mu_all_kernel_09.ps"
   writer Write
 
 wm withdraw .

More information about the CIG-COMMITS mailing list