s="comment">%fprintf('variable frame %s\n does not exist',vname)</span>
0398     <span class="keyword">else</span>
0399         <span class="comment">%compute wave induced bed stress using Soulsby formulas</span>
0400         <span class="comment">%friction factor (fw) is eq 62A</span>
0401         <span class="comment">%tau_w = 0.5*fw*Ubot^2</span>
0402         var1 = eval(vname1)';
0403         var2 = eval(vname2)';
0404         omega_wave = (2.0*pi./max(var1,.05))'; <span class="comment">%wave-orbital fr