function elemAreas = AreaCalc(H_sub,elementNodes,elSize)

    subdomainArea(:,1:4) = H_sub(elementNodes(:,2:5,1));
    elemAreas(:,1) = sum(subdomainArea(:,1:4),2)/4*((elSize)^2);
%     elemAreas = sum(subdomainArea(:,5));

end