;
; WRITTEN: 25 APR2003	Xuepu Zhao
; MODIFIED: 01APR2004   Xuepu for adding /dbl
;
pro zgfni,cmphi,smphi,gg,hh,Gfi,dGfi,DBL=dbl

    nmax=n_elements(cmphi)-1
    if keyword_set(dbl) then Gfi = dblarr(nmax+1,nmax+1) $
                        else Gfi = fltarr(nmax+1,nmax+1)
    dGfi = Gfi
    for n=0, nmax do begin
        for m=0, n do begin
            Gfi(n,m) = gg(n,m)*cmphi(m)+hh(n,m)*smphi(m)
            dGfi(n,m) = m*(-1.*gg(n,m)*smphi(m)+hh(n,m)*cmphi(m))
        endfor
    endfor

end