Porosity change submodel not running for Chen2020_composite parameter set

[SagarPal8894]

Hi all,
I am trying to simulate a _coupled degradation DFN model for a cell having composite (graphite + silicon) anode composition. For this, I am using the "Chen 2020 composite" as a base parameter set. The degradation mechanisms active are SEI growth (solvent-diffusion limited) and partially reversible lithium plating. Whenever I switch ON the porosity change due to SEI growth and lithium plating, it gives me the following error:

Important

ModelError: Missing variable for submodel 'porosity': 'Negative total SEI thickness [m]'.
Check the selected submodels provide all of the required variables.

Few code line snippets are provided for your kind reference:

param=pybamm.ParameterValues("Chen2020_composite")

model=pybamm.lithium_ion.DFN(options={
    "particle phases": ("2", "1"),
        "open-circuit potential": (("single", "current sigmoid"), "single"),
        # "thermal":"lumped",
        "cell geometry": "arbitrary",
        "SEI": "solvent-diffusion limited",
        "SEI porosity change": "true",
        "lithium plating": "partially reversible",
        "lithium plating porosity change": "true",
        # "particle mechanics": ("swelling and cracking", "swelling only"),
        # "loss of active material": "stress-driven",
        # "SEI on cracks": "true",
        "calculate discharge energy": "true"
        })

Further clarification is eagerly requested from the experience modellers who have performed such degradation simulations on composite electrode models.