[Nek5000-users] comp_vort3, conj_ht

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

I guess the routine is hardwired to to compute the curl of a field defined on the velocity mesh.

On 9 Mar 2018, at 15:11, "nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> " <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> > wrote:

Hello Nek team,

In conj_ht example, I compute curl 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 being computed correctly in the fluid domain, but in the solid domain it is showing zero.


Can you please correct it to compute the curl in the solid domain?


Thanks in advance.






_______________________________________________

Nek5000-users mailing list

Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov> 

https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180310/21f11cee/attachment.html>