<div dir="ltr">Hi Hu,<div><br></div><div>Is there any specific reason why you cant use the outpost routine.</div><div>call outpost(vx,vy,vz,pr,t,'his') might be a better option.</div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Best Regards,<br></div>Tanmoy<br></div></div></div>
<br><div class="gmail_quote">On Wed, Mar 8, 2017 at 4:32 AM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><p>hi, all</p>
<p>I'm trying to write down velocity in each processor into files. Below is an sample in the subroutine userchk:</p>
<p>open(unit=80,file='result.dat'<wbr>)<br>write(80,*) 'x, y, z, u, v, z, p'<br>do e=1,<font color="#ff0000">nelv<br></font>   do k=1,lz1<br>      do j=1,ly1<br>         do i=1,lx1<br>            write(80,1) xm1(i,j,k,e),ym1(i,j,k,e),zm1(<wbr>i,j,k,e),<br>&           vx(i,j,k,e),vy(i,j,k,e),vz(i,<wbr>j,k,e),pr(i,j,k,e)<br>         enddo<br>      enddo<br>   enddo<br>enddo<br>1 format(1p15e17.9)</p>
<p>But in this case ,I can get the data only in one processor. When I use more than one processor, how to get velocity in all processors?</p>
<p> </p>
<p>Best regards!</p>
<p>Hu</p><br>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>