[Nek5000-users] snapshots in ascii file

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

Hi Neks,

I would like to save a data set of snapshots in ascii, containing 
position and velocity.

I've tried to modify the userchk subroutine in the following way:

       real us(lx1,ly1,lz1,lelt), vs(lx1,ly1,lz1,lelt)
       integer i, j, k, l

       call opcopy(us,vs,ws,vx,vy,vz)

       ...

       open(223,file='example.dat')

       do i=1,lx
          do j=1,ly
            do k=1,lz
               do l=1,lelt
                  write(223,*) 
xm(i,j,k,l),ym(i,j,k,l),us(i,j,k,l),vs(i,j,k,l)
               end do
            end do
          end do
       end do

       close(223)


but the example.dat file is empty.

Could you help me?

SL