Raytracing though ionosphere

[lbbrhzn]

How would you do raytracing including ionosphere, for instance based on data from ITU P.531, NeQuick2, IRI
I guess it should be possible to add layers to the atmospheric model, correct?