[Nek5000-users] comp_vort3, nelt, conj_ht

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

Hi Sandeep,

This issue should now be resolved.

comp_vort3() will now compute the curl of the supplied vector field over all elements, e=1,...,nelt.

Moreover, prepost() has been modified to not put zeros in the velocity field when e > nelv.

So, you should be able to do something like:

       call comp_vort3(curlv,w1,w2,u,v,w)

       call outpost(curlv(1,1),curlv(1,2),culrv(1,3),pr,t,'   ')

and have the right thing happen.

Paul
________________________________
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, March 23, 2017 2:57 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] comp_vort3, nelt, conj_ht

Hello Neks,
I modified conj_ht example and defined three scalars. NEL>NELV. The domain of scalars is bigger than the fluid.
I create a vector by taking the curl of three scalars, using
  common /myvort/ b1(lx1,ly1,lz1,lelt,3),
  call comp_vort3(b1,work1,work2,t(1,1,1,1,1),
$     t(1,1,1,1,2),t(1,1,1,1,3))

b1 is computed correctly in the fluid domain, but is solid domain it showing zero.
Can someone tell me how to modify the comp_vort3 subroutine, to correctly compute b1 over the entire domain.
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.
I use the three components of computed b1 in userq
Thank you in advance.

Sincerely,
Sandeep


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170325/840abed2/attachment.html>