pro csmph,phd,nmax,cmph,smph

    iph=n_elements(phd)
    cmph=fltarr(nmax+1,iph) & smph=cmph
    for i=0,iph-1 do begin
        phri=phd(i)*!DTOR
        csmphi,nmax,phri,cmphi,smphi
        for m=0,nmax do begin
            cmph(m,i)=cmphi(m)
            smph(m,i)=smphi(m)
        endfor
    endfor

end