;
; MODIFIED:  01APR2004	Xuepu for adding /dbl
;
pro zcsmphi,nmax,phri,cmphi,smphi,DBL=dbl

    if keyword_set(dbl) then cmphi=dblarr(nmax+1) $
                        else cmphi=fltarr(nmax+1)
    smphi=cmphi
    for m=0,nmax do begin
        mph=m*phri
        cmphi(m)=COS(mph)
        smphi(m)=SIN(mph)
    endfor

end