levation (positive up) of observation point</span>
0085 <span class="comment">%-------------------------------------------------</span>
0086 zobs = -griddata(x,y,h,xobs,yobs); 
0087 
0088 <span class="comment">%-------------------------------------------------</span>
0089 <span class="comment">% loop over angles and depths, compute fetch</span>
0090 <span class="comment">%-------------------------------------------------</span>
0091 
0092 radline = 0:darc:max(max(x)-min(x),max(y)-min(y)); nRad = numel(radli