<div dir="ltr"><div class="gmail_default" style="font-size:small">Dear Neks,</div><div class="gmail_default" style="font-size:small">I solve</div><div class="gmail_default" style="font-size:small">\partial A/\partial t = \nabla^2 A + u x Curl A,</div><div class="gmail_default" style="font-size:small">where A = (s1,s2,s3) </div><div class="gmail_default" style="font-size:small">I solve three scalar equations with respective components of u x Curl A as forcing terms in qvol.</div><div class="gmail_default" style="font-size:small">I would like to compute the dealiased product of u x Curl A</div><div class="gmail_default" style="font-size:small">I build my code based on kov_stokes_test.usr and  subroutines in induct.f</div><div class="gmail_default" style="font-size:small">I compute local curl, project curl and U onto fine mesh, perform product and project the product back to coarse mesh using interp_rstd<br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default" style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial">However, I am not sure about the multiplication of mass matrices for the product in these subroutines.<br></div><br></div><div class="gmail_default" style="font-size:small">I attach my usr, rea and SIZE files.</div><div class="gmail_default" style="font-size:small">Your suggestions are welcome.</div><div class="gmail_default" style="font-size:small">Thanks,</div><div class="gmail_default" style="font-size:small">Sandeep</div><div><div class="gmail_signature" data-smartmail="gmail_signature"></div></div>
</div>