_RL locfc
locfc = 0.0
locfc = locfc + maskW(i,j,k,bi,bj)*
locfc = locfc + maskS(i,j,k,bi,bj)*
objf_eddytau(bi,bj) = locfc
print*,'objf_eddytau =',locfc
locfc = 0. _d 0
locfc = locfc + maskC(i,j,k,bi,bj)*
objf_tracer(bi,bj) = objf_tracer(bi,bj) + locfc
locfc = locfc + sum*dxG(i,j,bi,bj)