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

    if N_params() lt 1 then begin
       print,'pro pb_ss06,g,h,r,cthj,sthj,phri,Br,Bt,Bp,Rss=,Apar=,Nmax='
       return
    endif

    if not keyword_set(rss) then rss=2.5
    if not keyword_set(apar) then apar=0.0
    if not keyword_set(nmax) then nmax=9

    sshf06,apar,nmax,rss,r,Hr,dHr,d2Hr
    csmphi06,Nmax,phri,cmphi,smphi
    gfni06,nmax,cmphi,smphi,g,h,Gfi,dGfi
    pdpj06,cthj,Nmax,Pj,dPj
    brtp06,r,sthj,Hr,dHr,Pj,dPj,Gfi,dGfi,Br,Bt,Bp

end