disp(strcat('Maximum Courant number (Co=dt*U/dx): ', num2str(dt*max(max(U))/dxmin)))
disp(strcat('Maximum CFL number (CFL=dt*nu/dx^2): ', num2str(dt*nu/dxmin^2)))