[Nek5000-users] Pauls v2 Integration routine...

[nek5000-users at lists.mcs.anl.gov]

Yes, Murph.

Aleks
________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Thursday, October 27, 2016 1:02 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Pauls v2 Integration routine...


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<mailto:nek5000-users-request at lists.mcs.anl.gov> wrote:

Send Nek5000-users mailing list submissions to
        nek5000-users at lists.mcs.anl.gov<mailto: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<mailto: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<mailto: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<mailto:nek5000-users at lists.mcs.anl.gov>)
   2. Re: Integrate velocity field along a line
      (nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>)




_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto: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/924b60ca/attachment-0001.html>