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 change problem when using phase diagram #12

Closed NicolasRiel closed 8 months ago

NicolasRiel commented 8 months ago

Hello,

When trying to use phase change with phase diagram ("rho_ph" option in the Material structure) we get the following error message. Without the phase diagram option "rho_ph" everything works fine.


============================== INITIAL GUESS =============================


** NONLINEAR SOLVER FAILED TO CONVERGE! **

SNES Divergence Reason : residual norm is NAN [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: Severe divergence reason detected (see above) [0]PETSC ERROR: WARNING! There are option(s) set that were not used! Could be the program crashed before they were used or a spelling mistake, etc!

For instance: (I cannot copy the code, it put everything as a block)

![Uploading Untitled-1.png…]()

Nico