Answer to: R code for a constrained randomization in a three-dimensional array?

This is excellent! I think I was unclear in the specification of the constraints, though: ideally the number of samples per treatment should be either 5 or 6 per X,Y coordinate for the treatments with more reps, and 3 or 4 for the treatments with fewer. That is reflected in the code I edited in to the answer.