function [X, Y, DIM] =  readgrid(file)

	TMP = importdata(file);
	DIM = TMP(2:3);
	DATA = TMP(4:length(TMP));
	X = transpose(reshape(DATA(1:DIM(1)*DIM(2)), DIM(1), DIM(2)));
	a = DIM(1)*DIM(2)+1;
	Y = transpose(reshape(DATA(a:length(DATA)), DIM(1), DIM(2)));
end


function [] = plotgrid(X,Y)
	hold on;
	plot(X, Y);
	plot(tranpose(X), transpose(Y));
end