I am modeling a 3D convecting system, and trying some model ingredients. I am specifically trying to incorporate phase changes, just one for now, at a certain depth. This is easily done with the latent heat model..specifying Clapeyron slope, density jump, transition depth and width, etc. I am finding trouble when trying to combine that with an ambient viscosity structure and a compositional field. I am using the Depth dependent model, providing a viscosity layering. As a base model I am using the latent heat model. The model does not run. The core of the error message is this: ''Exception 'ExcMessage("Error: At least one list that gives input parameters for the phase " "transitions has the wrong size. Currently checking against transition depths'' This seems easy to fix by writing and specifying the correct number of viscosity layers, and the correct number of phase transition interfaces.
In my model there is only one phase transition (for now) and several viscosity layers, and one compositional field.
I have been unable to find the right combination of specified interfaces of each kind. I have tried yes/no including the phase transition as a viscosity interface.
Do I need to, in addition of providing viscosity and depth lists, write the viscosity layers inside the latent heat model to ghost-disguise them as zero-Clapeyron-slope interfaces?
I would appreciate suggestions from someone who has tried this kind of combinations.
cheers many !