<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>