function find_a, lookhere, filenames = thesefiles numparams = n_params() if numparams eq 0 then cd, current = lookhere cd, lookhere ;---------RHESSI--------- ;reads ascii file and returns it as a text array ;c = rd_ascii('flux_output_april21_12-25_c50.txt') c = rd_ascii(thesefiles) ;removes top line which was: ;box date stime etime flux area centr peak st.dev peak.val Elow Ehi detectors c = c[1:*] ;makes 1 line string two rows b = str2cols(c) ;changes the dims of array (does not change number of elements) ;box is all the box numbers box = fix(reform(b[0,*])) ;area in arcsec^2 areas = float(reform(b[5,*])) q = where(box eq 0) ;q=lines where box=0 (0,25...) b0 = b[*,q] ;b0=all full lines where box=0 ;all boxs are now 0 box = (box[q]) ;areas with box number 0 areas = (areas[q]) box = fix(reform(b0[0,*])) areas = float(reform(b0[5,*])) ;print,b0[*,0] ;prints all columns first row ; 0 21-Apr-2002 00:40:00.000 00:40:20.000 0.20812 1216.00 1051.61 -246.67 1048.77 -259.32 9.37 18.96 ; 0.00029369 12.000 25.000 4F 5F 6F 7F 8F return, areas end