Uncoupled temperature and compositional field in UW Geodynamics

[bush2014]

I am currently mentoring a high school student using UW Geodynamics, and I must say it's been great for educational purposes; we really appreciate the effort put into it.

In this project, where we've imposed a velocity (1 cm/year) from one side of the wall. However, we've noticed that the temperature of the slabs and the composition don't seem to match well when the slabs are being subducted. In my experience with UW, we usually address this issue with the following comment: dt = min(advDiff.get_max_dt(), advector.get_max_dt()). Is there a similar approach or workaround in UW Geodynamics to address this temperature and composition mismatch? We are using the version underworld2_2.13.1b, which combines UW and UW Geodynamics.

Any information or guidance on this matter would be greatly appreciated. Thanks!