; NOTE: jph here is 71, not 72
pro gfn,Nmax,jph,cmph,smph,gg,hh,Gf,dGf

    Gf = fltarr(Nmax+1,Nmax+1,jph+1)
    dGf = fltarr(Nmax+1,Nmax+1,jph+1)
    cmphj = fltarr(Nmax+1)
    smphj = cmphj
    for j=0,jph do begin
        cmphj(*) = cmph(0:Nmax,j)
        smphj(*) = smph(0:Nmax,j)
        Gfnj,Nmax,cmphj,smphj,gg,hh,Gfj,dGfj
        Gf(*,*,j) = Gfj(0:Nmax,0:Nmax)
        dGf(*,*,j) = dGfj(0:Nmax,0:Nmax)
    endfor
end