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