sumEmP = 0.
CALL GLOBAL_SUM_TILE_RL( tileEmP, sumEmP, myThid )
& freeSurfFac*deltaTFreeSurf*mass2rUnit*sumEmP/globalArea
_RL sumEmP, tileEmP(nSx,nSy)