[Nek5000-users] how to dumping velocity data in each processor to file

[nek5000-users at lists.mcs.anl.gov]

hi, all

I'm trying to write down velocity in each processor into files. Below is an sample in the subroutine userchk:

open(unit=80,file='result.dat')
write(80,*) 'x, y, z, u, v, z, p'
do e=1,nelv
   do k=1,lz1
      do j=1,ly1
         do i=1,lx1
            write(80,1) xm1(i,j,k,e),ym1(i,j,k,e),zm1(i,j,k,e),
&           vx(i,j,k,e),vy(i,j,k,e),vz(i,j,k,e),pr(i,j,k,e)
         enddo
      enddo
   enddo
enddo
1 format(1p15e17.9)

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?

 

Best regards!

Hu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170308/3d542c7f/attachment.html>