need to get rotation to obey virial theorem at low g, Maybe an error in kinetic term?