mask_cur(k)=1
mask_cur(k)=0
prof_mask1D_cur(kk,bi,bj)=mask_cur(1)
prof_mask1D_cur(kk,bi,bj)=mask_cur(nr)
if (mask_cur(kcur+1).EQ.1.) then
if (mask_cur(kcur).EQ.0.) then
prof_mask1D_cur(kk,bi,bj)=mask_cur(kcur)
_RL traj_cur(nR),mask_cur(nR)