ne);
0093 fetch = zeros(nTheta,nZeta);
0094 hc = nodes2elems(Mobj.h,Mobj);
0095 <span class="keyword">for</span> i=1:nTheta
0096   fprintf(<span class="string">'searching theta %f\n'</span>,theta(i)*180/pi);
0097   xline = xobs+radline*cos(theta(i)+pi);
0098   yline = yobs+radline*sin(theta(i)+pi);
0099   k = 1;
0100   <span class="keyword">for</span> j=1:nZeta
0101     elev = zeta(j); depth = hc+elev;
0102     found = 0;
0103     k = 1;
0104     <span class="keyword">while</span> ~found
0105       cell = i