pro gfn06,gg,hh,cmph,smph,Gf,dGf

    sz=SIZE(cmph) & nmax=sz(1)-1 & iph=sz(2)
    Gf = fltarr(nmax+1,nmax+1,iph)
    dGf = Gf
    for i=0,iph-1 do begin
        cmphi = cmph(*,i)
        smphi = smph(*,i)
        gfni06,nmax,cmphi,smphi,gg,hh,Gfi,dGfi
        Gf(*,*,i) = Gfi
        dGf(*,*,i) = dGfi
    endfor

end