<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Paul,</p>
<p><br>
</p>
<p>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/?</p>
<p><br>
</p>
<p><i class="moz-txt-slash"><span class="moz-txt-tag"></span><span
class="moz-txt-tag"></span></i> if (cbc(f,e,1).eq.'v
') v2surf=v2surf+facint_v(v,area,f,e)<br>
<br>
</p>
<p>Thanks,</p>
<p>Murph</p>
<p><br>
</p>
<p>-----------------------------------------------------------------------------------------------------------<br>
</p>
<p><br>
</p>
<p>common <i class="moz-txt-slash"><span class="moz-txt-tag">/</span>mystuff<span
class="moz-txt-tag">/</span></i> v2(lx1*ly1*lz1*lelt)
<br>
integer e,f
<br>
<br>
n = nx1*ny1*nz1*nelt
<br>
<br>
do i=1,n
<br>
v2(i)=vx(i,1,1,1)**2+vy(i,1,1,1)**2+vz(i,1,1,1)**2
<br>
enddo
<br>
<br>
v2surf = 0
<br>
do e=1,nelv ! Integrate over all surfaces with cbc='v '
<br>
do f=1,2*ndim
<br>
if (cbc(f,e,1).eq.'v ')
v2surf=v2surf+facint_v(v,area,f,e)
<br>
enddo
<br>
enddo
<br>
<br>
v2surf = glsum(v2surf,1) ! Sum across all processors
<br>
<br>
Note that "area" is passed in the TOTAL common block, which I
assume
<br>
you are including in the routine.
<br>
</p>
<br>
<div class="moz-cite-prefix">On 10/26/2016 1:00 PM,
<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a> wrote:<br>
</div>
<blockquote
cite="mid:mailman.43.1477501210.31010.nek5000-users@lists.mcs.anl.gov"
type="cite">
<pre wrap="">Send Nek5000-users mailing list submissions to
<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>
To subscribe or unsubscribe via the World Wide Web, visit
<a class="moz-txt-link-freetext" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>
or, via email, send a message with subject or body 'help' to
<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a>
You can reach the person managing the list at
<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users-owner@lists.mcs.anl.gov">nek5000-users-owner@lists.mcs.anl.gov</a>
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Nek5000-users digest..."
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">Today's Topics:
1. Order of appearance : USERBC and USERCHECK
(<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)
2. Re: Integrate velocity field along a line
(<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Nek5000-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a>
<a class="moz-txt-link-freetext" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>
</pre>
</blockquote>
<br>
</body>
</html>