traj_cur(k)=0. _d 0
traj_cur(k)=traj_cur(k)
traj_cur_out(kk)=traj_cur(1)
traj_cur_out(kk)=traj_cur(nr)
traj_cur_out(kk)=(1-tmp_coeff)*traj_cur(kcur)
& +tmp_coeff*traj_cur(kcur+1)
traj_cur_out(kk)=traj_cur(kcur)
_RL traj_cur(nR),mask_cur(nR)