[B(:).So]=deal(1e-4); % set the slope for the two reaches B(1).Chez=65; % set roughness for upstream reach B(2).Chez=55; % set roughness for downstream reach B(2).a0 ...
verifies the results against known expected values. Requires: Windows 10/11 with standard SDK (no extra installs needed).