keyword">end</span>;
0070 z0 = kn/30.;  <span class="comment">%hydraulic roughness</span>
0071 
0072 eval([<span class="string">'load '</span> matfile]);              <span class="comment">% load binary file containing SWAN results</span>
0073 <span class="comment">% obtained using BLOCK command with COMPGRID-set</span>
0074 <span class="comment">% load the connectivity</span>
0075 elefile=[basename <span class="string">'.ele'</span>];
0076 fid = fopen(elefile);                 <span class="comment">% load 