A study is made of an anisotropic Potts model in three dimensions where the coupling depends on both the Potts state on each site and also the direction of the bond between them using both analytical and numerical methods. The phase diagram is mapped out for all values of the exchange interactions. Six distinct phases are identified. Monte Carlo simulations have been used to obtain the order parameter and the values for the energy and entropy in the ground state and also the transition temperatures. Excellent agreement is found between the simulated and analytic results. We find one region where there are two phase transitions with the lines meeting in a triple point. The orbital ordering that occurs in LaMnO 3 occurs as one of the ordered phases.