Hi I’m discretizing a PDE system using MethodOfLines.jl with mixed boundary conditions. It works well. I would like to solve the discretized PDE (ODE) system with new initial conditions (at each discretized point). I do newprob = remake(prob). I would like to do something like solve(newprob,u0=newu0). However, the newprob.u0 stores the initial conditions in an unexpected order, so I do not how to pass the new initial conditions when solving the newprob. Note: the discretized PDE is actually a DAE system.