[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
Thu Mar 9 04:46:59 CST 2017
Thanks Paul! I'll try.
Best regards!
Hu
-----Original Messages-----
From: nek5000-users at lists.mcs.anl.gov
Sent Time: Thursday, March 9, 2017
To: "nek5000-users at lists.mcs.anl.gov" <nek5000-users at lists.mcs.anl.gov>
Cc:
Subject: Re: [Nek5000-users] how to dumping velocity data in each processor to file
You can set param(66)=0 in usrdat2:
param(66)=0
this will cause the output to be in ascii, in a .fld file...
I wouldn't recommend this for large parallel cases, however, as the ascii output
is not fast.
Paul
From:nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Wednesday, March 08, 2017 8:27 PM
To:nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] how to dumping velocity data in each processor to file
Hi Tanmoy,
Thanks! In fact I can use the outpost routine, but I don't know how to read the hisxxxx.f0001 file except by VisIt . I want to write down the data in a familar format so that I can handle it very easy. Or can I transform the outpost format into ASCII that I can also read it easy?
Best regards!
Hu
-----Original Messages-----
From:nek5000-users at lists.mcs.anl.gov
Sent Time: Thursday, March 9, 2017
To:nek5000-users at lists.mcs.anl.gov
Cc:
Subject: Re: [Nek5000-users] how to dumping velocity data in each processor to file
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/20170309/2f0e5f12/attachment-0001.html>
More information about the Nek5000-users
mailing list