Stream function for optical energy flow

A formula is presented for the stream function S corresponding to optical fields in the plane that can be represented as a superposition of plane waves with the same transverse wavenumbers, including evanescent waves. Energy flow lines (parallel to the Poynting vector) are easily computed as contour lines of S. Optical vortices (including higher-order vortices) correspond to extrema of S.