;
; PURPOSE:    calculate vector field at a point (r,thri(ASIN(sthi)),phrj)
; MODIFIED:   14JUN2000 Xuepu Zhao
;
pro pb_hc,g,h,r,cthi,sthi,phrj,Br,Bt,Bp,Apar=apar

    if N_params() lt 1 then begin
       print,'pb_hc,g,h,r,cthi,sthi,phrj,Br,Bt,Bp,apar='
       return
    endif

    sz=SIZE(g)
    nmax=sz(1)-1
    if not keyword_set(apar) then apar=0.2

    hchf,Nmax,apar,r,Hr,dHr,d2Hr
    csmphj,nmax,phrj,cmphj,smphj
    gfnj,nmax,cmphj,smphj,g,h,Gfj,dGfj
    pdpi,cthi,nmax,Pi,dPi
    brtp,Nmax,apar,r,sthi,Hr,dHr,Pi,dPi,Gfj,dGfj,Br,Bt,Bp

end