pro gfn,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)
        gfni,gg,hh,cmphi,smphi,Gfi,dGfi
        Gf(*,*,i) = Gfi
        dGf(*,*,i) = dGfi
    endfor
end