UniMainzGeo / LaMEM

LaMEM (Lithosphere and Mantle Evolution Model) - a 3D parallel code to simulate geodynamic and geomechanical processes
MIT License
50 stars 16 forks source link

Phase Transitions don't work when one of the Phases==0 #14

Closed boriskaus closed 8 months ago

boriskaus commented 8 months ago

Phase Transitions don't work if the involved phases have phase==0. This can be reproduced with t16_PhaseTransitions/Plume_PhaseTransitions.dat by making the air phase 3 instead of 0, and the mantle phase 0 instead of 3.

This results in:

Screenshot 2024-02-04 at 18 48 51

whereas it should look like:

Screenshot 2024-02-04 at 18 50 23

A test file is here: [Plume_PhaseTransitions_SwappedPhases.dat.zip] (https://github.com/UniMainzGeo/LaMEM/files/14157744/Plume_PhaseTransitions_SwappedPhases.dat.zip)

boriskaus commented 8 months ago

this has been fixed upstream in LaMEM and should become available in the next release of LaMEM