Computational recovery of the homoclinic orbit in porous media convection