;-------------------------------------------------------------------------

pro fits::read,file

 if n_elements(file) ne 1 then return ;-- enter at least one file name
 check=findfile(file,count=count)     ;-- check that file exists
 if count eq 0 then return
 data=mrdfits(file)                   ;-- call Astronomy library FITS reader
 self->set,data                       ;-- insert data into itself
 self.filename=file                   ;-- save filename in property 
 return

end

;-------------------------------------------------------------------------

pro fits__define
 
 void={fits,filename:'', inherits data}
 return 

end