Removed passing of call_cnt for debugging, extra debugging output, added option of using 2D strain rate invariants and fixed Mohr-Coulomb plasticity formatting to be similar to the other options