Документ взят из кэша поисковой машины. Адрес оригинального документа : http://star.arm.ac.uk/~csj/idl/UCLES/ucles_read.pro Дата изменения: Wed Aug 22 19:06:24 2007 Дата индексирования: Tue Oct 2 07:54:45 2012 Кодировка: Поисковые слова: р р р р р р рер р

PRO ucles_read, file, nw, no, w, f

;+
; Program:
; ucles_read

; Purpose:
; reads a 2d dataframe into IDL where the
; original data comes in a simple 2-column file.
; the user MUST tell the programme what the 2 dimensions are.

; Example:
; ucles_read,'aat_fram_0001.dat',1024,46,w,f

; INPUT

; file : string variable containing file name
; file contains 2 columns of wl and fl
; and EXACTLY nw * no lines.

; nw : number of wavelengths
; no : number of orders

; OUTPUT

; w ; 2d real array (nw,no) containing wavelengths
; f : 2d real array (nw,no) containing fluxes

;-

PRINT,'UCLES_READ: Reading echelle spectrum:',file

w = fltarr(nw,no)
f = fltarr(nw,no)
title = ' dummy '

; Open the data file
on_ioerror,endofdata
openr, lun, file, /get_lun

; read the header
ReadF, lun, title
print, title
ReadF, lun, title
ReadF, lun, nw_tot
print, nw_tot, ' ? ', nw,no,'=',nw*no

FOR j=0,no-1 DO BEGIN
FOR i=0,nw-1 DO BEGIN

ReadF, lun, ww, ff
w(i,j) = ww
f(i,j) = ff

ENDFOR
ENDFOR

close,lun
free_lun, lun

RETURN

endofdata:
close,lun
free_lun, lun
print,'bad input file ',file

END