Dear all,
I have a very large set of data in an array of 6 columns and 32400
rows. The columns are the coordinates x,y,z and the velocity
components Vx,Vy,Vz. There is a large number of these points that are
zero for vx, vy, vz. However, I my macro is plotting these points.
How can I avoid this?
Regards,
Adhara
P.D. Also...If I decide to get rid of the zeros, how can I define the
array? I say this because, every "z" will have different points (x,y),
at which (Vx, Vy, Vz) is not equal to zero.
Thanks for any advice, I know that you all are very helpfull !
Adhara
openr,q,'VOa.dat',/Get_Lun ;open files for reading
A=fltarr(6,90*360L) ;create arrays to hold data in files
readf,q,A ;read files into arrays
Free_Lun,q ;close files and clear logical unit assigned to each
file
Vx=fltarr(360,360,90) ;create 3 dimensional arrays to hold vector
data.
Vy=fltarr(360,360,90) ;the indices will be Vx(x,y,z)
Vz=fltarr(360,360,90)
factor = 1
For i = 0, 32399L Do Begin ;change here the number of records
ifneeded.
;(0) first column of the array is z (0<z<90)
;(1) second column of the array is y (constant)
;(2) third column of the array is x (0<x<360)
;(3) fourth column of the array is Vx
;(4) fifth column of the array is Vy
;(5) sixth column of the array is Vz
Vx(A(2,i)-1,A(1,i)-1,A(0,i)-1) = A(3,i)*10.0^factor
Vy(A(2,i)-1,A(1,i)-1,A(0,i)-1) = A(4,i)*10.0^factor
Vz(A(2,i)-1,A(1,i)-1,A(0,i)-1) = A(5,i)*10.0^factor
Endfor
sx=A[2,*]
sy=A[1,*]
sz=A[0,*]
vol = FLTARR(360, 360, 90)
WINDOW, XSIZE = 390, YSIZE = 400
;CREATE_VIEW, XMAX = 360, YMAX = 360, ZMAX = 100
Scale3, zr=[0,90], yr=[0,360], xr = [0,360] ;create 3D scaling
system
Flow3, Vx, Vy, Vz, SX=Reform(sx), SY=Reform(sy), SZ=Reform(sz)
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
