EnsembleSummary support for DiffEqArray

[ChrisRackauckas]

It would be nice to have EnsembleSummary support DiffEqArray outputs, i.e.:

function save_mortality(sol, ind)
    # Save x[5]=flow and N*x[8]=total from solution
    flow_and_total = [[x[5],p.N .* x[8]] for x in sol]
    return (DiffEqArray(flow_and_total,sol.t), false)
end
ensembleprob = EnsembleProblem(prob,  output_func = save_mortality)
sol = solve(ensembleprob, SOSRI(), EnsembleThreads(), saveat = 0.5, trajectories = 1000)

The issue is that the current codepath assumes an array output instead of an array of arrays, so an automated conversion would be the quickest way of getting this supported.