;
; Determin the three components (without sign) of B between Rcp and Rss
; NOTE: Br, .. in mT
;
pro pb0_csss,gc,hc,apar,Rcp,Rss,r,cthi,sthi,phrj,Br,Bt,Bp

    sz=size(gc)
    Nmax=sz(1)-1
    cssshf,Nmax,apar,Rcp,Rss,r,Hr,dHr,d2Hr
    csmphj,Nmax,phrj,cmphj,smphj
    gfnj,Nmax,cmphj,smphj,gc,hc,Gfj,dGfj
    pdpi,cthi,Nmax,Pi,dPi
    brtp_csss,Nmax,r,sthi,Hr,dHr,Pi,dPi,Gfj,dGfj,Br,Bt,Bp

end