pro csmph06,phd,nmax,cmph,smph

    if n_params( ) lt 1 then begin
       print,'pro csmph06,phd,nmax,cmph,smph'
       return
    endif
    iph=n_elements(phd)
    cmph=fltarr(nmax+1,iph) & smph=cmph
    for i=0,iph-1 do begin
        phri=phd(i)*!DTOR
        csmphi06,nmax,phri,cmphi,smphi
        for m=0,nmax do begin
            cmph(m,i)=cmphi(m)
            smph(m,i)=smphi(m)
        endfor
    endfor

end