<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
<div><br>
</div>
Hi Sandeep,
<div><br>
</div>
<div>This issue should now be resolved.</div>
<div><br>
</div>
<div>comp_vort3() will now compute the curl of the supplied vector field over all elements, e=1,...,nelt.</div>
<div><br>
</div>
<div>Moreover, prepost() has been modified to not put zeros in the velocity field when e > nelv.</div>
<div><br>
</div>
<div>So, you should be able to do something like:</div>
<div><br>
</div>
<div> call comp_vort3(curlv,w1,w2,u,v,w)</div>
<div><br>
</div>
<div> call outpost(curlv(1,1),curlv(1,2),culrv(1,3),pr,t,' ')</div>
<div><br>
</div>
<div>and have the right thing happen.</div>
<div><br>
</div>
<div>Paul</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF405517" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov [nek5000-users-bounces@lists.mcs.anl.gov] on behalf of nek5000-users@lists.mcs.anl.gov [nek5000-users@lists.mcs.anl.gov]<br>
<b>Sent:</b> Thursday, March 23, 2017 2:57 PM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] comp_vort3, nelt, conj_ht<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_default" style="font-size:small">Hello Neks,</div>
<div class="gmail_default" style="font-size:small">I modified conj_ht example and defined three scalars. NEL>NELV. The domain of scalars is bigger than the fluid.</div>
<div class="gmail_default" style="font-size:small">I create a vector by taking the curl of three scalars, using</div>
<div class="gmail_default" style="font-size:small">
<div class="gmail_default"> common /myvort/ b1(lx1,ly1,lz1,lelt,3),</div>
<div class="gmail_default"> call comp_vort3(b1,work1,work2,t(1,1,1,1,1),<br>
</div>
</div>
<div class="gmail_default" style="font-size:small">
<div class="gmail_default">$ t(1,1,1,1,2),t(1,1,1,1,3))</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default">b1 is computed correctly in the fluid domain, but is solid domain it showing zero. </div>
<div class="gmail_default">Can someone tell me how to modify the comp_vort3 subroutine, to correctly compute b1 over the entire domain.</div>
<div class="gmail_default">I tried to modify call dudxyz(work1,w,rym1,sym1,tym1,jacm1,2,2), changing the imsh=2 (the last second parameter in the above call), but the the result remains the same.</div>
<div class="gmail_default">I use the three components of computed b1 in userq</div>
<div class="gmail_default">Thank you in advance.</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default">Sincerely,</div>
<div class="gmail_default">Sandeep</div>
<div class="gmail_default"><br>
</div>
<div><br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>