Dear Sandeep, Please take a look at hemi.usr in the nekexamples repo. There is a routine interp_v() that returns the interpolated velocity field at a user-prescribed set of points. You could easily modify this routine to pass in vorticity and interpolate that (or any other scalar or vector field). Please note that you should not duplicate interrogation points on multiple processors (i.e., each processor should have a unique set of points, xyz(); one possible configuration is that MPI rank 0 has a list of points of length n and the other ranks have a list of length 0; this approach is valid if the length of the list is not too large). hth, Paul ________________________________ From: Nek5000-users <[email protected]> on behalf of [email protected] <[email protected]> Sent: Tuesday, December 5, 2017 8:57:20 AM To: [email protected] Subject: Re: [Nek5000-users] hpts for vorticity field poor man`s solution: Store your vorticity as an additional scalar field. -----Original message-----

From:[email protected] <[email protected]> Sent: Tuesday 5th December 2017 15:55 To: [email protected] Subject: [Nek5000-users] hpts for vorticity field

Hello Neks, Does anyone knows how to output the vorticity probes using call hpts() When I use the function "call hpts()", it outputs only the probes for native arrays like velocity, prssure (and scalar fields, if they exist). 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? i will be grateful for your answers.

Thank you in advance.

Cheers, Sandeep _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users