The numerical modeling presented in this work shows the importance of taking into account an effect of the LH power on the SOL heating, and hence on the ionization, when modeling the SOL plasma. The modeled density growth due to SOL heating and gas puff is consistent with the modeled SOL ionization source profiles, which for puff and heating are strongly enhanced and extend into the far SOL, contrary to the case without heating and/or without the gas puff. This could explain the observed improvement of the LH wave coupling at the near grill gas puffing.