;
; calculate Br,Bt,Bp,Dens,Ptot,curtt,curtp at r,thi,phj using HCCSSS model
; r > Rcp
;
pro Pbp_csss,Rcp,Rss,c1,d1,r,cthi,sthi,phrj,Br,Bt,Bp,Dens,Ptot,curtt,curtp

    common share1, apar,Nmax,g,h,gc,hc

    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
    bprtp,Nmax,apar,c1,d1,r,sthi,Hr,dHr,d2Hr,Pi,dPi,Gfj,dGfj,$
          Br,Bt,Bp,Dens,Ptot,curtt,curtp

end