; 
; PURPOSE:  calculate vector field at a point (r,thrj(ASIN(sthj)),phri)
; MODIFIED:   14JUN2000	Xuepu Zhao
;
pro zpb_ss,g,h,r,cthj,sthj,phri,Br,Bt,Bp,RSS=rss,NMAX=nmax,N0=n0

    if N_params() lt 1 then begin
       print,'zpb_ss,g,h,r,cthj,sthj,phri,Br,Bt,Bp,rss=,nmax=,/n0'
       return
    endif

    if not keyword_set(rss) then rss=2.5
    if not keyword_set(nmax) then begin
       sz=size(g)
       nmax=sz(1)-1
    endif

    zsshf,nmax,rss,r,Hr,dHr,d2Hr
    zcsmphi,nmax,phri,cmphi,smphi
    zpdpj,cthj,nmax,Pj,dPj
    zgfni,cmphi,smphi,g,h,Gfi,dGfi
    if keyword_set(n0) then $
       zbrtp,r,sthj,Hr,dHr,Pj,dPj,Gfi,dGfi,Br,Bt,Bp,/n0 else $
       zbrtp,r,sthj,Hr,dHr,Pj,dPj,Gfi,dGfi,Br,Bt,Bp

end