Hello,
another newbie question: I am wondering about the coordinate directions
in DX, especially when handling a vector field. My data set consists of
a series of 3D velocities in a 3D volume; each step of the series is
stored as follows:
write (30) (((field1(i,j,k),field2(i,j,k),field3(i,j,k),
& k=1,nz),j=1,ny),i=1,nx)
i.e. vx(1,1,1),vy(1,1,1),vz(1,1,1),vx(1,1,2),... the corresponding
header is shown at the end. The problem is that the vector components
are shifted: rather than being arranged in the order vx,vy,vz, they seem
to be in the order vy,vz,vx! Of course, I am sure that I didn't mix up
the components in the program itself. I wonder if it has to do with the
data format, but I can't figure out what it could have to do with
row/column majority and interleaving. My understanding is that the above
is record-vector, which is the default anyway; it didn't become clear to
me, however, how I could have set that explicitly in the header file.
For a better understanding I have put 2 images under
http://www.geophysik.uni-frankfurt.de/~ruedas/correct.jpg
http://www.geophysik.uni-frankfurt.de/~ruedas/false.jpg
In connection with this, I also wonder why I have to set "(Off) Bottom"
to see the dataset with x horizontal, y away from me and z vertical;
this is what I would expect when setting "(Off) Front", but then I get
actually the view from above. The abovementioned images are "Off
Bottom", BTW.
Somehow it seems that I write the data in a different way from what I
intend, but then, how can I resolve this?
Any comments and hints are welcome.
# grid
object 1 class gridpositions counts 60 14 100
origin 0 0 660.0
delta 25.4237 0 0
delta 0 26.9231 0
delta 0 0 -6.6667
# connections
object 2 class gridconnections counts 60 14 100
attribute "element type" string "cubes"
attribute "ref" string "positions"
# data
object 3 class array type float rank 1 shape 3 items 84000 msb ieee
data file vplt.dat , 0
attribute "dep" string "positions"
# field
object 4 class field
component "positions" value 1
component "connections" value 2
component "data" value 3
# data
object 5 class array type float rank 1 shape 3 items 84000 msb ieee
data file vplt.dat , 1008008
attribute "dep" string "positions"
# field
object 6 class field
component "positions" value 1
component "connections" value 2
component "data" value 5
# series
object 7 class series
member 0 value 4
member 1 value 6
#
end
--
------------------------------------------------------------------------
Thomas Ruedas
Institute of Meteorology and Geophysics, J.W.Goethe University Frankfurt
e-mail: [EMAIL PROTECTED]
http://www.geophysik.uni-frankfurt.de/~ruedas/
------------------------------------------------------------------------
