load uhov.out
a=fft2(uhov);
p=a.*conj(a);
p=log(p+1);
p2=fliplr(p);
p3(1:120,1:22)=p2(1:120,22:43);
p3(1:120,23:43)=p2(1:120,1:21);
x=(-21:21);
y=(1:120)./120;
contourf(x,y,p3);
