surface();
surface("surface.dat", center=true);