A drift instability of a collisional magnetized plasma, unstable due to the Braginskii thermal force but not requiring any direct dissipation such as resistivity or electron inertia, is examined. Unlike conventional drift-modes, the maximum growth rate of the thermal force drift wave (TFDW) is of order the drift frequency, making for a strongly turbulent nonlinear state. A 3D, magnetized two-fluid code is developed to allow the study of both ideal MHD modes as well as lower frequency drift modes. The governing equations are essentially the ideal MHD equations with the inclusion of Hall and thermal force terms in Ohm’s law. This set of equations is reduced in a finite β, long parallel wavelength, and small but significant Larmor radius ordering and tested for shear Alfven waves, parallel sound waves, and drift modes. The code is employed to recover the TFDW instability, to verify the code against the mode’s analytic linear characteristics, and to study the nonlinear behavior of the TFDW. The TFDW growth is strongly suppressed by parallel thermal conduction and thus this mode is more likely to be observed in low temperature plasmas.