def surf_to_phi(surf, phi, x):
    """
    move the surf density back to the full phi density function
    """
    for ii in range(len(x)):
        for jj in range(len(x)):
            kk = len(x) - 1 - ii - jj
            phi[ii, jj, kk] = surf[ii, jj]
    return phi
