<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<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:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Dear Sandeep,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Please take a look at hemi.usr in the nekexamples repo.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">There is a routine interp_v()  that returns the interpolated</p>
<p style="margin-top:0;margin-bottom:0">velocity field at a user-prescribed set of points.  You could</p>
<p style="margin-top:0;margin-bottom:0">easily modify this routine to pass in vorticity and interpolate</p>
<p style="margin-top:0;margin-bottom:0">that (or any other scalar or vector field).</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Please note that you should not duplicate interrogation points</p>
<p style="margin-top:0;margin-bottom:0">on multiple processors (i.e., each processor should have a </p>
<p style="margin-top:0;margin-bottom:0">unique set of points, xyz();  one possible configuration is that</p>
<p style="margin-top:0;margin-bottom:0">MPI rank 0 has a list of points of length n and the other ranks</p>
<p style="margin-top:0;margin-bottom:0">have a list of length 0;  this approach is valid if the length of the</p>
<p style="margin-top:0;margin-bottom:0">list is not too large).</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">hth,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Paul</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Sent:</b> Tuesday, December 5, 2017 8:57:20 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] hpts for vorticity field</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:10pt;">
<div class="PlainText">poor man`s solution: Store your vorticity as an additional scalar field.
<br>
<br>
-----Original message-----<br>
> From:nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
> Sent: Tuesday 5th December 2017 15:55<br>
> To: nek5000-users@lists.mcs.anl.gov<br>
> Subject: [Nek5000-users] hpts for vorticity field<br>
> <br>
> Hello Neks,<br>
> Does anyone knows how to output the vorticity probes using call hpts()<br>
> When I use the function "call hpts()", it outputs only the probes for native arrays like velocity, prssure (and scalar fields, if they exist).<br>
> I compute the vorticity using the function comp_vort3(vort,w1,w2,vx,vy,vz). How can I output the values of the vorticity probes in the flow field using hpts?<br>
> i will be grateful for your answers.<br>
> <br>
> Thank you in advance.<br>
> <br>
> Cheers,<br>
> Sandeep<br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div>
</span></font></div>
</body>
</html>