0001 function launch = imake_tail_yp(launch)
0002
0003
0004
0005 if ~isfield(launch,'tail')
0006 if length(launch.bNpar0) <= 10,
0007 tail = struct;
0008 tail.dtn = input_dke_yp('Tail dtn ',Inf);
0009 tail.mode = input_dke_yp('Tail mode spectrum: (0) no tail model, (1) Gaussian, (2) Linear, (3) calc_tail_jd ',3,0:3,'',[1,1]);
0010 tail.bNparmax_tail = input_dke_yp('Nparmax_tail: npar0 extend from npar_antenna to Nparmax_tail. If Nparmax_tail is 0, it is calculated by the code in main_rayinit_launch_jd.m -> 6.5/sqrt(Te0)',0);
0011 tail.bn_tail = input_dke_yp('n_tail: number of lobes to describe the tail. If n_tail = 0, no tail up to Nparmax_tail (no used if linear_tail is chosen)',5);
0012 tail.bP_tail = input_dke_yp('P_tail: power fraction in the tail. If NaN, it is calculated from n_tail and opt_tail (no used if linear_tail is chosen))',0.5);
0013 tail.bopt_tail = input_dke_yp('opt_tail: final power fraction of the tail @ Nparmax_tail (+1i) Adjust initial dNpar0 such that initial and tail lobes have the same width',1 + 1i);
0014
0015 launch.tail = tail;
0016
0017 if isfield(launch,'tail')
0018 if launch.tail.bn_tail > 0 && (launch.tail.bP_tail > 0 || isnan(launch.tail.bP_tail)),
0019 launch.id = [launch.id,'_TAIL_',make_tail_simulid_jd(launch.tail)];
0020 end
0021 end
0022 else
0023 disp('WARNING: the launch structure has more than ten lobes')
0024 end
0025 else
0026 disp('WARNING: a structure for the tail LH model already exists')
0027 end