;
; WRITTEN: 25APR2003	Xuepu Zhao 
; MODIFIED: 01APR2004   Xuepu for adding /dbl
;
pro zgfn,cmph,smph,gg,hh,Gf,dGf,NMAX=nmax,DBL=dbl

    sz=size(cmph)
    iph=sz(2)
    if not keyword_set(nmax) then begin
       nmax=sz(1)-1
    endif
    if keyword_set(dbl) then begin
       Gf = dblarr(nmax+1,nmax+1,iph)
       cmphj = dblarr(nmax+1)
    endif else begin
       Gf = fltarr(nmax+1,nmax+1,iph)
       cmphj = fltarr(nmax+1)
    endelse
    dGf = Gf
    smphj = cmphj
    for j=0,iph-1 do begin
        cmphj(*) = cmph(0:nmax,j)
        smphj(*) = smph(0:nmax,j)
        if keyword_set(dbl) then zgfni,cmphj,smphj,gg,hh,Gfj,dGfj,/dbl $
                            else zgfni,cmphj,smphj,gg,hh,Gfj,dGfj
        Gf(*,*,j) = Gfj(0:nmax,0:nmax)
        dGf(*,*,j) = dGfj(0:nmax,0:nmax)
    endfor
end