@bruhwiler Two notes before creating the pull request:
Because of the windowing applied to the file-initialized fields, we now get fewer photons in file-init than in function-init (where previously they were the same)
We assume electric_field_y is zero in the file-init and other places. Should we change this?
Currently not getting good phase output of propagators + changes to mesh dimensions and extent.