; WRITTEN: 	28Jun2006, Xuepu
pro opnl_brf06,brf,px,py,SINL=sinl,CLR=clr,BRA=bra

    if N_params() lt 1 then begin
       print,'pnl_brf06,brf,px,py,/sinl,/clr,/bra'
       return
    endif
;
    if keyword_set(bra) then bbr=brf else bbr=readfits(brf,hd)
    sz=SIZE(bbr)
    zgrid06,sz(1),sz(2),phd,thd,lad
    lv=0.025
    if keyword_set(sinl) then begin
       lad=SIN(lad*!DTOR)
       yr=[-1,1]
    endif else begin
       yr=[-90,90]
    endelse
    xy_pt,phd,lad,xx,yy
    if keyword_set(clr) then begin
       tek_color
       contour,bbr,xx,yy,position=[px(0),py(0),px(1),py(1)],xrange=[0,360],$
               yrange=yr,xstyle=5,ystyle=5,c_colors=[3], $
               level=[0],thick=3,/device,/noerase
    endif else begin
       contour,bbr,xx,yy,position=[px(0),py(0),px(1),py(1)],xrange=[0,360],$
               yrange=yr,xstyle=5,ystyle=5,$
               level=[0],thick=3,/device,/noerase
    endelse

end