We previously showed that the adult ocellated lizard skin colour pattern is effectively generated by a stochastic cellular automaton (CA) of skin scales. We additionally suggested that the canonical continuous 2D reaction-diffusion (RD) process of colour pattern development is transformed into this discrete CA by reduced diffusion coefficients at the borders of scales (justified by the corresponding thinning of the skin). Here, we use RD numerical simulations in 3D on realistic lizard skin geometries and demonstrate that skin thickness variation on its own is sufficient to cause scale-by-scale coloration and CA dynamics during RD patterning. In addition, we show that this phenomenon is robust to RD model variation. Finally, using dimensionality-reduction approaches on large networks of skin scales, we show that animal growth affects the scale-colour flipping dynamics by causing a substantial decrease of the relative length scale of the labyrinthine colour pattern of the lizard skin.