Nx = 2*102; Ny = 2*40; Lx = 0.05; Ly = 0.004;


inx = 2:(Nx-1); iny = 2:(Ny-1); 
north = iny+1;
east = inx+1; 
south = iny-1; 
west = inx-1;

dx = Lx/(Nx-2); dy = Ly/(Ny-2); 

x=linspace(0,Lx,Nx-2); 
y=linspace(0,Ly,Ny-2); 

[X,Y] = meshgrid(x,y); 

T = 293*ones(Ny,Nx);
CFLmax = 0.1;