Finally debugged nonconservation of energy and momentum. Missing factor of 4pi in momneutm update. adiabatic-1 in TE update