| Re: Reading in to array elements from a file. [message #12272] |
Wed, 22 July 1998 00:00 |
davidf
Messages: 2866
Registered: September 1996
|
Senior Member |
|
|
Neil Rogers (nrogers@dera.gov.uk) writes:
> I wonder if anyone could solve the following mystery for me.
> I am trying to read in *individual* array elements from a file using
> the readf procedure. I have an input file called test.txt containing
> (say)
> 1 2 3 4 5 6
> and the following IDL procedure to read them in and print them out
> again.
> ;------------------
> pro test
> array=intarr(6)
> infile='test.txt'
> openr,inlun,infile,/get_lun
> for cnt=0,5 do begin
> readf,inlun,array[cnt]
> endfor
> close,inlun
> free_lun,inlun
> print,array
> return
> end
> ;-----------------
> Unfortunately, although no compile or run-time errors are encountered,
> nothing is assigned to the array elements and the file pointer stays
> at the beginning of the file. The output is then,
> 0 0 0 0 0 0
> Why is this?
The reason you can't read into a subscripted array has to
do with the way IDL passes information into and out of
procedures and functions. Variables are passed by
*reference*, meaning that it is actually the pointer
(a _real_ pointer here) or address that is passed into the
procedure or function. Changes made to the variable inside
the procedure or function are made to the "real thing".
For example, consider this procedure:
PRO JUNK, variable
variable = 10
END
If I call it like this:
thisVar = 5
Junk, thisVar
Print, thisVar
10
We see that thisVar has been changed by the procedure.
But variables are the *only* things that are passed by
reference. Everything else (and this includes things
like structure defererences, expressions, constants,
and--yes--subscripted variables) are passed by *value*.
This means that a COPY of the thing itself is passed
into the procedure or function.
So, for example, in our little example if we subscript
the variable like this then the variable outside the
procedure doesn't change:
thisVar = 5
Junk, thisVar[0]
Print, thisVar
5
Notice that this code doesn't cause any error messages.
In fact it is perfectly valid IDL code to pass a parameter
by value and not by reference.
As it happens, this is *exactly* what happens when you
do this:
array = IntArr(6)
For j=0,5 DO ReadF, 1, array[j]
The value array[j] is passed by VALUE, not by REFERENCE
and thus it's value outside the ReadF procedure can never
change, although internally to the ReadF procedure it is
doing exactly what it is suppose to be doing (i.e., reading
a value from the file into a *copy* of a subscripted
variable).
This behavior is so fundamental to the way IDL works that
it is unlikely to change in our lifetime. :-)
Cheers,
David
-----------------------------------------------------------
David Fanning, Ph.D.
Fanning Software Consulting
E-Mail: davidf@dfanning.com
Phone: 970-221-0438
Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
|
| Re: Reading in to array elements from a file. [message #12273 is a reply to message #12272] |
Wed, 22 July 1998 00:00  |
Kevin Ivory
Messages: 71
Registered: January 1997
|
Member |
|
|
[original message cited below]
Hi Neil,
There are two reasons the readf doesn't work as expected:
1. You cannot read into an array element. The only way to put a value
into a single array element is to read it and assign it in the next
statement:
value = 0
for i = 0, 5 do begin
readf, inlun, value
array[i] = value
endfor
2. Each readf statement reads a new line, so if there is more than one
value per line, have to get all values from that line at once.
Instead of using the for loop, you should get the complete array in
one read. This way, you don't have to worry about the items mentioned
above.
readf, inlun, array
Hope this helps,
Kevin
--
Kevin Ivory Tel: +49 5556 979 434
Max-Planck-Institut fuer Aeronomie Fax: +49 5556 979 240
Max-Planck-Str. 2 mailto:Kevin.Ivory@linmpi.mpg.de
D-37191 Katlenburg-Lindau, GERMANY http://www.gwdg.de/~kivory2/
Neil Rogers wrote:
>
> I wonder if anyone could solve the following mystery for me.
> I am trying to read in *individual* array elements from a file using
> the readf procedure. I have an input file called test.txt containing
> (say)
> 1 2 3 4 5 6
> and the following IDL procedure to read them in and print them out
> again.
> ;------------------
> pro test
> array=intarr(6)
> infile='test.txt'
> openr,inlun,infile,/get_lun
> for cnt=0,5 do begin
> readf,inlun,array[cnt]
> endfor
> close,inlun
> free_lun,inlun
> print,array
> return
> end
> ;-----------------
> Unfortunately, although no compile or run-time errors are encountered,
> nothing is assigned to the array elements and the file pointer stays
> at the beginning of the file. The output is then,
> 0 0 0 0 0 0
> Why is this?
> Thanks for your help.
>
> Neil.
--
Kevin Ivory Tel: +49 5556 979 434
Max-Planck-Institut fuer Aeronomie Fax: +49 5556 979 240
Max-Planck-Str. 2 mailto:Kevin.Ivory@linmpi.mpg.de
D-37191 Katlenburg-Lindau, GERMANY http://www.gwdg.de/~kivory2/
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
