;
; PURPOSE: produce arrays for P and dP
; Nmax le 63 OK for the case of ith=30    04Dec1999
;
pro pdp,cth,Nmax,P,dP

    if N_params() LT 1 then begin
       print,'SYNTAX - pdp,cth,Nmax,P,dP'
       return
    endif
    ith=n_elements(cth)
    P = fltarr(Nmax+1,Nmax+1,ith)
    dP = fltarr(Nmax+1,Nmax+1,ith)
    for i=0,ith-1 do begin
        x = cth(i)
        pdpi,x,Nmax,Pi,dPi
        P(*,*,i) = Pi(*,*)
        dP(*,*,i) = dPi(*,*)
    endfor

end