[Nek5000-users] Pauls v2 Integration routine...
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Thu Oct 27 13:02:46 CDT 2016
Paul,
Hi... I'd like to use this also; however, was your intent for the
argument of the facint_v() sub to be v2, and not v? Also, would you
want to pass v2surf back in /mystuff/?
if (cbc(f,e,1).eq.'v ') v2surf=v2surf+facint_v(v,area,f,e)
Thanks,
Murph
-----------------------------------------------------------------------------------------------------------
common /mystuff/ v2(lx1*ly1*lz1*lelt)
integer e,f
n = nx1*ny1*nz1*nelt
do i=1,n
v2(i)=vx(i,1,1,1)**2+vy(i,1,1,1)**2+vz(i,1,1,1)**2
enddo
v2surf = 0
do e=1,nelv ! Integrate over all surfaces with cbc='v '
do f=1,2*ndim
if (cbc(f,e,1).eq.'v ') v2surf=v2surf+facint_v(v,area,f,e)
enddo
enddo
v2surf = glsum(v2surf,1) ! Sum across all processors
Note that "area" is passed in the TOTAL common block, which I assume
you are including in the routine.
On 10/26/2016 1:00 PM, nek5000-users-request at lists.mcs.anl.gov wrote:
> Send Nek5000-users mailing list submissions to
> nek5000-users at lists.mcs.anl.gov
>
> To subscribe or unsubscribe via the World Wide Web, visit
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> or, via email, send a message with subject or body 'help' to
> nek5000-users-request at lists.mcs.anl.gov
>
> You can reach the person managing the list at
> nek5000-users-owner at lists.mcs.anl.gov
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Nek5000-users digest..."
>
>
> Today's Topics:
>
> 1. Order of appearance : USERBC and USERCHECK
> (nek5000-users at lists.mcs.anl.gov)
> 2. Re: Integrate velocity field along a line
> (nek5000-users at lists.mcs.anl.gov)
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161027/d739e2fb/attachment.html>
More information about the Nek5000-users
mailing list