hello, So I followed the advice.... from the 32400 records that I
initially had, now all my coordinate arrays (sx,sy,sz) and velocity
arrays (vx,vy,vz) are of 521 records each (i.e. sxx, syy, szz and vxx,
vyy and vzz).
However, I am getting and error message that I have not been able to
fix...
IDL> .COMPILE "C:\RSI\IDL55\lib\mycount.pro"
% Compiled module: $MAIN$.
IDL> .CONTINUE
% Type of FOR statement index variable <No name> may not be changed.
% Execution halted at: $MAIN$ 309
C:\RSI\IDL55\lib\mycount.pro
How can I fix this?
Thank you in advance for any asistance that you can give me.
Regards,
Adhara
Here is the new version of my code:
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 = 0
For i = 0, 32399L Do Begin ;change here the number of records
ifneeded.
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
Vxx=WHERE(Vx NE 0,count1)
Vyy=WHERE(Vy NE 0,count2)
Vzz=WHERE(Vz NE 0,count3)
sx=A[2,*]
sy=A[1,*]
sz=A[0,*]
Sxx=where(Vxx)
Syy=where(Vyy)
Szz=where(Vzz)
vol = FLTARR(360, 360, 90)
WINDOW, XSIZE = 390, YSIZE = 400
CREATE_VIEW, XMAX = 360, YMAX = 360, ZMAX = 100
Scale3, zr=[0,10], yr=[0,120], xr = [0,120] ;create 3D scaling
system
Scale3, zr=[0,90], yr=[0,360], xr = [0,360]
Flow3, Vxx, Vyy, Vzz, SXX=Reform(sxx), SYY=Reform(syy),
SZZ=Reform(szz)
end
eoraptor@aol.com (Richard Hoffpauir) wrote in message news:<77bd9c11.0405151543.e20d23b@posting.google.com>...
> I'm trying to help Adhara get her flow velocity data to plot using
> FLOW3.
>
> FLOW3, Vx, Vy, Vz [,ARROWSIZE=value] [,/BLOB] [,LEN=value]
> [,NSTEPS=value] [,NVECS=value] [,SX=vector,SY=vector,SZ=vector]
>
> She has 3 float arrays named Vx, Vy, and Vz that are of dimensioned
> (90,360,360). Vx contains the magnitude of the flow vector in the x
> direction, Vy is the magnitude in the y direction, etc.
>
> The actual values in the arrays Vx, Vy, and Vz are almost all zero.
> The few values that are non-zero are +/- in the range of e-005.
>
> Since most of the values in her arrays are zero, should she be using
> WHERE to extract out the non-zero magnitudes before calling FLOW3?
> The arrays of SX, SY, and SZ in FLOW3 could be used to give the x,y,z
> coordinates of the non-zero magnitudes, right?
>
> Thanks,
> Richard
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
