Setting initial H concentration

User Support
2

Hi @HelenS!

You can set the initial concentration in a specific domain with sympy.Piecewise(...) function:

my_model.initial_conditions = [
    F.InitialCondition(field=0, value=10*sympy.Piecewise((1, F.x <= L), (0, True))),
]

where L is the length of the first domain, for example. See this topic for additional details.

Please, let me know if it is what you’ve asked about.