[en] A computer code LHRAY has been developed for modelling lower hybrid wave propagation, absorption and current drive in tokamaks. LHRAY has been designed to run either as a stand-alone program with fixed plasma profiles or as an addition to an OLYMPUS style 1-D tokamak transport code. We assume an input rf power spectrum given by Brambilla theory. The spectrum is broken up into a number of discrete modes, which are followed as they propagate in the tokamak using the toroidal ray equations. The waves interact non-linearly with the electrons as they travel through the plasma, flattening the electron velocity space distribution and driving current. The electron distribution function is computed using a 1-D quasilinear model similar to that of Fisch. Relativistic effects and the presence of a background Ohmic electric field are also included. From the electron distribution function the quasilinear wave damping and electron current are calculated self-consistently. Simulation results give a fairly good prediction (to within approx. = 30%) of the steady state lower hybrid current driven in PLT experiments