Mani, A variable temp can be used only in the subroutines which are called with the arguments of collocation point coordinates and global element number: ix,iy,iz,ieg These routines are userf, userq, userbc and useric where temp has the following meaning (updated through NEKUSE) temp = t(ix,iy,iz,gllel(ieg)) Note that gllel(ieg) returns the local element number on a processor for a global element number ieg. Similar meaning there have ux,uy,uz and ps(1),ps(2)...,ps(ldimt), namely, velocity components and passive scalars at the point ix,iy,iz,ieg Best, Aleks ----- Original Message ----- From: [email protected] To: [email protected] Sent: Thursday, September 2, 2010 7:51:55 AM GMT -06:00 US/Canada Central Subject: Re: [Nek5000-users] [*] Re: Function to calculate volume intergral Thanks! I tried that earlier but it was seg faulting cause I put in temperature as 'temp' and not 't'. Mani On 09/02/2010 05:56 PM, [email protected] wrote:

Hi Mani,

You can use

n = nx1*ny1*nz*nelv vytvol = glsc3(vy,t,bm1,n)

for the volume integral of Vy and temperature

Best, Aleks

----- Original Message ----- From: [email protected] To: [email protected] Sent: Thursday, September 2, 2010 12:38:10 AM GMT -06:00 US/Canada Central Subject: [Nek5000-users] Function to calculate volume intergral

Hi,

Is there a function one can use in Nek that computes volume integrals. I'd like to compute the Nusselt number for Rayleigh-Benard convection (=1 +<Vy*Temp>), where Vy is the vertical velocity and<> is the volume integral.

Thanks, Mani _______________________________________________ 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

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