[Nek5000-users] Instantaneous velocity field

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

Hello,

Thank you for the reply. Yes, outpost() in fact does yield the information I am looking for.  But I hope I could generate the field information in .dat file similar to following pattern:

         open(unit=56,file='vel_fluc_prof.dat')
         write(56,'(A,1pe14.7)') '#time = ', time
         write(56,'(A)')
     &   '#  y     y+     uu     vv     ww     uv'

Could you please teach me how to do this?

Regards,

Tim
________________________________

I forgot to add 'call' before outpost. Please read the block as:



  if (istep .eq. 10) then

       call outpost(vx,vy,vz,pr,t,'vel')

   endif





On 02/21/2018 04:33 PM, Navrose Navrose wrote:

>

> Doesn't outpost() (userchk subroutine) inside a conditional statement

> address your issue?

> Something like (if you have to flow information at istep=10):

>

>  if (istep .eq. 10) then

>       outpost(vx,vy,vz,pr,t,'vel')

>   endif

>

>

> On 02/21/2018 04:22 PM, nek5000-users at lists.mcs.anl.gov<https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users> wrote:

>>

>> Hello,

>>

>> Thank you for the reply. I meant the instantaneous velocity field for

>> the domain (u(x,y,z),v(x,y,z),w(x,y,z)) at some certain time step

>> (istep).

>>

>> Regards,

>>

>> Tim

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180221/937c5d0b/attachment.html>