[cig-commits] r22761 - seismo/3D/SPECFEM3D/trunk/utils/Visualization/vtk_scripts

[carltape at geodynamics.org]

Author: carltape
Date: 2013-09-03 12:58:10 -0700 (Tue, 03 Sep 2013)
New Revision: 22761

Modified:
   seismo/3D/SPECFEM3D/trunk/utils/Visualization/vtk_scripts/split_sr_vtk.pl
Log:
modified split_sr_vtk.pl to handle epicenter.vtk for the case of a finite fault; eventually we should have SPECFEM3D output source.vtk, epicenter.vtk, receiver.vtk, rather than using split_sr_vtk.pl with sr.vtk


Modified: seismo/3D/SPECFEM3D/trunk/utils/Visualization/vtk_scripts/split_sr_vtk.pl
===================================================================
--- seismo/3D/SPECFEM3D/trunk/utils/Visualization/vtk_scripts/split_sr_vtk.pl	2013-09-03 18:59:37 UTC (rev 22760)
+++ seismo/3D/SPECFEM3D/trunk/utils/Visualization/vtk_scripts/split_sr_vtk.pl	2013-09-03 19:58:10 UTC (rev 22761)
@@ -39,7 +39,7 @@
 
 #print "\n -- $lines[0] -- $lines[1] -- $lines[2] -- $lines[3] -- $lines[4] -- $lines[5] \n";
 
-($xsrc,$ysrc,undef) = split(" ",$lines[5]);
+#($xsrc,$ysrc,undef) = split(" ",$lines[5]);
 $nlines = @lines;
 $nrec = $nlines - $nsrc - 6;
 
@@ -73,8 +73,11 @@
 print FILE3 "Epicenter VTK file\n";
 print FILE3 "$lines[2]";
 print FILE3 "$lines[3]";
-print FILE3 "POINTS      1 float\n";
-print FILE3 "$xsrc $ysrc $zepi";
+print FILE3 "POINTS     $nsrc float\n";
+for ($k=5; $k <= $nsrc+4; $k=$k+1) {
+   ($x,$y,undef) = split("\ ",$lines[$k]);
+   print FILE3 "$x $y $zepi\n";
+}
 close(FILE3);
 
 print "\n-- writing out source.vtk, epicenter.vtk, receiver.vtk in directory $odir --\n";