Recent experimental data at the Large Hadron Collider (LHC) confirmed that the soft sector azimuthal anisotropies in central XeXe collisions are sensitive to deformations in the wave function of the Xenon nucleus. Additionally, the CMS experiment found that D meson flow is slightly suppressed compared to other particle species when considering quark number scaling in small systems of pPb compared to PbPb. In this talk we used the D and B Mesons Modular (DAB-MOD) code coupled to Trento+v-USPhydro to calculate the $R_{AA}$ and $v_n$’s of D mesons using different energy loss models and Langevin techniques. Comparing D mesons $R_{AA}$ and $v_n$ in PbPb to XeXe collisions it is found that in central collisions D meson azimuthal anisotropies are sensitive to details in the nuclear structure wave function of Xenon. Additionally, we find that in mid-central collisions the smaller system size of XeXe suppresses D meson flow.