added a solver for equilibrium crusts

still needs to be automated
moved deallocate after write statement