<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;" dir="ltr">
<p></p>
<p>Hi Everyone, </p>
<p><br>
</p>
<p>I'm very new to Paraview and python scripting. I have a series of time dependent .vtu files and would like to </p>
<p></p>
<ol style="margin-bottom: 0px; margin-top: 0px;">
<li><span style="font-size: 12pt;">Read the vtu files for each time step</span></li><li><span style="font-size: 12pt;"><span style="font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 12pt;">c</span><span style="font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 12pt;">alculate
the magnitude of shear stress</span><br>
</span></li><li><span style="font-size: 12pt;"></span><span style="font-size: 12pt;">c</span><span style="font-size: 12pt;">alculate time-averaged shear stress for all the elements in my geometry </span><br>
</li><li><span style="font-size: 12pt;">output/write a vtu file and load it back to paraview. </span></li></ol>
<div><br>
</div>
<div>I think I have been able to successfully implement steps 1 and 2 (please see bellow) but am completely clueless on how to do 3 and 4. Could anyone help me with this? I also should mention that I have tried the Temporal Statistics Filter and it seems to
be working fine, but would like to have a way to do it with pythons scripting as well.</div>
<div><br>
</div>
<div><br>
</div>
<div>Thanks,</div>
<div>Parastou</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>#### import the simple module from the paraview<br>
from paraview.simple import *<br>
import numpy as np<br>
#### disable automatic camera reset on 'Show'<br>
paraview.simple._DisableFirstRenderCameraReset()<br>
<br>
dumm_vWSS = []<br>
<br>
for i in range(4500, 4550, 50):<br>
f = "\...\.../all_results_0" + str(i) + ".vtu"<br>
print i<br>
# create a new 'XML Unstructured Grid Reader'<br>
all_results = XMLUnstructuredGridReader(FileName=f)<br>
all_results.CellArrayStatus = ['GlobalElementID']<br>
all_results.PointArrayStatus = ['GlobalNodeID', 'vinplane_traction', 'vWSS', 'average_speed', 'average_pressure', 'pressure', 'velocity', 'timeDeriv']<br>
# create a new 'Calculator'<br>
calculator1 = Calculator(Input=all_results)<br>
calculator1.Function = ''<br>
# Properties modified on calculator1<br>
calculator1.ResultArrayName = 'WSS mag'<br>
calculator1.Function = 'mag(vWSS)'<br>
# get the data points<br>
SetActiveSource(calculator1)<br>
calculator1.UpdatePipeline()<br>
vWSSdata=servermanager.Fetch(calculator1)<br>
vWSSPoints =vWSSdata.GetPoints()<br>
vWSSPointData =vWSSdata.GetPointData()<br>
vWSSArrayData=vWSSPointData.GetArray('WSS mag')<br>
dumm_vWSS.append(vWSSArrayData)<br>
print dumm_vWSS</div>
<br>
</div>
<p></p>
<p><br>
</p>
<div id="Signature">
<div style="font-family: Tahoma; font-size: 13px;"><font face="Times New Roman" size="3">Parastou Eslami, <i>PhD</i></font>
<div><font face="Times New Roman" size="3">Post-doctoral Fellow, Harvard University</font></div>
<div><font face="Times New Roman" size="3">Department of Radiology, Massachusetts General Hospital</font></div>
<div><font face="Times New Roman" size="3"><br>
</font></div>
</div>
</div>
<br>
<p></p>
<p><br>
</p>
<div style="color: rgb(0, 0, 0);">
<div id="divRplyFwdMsg" dir="ltr">
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;">
<p><br>
</p>
</div>
</div>
</div>
</div>
<p></p>
<p>The information in this e-mail is intended only for the person to whom it is<br>
addressed. If you believe this e-mail was sent to you in error and the e-mail<br>
contains patient information, please contact the Partners Compliance HelpLine at<br>
http://www.partners.org/complianceline . If the e-mail was sent to you in error<br>
but does not contain patient information, please contact the sender and properly<br>
dispose of the e-mail.</p></body>
</html>