There is a problem in the reconstruction using AOtools KL modes. I have generated 1000 atmospheric phase screens and reconstructed them using a 'perfect' wavefront sensor that just flattens the phase. I looked at the variances of the reconstruction modes and compared them with KL modes generated by the OOMAO simulation :

The variances should be both smooth. The peaks in the AOtools plot correspond to radial modes.