Constructing non-cyclic pandiagonal Latin squares of prime orders

Let p   be a prime of the form p=6m+1p=6m+1, where m is a positive integer. In this paper we propose an algorithm for constructing non-cyclic pandiagonal Latin squares of order p   for p>7p>7.