We present results from a continuing effort to simulate equatorial spread F (ESF) using observations from the Jicamarca Radio Observatory near Lima, Peru. Jicamarca measures vertical and zonal plasma drifts along with plasma number density profiles overhead. The number density profiles are used to initialize a three-dimensional regional model of the ionosphere capable of simulating plasma density irregularities produced during ESF conditions. The vertical drifts measurements are used to drive the numerical simulation continuously. Neutral winds are derived from the new Horizontal Wind Model '14 (HWM-14) model, and the zonal winds are scaled so as to make the zonal plasma flows at the start of the simulation agree with the ISR profile measurements. Coherent scatter radar imagery from Jicamarca is used to validate the simulation results. Campaign data were collected in April and December 2014, and a few events representative of low and high ESF activity were selected for analysis. The numerical simulations are able to reproduce the level of activity observed along with the gross features of the ESF irregularities and radar plumes. Data from a network of HF beacons are being incorporated into the forecast analysis in order to elucidate radar plumes which sometimes appear even when the simulation fails to predict them.