<div dir="ltr"><div>Hi, neks.<br></div>I am trying to get the 1d profile. vx(y) in turbChannel, for or example. Here the code below:<br><br>real myvx(ny1,nely)<br><br>do e=1,nelv<br><br> eg = gllel(e)<br> call get_exyz(ex,ey,ez,eg,nelx,nely,nelz)<br><br> if (ex.eq.1 .and. ez.eq.1) then<br> do i=1,ny1<br> myvx(i,ey) = vx(1,i,1,e)<br> enddo<br> endif<br><br>enddo<br><div><br></div><div>But after this, myvx contains incorrect value, that not looks like really vx(y).<br></div><div>What's wrong?<br></div><div><br clear="all"><div><div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">Vatslav</div></div></div></div></div>
</div></div></div></div>