Nondepleted Time-Domain Method for Simulating Difference Frequency Generation From Metallic Metasurfaces