plumekb(1)=ABS(rF(k))
plumekb(2)=ABS(rF(k+1))
U plumekb,
plumefrac(I,J) = (plumekb(2)-plumekb(1))*maskC(i,j,k,bi,bj)
_RL plumekb(two2), SPdepth(two2)