%
exx=grdread('exx.grd');
eyy=grdread('eyy.grd');
exy=grdread('exy.grd');
edot=log10(1.e9*sqrt(exx.*exx+eyy.*eyy+2.*exy.*exy));
%
imagesc(flipud(edot),[1,3]);
colorbar
axis('equal')
