Dear NEKs I think I could fix the problem regarding mid-plane flux computation somehow. I have attached the routine in case any users run into the same issue in the future! Thanks Azad P.S. Just wanted to mention I add the plane number manually (if (ex.eq.4) then) since I had 6 element in the x-direction. One could change it with respect to the domain.

Date: Thu, 01 Mar 2012 22:09:57 +0100 From: [email protected] Subject: [Nek5000-users] middle plane flux To: [email protected] Message-ID: <[email protected]> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"; DelSp="Yes"

Dear NEKs

I am using NEK5000 to simulate turbulent flow in duct (channel with side walls). I would like to calculate flux at the middle slab of the domain (a zy plane) in each time step in order to compare it with my bulk (Z is my streamwise direction). Attached file is the .usr file containing the routine I wrote to compute the mid-flux. Seems works fine for the primitive tests; however, when I increase number of processors to 256 or upper it gives me NaNs. I was wondering if you could help me with that.

With many thanks in advance. Azad

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

End of Nek5000-users Digest, Vol 37, Issue 1 ********************************************