This document examines the impact of source terms in the variational representation of traffic flow. It finds that variational theory solutions are valid in Eulerian coordinates for exogenous source terms but not for endogenous source terms that depend on traffic conditions. However, in discrete time steps, endogenous source terms from the previous step become exogenous for the current step, allowing improved numerical solution methods. Specifically, it presents Godunov's method and variational networks methods for solving the kinematic wave and Hamilton-Jacobi models of traffic flow in discrete time steps.