;
; WRITTEN:  30JAN2006   Xuepu
;
pro opfp_iflf06,iflf,fpsz,SINL=sinl,CLR=clr,OPN=opn

    phfp=fltarr(2)
    thfp=phfp
    get_lun,unit
    openr,unit,iflf
    WHILE ~ EOF(unit) DO BEGIN
      readf,unit,np,fte
      anp=ABS(np)
    ; if anp gt 1 and anp lt 3 then print,'np,fte: ',np,fte
      if anp gt 5 then begin
         for n=0,anp-1 do begin
             readf,unit,x0,y0,z0,r0,t0,p0
             if n eq anp-1 then begin
             ;  print,'t0,p0: ',t0,p0
                phfp(0)=p0
                thfp(0)=90-t0
                if keyword_set(sinl) then thfp(0)=SIN(thfp(0)*!DTOR)
                if keyword_set(norm) then begin
                   thfp(0)=SIN(thfp(0)*!DTOR)
                   thfp(0)=(thfp(0)+1)/2.0
                   phfp(0)=phfp(0)/360.0
                endif
             endif
         endfor
         if keyword_set(clr) then begin
            tek_color
            if keyword_set(opn) then mkpsym8,fpsz,/open else mkpsym8,fpsz
            if fte gt 0.0 then begin
               oplot,phfp(0:0),thfp(0:0),psym=8,color=4
            endif else begin
               oplot,phfp(0:0),thfp(0:0),psym=8,color=2
            endelse
         endif else begin
            if fte ge 0.0 then begin
               mkpsym8,fpsz
               oplot,phfp(0:0),thfp(0:0),psym=8
            endif else begin
               mkpsym8,fpsz,/open
               oplot,phfp(0:0),thfp(0:0),psym=8
            endelse
         endelse
      endif
    ENDWHILE
    free_lun,unit

end