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

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Mar 8 18:36:24 CST 2017


Hi Hu,

Is there any specific reason why you cant use the outpost routine.
call outpost(vx,vy,vz,pr,t,'his') might be a better option.

Best Regards,
Tanmoy

On Wed, Mar 8, 2017 at 4:32 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> 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
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170308/e9f409bb/attachment.html>


More information about the Nek5000-users mailing list