[Nek5000-users] comp_vort3, nelt, conj_ht

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

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/20170323/4e077852/attachment.html>