We have modeled the individual harmonic frequencies of the Ionospheric Alfvén Resonator (IAR) at Eskdalemuir by solving a one-dimensional wave equation and using non-uniform modeled Alfvén velocity profiles. By comparing the results of the modeling alongside harmonics obtained from the Eskdalemuir, UK, data set from 2013 to 2021, the effects of the non-uniformity of the Alfvén velocity profile on the IAR are considered. We calculated the offset between the fundamental frequency and the harmonic frequency separation and found that this is not constant. From this parameter, we infer that the lower boundary condition of the electric field of the IAR is closest to a node, which agrees with previous studies. We compare the results of the non-uniform model with previous uniform models and evaluate their interpretations and the implications for the lower boundary condition
