add the code that calculates the motion of a wave packet
in a 2D lattice while the orientation of the DC field
is changing adiabatically