pro fitsw,ofile,data,head sh=size(head) case sh(sh(0)+1) of 1: 8:head=fitsheader(head) else: head=fitsbasic(data) endcase openw,unit,ofile,/GET_LUN ; open file, get a lun for it. writeu,unit,head writeu,unit,data free_lun,unit en