<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div></div><div><br></div><div><div><div>On 12/04/2011, at 4:47 AM, Brad Aagaard wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>On 04/11/2011 09:01 AM, luigi vadacca wrote:<br><blockquote type="cite">Dear all,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I have two questions about the plot of the fault slip vector and principal<br></blockquote><blockquote type="cite">stress axes in paraview.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">In figure 7.33 of the example "step11" of the pylith manual is written:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">"Note that PyLith outputs slip in the fault coordinate system, so we transform<br></blockquote><blockquote type="cite">them to the global coordinate<br></blockquote><blockquote type="cite">system using the Calculator in ParaView. A more general approach involves<br></blockquote><blockquote type="cite">outputing the fault coordinate<br></blockquote><blockquote type="cite">system information and using these fields in the Calculator"<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">1 - How can I do this?<br></blockquote><br>For the fault output, request that the output include the strike, dip, <br>and normal directions in the info file using something similar to the <br>following (you may need to adjust the name of the fault output component <br>if you have more than one fault or name your fault component something <br>else; you may also want to add other info fields to the output such as <br>final slip and slip time):<br><br>[pylithapp.timedependent.interfaces.fault.output]<br>vertex_info_fields = [strike_dir,dip_dir,normal_dir]<br><br>In ParaView you can use the calculator to take the dot product of the <br>slip with the three directions to get the slip vector in global <br>coordinates. For example,<br><br>(slip[0]*strike_dir[0] +<br> slip[1]*dip_dir[0] +<br> slip[2]*normal_dir[0]) * iHat +<br>(slip[0]*strike_dir[1] +<br> slip[1]*dip_dir[1] +<br> slip[2]*normal_dir[1]) * jHat +<br>(slip[0]*strike_dir[2] +<br> slip[1]*dip_dir[2] +<br> slip[2]*normal_dir[2]) * kHat<br><br><blockquote type="cite"><br></blockquote><blockquote type="cite">2 - Is possible to plot the three principal stress axes directly in paraview?<br></blockquote><br>ParaView does not supply a filter to calculate the principal stresses, <br>so you will need to write your own post-processing script.<br><br><br>Brad<br>_______________________________________________<br>CIG-SHORT mailing list<br><a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short<br></div></blockquote></div><br><div>
<span class="Apple-style-span" style="font-size: 12px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" style="font: normal normal normal 12px/normal Helvetica; ">Charles A. Williams</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Scientist<br>GNS Science<br>1 Fairway Drive, Avalon<br>PO Box 30368<br>Lower Hutt 5040<br>New Zealand<br>ph (office): 0064-4570-4566<br>fax (office): 0064-4570-4600</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><a href="mailto:C.Williams@gns.cri.nz">C.Williams@gns.cri.nz</a></div></div></div></span>
</div>
<br></div></body></html>