  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.
