[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