  We construct a tensor network representation of the partition function for
the massless Schwinger model on a two dimensional lattice using staggered
fermions. The tensor network representation allows us to include a topological
term. Using a particular implementation of the tensor renormalization group
(HOTRG) we calculate the phase diagram of the theory. For a range of values of
the coupling to the topological term $\theta$ and the gauge coupling $\beta$ we
compare with results from hybrid Monte Carlo when possible and find good
agreement.
