Changed the E spin orbit constant from the product of BzzE and zeta to just zeta. Also, changed DHamilt back to the default SV from the special version I tried. This seems to fix problems with fits and the model no longer crashes when putting nonzero SO coupling.