A nonequilibrium molecular dynamics method for calculating the thermal conductivity is presented. It reverses the usual cause and effect picture. The "effect," the heat flux, is imposed on the system and the "cause," the temperature gradient is obtained from the simulation. Besides being very simple to implement, the scheme offers several advantages such as compatibility with periodic boundary conditions, conservation of total energy and total linear momentum, and the sampling of a rapidly converging quantity (temperature gradient) rather than a slowly converging one (heat flux). The scheme is tested on the Lennard-Jones fluid.