A wind flow sensing system for determining velocity fields of a wind flow at a set of different altitudes from a set of measurements of radial velocities at each of the altitudes, comprising: an input interface to accept the set of measurements of radial velocities at line-of-site points for each of the altitudes; a processor configured to determine a first approximation of the velocity field at each of the altitudes by simulating computational fluid dynamics, CFD, of the wind flow using operating parameters that minimize a cost function of a weighted combination of errors, wherein each error corresponds to one of the altitudes and includes a difference between measured velocities at the measurement points at the corresponding altitude and simulated line-of-sight velocities at the measurement points simulated by the CFD for the corresponding altitude with current values of the operating parameters, wherein the errors include a first error corresponding to a first altitude and a second error corresponding to a second altitude, wherein a weight for the first error in the weighted combination of errors is different from a weight of the second error in the weighted combination of errors; determine a horizontal derivative of vertical velocity at each of the altitudes from the first approximation of the velocity field at each of the altitudes; and determine a second approximation of the velocity field at each of the altitudes using geometric relationships between a velocity field for each of the altitudes, projections of the measurements of radial velocities on the three-dimensional axes, and the horizontal derivative of vertical velocity for the corresponding velocity field; and an output interface to render the second approximation of velocity fields of the wind flow. an input interface to accept the set of measurements of radial velocities at line-of-site points for each of the altitudes; a processor configured to determine a first approximation of the velocity field at each of the altitudes by simulating computational fluid dynamics, CFD, of the wind flow using operating parameters that minimize a cost function of a weighted combination of errors, wherein each error corresponds to one of the altitudes and includes a difference between measured velocities at the measurement points at the corresponding altitude and simulated line-of-sight velocities at the measurement points simulated by the CFD for the corresponding altitude with current values of the operating parameters, wherein the errors include a first error corresponding to a first altitude and a second error corresponding to a second altitude, wherein a weight for the first error in the weighted combination of errors is different from a weight of the second error in the weighted combination of errors; determine a horizontal derivative of vertical velocity at each of the altitudes from the first approximation of the velocity field at each of the altitudes; and determine a second approximation of the velocity field at each of the altitudes using geometric relationships between a velocity field for each of the altitudes, projections of the measurements of radial velocities on the three-dimensional axes, and the horizontal derivative of vertical velocity for the corresponding velocity field; and determine a horizontal derivative of vertical velocity at each of the altitudes from the first approximation of the velocity field at each of the altitudes; and determine a second approximation of the velocity field at each of the altitudes using geometric relationships between a velocity field for each of the altitudes, projections of the measurements of radial velocities on the three-dimensional axes, and the horizontal derivative of vertical velocity for the corresponding velocity field; and an output interface to render the second approximation of velocity fields of the wind flow. <EOS>