We would like to thank the referee for their kind feedback and helpful comments. We will correct the typo errors pointed out.

Regarding the comment on line 185: 

To compute the electric field amplitudes at height z, we use a simplified version of equation 16 from Shoucri and Morales 1984, given by: 

$|E(z)|= 0.25*(beta*W*G)^{1/2}/D$

where beta is the absorption at low altitude, assuming no loss in the D region so that beta = 1, W is the average radiated power, G is the gain of the antenna, and D is the altitude from the ground. Here (W*G) is the ERP, the effective radiated power. In principle, the actual electric field will depend on the electron density at altitudes below, which controls how much of the Poynting flux is absorbed. For this paper, we rather choose some representative, but somewhat realistic, electric field amplitude to be able to compare the shape of the phase space density for the Monte Carlo simulation and the solution based on Boltzmann equation, and in order to compare the response at different altitudes, we do not take into account the absorption below. 

References:

Shoucri, M. M. and Morales, G. J.: Ohmic Heating of the Polar F Region by HF Pulses, Journal of geophysical research, 89, 2907–2917, 1984.