pro norh_rd_img,infil,index,data

;+ 
; NAME:
;       NORH_RD_IMG
;
; PURPOSE:
;       This procedure reads the NORH image file.
;
; CALLING SEQUENCE:
;       norh_rd_img,infil,index,data
;
; INPUTS:
;       infil: list of image files to be read
; OUTPUTS:
;       index: index structure
;       data: Image data
; HISTORY:
;       Writtern 1998-05-29 JS and TY
;       Change name 1998-01-03 TY from norh_rd_maps
;-
  
mx=n_elements(infil)

if n_params() eq 3 then begin
if (mx eq 1) then begin
  file=infil(0)
  data_el=readfits(file,header,/silent)
  if (n_elements(data_el) eq 1) then begin
    if (data_el eq -1) then begin
      print,'No data file exist'
      return
    endif
  endif
  index_el=norh_img2idx(header)

  data=data_el
  index=index_el
endif else begin
  m=0
    file=infil(m)
    data_el=readfits(file,header,/silent)
    index_el=norh_img2idx(header)

    sz=size(data_el)
    data=replicate(data_el(0,0),sz(1),sz(2),mx)

    data(*,*,m)=data_el
    index=index_el
  for m=1,mx-1 do begin
    file=infil(m)
    data_el=readfits(file,header,/silent)
    index_el=norh_img2idx(header)

    data(*,*,m)=data_el
    index=[index,index_el]
  endfor
endelse

endif else begin  ; n_params() = 2
if (mx eq 1) then begin
  file=infil(0)
  header=headfits(file)
  index_el=norh_img2idx(header)

  index=index_el
endif else begin
  m=0
    file=infil(m)
    header=headfits(file)
    index_el=norh_img2idx(header)

    index=index_el
  for m=1,mx-1 do begin
    file=infil(m)
    header=headfits(file)
    index_el=norh_img2idx(header)

    index=[index,index_el]
  endfor
endelse
endelse

end