or</span> i = I1+1:I2
0135                         k = k+1;
0136                         E(i,NRR-j) = elev(k);
0137                     <span class="keyword">end</span>
0138                 <span class="keyword">end</span>
0139                 
0140                 <span class="comment">% Interpolate the elevation onto the open boundary nodes</span>
0141                 <span class="keyword">if</span> ~exist(<span class="string">'SurfaceElevation'</span>,<span class="string">'var'</span>)
0142              