This repository contains python scripts used to create visualizations for the paper "Multipolar origin of shape birefringence in metasurfaces", submitted to Optics Express. The output of this script is an interactive html file and/or static image in png format. You can select your preference by commenting out save_in_2D/3D function at the end of the script. To choose what multipoles you want to use for the visualization, please navigate to line 634 in get_E() function where you can select them by editing this line of code. Many examples are already present in the form of comments. The general formula is: total_E = total_E + get_E_p(element_position, dipole_px, pi) field of which multipole to calculate ↑ ↑ select multipole orientation ↑ ↑ always set "element position" ↑ select phase shift of multipole ↑ , where p stands for electric dipole, m for magnetic dipole, eq for electric quadrupole and mq for magnetic quadrupole.