;
pro psb_hc,Nmax,apar,r,ith,jph,sth,cmph,smph,g,h,P,dP,bbr,bbt,bbp

    if N_params() LT 1 then begin
       print,'SYNTAX - pro sb_hc,Nmax,apar,r,ith,jph,sth,'
       print,'                   cmph,smph,g,h,P,dP,bbr,bbt,bbp'
       return
    endif

    Gfj = fltarr(Nmax+1,Nmax+1)
    dGfj = Gfj
    pi = Gfj
    dpi = pi
    iith = fix(ith + 1 + 0.5)
    jjph = fix(jph + 1 + 0.5)
    bbr = fltarr(jjph,iith)
    bbt = bbr
    bbp = bbr
    hcHf,Nmax,apar,r,Hr,dHr,d2Hr
    Gfn,Nmax,jph,cmph,smph,g,h,Gf,dGf
    FOR j = 0, jph do begin
        Gfj(*,*) = Gf(*,*,j)
        dGfj(*,*) = dGf(*,*,j)
        FOR i = 0, ith do begin
            sthi = sth(i)
            pi(*,*) = P(*,*,i)
            dpi(*,*) = dP(*,*,i)
            brtp,Nmax,apar,r,sthi,Hr,dHr,pi,dpi,Gfj,dGfj,Br,Bt,Bp
            bbr(j,i) = Br
            bbt(j,i) = Bt
            bbp(j,i) = Bp
        ENDFOR
    ENDFOR

end