Simulation of many-body localization and time crystals in two dimensions with the neighborhood tensor update