raypath_jd

0001 function zF = raypath_jd(xyF,xyprop_dke,yb,r_dke);
0002 %
0003 [nr_dke,ny] = size(xyprop_dke);
0004 %
0005 xymask_prop_dke = xyprop_dke ~= 0;%DKE points visited by the ray
0006 %
0007 ib = 0;
0008 %
0009 for iy = 1:ny 
0010     if max(xyprop_dke(:,iy)) == 1 %ray portion is propagating outward
0011         ir_dke = min(find(xyprop_dke(:,iy) == 1)); %position of ray start
0012     %
0013     elseif min(xyprop_dke(:,iy)) == -1 %ray portion is propagating inward
0014         ir_dke = max(find(xyprop_dke(:,iy) == -1)); %position of ray start
0015         %
0016     else %no RT-DKE consistency seeked for current ray
0017         continue
0018     end
0019     ir = r_dke(ir_dke);
0020     %
0021     % check if new ray or only new portion
0022     %
0023     if ~isempty(find(iy == yb)); %new ray
0024         ib = ib + 1;
0025         zF{ib} = [];
0026     end
0027     %
0028     while ir_dke <= nr_dke & ir_dke >= 1 & xyprop_dke(ir_dke,iy) ~= 0 % stay on current ray
0029         ir = r_dke(ir_dke);
0030         %
0031         zF{ib} = [zF{ib},xyF(ir,iy)];
0032         %
0033         % next FS
0034         %
0035         ir_dke = ir_dke + xyprop_dke(ir_dke,iy);
0036     end
0037     %
0038     % interpolation for the FP positions (only necessary if dke_mode == 1)
0039     %
0040     %if dke_mode == 1
0041         %
0042     %end
0043 end
0044         
0045