# Solve discretized PDE with new initial conditions

**URL:** https://discourse.julialang.org/t/solve-discretized-pde-with-new-initial-conditions/125319
**Category:** Modelling & Simulations
**Tags:** package
**Created:** [January 28, 2025, 7:40pm UTC](https://discourse.julialang.org/t/solve-discretized-pde-with-new-initial-conditions/125319 "2025-01-28T19:40:45Z")
**Posts on this page:** 1
**Page:** 1

<div class="post-metadata">

### Author: ![oswaldo](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/oswaldo/32/213940_2.png) [@oswaldo](https://discourse.julialang.org/u/oswaldo)
#### Post date: [January 28, 2025, 7:40pm UTC](https://discourse.julialang.org/t/solve-discretized-pde-with-new-initial-conditions/125319/1 "2025-01-28T19:40:46Z")

</div>

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.