RHOsumGlob_0=RHOsumGlob
& *(1. _d 0 - RHOsumGlob/RHOsumGlob_0)
& RHOsumGlob,
_RL VOLsumGlob_0, VOLsumGlob, RHOsumGlob_0, RHOsumGlob
RHOsumGlob=0. _d 0
CALL GLOBAL_SUM_TILE_RL( RHOsumTile, RHOsumGlob, myThid )
RHOsumGlob=RHOsumGlob/VOLsumGlob