<HTML><BODY>Hi, Neks!<br><br>1) I am doing a DNS study of turbulent pipe jets and looking at turbulent kinetic energy budgets. <br><span id="result_box" class="short_text" lang="en"><span><br>I take derivatives via gradm1 subroutine, for example:<br> <br></span></span><em>call gradm1(uxdx,uxdy,uxdz,ux)<br><br></em><span>After all, if I look at my dissipation field (via outpost), I see some wrong values on the surface between spectral elements. <br>I want to filter them using:<br><br><em>call col2(<em>uxdx</em>,bm1,ntot) </em><br><em>call dssum(<em>uxdx</em>,lx1,ly1,lz1) </em><br><em>call col2(<em>uxdx</em>,binvm1,ntot)</em><br><br><em>call dsavg(dxux)</em><br><br>for every array. But <span id="result_box" lang="en"><span>conversely it leads <span id="result_box" class="short_text" lang="en"><span>to an increase in gaps. What am I doing wrong?<br><br></span></span></span></span>2) How does a filter-function in <span><em>dssum </em>look and how much points does it use for filtering? <span></span></span><br></span><br>Best regards,<br>Vlad.<br></BODY></HTML>