;
; PURPOSE:      Overplot foot-points of IFLs at rr
; NOTE:         bottom=1.001, dstep=0.02
; WRITTEN:      16Aug2000    refer to ~/ZIDLLIB/JOB/opfp_brfiflf.pro
;
;pro opfl_oflf06,oflf,r1,r2,flsz,SINL=sinl,CLR=clr,OPN=opn
pro opfl_oflf06,oflf,r1,r2,SINL=sinl,CLR=clr,OPN=opn

    if n_params( ) lt 1 then begin
       print,'pro opfl_oflf06,oflf,r1,r2,/sinl,/clr,/opn'
       return
    endif
    if keyword_set(clr) then tek_color
    get_lun,unit
    openr,unit,oflf
    while eof(unit) ne 1 do begin
       readf,unit,np
       anp=abs(np)
       phfl=fltarr(anp)
       lafl=phfl
       rrfl=phfl
       for n=0,anp-1 do begin
           readf,unit,x0,y0,z0,r0,t0,p0
           phfl(n)=p0
           lafl(n)=90-t0
           rrfl(n)=r0
       endfor
;      mkpsym8,0.08
;      mkpsym8,flsz
       if keyword_set(sinl) then lafl=SIN(lafl*!DTOR)
       if anp ge 3 then begin
          rmax=max(rrfl)
          ind=where(rrfl eq rmax)
          in_rmax=ind(0)
          if in_rmax lt (anp-1) then begin   ; closed field lines
             if rmax gt r1 and rmax lt r2 then begin
                if in_rmax gt 0 AND in_rmax lt (anp-1) then begin
                   if np gt 0 then begin
                      phflp=phfl(0:(in_rmax)) & laflp=lafl(0:(in_rmax))
                      phfln=phfl(in_rmax:(anp-1)) & lafln=lafl(in_rmax:(anp-1))
                      oplot,phflp,laflp,color=4
                      oplot,phfln,lafln,color=2
           ; for k=0,in_rmax-1 do oplot,phfl(k:k),lafl(k:k),psym=8,color=4
           ; for k=in_rmax,anp-1 do oplot,phfl(k:k),lafl(k:k),psym=8,color=2
                   endif else begin
                      phfln=phfl(0:(in_rmax)) & lafln=lafl(0:(in_rmax))
                      phflp=phfl(in_rmax:(anp-1)) & laflp=lafl(in_rmax:(anp-1))
                      oplot,phfln,lafln,color=2
                      oplot,phflp,laflp,color=4
           ; for k=0,in_rmax-1 do oplot,phfl(k:k),lafl(k:k),psym=8,color=2
           ; for k=in_rmax,anp-1 do oplot,phfl(k:k),lafl(k:k),psym=8,color=4
                   endelse
                 endif
             endif
          endif else begin
             if keyword_set(opn) then oplot,phfl,lafl,color=3
          ;  if keyword_set(opn) then $
          ;     for k=0,anp-1 do oplot,phfl(k:k),lafl(k:k),psym=8,color=3
          ;  if np gt 0 then $
          ;     for k=0,anp-1 do oplot,phfl(k:k),lafl(k:k),psym=8,color=4 else $
          ;     for k=0,anp-1 do oplot,phfl(k:k),lafl(k:k),psym=8,color=2
          endelse 
       endif
    endwhile
    free_lun,unit

end