pro pdp,cth,Nmax,P,dP

;+
; NAME:
;	pdp
; INCLUDE:
	@compile_opt.pro		; On error, return to caller
; PURPOSE:
;	produce arrays for P and dP
;	Nmax le 63 OK for the case of ith=30    04Dec1999
;-

if N_params() LT 1 then message,'SYNTAX - pdp,cth,Nmax,P,dP'

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
	pdpi,cth[i],Nmax,Pi,dPi
	P [*,*,i] = Pi [*,*]
	dP[*,*,i] = dPi[*,*]
endfor

return  &  end