PhaseTransitions were inactive because a keyword was not set.
On the occasion, a few additional missing keywords were added & I made "rk2" the default advection scheme (which was reported to work better by @NicolasRiel).
Fixes issue #26
NOTE: there was also an issue in LaMEM which made phase transitions to be inactive when Phase==0 was selected. See here. This PR fixes that, which should be integrated in version 2.1.4. As a workaround, avoid using phase=0 for the mantle etc.
PhaseTransitions were inactive because a keyword was not set. On the occasion, a few additional missing keywords were added & I made "rk2" the default advection scheme (which was reported to work better by @NicolasRiel).
Fixes issue #26
NOTE: there was also an issue in LaMEM which made phase transitions to be inactive when
Phase==0was selected. See here. This PR fixes that, which should be integrated in version 2.1.4. As a workaround, avoid usingphase=0for the mantle etc.