Fix quadratic summary stats with delta models

[ericward-noaa]

This is the fix for the quadratics roots code -- using this to finish the spatial indicators paper / ECCWO project. Not sure whether this is better left on its own branch, merged in, or moved to the experimental. Test code exists but currently commented out