_RL prof_traj1D(NLEVELMAX), prof_traj1D_mean(NLEVELMAX)
prof_traj1D_mean(k)=0.
call PROFILES_INTERP(prof_traj1D_mean,
prof_traj1D_mean(K)=0.
& *(prof_traj1D(K)-prof_data1D(K)-prof_traj1D_mean(K))