
% Graph of 
% grid

    function meshgen(nf,nc)
    
    load mallaTTM.dat;
	 for ii=1:nf
	 	for j=1:nc
		   X(ii,j) = mallaTTM(j+(ii-1)*nc,1);
		   Y(ii,j) = mallaTTM(j+(ii-1)*nc,2);
		   U(ii,j) = 0;
		end
   end
   figure
   orient landscape
   mesh(X,Y,U)
   grid on
   xlabel('x')
	  ylabel('y')
    view(0,90)

    load mallainicial.dat;
	 for ii=1:nf
	 	for j=1:nc
		   X(ii,j) = mallainicial(j+(ii-1)*nc,1);
		   Y(ii,j) = mallainicial(j+(ii-1)*nc,2);
		   U(ii,j) = 0;
		end
   end
   figure
   orient landscape
   mesh(X,Y,U,'EdgeColor','black')
   grid on
   xlabel('x')
	 ylabel('y')
   view(0,90)
    
    load mallac1.dat;
	 for ii=1:nf
	 	for j=1:nc
		   X(ii,j) = mallac1(j+(ii-1)*nc,1);
		   Y(ii,j) = mallac1(j+(ii-1)*nc,2);
		   U(ii,j) = 0;
		end
   end
   figure
   orient landscape
   mesh(X,Y,U,'EdgeColor','black')
   grid on
   xlabel('x')
	 ylabel('y')
   view(0,90)
 
