Discrete phase-space approach to mutually orthogonal Latin squares