;
; PURPOSE: produce arrays for P and dP
; WRITTEN: Xuepu Zhao  28Feb1996
;
pro pdp,cth,Nmax,P,dP

    if N_params() LT 1 then begin
       print,'SYNTAX - pdp,cth,Nmax,P,dP'
       return
    endif
    jth=n_elements(cth)
    P = fltarr(Nmax+1,Nmax+1,jth)
    dP = fltarr(Nmax+1,Nmax+1,jth)
    for j=0,jth-1 do begin
        x = cth(j)
        pdpj,x,Nmax,Pj,dPj
        P(*,*,j) = Pj(*,*)
        dP(*,*,j) = dPj(*,*)
    endfor

end