[cig-commits] r22761 - seismo/3D/SPECFEM3D/trunk/utils/Visualization/vtk_scripts
carltape at geodynamics.org
carltape at geodynamics.org
Tue Sep 3 12:58:11 PDT 2013
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";
More information about the CIG-COMMITS
mailing list