Future satellite-to-ground optical communication systems will benefit from accurate forecasts of atmospheric optical turbulence; namely for site selection, for the routing and the operation of optical links, and for the design of optical communication terminals. This work presents a numerical approach based on the Weather Research and Forecasting software that enables continuous forecast of the refractive index structure parameter, C n2, vertical profiles. Two different C n2 models are presented and compared. One is based on monitoring the turbulent kinetic energy, while the other is a hybrid model using the Tatarskii equation to depict the free atmosphere region, and the Monin-Obukhov similarity theory for describing the boundary layer. The validity of both models is assessed by using thermosonde measurements from the Terrain-induced Rotor Experiment campaign, and from day and night measurements of the coherence length collected during a six-day campaign at Paranal observatory by a Shack-Hartmann Image Motion Monitor. The novelty of this work is the ability of the presented approach to continuously predict optical turbulence both during daytime and nighttime, and its validation with measurements in day and night conditions.