function [xElCount,yElCount,nodeCoords,elementNodes,x,y]=CreateMesh(x0,y0,xDim,yDim,elSize) 

[x,y]=meshgrid(x0:elSize:xDim,y0:elSize:yDim);
xElCount = size(x,2)-1;
yElCount = size(x,1)-1;
nodeCoords(:,1)=x(:); % x-coordinate of nodes 
nodeCoords(:,2)=y(:); % y-coordinate of nodes 
%Grid ID associated with the nodes for use in TDF function
elementNodes = zeros((xElCount*yElCount),5);
elementNodes(:,1,1) = 1:(xElCount*yElCount); %TDF Elem ID
elementNodes(:,2,1) = fix((elementNodes(:,1)-1)/yElCount)+elementNodes(:,1); %Node 1 ID
elementNodes(:,3,1) = elementNodes(:,2)+1; %Node 2 ID
elementNodes(:,4,1) = elementNodes(:,2)+yElCount+1; %Node 3 ID
elementNodes(:,5,1) = elementNodes(:,4)+1; %Node 4 ID

end