_RL plumefrac(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
plumefrac(I,J) = 0. _d 0
plumefrac(I,J) = plumekb(1)*maskC(i,j,klev,bi,bj)
& +saltPlumeFlux(i,j,bi,bj)*plumefrac(I,J)*mass2rUnit
& plumefrac,'PLUMEKB ',kLev,1,2,bi,bj,myThid )