Module-level functions for size of LM and LMpM arrays
Add Wigner_D_matrices function to create and return the D array