_RL locfc,tmpC
tmpC = 0. _d 0
tmpC = tmpC + maskC(i,j,1,bi,bj)
_GLOBAL_SUM_RL( tmpC , myThid )
IF ( tmpC.GT.0. ) tmpC = 1. _d 0 / tmpC
locfc = locfc + tmpC*maskC(i,j,1,bi,bj)*