;
; PURPOSE:	Calculate Br,Bt,Bp at any r between r=1 and r=2.5 using
;               the PFSS model with Nmax
; INPUT: 	br1f - the data file for global distribution of the radial
;               photospheric field; Nmax - The maximum order number of
;               the spherical harmonic coefficient; r
; OUTPUT:	br1a - the array for field distribution at r=1 solar radia;
;               brra, btra, bpra - the arraies for surface distribution of three
;                                  components at r=r solar radii
; WRITTEN:	Jan. 3, 2005	Xuepu Zhao
;
pro cbr_brf,brf,Nmax,r,br1a,brra,btra,bpra,FTS=fts

    if keyword_set(fts) then br1a=readfits(brf,hd) $
                        else br1a=readwsosc(brf,/cali,/br)
    sz=SIZE(br1a)
    xsz=sz(1) & ysz=sz(2)
    zgrid,xsz,ysz,phd,thd,lad,cth,sth
    pdp,cth,Nmax,P,dP
    csmph,phd,Nmax,cmph,smph
    ghf_bra,br1a,Nmax,ga,ha,/cmpt,P,cmph,smph 
    zsb_ss,ga,ha,r,sth,P,dP,cmph,smph,brra,btra,bpra

end