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

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

    if not keyword_set(apar) then apar=0.2
    if not keyword_set(nmax) then nmax=9

    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