<div dir="ltr"><div class="gmail_default" style="font-size:small">Thank you Paul. Now comp_vort3 doesnt give zeros in the solid domain.</div><div class="gmail_default" style="font-size:small">What is the difference between comp_vort3 and the calculation of curl using </div><div class="gmail_default"> call gradm1(dAxdx,dAxdy,dAxdz,t(1,1,1,1,1)) </div><div class="gmail_default"> call gradm1(dAydx,dAydy,dAydz,t(1,1,1,1,2)) </div><div class="gmail_default"> call gradm1(dAzdx,dAzdy,dAzdz,t(1,1,1,1,3))</div><div class="gmail_default"> call sub3(bfx1,dAzdy,dAydz,m) </div><div class="gmail_default"> call sub3(bfy1,dAxdz,dAzdx,m) </div><div class="gmail_default"> call sub3(bfz1,dAydx,dAxdy,m) </div><div class="gmail_extra"><div class="gmail_default" style="font-size:small">When I compare the glsc3 of a curls computed by comp_vort3 and gradm1, there is a difference of 0.5% to 1%.</div><div class="gmail_default" style="font-size:small">Which procedure suits well, if I want to use u x curl in the userq</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Sincerely,</div><div class="gmail_default" style="font-size:small">Sandeep</div><br></div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
<br>
<br>
Hi Sandeep,<br>
<br>
This issue should now be resolved.<br>
<br>
comp_vort3() will now compute the curl of the supplied vector field over all elements, e=1,...,nelt.<br>
<br>
Moreover, prepost() has been modified to not put zeros in the velocity field when e > nelv.<br>
<br>
So, you should be able to do something like:<br>
<br>
call comp_vort3(curlv,w1,w2,u,v,w)<br>
<br>
call outpost(curlv(1,1),curlv(1,2),<wbr>culrv(1,3),pr,t,' ')<br>
<br>
and have the right thing happen.<br>
<br>
Paul<br>
______________________________<wbr>__<br>
From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.mc<wbr>s.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.m<wbr>cs.anl.gov</a>] on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>]<br>
Sent: Thursday, March 23, 2017 2:57 PM<br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
Subject: [Nek5000-users] comp_vort3, nelt, conj_ht<br>
<br>
Hello Neks,<br>
I modified conj_ht example and defined three scalars. NEL>NELV. The domain of scalars is bigger than the fluid.<br>
I create a vector by taking the curl of three scalars, using<br>
common /myvort/ b1(lx1,ly1,lz1,lelt,3),<br>
call comp_vort3(b1,work1,work2,t(1,<wbr>1,1,1,1),<br>
$ t(1,1,1,1,2),t(1,1,1,1,3))<br>
<br>
b1 is computed correctly in the fluid domain, but is solid domain it showing zero.<br>
Can someone tell me how to modify the comp_vort3 subroutine, to correctly compute b1 over the entire domain.<br>
I tried to modify call dudxyz(work1,w,rym1,sym1,tym1,<wbr>jacm1,2,2), changing the imsh=2 (the last second parameter in the above call), but the the result remains the same.<br>
I use the three components of computed b1 in userq<br>
Thank you in advance.<br>
<br>
Sincerely,<br>
Sandeep<br>
<br><br></blockquote></div></div></div>