_RL delTrac
delTrac=Tracer(i,j,k,bi,bj)-Tracer(i,j,k-1,bi,bj)
& +weightA(i,j)*delTrac
& -weightB(i,j)*delTrac