[Nek5000-users] output filed

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Tue Feb 23 12:16:49 CST 2016


Hi neks,
i need to print the output field u,v,w for every x coordinate.
Now in userchk i wrote:

         do i=1,ntot

          write(16,*) istep,time, xm1(i,1,1,1),
      & vx(i,1,1,1),vy(i,1,1,1),vz(i,1,1,1)

         enddo

and i receive the following output:

            1   7.1530856036040695E-003   0.0000000000000000 
0.37293630591777449      -0.32084904392637759 0.33230910496406607
            1   7.1530856036040695E-003  0.17357493101301319 
0.21130775186288545      -0.27768404661864132 0.20739454454471290
            1   7.1530856036040695E-003  0.45442506898698676 
0.46084398518785818      -0.27156542355496571 0.18149139586311763
            1   7.1530856036040695E-003  0.62800000000000000 
0.46110890335130450       -3.9277174842836750E-002 8.5610114286700789E-002
            1   7.1530856036040695E-003   0.0000000000000000 
0.74352058860894332       -7.5686819189858645E-002 1.3289264769474618E-002
            1   7.1530856036040695E-003  0.17357493101301319 
0.40009961054319920       0.13480828683699964 0.44478023194716509
            1   7.1530856036040695E-003  0.45442506898698676 
0.81382696494932794      -0.11005654233044537 0.17652132504468790
            1   7.1530856036040695E-003  0.62799999999999989 
0.67943618184656285       -9.6081837856840727E-002 -0.14185885411564772
            1   7.1530856036040695E-003   0.0000000000000000 
0.32147262016787648      -0.24804263078583569 -0.29171129659235678
            1   7.1530856036040695E-003  0.17357493101301319 
0.37769256515663086      -0.15125228365962307 -3.2429247313227494E-002
            1   7.1530856036040695E-003  0.45442506898698676 
0.24361827129204860       -3.1253044882220637E-002 0.19706743377184774
            1   7.1530856036040695E-003  0.62799999999999989 
0.63567408933529579       0.14261690852412279 -0.10888850075602638
            1   7.1530856036040695E-003   0.0000000000000000 
0.97344927314055640       0.24670700538197837 0.11490941870543755
            1   7.1530856036040695E-003  0.17357493101301319 
0.70722038909668339       0.30861285610239220 0.13262927045873465
            1   7.1530856036040695E-003  0.45442506898698676 
-0.26663222342418402        8.8968243007683773E-002 -5.8482352542457344E-002
            1   7.1530856036040695E-003  0.62800000000000000 
0.68703349900478228        5.2212760358220950E-002 -0.20495777585812078

I see that in x i have, every 4 rows, the numbers are repeated.
Can someone tell me how should i read this numbers? I think they are 
vertices of every cell... it's all right?
How can i have only x vx vy vz ordered by x from low to high?


More information about the Nek5000-users mailing list