Several points in the code require conversion from specific to conserved quantities. This generalizes the selection of fields to a single function call.
This changes the answer on the Photon fields, as there was a missing flag to prevent RadiationFields from being multiplied by density in ProjectSolutionToParentGrid.
The results for the PhotonTest still agrees with the analytic solution. I'm accepting. This makes things a lot nicer, btw!