function domainplot(panelsL,panelsR,panels0,panelsF)

linewidth = 2;
markerize = 20;
fontsize = 18;

% Plots the domain with the structure nodes in red
set(gca,'FontSize',18,'LineWidth',linewidth)

for i=1:size(panelsL,1),
    plot(panelsL(i,1),panelsL(i,2),'.r','MarkerSize',markerize)
    plot(panelsL(i,3),panelsL(i,4),'.r','MarkerSize',markerize)
    plot([panelsL(i,1),panelsL(i,3)],[panelsL(i,2),panelsL(i,4)],'r','LineWidth',linewidth)
    hold on
end
% dsafsdaf
for i=1:size(panelsR,1),
       plot(panelsR(i,1),panelsR(i,2),'.r','MarkerSize',markerize)
    plot(panelsR(i,3),panelsR(i,4),'.r','MarkerSize',markerize)
    plot([panelsR(i,1),panelsR(i,3)],[panelsR(i,2),panelsR(i,4)],'r','LineWidth',linewidth)
end

for i=1:size(panels0,1),
       plot(panels0(i,1),panels0(i,2),'.b','MarkerSize',markerize)
    plot(panels0(i,3),panels0(i,4),'.b','MarkerSize',markerize)
    plot([panels0(i,1),panels0(i,3)],[panels0(i,2),panels0(i,4)],'b','LineWidth',linewidth)
end

for i=1:size(panelsF,1),
       plot(panelsF(i,1),panelsF(i,2),'.g','MarkerSize',markerize)
    plot(panelsF(i,3),panelsF(i,4),'.g','MarkerSize',markerize)
    plot([panelsF(i,1),panelsF(i,3)],[panelsF(i,2),panelsF(i,4)],'g','LineWidth',linewidth)
end

xlabel('x','FontSize',fontsize)
ylabel('z','FontSize',fontsize)
hold off
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%