proc_wave_alphaphi_luke_jd

0001 function waves_out = proc_wave_alphaphi_luke_jd(output,waves,opt_abs)
0002 %
0003 if nargin < 3 || isempty(output),
0004     opt_abs = 0;% only correct data
0005 end
0006 if nargin < 2,
0007     if ~isempty(output),
0008         waves = output.dke_out.waves;
0009     else
0010         waves = {};
0011     end
0012 end
0013 %
0014 for iw = 1:length(waves),
0015     for ib = 1:length(waves{iw}.rays),
0016         ray = waves{iw}.rays{ib};
0017         %
0018         if isfield(ray,'stau') && ~isfield(ray,'stau_lin'),
0019             ray.stau_lin = ray.stau;
0020             ray = rmfield(ray,'stau');
0021         end
0022         %
0023         if isfield(ray,'sP_2piRp') && ~isfield(ray,'sP_2piRp_lin'),
0024             ray.sP_2piRp_lin = ray.sP_2piRp;
0025             ray = rmfield(ray,'sP_2piRp');
0026         end
0027         %
0028         if isfield(ray,'sdP_2piRp_ds') && ~isfield(ray,'sdP_2piRp_ds_lin'),
0029             ray.sdP_2piRp_ds_lin = ray.sdP_2piRp_ds;
0030             ray = rmfield(ray,'sdP_2piRp_ds');
0031         end
0032         %
0033         waves{iw}.rays{ib} = ray;
0034     end
0035 end
0036 %
0037 if opt_abs == 1,
0038     waves_out = alphaphi_luke_jd(waves,output.dke_out,output.ZP0,output.equilDKE,output.mksa);
0039 else
0040     waves_out = waves;
0041 end
0042 %