[cig-commits] commit: add strike and dip angles, and rake unit vector to vtp files created with flt2vtk.sh

[Mercurial]

changeset:   217:21cb95fabea0
tag:         tip
user:        Sylvain Barbot <sbarbot at ntu.edu.sg>
date:        Sat Sep 07 15:27:22 2013 +0800
files:       util/flt2vtk.sh util/seg2flt.py
description:
add strike and dip angles, and rake unit vector to vtp files created with flt2vtk.sh


diff -r 6d37c413e667 -r 21cb95fabea0 util/flt2vtk.sh
--- a/util/flt2vtk.sh	Fri Aug 30 16:48:45 2013 +0800
+++ b/util/flt2vtk.sh	Sat Sep 07 15:27:22 2013 +0800
@@ -22,6 +22,9 @@
 		W=$6;
 		strike=$7*pi/180;
 		dip=$8*pi/180;
+		rake=$9*pi/180;
+		nstrike=7;
+		ndip=8;
 		} else {
 		slip=$2
 		x1=$3;
@@ -32,9 +35,8 @@
 		strike=$8*pi/180;
 		dip=$9*pi/180;
 		rake=$10*pi/180;
-		b[0]=slip*cos(strike)*cos(rake)+sin(strike)*cos(dip)*sin(rake);
-		b[1]=slip*sin(strike)*cos(rake)-cos(strike)*cos(dip)*sin(rake);
-		b[2]=slip*(-1)*sin(dip)*sin(rake);
+		nstrike=8;
+		ndip=9;
 		}
 		s[0]=cos(strike);
 		s[1]=sin(strike);
@@ -42,6 +44,9 @@
 		d[0]=+sin(strike)*cos(dip);
 		d[1]=-cos(strike)*cos(dip);
 		d[2]=-sin(dip);
+		b[0]=+cos(strike)*cos(rake)+sin(strike)*cos(dip)*sin(rake);
+		b[1]=+sin(strike)*cos(rake)-cos(strike)*cos(dip)*sin(rake);
+		b[2]=-sin(dip)*sin(rake);
 		printf("    <Piece NumberOfPoints=\"4\" NumberOfPolys=\"1\">\n");
 		printf("      <Points>\n");
 		printf("        <DataArray type=\"Float32\" Name=\"Fault Patch\" NumberOfComponents=\"3\" format=\"ascii\">\n");
@@ -59,13 +64,23 @@
 		printf("          4\n");
 		printf("        </DataArray>\n");
 		printf("      </Polys>\n");
+		printf("      <CellData Scalar=\"geometry\">\n");
+		printf("        <DataArray type=\"Float32\" Name=\"strike\" NumberOfComponents=\"1\" format=\"ascii\">\n");
+		printf("%f\n",((($nstrike+180)%360-360)%360+180));
+		printf("        </DataArray>\n");
+		printf("        <DataArray type=\"Float32\" Name=\"dip\" NumberOfComponents=\"1\" format=\"ascii\">\n");
+		printf("%f\n",$ndip);
+		printf("        </DataArray>\n");
 		if (10==NF){
-		printf("      <CellData Normals=\"slip\">\n");
 		printf("        <DataArray type=\"Float32\" Name=\"slip\" NumberOfComponents=\"3\" format=\"ascii\">\n");
+		printf("%f %f %f\n",slip*b[0],slip*b[1],slip*b[2]);
+		printf("        </DataArray>\n");
+		} else {
+		printf("        <DataArray type=\"Float32\" Name=\"unit slip\" NumberOfComponents=\"3\" format=\"ascii\">\n");
 		printf("%f %f %f\n",b[0],b[1],b[2]);
 		printf("        </DataArray>\n");
+		}
 		printf("      </CellData>\n");
-		}
 		printf("    </Piece>\n");
 	}
 	END{
diff -r 6d37c413e667 -r 21cb95fabea0 util/seg2flt.py
--- a/util/seg2flt.py	Fri Aug 30 16:48:45 2013 +0800
+++ b/util/seg2flt.py	Sat Sep 07 15:27:22 2013 +0800
@@ -1,4 +1,4 @@
-#!/opt/local/bin/python2
+#!/opt/local/bin/python
 
 """
 Created on Mon Jul  2 06:15:53 2012
@@ -104,7 +104,7 @@
         fid=stdin
     else:
         fname=argv[1]
-        print fname, len(argv)
+        #print fname, len(argv)
         #if not path.isfile(fname):
         #    raise ValueError("invalid file name: " + fname)
         fid=open(fname, 'r')