add a function that calculates the phase
change check for crossed phase transition
make width of phase transitions an input parameter
Disabling full_pressure and use adiabatic pressure for viscosity.
Introduce minimum and maximum viscosity.