pro cbriflofl_bra,bra,Nmax,r,ghf=ghf,scrn=scrn,brf=brf,iflfpf=iflfpf,oflf=oflf

    writefits,'br1_2000.fits',bra
    sz=SIZE(bra)
    xsz=sz(1) & ysz=sz(2)
    zgrid,xsz,ysz,phd,thd,lad,cth,sth
    pdp,cth,Nmax,P,dP
    csmph,phd,Nmax,cmph,smph
    if keyword_set(ghf) then $ 
       ghf_bra,bra,Nmax,ga,ha,/cmpt,P,cmph,smph,ghf=ghf,scrn=scrn $
                        else ghf_bra,bra,Nmax,ga,ha,/cmpt,P,cmph,smph

    if keyword_set(brf) then begin
       zsb_ss,ga,ha,r,sth,P,dP,cmph,smph,bbr,bbt,bbp
       writefits,brf,bbr
    endif 
    if keyword_set(iflfpf) then begin
       get_lun,unit
       openw,unit,iflfpf
       printf,unit,'xsz,ysz: ',xsz,ysz
       for i=0,xsz-1 do begin
           phri=phd(i)*!DTOR
           for j=0,ysz-1 do begin
               thrj=thd(j)*!DTOR
               ziflij_ss,ga,ha,phri,thrj,np,xv,yv,zv,rv,tv,pv,fte,nmax=Nmax
               printf,unit,'phi,thj,Np,FTE: ',phd(i),thd(j),np,fte
               if np gt 5 then $
                  for k=0,np-1 do printf,unit,xv(k),yv(k),zv(k),rv(k),tv(k),pv(K)
           endfor
       endfor
       free_lun,unit
    endif
    if keyword_set(oflf) then begin
       get_lun,unit
       openw,unit,oflf
       for j=0,ysz-1 do begin
           thrj=thd(j)*!DTOR
           for i=0,xsz-1 do begin
               phri=phd(i)*!DTOR
               if keyword_set(n0) then begin
                  zoflij_ss,ga,ha,phri,thrj,np,x,y,z,rr,th,ph,/n0
               endif else begin
                  zoflij_ss,ga,ha,phri,thrj,np,x,y,z,rr,th,ph
               endelse
               printf,unit,np
               anp=abs(np)
               for k=0,anp-1 do printf,unit,x(k),y(k),z(k),rr(k),th(k),ph(k)
           endfor
       endfor
       free_lun,unit
    endif

end